Nothing Special   »   [go: up one dir, main page]

CN112116922A - 一种噪声盲源信号分离方法、终端设备及存储介质 - Google Patents

一种噪声盲源信号分离方法、终端设备及存储介质 Download PDF

Info

Publication number
CN112116922A
CN112116922A CN202010982357.XA CN202010982357A CN112116922A CN 112116922 A CN112116922 A CN 112116922A CN 202010982357 A CN202010982357 A CN 202010982357A CN 112116922 A CN112116922 A CN 112116922A
Authority
CN
China
Prior art keywords
signal
phase
optimization interval
optimal
correlation coefficient
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202010982357.XA
Other languages
English (en)
Other versions
CN112116922B (zh
Inventor
周海峰
方艺鹏
张恩来
林开荣
郑东强
林忠华
杨志荣
焦健
王新乡
陈虹宇
汪海志
包素钦
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Jimei University
Original Assignee
Jimei University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Jimei University filed Critical Jimei University
Priority to CN202010982357.XA priority Critical patent/CN112116922B/zh
Publication of CN112116922A publication Critical patent/CN112116922A/zh
Application granted granted Critical
Publication of CN112116922B publication Critical patent/CN112116922B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G10—MUSICAL INSTRUMENTS; ACOUSTICS
    • G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
    • G10L21/00—Speech or voice signal processing techniques to produce another audible or non-audible signal, e.g. visual or tactile, in order to modify its quality or its intelligibility
    • G10L21/02—Speech enhancement, e.g. noise reduction or echo cancellation
    • G10L21/0272—Voice signal separating
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00—Pattern recognition
    • G06F18/20—Analysing
    • G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
    • G06F18/213—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods
    • G06F18/2134—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods based on separation criteria, e.g. independent component analysis
    • G—PHYSICS
    • G10—MUSICAL INSTRUMENTS; ACOUSTICS
    • G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
    • G10L21/00—Speech or voice signal processing techniques to produce another audible or non-audible signal, e.g. visual or tactile, in order to modify its quality or its intelligibility
    • G10L21/02—Speech enhancement, e.g. noise reduction or echo cancellation
    • G10L21/0272—Voice signal separating
    • G10L21/0308—Voice signal separating characterised by the type of parameter measurement, e.g. correlation techniques, zero crossing techniques or predictive techniques

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Signal Processing (AREA)
  • Health & Medical Sciences (AREA)
  • Audiology, Speech & Language Pathology (AREA)
  • Human Computer Interaction (AREA)
  • Quality & Reliability (AREA)
  • Acoustics & Sound (AREA)
  • Multimedia (AREA)
  • Computational Linguistics (AREA)
  • Data Mining & Analysis (AREA)
  • Theoretical Computer Science (AREA)
  • Artificial Intelligence (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Evolutionary Biology (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Complex Calculations (AREA)
  • Noise Elimination (AREA)

Abstract

本发明涉及一种噪声盲源信号分离方法、终端设备及存储介质，该方法中包括：S1：将观测信号通过Fastica算法得到解混矩阵和分离后的独立分量；S2：对独立分量做多次相位调整，并根据相位调整前后对应的观测信号之间的总相关系数的大小确定最佳调整相位；S3：查找权值调整函数的最佳寻优区间，从最佳寻优区间内查找最佳权值调整函数；S4：将独立分量进行相位和幅值的调整得到调整独立分量，将其与观测信号均进行傅里叶变换得到对应的频谱，计算每个调整独立分量频谱和其对应的观测信号频谱的总相关系数，按照总相关系数的大小将所有调整独立分量进行排序，得到分离后的信号源。该解决方案对波形的不确定性进行消除，达到较为完美的分离结果。

Description

一种噪声盲源信号分离方法、终端设备及存储介质
技术领域
本发明涉及噪声信号分离领域，尤其涉及一种噪声盲源信号分离方法、终端设备及存储介质。
背景技术
随着经济的发展，越来越多的大功率机械出现，使得噪声污染变得日益的严重。噪声污染对人体的伤害具有不可逆性，长期处于噪声下会对人的听力噪声损害以及诱发心血管疾病等。近年来，人们对于噪声污染的重视逐渐的加强，采取了很多控制噪声的措施。针对于这些大功率机械及动力设备的噪声控制，最有效的手段就是从噪声源入手。如何准确的识别噪声源以及确定噪声源的主次，是最关键的一个步骤。经研究发现，对于大功率机械及动力设备的噪声源，由于其结构复杂，激励源众多且不可预测，是属于典型的盲源信号。所以，识别噪声源的前提是将盲源信号进行分离。
独立分量分析(Independent component analysis，ICA)算法是近几十年来在盲源分离方面运用最为广泛而且有效的方法。这个算法最早是由Jutten C,Herault J在1986年提出来的，被广泛的应用语音信号处理，信号特征处理等领域。1997年Hyvarinen基于ICA算法的基础，研究并提出一种新的ICA方法，即快速独立分量分析(Fast independentcomponent analysis，FastICA)。该算法相比于ICA而言，收敛速度更快。但是，独立分量分析在信号分析当中存在一些不足的地方。就是经由该算法进行分离的信号在幅值，相位以及排序等方面存在不确定性。由于噪声的贡献度与能量大小有着密切的关系，所以对分离后的独立分量进行不确定性校正对噪声源的识别以及确定噪声源的主次有着重大的意义。
发明内容
为了解决上述问题，本发明提出了一种噪声盲源信号分离方法、终端设备及存储介质。
具体方案如下：
一种噪声盲源信号分离方法，包括以下步骤：
S1：将由多个信号源发出的包含噪声的观测信号通过Fastica算法得到解混矩阵和分离后的独立分量；
S2：对独立分量做多次相位调整，并根据相位调整前后对应的观测信号之间的总相关系数的大小确定最佳调整相位；
S3：查找权值调整函数的最佳寻优区间，从最佳寻优区间内查找最佳权值调整函数，使得通过最佳权值调整函数加权得到的混合矩阵计算得到的观测信号与通过Fastica算法处理前的观测信号的总相关系数小于阈值；
S4：将独立分量通过步骤S2得到的最佳调整相位和步骤S3得到的最佳权值调整函数进行相位和幅值的调整得到调整独立分量，将调整独立分量和观测信号均进行傅里叶变换得到对应的频谱，计算每个调整独立分量频谱和其对应的观测信号频谱的总相关系数，按照总相关系数的大小将所有调整独立分量进行排序，得到最终分离后的信号源。
进一步的，步骤S2具体包括以下步骤：
S21：对分离出来的独立分量做多次相位调整；
Figure BDA0002688014900000031
其中，yj (M)为第j个独立分量在第M次相位调整后的值，j为独立分量的序号，yj (0)为通过Fastica算法得到的第j个独立分量的值，εM为第M次相位调整函数；
S22：根据相位调整后的独立分量，计算相位调整后对应的观测信号：
Figure BDA0002688014900000032
其中，xi (M)为第i个观测信号在第M次相位调整后的值，i为观测信号的序号，m为独立分量的总数，aij (0)为的初始混合矩阵中第i行第j列的值，初始混合矩阵为通过Fastica算法得到的解混矩阵求逆后得到的矩阵；
S23：针对每次相位调整函数εM，计算相位调整后对应的观测信号xi (M)与相位调整前对应的观测信号xi之间的总相关系数ACεM(xi,xi (M))：
Figure BDA0002688014900000033
其中，Xcor(·)为自相关系数，n为观测信号的总数；
S24：根据最大相似准则，将总相关系数的最大值所对应相位调整函数作为最佳调整相位。
进一步的，权值调整函数的最佳寻优区间的查找方法为：
S31：初始设定参数k＝1，设定权值调整函数的第k个寻优区间；
S32：将混合矩阵乘以第k个寻优区间的最大值和最小值得到对应的最大混合矩阵和最小混合矩阵；
S33：根据最大混合矩阵、最小混合矩阵和独立分量，计算第k个最大混合矩阵对应的最大观测信号与观测信号的总相关系数和最小混合矩阵对应的最小观测信号与观测信号的总相关系数的差值；判断第k个差值的符号与第k-1个差值的符号是否相同，如果是，进入S34；否则，进入S38；
S34：判断第k个差值的符号与第k个寻优区间的最大值与最小值的差值的符号是否相同，如果是，进入S35；否则，进入S36；
S35：设定第k+1个寻优区间，且令第k+1个寻优区间的最小值为第k个寻优区间的最大值，进入S37；
S36：设定第k+1个寻优区间，且令第k+1个寻优区间的最大值为第k个寻优区间的最小值，进入S37；
S37：令k＝k+1，返回S32；
S38：将第k-1个和第k个寻优区间组成最佳寻优区间。
进一步的，从最佳寻优区间内查找最佳权值调整函数的方法为在最佳寻优区间内采用黄金分割算法进行逐级查找。
一种噪声盲源信号分离终端设备，包括处理器、存储器以及存储在所述存储器中并可在所述处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现本发明实施例上述的方法的步骤。
一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，所述计算机程序被处理器执行时实现本发明实施例上述的方法的步骤。
本发明采用如上技术方案，实现将混合的波形信号分离出源信号波形，同时将波形的不确定性进行消除，以达到较为完美的分离结果。
附图说明
图1所示为本发明实施例一的流程图。
图2所示为该实施例中信号混合与ICA分离过程示意图。
图3所示为该实施例中不确定性校正过程的流程图。
图4所示为该实施例中四个典型的源信号波形图。
图5所示为该实施例中混合后的观测信号波形图。
图6所示为该实施例中通过Fastica算法处理得到的信号波形图。
图7所示为该实施例中经过不确定性校正后得到的信号波形图。
图8所示为该实施例中语音信号的源信号波形图。
图9所示为该实施例中语音信号的观测信号波形图。
图10所示为该实施例中语音信号的分离信号波形图。
图11所示为该实施例中语音信号经过不确定性校正后的信号波形图。
具体实施方式
为进一步说明各实施例，本发明提供有附图。这些附图为本发明揭露内容的一部分，其主要用以说明实施例，并可配合说明书的相关描述来解释实施例的运作原理。配合参考这些内容，本领域普通技术人员应能理解其他可能的实施方式以及本发明的优点。
现结合附图和具体实施方式对本发明进一步说明。
实施例一：
本发明实施例提供了一种噪声盲源信号分离方法，如图1所示，所述方法包括以下步骤：
S1：将由多个信号源发出的包含噪声的观测信号通过Fastica算法得到解混矩阵和去除噪声后的独立分量。
1.ICA原理与不确定性
对于独立分量分析的描述，最经典的就是“鸡尾酒宴会”问题。假设宴会中装有n个声音收录器，宴会中的人数记为m个且他们同时说话。那么在宴会结束之后，n个声音收录器里面收集了一组数据，这组数据记为观测信号X，X＝[x1(t),x2(t),x3(t),...,xn(t)]T；t＝(1,...,m)。独立分量分析的作用就是分辨出m维数组的源，也就是每个人发出的声音源。
其中m个信号源都是独立的，对应的有n个传感器。这里不考虑环境中的其他噪音影响，传感器所检测到的观测信号与源信号通过线性混合，则有如下的公式：
Figure BDA0002688014900000061
式中，S(t)＝[s1(t),s2(t),...,sm(t)]T为源信号，源信号无法直接识别。该源信号是由m个互相独立的源信号构成。aij∈A是n×m维的混合矩阵。其中需要满足n>m这个条件。式中有A，S(t)均是未知的，仅有X(t)为已知的信号。独立分量分析就是需要找到一个解混矩阵W(W＝A-1)使得解混后得到一个最接近源信号S(t)的矩阵Y(t)，于是可以得到下式：
Y(t)＝WAS(t)＝A-1AS(t) (2)
式中，Y(t)为分离出的独立分量，其信号特征接近于源信号。图2所示为解混的过程。仅有当W为A的逆矩阵的时候才能实现源信号的分离。
ICA模型在应用的时候存在有几点问题。
首先，因为源信号所经历的混合过程无法进行预测，要直接进行信号的分离是一件比较困难的事情。所以在应用ICA模型求分离的唯一解时候，需要有以下三个基本前提条件：
(1)每个源信号必须的互相统计独立的信号，即其中一个源信号出现的概率不受其它源信号的影响；
(2)最多只能存在一个源信号服从高斯分布。因为高斯信号在线性叠加的结果不会改变其高斯性，源信号若有两个以上服从高斯分布，那么将无法使用ICA模型进行分离。
(3)传感器的数量n不能够小于源信号的数量m。在实际应用当中，噪声源的数量无法确定，所以通常做出二者相等的假设情况。
其次就是ICA的不确定性。虽然分离结果能够保证唯一的，但结果仍存在一些不确定性。比如幅值和相位会出现偏差，这是由于A与S均为未知的。如式(3)所示，当变换过程当中出现S增大G倍，则相应的A只需要除以G倍，得出的结果不被影响，但是估计的源信号却发生了变化。
Figure BDA0002688014900000071
估计的源信号分量在分离之后有出现排序方面的问题。这是由于分离过程当中会插入一个置换矩阵P，使得X＝APP-1S(t)，排序会随之发生改变。
2.Fastica算法
Fastica算法，也被称为基于固定点(Fixed-point)的算法。该算法包含基于负熵最大、基于最大似然和基于峭度等三种不同的形式。该实施例中采用的是基于负熵最大的Fastica算法作为基础。其负熵的近似值函数表达式如下：
Figure BDA0002688014900000081
式中，ki为常量，E[·]为计算数学期望，yi和vi为两个高斯随机变量，两个变量的协方差值是相等的，在使用变量之前，需要先进行归一化处理。G为非二次函数，选择合适的G，得到的估计量会相对稳定。该实施例中设定G(y)为：
Figure BDA0002688014900000082
式中，a1的取值一般是[1,2]之间，该实施例中取值为1，a2取1。其中，函数G1适合应用在超高斯信号与亚高斯信号都存在的情形，G2适用于超高斯信号源，G3适用于亚高斯信号源。因为有y＝wi Tx是一个独立分量，其中wi是解混矩阵W中的一个n维向量。由K-T条件可知，要让式(4)的负熵达到最大化，仅有当E{(wi Tx)2}＝||Wi||2＝1时，才能达到最优值满足：
Figure BDA0002688014900000083
式中，g(·)是G(·)的一阶导，系数β为拉格朗日乘子，这里假设左侧为函数F，则可以得到雅克比(Jacobi)矩阵：
JF(wi)＝E[x·xTg'(wi Tx)]-βI (7)
式(7)的右边第一项近似为：
E[x·xTg'(wi Tx)]≈E[x·xT]E[g'(wi Tx)]＝E[g'(wi Tx)]I (8)
由上式的近似结果，可以使得雅可比矩阵转化为对角阵，从而能更容易的求得相应的逆矩阵。由此可以易得牛顿迭代式：
Figure BDA0002688014900000091
式(9)两边同时乘以β＝E[xg(wT(i)x)]-E[xg'(wT(i)x)]w(i)，经过化简之后能够得到相应的FastICA算法的迭代公式如下式所示：
w(i+1)＝E[xg(wT(i)x)]-E[xg'(wT(i)x)]w(i) (10)
Figure BDA0002688014900000092
通过式(11)的迭代式可以得到FastICA的算法过程如下；
1)对观测信号X计算期望，并减去相应的期望值，此步也叫中心化处理；
2)对X做白化处理，做正交变换使得E[XTX]＝I；
3)选取待估计的数量m，并设定迭代次数i；
4)任取一个初始向量wi，该向量必须存在单位范数；
5)令w(i+1)＝E[xg(wT(i)x)]-E[xg'(wT(i)x)]w(i)，其中，g(·)原函数选取参见式(5)；
6)wi标准化；
7)迭代后判断wi是否收敛，若出现不收敛的情况则返回步骤5)，若收敛则继续执行接下去的程序；
8)令i＝i+1，当i≤m时，跳回第4)步，否则结束算法程序。
由于经过Fastica算法所分离出的独立分量yi存在幅值，相位和排序等方面的不确定性问题，特别是幅值和相位的问题，将会对噪声源的识别产生直接的影响，因此下面步骤中为对其结果进行处理，达到校正的目的，如图3所示。
S2：对独立分量做多次相位调整，并根据相位调整前后对应的观测信号之间的总相关系数的大小确定最佳调整相位。
步骤S2具体包括以下步骤：
S21：对分离出来的独立分量做多次相位调整；
Figure BDA0002688014900000103
其中，yj (M)为第j个独立分量在第M次相位调整后的值，j为独立分量的序号，yj (0)为通过Fastica算法得到的第j个独立分量的值，εM为第M次相位调整函数。
S22：根据相位调整后的独立分量，计算相位调整后对应的观测信号：
Figure BDA0002688014900000101
其中，xi (M)为第i个观测信号在第M次相位调整后的值，i为观测信号的序号，m为独立分量的总数，aij (0)为的初始混合矩阵中第i行第j列的值，初始混合矩阵为通过Fastica算法得到的解混矩阵求逆后得到的矩阵；
S23：针对每次相位调整函数εM，计算相位调整后对应的观测信号xi (M)与相位调整前对应的观测信号xi之间的总相关系数
Figure BDA0002688014900000104
Figure BDA0002688014900000102
其中，Xcor(·)为自相关系数，n为观测信号的总数；
S24：根据最大相似准则，将总相关系数的最大值所对应相位调整函数作为最佳调整相位。
当总相关系数最大时，认为此时的相位调整至正常。
S3：查找权值调整函数的最佳寻优区间，从最佳寻优区间内查找最佳权值调整函数，使得通过最佳权值调整函数加权得到的混合矩阵计算得到的观测信号与通过Fastica算法处理前的观测信号的总相关系数小于阈值。
该实施例中从最佳寻优区间内查找最佳权值调整函数的方法优选采用在最佳寻优区间内采用黄金分割算法进行逐级查找。
经由Fastica算法的白化过程，使得分离后的独立分量在幅值上呈现的不确定性包含放大和缩小两种可能。故需要对混合矩阵乘以一个最佳权值调整函数进行调整。为了加快权值调整函数的查找速度，该实施例中采用以下步骤查找权值调整函数的最佳寻优区间。
S31：初始设定参数k＝1，设定权值调整函数的第k个寻优区间。
S32：将混合矩阵乘以第k个寻优区间的最大值和最小值得到对应的最大混合矩阵和最小混合矩阵。
S33：根据最大混合矩阵、最小混合矩阵和独立分量，计算第k个最大混合矩阵对应的最大观测信号与观测信号的总相关系数和最小混合矩阵对应的最小观测信号与观测信号的总相关系数的差值；判断第k个差值的符号与第k-1个差值的符号是否相同，如果是，进入S34；否则，进入S38。
两个观测信号之间的总相关系数的计算公式为：
Figure BDA0002688014900000111
其中，xi,x′i分别表示两个观测信号。
S34：判断第k个差值的符号与第k个寻优区间的最大值与最小值的差值的符号是否相同，如果是，进入S35；否则，进入S36。
S35：设定第k+1个寻优区间，且令第k+1个寻优区间的最小值为第k个寻优区间的最大值，进入S37。
S36：设定第k+1个寻优区间，且令第k+1个寻优区间的最大值为第k个寻优区间的最小值，进入S37。
S37：令k＝k+1，返回S32。
S38：将第k-1个和第k个寻优区间组成最佳寻优区间。
S4：将独立分量通过步骤S2得到的最佳调整相位和步骤S3得到的最佳权值调整函数进行相位和幅值的调整得到调整独立分量，将调整独立分量和观测信号均进行傅里叶变换得到对应的频谱，计算每个调整独立分量频谱和其对应的观测信号频谱的总相关系数，按照总相关系数的大小将所有调整独立分量进行排序，得到最终分离后的信号源。
本发明实施例通过上述步骤对独立分量进行校正，可以准确的得到源信号的信息。而在处理噪声的过程当中，校正信号成为了关键的一步。这是由于噪声信号与能量之间有着直接的关系，源信号估计的准确度，会直接影响后续噪声的处理。
试验验证：
该实施例中采用四个较为典型的源信号Sn(t)(n＝1,2,3,4)，通过本实施例方法进行实验验证。其中，如图4所示，典型信号分别是S1余弦信号，S2方波信号，S3正弦波信号，S4锯齿波信号。每个信号波形均设置200个样本点。混合矩阵A采用的是4×4的随机方阵。将源信号Sn(t)与矩阵A进行混合，可以得出观测信号Xn(t)。再通过Fastica算法进行迭代计算，得到相应的分离后的独立分量Yn(t)。
将源信号与矩阵A混合后，得到的观测信号如图5所示，可以看出每一个源信号在经过与随机矩阵A的混合之后，波形明显发生变化。
通过Fastica算法处理信号，分离后得到的信号如图6所示。由于Fastica算法处理信号的过程中对信号做了白化处理，所以在分离之后出现幅值，相位以及排序发生了改变。
经过本实施例方法上述S2-S4步骤的计算，将源信号的不确定性加以消除，包括幅值，相位以及排序等。如图7所示，是经过校正后得到的信号波形。校正后波形基本上恢复至与源信号相同。
为了能够更进一步的对本实施例方法进行验证，该实施例中还采用了三段语音信号进行混合，分离。其中语音信号的采样频率为44.1KHz，采样点数为441000点。三段语音信号与3维随机矩阵A进行混合，得出相应的观测信号。图8所示为语音信号的源信号波形，图9所示为语音信号的观测信号波形。图10所示，为语音信号的分离信号波形，由波形可以看出经过Fastica的白化处理过程，分离后所得出的源信号波形，幅值，相位，排序均发生了变化。图11所示为语音信号经过不确定性校正后的信号波形。由波形可以看出，该源信号波形的幅值，相位以及排序问题基本得到了校正。
通过上述的仿真实验证实了该实施例方法可以实现将混合的波形信号分离出源信号波形，同时将波形的不确定性进行消除，以达到较为完美的分离结果。这一结论为噪声源的识别提供了更加准确的分离信号，使得噪声源识别中可能存在的能量减弱等问题能够被降到最低，以达到源信号最高程度的保真，为噪声源识别研究提供更加可靠的技术支持。
实施例二：
本发明还提供一种噪声盲源信号分离终端设备，包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现本发明实施例一的上述方法实施例中的步骤。
进一步地，作为一个可执行方案，所述噪声盲源信号分离终端设备可以是桌上型计算机、笔记本、掌上电脑及云端服务器等计算设备。所述噪声盲源信号分离终端设备可包括，但不仅限于，处理器、存储器。本领域技术人员可以理解，上述噪声盲源信号分离终端设备的组成结构仅仅是噪声盲源信号分离终端设备的示例，并不构成对噪声盲源信号分离终端设备的限定，可以包括比上述更多或更少的部件，或者组合某些部件，或者不同的部件，例如所述噪声盲源信号分离终端设备还可以包括输入输出设备、网络接入设备、总线等，本发明实施例对此不做限定。
进一步地，作为一个可执行方案，所称处理器可以是中央处理单元(CentralProcessing Unit，CPU)，还可以是其他通用处理器、数字信号处理器(Digital SignalProcessor，DSP)、专用集成电路(Application Specific Integrated Circuit，ASIC)、现场可编程门阵列(Field-Programmable Gate Array，FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件等。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等，所述处理器是所述噪声盲源信号分离终端设备的控制中心，利用各种接口和线路连接整个噪声盲源信号分离终端设备的各个部分。
所述存储器可用于存储所述计算机程序和/或模块，所述处理器通过运行或执行存储在所述存储器内的计算机程序和/或模块，以及调用存储在存储器内的数据，实现所述噪声盲源信号分离终端设备的各种功能。所述存储器可主要包括存储程序区和存储数据区，其中，存储程序区可存储操作系统、至少一个功能所需的应用程序；存储数据区可存储根据手机的使用所创建的数据等。此外，存储器可以包括高速随机存取存储器，还可以包括非易失性存储器，例如硬盘、内存、插接式硬盘，智能存储卡(Smart Media Card,SMC)，安全数字(Secure Digital,SD)卡，闪存卡(Flash Card)、至少一个磁盘存储器件、闪存器件、或其他易失性固态存储器件。
本发明还提供一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，所述计算机程序被处理器执行时实现本发明实施例上述方法的步骤。
所述噪声盲源信号分离终端设备集成的模块/单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明实现上述实施例方法中的全部或部分流程，也可以通过计算机程序来指令相关的硬件来完成，所述的计算机程序可存储于一计算机可读存储介质中，该计算机程序在被处理器执行时，可实现上述各个方法实施例的步骤。其中，所述计算机程序包括计算机程序代码，所述计算机程序代码可以为源代码形式、对象代码形式、可执行文件或某些中间形式等。所述计算机可读介质可以包括：能够携带所述计算机程序代码的任何实体或装置、记录介质、U盘、移动硬盘、磁碟、光盘、计算机存储器、只读存储器(ROM，Read-OnlyMemory)、随机存取存储器(RAM，Random Access Memory)以及软件分发介质等。
尽管结合优选实施方案具体展示和介绍了本发明，但所属领域的技术人员应该明白，在不脱离所附权利要求书所限定的本发明的精神和范围内，在形式上和细节上可以对本发明做出各种变化，均为本发明的保护范围。

Claims (6)

1.一种噪声盲源信号分离方法，其特征在于，包括以下步骤：
S1：将由多个信号源发出的包含噪声的观测信号通过Fastica算法得到解混矩阵和分离后的独立分量；
S2：对独立分量做多次相位调整，并根据相位调整前后对应的观测信号之间的总相关系数的大小确定最佳调整相位；
S3：查找权值调整函数的最佳寻优区间，从最佳寻优区间内查找最佳权值调整函数，使得通过最佳权值调整函数加权得到的混合矩阵计算得到的观测信号与通过Fastica算法处理前的观测信号的总相关系数小于阈值；
S4：将独立分量通过步骤S2得到的最佳调整相位和步骤S3得到的最佳权值调整函数进行相位和幅值的调整得到调整独立分量，将调整独立分量和观测信号均进行傅里叶变换得到对应的频谱，计算每个调整独立分量频谱和其对应的观测信号频谱的总相关系数，按照总相关系数的大小将所有调整独立分量进行排序，得到最终分离后的信号源。
2.根据权利要求1所述的噪声盲源信号分离方法，其特征在于：步骤S2具体包括以下步骤：
S21：对分离出来的独立分量做多次相位调整；
yj (M)＝yj (0)|εM
其中，yj (M)为第j个独立分量在第M次相位调整后的值，j为独立分量的序号，yj (0)为通过Fastica算法得到的第j个独立分量的值，εM为第M次相位调整函数；
S22：根据相位调整后的独立分量，计算相位调整后对应的观测信号：
Figure FDA0002688014890000021
其中，xi (M)为第i个观测信号在第M次相位调整后的值，i为观测信号的序号，m为独立分量的总数，aij (0)为的初始混合矩阵中第i行第j列的值，初始混合矩阵为通过Fastica算法得到的解混矩阵求逆后得到的矩阵；
S23：针对每次相位调整函数εM，计算相位调整后对应的观测信号xi (M)与相位调整前对应的观测信号xi之间的总相关系数
Figure FDA0002688014890000022
Figure FDA0002688014890000023
其中，Xcor(·)为自相关系数，n为观测信号的总数；
S24：根据最大相似准则，将总相关系数的最大值所对应相位调整函数作为最佳调整相位。
3.根据权利要求1所述的噪声盲源信号分离方法，其特征在于：权值调整函数的最佳寻优区间的查找方法为：
S31：初始设定参数k＝1，设定权值调整函数的第k个寻优区间；
S32：将混合矩阵乘以第k个寻优区间的最大值和最小值得到对应的最大混合矩阵和最小混合矩阵；
S33：根据最大混合矩阵、最小混合矩阵和独立分量，计算第k个最大混合矩阵对应的最大观测信号与观测信号的总相关系数和最小混合矩阵对应的最小观测信号与观测信号的总相关系数的差值；判断第k个差值的符号与第k-1个差值的符号是否相同，如果是，进入S34；否则，进入S38；
S34：判断第k个差值的符号与第k个寻优区间的最大值与最小值的差值的符号是否相同，如果是，进入S35；否则，进入S36；
S35：设定第k+1个寻优区间，且令第k+1个寻优区间的最小值为第k个寻优区间的最大值，进入S37；
S36：设定第k+1个寻优区间，且令第k+1个寻优区间的最大值为第k个寻优区间的最小值，进入S37；
S37：令k＝k+1，返回S32；
S38：将第k-1个和第k个寻优区间组成最佳寻优区间。
4.根据权利要求3所述的噪声盲源信号分离方法，其特征在于：从最佳寻优区间内查找最佳权值调整函数的方法为在最佳寻优区间内采用黄金分割算法进行逐级查找。
5.一种噪声盲源信号分离终端设备，其特征在于：包括处理器、存储器以及存储在所述存储器中并在所述处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现如权利要求1～4中任一所述方法的步骤。
6.一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，其特征在于：所述计算机程序被处理器执行时实现如权利要求1～4中任一所述方法的步骤。
CN202010982357.XA 2020-09-17 2020-09-17 一种噪声盲源信号分离方法、终端设备及存储介质 Active CN112116922B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010982357.XA CN112116922B (zh) 2020-09-17 2020-09-17 一种噪声盲源信号分离方法、终端设备及存储介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010982357.XA CN112116922B (zh) 2020-09-17 2020-09-17 一种噪声盲源信号分离方法、终端设备及存储介质

Publications (2)

Publication Number Publication Date
CN112116922A true CN112116922A (zh) 2020-12-22
CN112116922B CN112116922B (zh) 2024-04-12

Family

ID=73799872

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010982357.XA Active CN112116922B (zh) 2020-09-17 2020-09-17 一种噪声盲源信号分离方法、终端设备及存储介质

Country Status (1)

Country Link
CN (1) CN112116922B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112869737A (zh) * 2021-02-01 2021-06-01 浙江大学山东工业技术研究院 一种非接触式人体血氧饱和度检测方法
CN113223548A (zh) * 2021-05-07 2021-08-06 北京小米移动软件有限公司 声源定位方法及装置
CN113948104A (zh) * 2021-12-07 2022-01-18 思必驰科技股份有限公司 盲语音信号分离方法、电子设备和存储介质

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003078423A (ja) * 2001-09-03 2003-03-14 Clarion Co Ltd ブラインド信号分離処理装置
US20090086998A1 (en) * 2007-10-01 2009-04-02 Samsung Electronics Co., Ltd. Method and apparatus for identifying sound sources from mixed sound signal
CN103679199A (zh) * 2013-12-11 2014-03-26 中国石油大学(华东) 基于入侵性杂草算法的有噪独立分量分析方法
CN104485113A (zh) * 2014-12-23 2015-04-01 长沙理工大学 一种多故障源声发射信号分离方法
CN105355212A (zh) * 2015-10-14 2016-02-24 天津大学 一种稳健的欠定盲分离源数及混合矩阵估计方法及装置
WO2018059406A1 (zh) * 2016-09-29 2018-04-05 合肥华凌股份有限公司 盲音分离方法、结构及语音控制系统和电器总成
CN109243483A (zh) * 2018-10-17 2019-01-18 西安交通大学 一种含噪频域卷积盲源分离方法
CN109973251A (zh) * 2019-03-24 2019-07-05 西安航天动力测控技术研究所 一种用于固体火箭发动机振动信号降噪的盲源分离方法
CN110263876A (zh) * 2019-06-27 2019-09-20 南京信息工程大学 一种雨声信号降噪处理方法和系统
CN110807524A (zh) * 2019-11-13 2020-02-18 大连民族大学 单通道信号盲源分离幅度校正方法

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003078423A (ja) * 2001-09-03 2003-03-14 Clarion Co Ltd ブラインド信号分離処理装置
US20090086998A1 (en) * 2007-10-01 2009-04-02 Samsung Electronics Co., Ltd. Method and apparatus for identifying sound sources from mixed sound signal
CN103679199A (zh) * 2013-12-11 2014-03-26 中国石油大学(华东) 基于入侵性杂草算法的有噪独立分量分析方法
CN104485113A (zh) * 2014-12-23 2015-04-01 长沙理工大学 一种多故障源声发射信号分离方法
CN105355212A (zh) * 2015-10-14 2016-02-24 天津大学 一种稳健的欠定盲分离源数及混合矩阵估计方法及装置
WO2018059406A1 (zh) * 2016-09-29 2018-04-05 合肥华凌股份有限公司 盲音分离方法、结构及语音控制系统和电器总成
CN109243483A (zh) * 2018-10-17 2019-01-18 西安交通大学 一种含噪频域卷积盲源分离方法
CN109973251A (zh) * 2019-03-24 2019-07-05 西安航天动力测控技术研究所 一种用于固体火箭发动机振动信号降噪的盲源分离方法
CN110263876A (zh) * 2019-06-27 2019-09-20 南京信息工程大学 一种雨声信号降噪处理方法和系统
CN110807524A (zh) * 2019-11-13 2020-02-18 大连民族大学 单通道信号盲源分离幅度校正方法

Non-Patent Citations (7)

* Cited by examiner, † Cited by third party
Title
焦卫东;杨世锡;钱苏翔;严拱标;: "基于FFT-MCC分析的ICA(BSS)盲不确定性消除", 中国机械工程, no. 07, 15 April 2006 (2006-04-15) *
王卓 等: "用ICA方法消除无线通信与雷达信号中的多频混合噪声", 现代电子技术, no. 13, 12 March 2010 (2010-03-12) *
王彬 等: "基于FastICA的低信噪比雷达信号分选算法", 东北大学学报（自然科学版）, vol. 40, no. 11, 7 November 2019 (2019-11-07) *
王雁涛;吉磊;: "基于独立分量分析的混合信号盲分离技术", 通信技术, no. 05, 10 May 2018 (2018-05-10) *
贾银洁;许鹏飞;: "基于FastICA的混合音频信号盲分离", 信息与电子工程, no. 04, 25 August 2009 (2009-08-25) *
郁芸;张昕;胡杰;万程;周高信;王捷;殷悦楚楚;冯伟;王娟;: "快速独立分量分析算法解析及其仿真", 教育教学论坛, no. 28, 8 July 2020 (2020-07-08) *
闫姝;徐岩;: "负熵最大化条件下的语音信号盲源分离技术研究", 自动化仪表, no. 10, 20 October 2017 (2017-10-20) *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112869737A (zh) * 2021-02-01 2021-06-01 浙江大学山东工业技术研究院 一种非接触式人体血氧饱和度检测方法
CN113223548A (zh) * 2021-05-07 2021-08-06 北京小米移动软件有限公司 声源定位方法及装置
CN113948104A (zh) * 2021-12-07 2022-01-18 思必驰科技股份有限公司 盲语音信号分离方法、电子设备和存储介质

Also Published As

Publication number Publication date
CN112116922B (zh) 2024-04-12

Similar Documents

Publication Publication Date Title
US11996091B2 (en) Mixed speech recognition method and apparatus, and computer-readable storage medium
AU2021245165B2 (en) Method and device for processing quantum data
CN112116922A (zh) 一种噪声盲源信号分离方法、终端设备及存储介质
CN110164465B (zh) 一种基于深层循环神经网络的语音增强方法及装置
CN112509600A (zh) 模型的训练方法、装置、语音转换方法、设备及存储介质
CN109360581A (zh) 基于神经网络的语音增强方法、可读存储介质及终端设备
CN108732559B (zh) 一种定位方法、装置、电子设备及可读存储介质
US20110213587A1 (en) Method and computer program product for finding statistical bounds, corresponding parameter corners, and a probability density function of a performance target for a circuit
CN112183620A (zh) 基于图卷积神经网络的小样本分类模型的发育方法及系统
Chen et al. Testing for trends in high-dimensional time series
CN113646833A (zh) 语音对抗样本检测方法、装置、设备及计算机可读存储介质
KR20220048941A (ko) 인공 신경망을 이용한 트랜지스터 컴팩트 모델링 시스템, 방법 및 컴퓨터 프로그램 제품
CN105355198A (zh) 一种基于多重自适应的模型补偿语音识别方法
CN113095353B (zh) 一种基于ap聚类的欠定盲源分离方法
Guo et al. Underdetermined mixing matrix estimation algorithm based on single source points
CN109598175A (zh) 一种基于多小波基函数和超正交前向回归的时频分析方法
CN109657693B (zh) 一种基于相关熵和迁移学习的分类方法
CN110334716B (zh) 特征图处理方法、图像处理方法及装置
Yan et al. Computationally efficient direction of arrival estimation with unknown number of signals
Zhang et al. K-plane clustering algorithm for analysis dictionary learning
Wei et al. ATA: Attentional non-linear activation function approximation for VLSI-based neural networks
CN114661940B (zh) 一种适用于黑盒攻击下快速获取语音对抗样本的方法
Sahnoun et al. Sparse modal estimation of 2-D NMR signals
CN111508525B (zh) 一种全参考音频质量评价方法及装置
Liu et al. Autoregressive moving average graph filter design

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant