一种基于阵列信号的多目标定位识别方法
阅读说明：本技术 一种基于阵列信号的多目标定位识别方法 (Multi-target positioning identification method based on array signals ) 是由 理华 孟晓辉 王耀辉 于 2018-09-11 设计创作，主要内容包括：本发明公开了一种基于阵列信号的多目标定位识别方法,所述方法包括：一种基于阵列信号的多目标定位识别方法,所述方法包括：步骤1)利用阵列盲信号处理将N个独立信源从混叠信号中分离出来；步骤2)对多通道的阵列信号进行频带分解,分离出M个不同频带的信号；将空间平面分为P*Q个网格,P为网格的行总数,Q为网格的列总数；步骤3)利用阵列信号处理的声源定位算法分别计算M个频带在每个网格位置上输出的功率；步骤4)基于步骤3)的每个网格的M个频带的功率,将每一个网格的信号恢复到时域,获得每个网格上的时域信号；步骤5)将步骤1)的每个独立信源分别与P*Q个时域信号进行匹配,匹配到的网格位置为独立信源的位置,共匹配P*Q*N次。(The invention discloses a multi-target positioning identification method based on array signals, which comprises the following steps: a multi-target positioning identification method based on array signals comprises the following steps: step 1) separating N independent information sources from aliasing signals by array blind signal processing; step 2) carrying out frequency band decomposition on the multi-channel array signal to separate out M signals with different frequency bands; dividing the space plane into P × Q grids, wherein P is the total number of rows of the grids, and Q is the total number of columns of the grids; step 3) respectively calculating the output power of the M frequency bands at each grid position by using a sound source positioning algorithm of array signal processing; step 4) recovering the signal of each grid to the time domain based on the power of the M frequency bands of each grid in the step 3), and obtaining a time domain signal on each grid; and 5) respectively matching each independent signal source in the step 1) with P x Q time domain signals, wherein the matched grid position is the position of the independent signal source and is matched for P x Q x N times.)
技术领域
本发明涉及语音信号处理领域，具体涉及一种基于阵列信号的多目标定位识别方法。
背景技术
传声器阵列信号处理技术是源于雷达、声呐等系统。传声器阵列声源定位在通信、移动机器人和助听器装置等许多领域具有广泛应用价值。在这些应用中，最基本的就是要估计声源位置，并且需要兼顾定位精度和算法的实时性。传统的DOA估计方法，比如可控功率响应的声源定位方法(SRP-PHAT)，只能估计声源的位置，但是对于频率相似的声源不能分辨其位置，即无法建立声源和位置的一一对应的关系。
发明内容
本发明的目的在于克服上述技术缺陷，提出了一种基于阵列信号的多目标定位识别方法，该方法可以将频率相似的独立信源或者稀疏信源与其声源位置建立一一对应的关系，不但实现了声源的位置定位，也实现了声源的种类定位。
为了实现上述目的，本发明提出了一种基于阵列信号的多目标定位识别方法，所述方法包括：
步骤1)利用阵列盲信号处理将N个独立信源从混叠信号中分离出来；
步骤2)对多通道的阵列信号进行频带分解，分离出M个不同频带的信号；将空间平面分为P*Q个网格，P为网格的行总数，Q为网格的列总数；
步骤3)利用阵列信号处理的声源定位算法分别计算M个频带在每个网格位置上输出的功率；
步骤4)基于步骤3)的每个网格的M个频带的功率，将每一个网格的信号恢复到时域，获得每个网格上的时域信号；
步骤5)将步骤1)的每个独立信源分别与P*Q个时域信号进行匹配，匹配到的网格位置为独立信源的位置，共匹配P*Q*N次。
作为上述方法的一种改进，所述步骤3)的具体实现过程为：
利用SRP-PHAT方法计算每个网格位置上输出的第m，1≤m≤M个频带的功率Pm(s)：
其中，L为的阵列信号的通道个数，Xk(ω)为第k个通道信号Xk(t)的加窗傅里叶变换，τk为第k个通道指向网格(p，q)处的可控时延；
M个频带在P*Q个网格上输出P*Q*M个功率。
作为上述方法的一种改进，所述步骤4)的具体实现过程为：网格(p，q)处的时域信号x(m，s)为：
其中，m＝0，1，2，...M-1，m表示时域信号的时间点，
作为上述方法的一种改进，所述步骤5)的具体实现过程为：
对于步骤1)得到的每个独立信源X和步骤4)得到的网格(p，q)处的时域信号Y(p，q)＝x(m，s)，使用互相关法计算相关系数：
其中，Cov(X，Y(p，q))为X和Y(p，q)的协方差，Var[X]为X的方差，Var[Y(p，q)]为Y(p，q)的方差；
则相关系数最大的网格为：
本发明的优势在于：
1、本发明首先利用盲源分离的方法，将混叠在一起的独立声源或者稀疏声源还原出来，可以得到多个声源的时域和频域特征，能够实现声源的种类区分；
2、本发明的方法利用多频段的SRP算法将声源按照不同频率进行声源定位，利用不同声源的特征进行匹配最终建立声源和位置的一一对应关系。
附图说明
图1为本发明的不同信源的频谱特征的计算示意图；
图2为本发明的不同位置的频谱特征的计算示意图；
图3为本发明的利用匹配算法确定不同信源位置的示意图。
具体实施方式
下面结合附图和具体实施例对本发明进行详细的说明。
本发明提供了一种基于阵列信号的多目标定位识别方法，所述方法包括：
步骤1)利用阵列盲信号处理将N个独立信源s1(t)，…sN(t)从混叠信号中分离出来；
采用自然梯度法进行迭代计算，算法做出了如下的假定：具有N个独立信号源的信号s1(t)，…sN(t)和L个独立通道的观察量x1(t)，…xL(t)(已知)，观察量和信号源具有如下的关系：
x(t)＝As(t)
其中x(t)＝[x1(t)，…xL(t)]T，s(t)＝[s1(t)，…sN(t)]T，A是一个N×L的系数矩阵，原问题变成了已知x(t)和s(t)的独立性，求对s(t)的估计问题。假定有如下公式：
y(t)＝Wx(t)
其中y(t)是对s(t)的估计，W是一个L×N系数矩阵，问题变成了如何有效的对矩阵W做出估计。估计的迭代过程为：
1)初始化：W(0)为单位矩阵；
2)循环执行如下的步骤，直到W(n+1)与W(n)差异小于规定值ε(计算矩阵差异的方法可以设置)，也可以规定迭代次数。
3)利用公式
y(n)＝W(n)y(n-1)，其中，y(-1)＝x
4)利用计算公式如下：
W(n+1)＝W(n)+η(n)[I-φ(y(n))yT(n)]W(n)
其中，W(n)为需要估计的矩阵，η(n)为步长，φ(y)是一个非线性变换，比如φ(y)＝φ(y3)，实际计算时y为一个L×k矩阵，k为采样点个数。
5)W(n)收敛后为估计的W
利用y(t)＝Wx(t)得到N个独立信源的信号s1(t)，…sN(t)。
步骤2)将多通道的阵列信号使用快速傅里叶变换(FFT)进行频带分解，分离出不同频带的信号；假设共有M个频带：分别为M1，M2，…MM；将空间平面分为P*Q个网格，P为网格的行总数，Q为网格的列总数；
步骤3)利用阵列信号处理的PHAT算法分别计算M个频带在每个网格位置上输出的功率；
利用SRP-PHAT方法计算每个网格位置上输出的第m，1≤m≤M个频带的功率Pm(s)：
其中，L为的阵列信号的通道个数，Xk(ω)为第k个通道信号Xk(t)的加窗傅里叶变换，τk为第k个通道指向网格(p，q)处的可控时延；
M个频带在P*Q个网格上输出P*Q*M个功率。
经过SRP-PHAT算法后，可以得到每个网格的输出功率的大小；因此共有P*Q组数据，M个频带会有P*Q*M个数据；
步骤4)使用逆快速傅里叶变换(ifft)将每一个网格的信号进行恢复，恢复到时域；共有P*Q组波形；如图2所示；
在某一网格中，有M个频段上有M个数据，利用这些数据，使用反傅里叶变换，将每个网格中的数据恢复到时域，得到P*Q个时域信号。
网格(p，q)处的时域信x(m，s)为：
其中，m＝0，1，2，...M-1，m表示时域信号的时间点，
步骤5)对于步骤1)的每个独立信源X，与P*Q个时域信号Y(p，q)进行匹配：使用互相关法计算相关系数：
其中，Y(p，q)＝x(m，s)，Cov(X，Y)为X，Y的协方差，Var[X]为X的方差，Var[Y(p，q)]为Y(p，q)的方差；
则相关系数最高的网格为：
针对N个独立信源上述过程需要进行N次，则可以确定N个独立信源在为P*Q网格上的位置；如图3所示。
最后所应说明的是，以上实施例仅用以说明本发明的技术方案而非限制。尽管参照实施例对本发明进行了详细说明，本领域的普通技术人员应当理解，对本发明的技术方案进行修改或者等同替换，都不脱离本发明技术方案的精神和范围，其均应涵盖在本发明的权利要求范围当中。
- 上一篇：一种医用注射器针头装配设备
- 下一篇：车载雷达控制方法、装置、设备和存储介质