一种基于去噪核范数最小化的互质相干信源估计方法

1.本发明涉及互质阵列信号处理技术，尤其涉及一种基于去噪核范数最小化的互质相干信源估计方法。


背景技术：

2.随着空间电磁环境的日趋复杂，夺取制空权和制电磁权在现代信息化战争中的地位越发重要，因此军事行动的效果取决于拥有制电磁权的电子侦察设备的探测能力。然而，面对敌/我双方诸多电磁武器装备所释放的电磁辐射、民用电磁设备的辐射和自然界本身所产生的电磁波，电子侦察设备迫切需要解决严重的相干信源估计问题以提高其的多目标测向能力。
3.相较于之前电子侦察设备所配置的均匀线阵，互质阵具有系统化的阵列结构、增广的阵列孔径和高自由度这三大特点，使得在提高电子侦察设备的测向性能的同时增加了可探测目标数。基于上述互质阵的优点，已经有很多基于互质阵的相干信源估计算法和专利被提出，例如基于核范数最小化(nuclear norm minimization,nnm)的互质阵内插算法，基于原子范数最小化(atom norm minimization,anm)的互质阵内插算法和doa估计方法。然而，这些方法在低信噪比条件下测向能力急剧下降，并且在现代电子战下的信噪比大多数都是很低的，因此针对互质阵列迫切需要提高在低信噪比条件下相干信源的测向能力。


技术实现要素：

4.发明目的：本发明的目的是提供一种基于去噪核范数最小化的互质相干信源估计方法，该方法在基于核范数最小化的互质阵内插算法的基础上，引入去噪操作，通过凸优化求解由接收信号所构成的toeplitz矩阵的核范数最小化问题，实现对多个相干信源的高精度超分辨率doa(direction-of-arrival)估计。
5.技术方案：本发明的一种基于去噪核范数最小化的互质相干信源估计方法，包括以下步骤：
6.s1、假设现有一对稀疏均匀子阵所构成互质阵列，其中第一个子阵由间距为nd的m个阵元构成，第二个子阵由间距为md的n个阵元构成，d为半波长。因此该互质阵列位置集合标记为同时k个远场相干信源按不同的doa均向该互质阵发射相同的电磁波信号s1(t)，因此得到该互质阵列的接收信号为
7.s2、借助内插技巧和上述互质阵列构造内插均匀线阵(ula)，获取到内插ula的接收信号其中内插ula阵元位置集合标记为l为中间变量；
8.s3、根据协方差公式计算内插ula的采样协方差矩阵并取出第m行中所有元素用以定义具有toeplitz结构的矩阵其中，同时，根据进一步构造二值标识矩阵c(m)；
9.s4、由于借助插值技巧必然导致由构成的toeplitz结构矩阵缺失部分元素，即的部分元素值为零；为了将的零元素补全得到具有完整元素的toeplitz矩阵f(m)，首先将中m＝1的第一行和的最后一行取出，分别得到和然后利用矩阵补全和低秩理论分别重构出具有完整元素的toeplitz矩阵f(1)和其中f(1)是在基础上补全缺失元素后所得的矩阵，是在基础上补全缺失元素后所得的矩阵；
10.s5、分别对根据cvx工具包所求解得出的f(1)和进行特征值分解，并分别通过谱峰搜索得到第k个相干信源的角度参数和其中k＝1,
…
,k；
11.s6、根据f(1)和分别对应的噪声子空间的第k个相干信源的角度参数和最终得到第k个相干信源的估计doa参数其中k＝1,
…
,k。
12.进一步的，步骤s1中互质阵列位置集合互质阵列的接收信号为：
[0013][0014]
其中，为噪声，为阵列流型，和分别表示为第1个和第k个相干信源的导行矢量，θ1和θk分别为第1个和第k个相干信源的真实doa；第k个信源的导行矢量以及θk为第k个相干信源的真实doa，l2和分别为第2个和第个阵元所在位置，为该互质阵列的势，α∈ck为囊括了第k个元素为非零复数αk的列向量，其中k＝1,
…
,k,(
·
)h表示为共轭转置操作算子。
[0015]
进一步的，步骤s2中内插ula的接收信号的第l个元素也就是内插ula阵元位置集合中位于ld位置的阵元所对应的接收信号为：
[0016][0017]
其中，表示内插ula阵元位置集合位于ld位置的阵元所对应的接收信号，表示互质阵列阵元位置集合中位于ld位置的阵元所对应的接收信号，为在上的补集。
[0018]
进一步的，步骤s3具体为：
[0019]
根据公式(3)计算内插ula的采样协方差矩阵表达式为：
[0020][0021]
其中，(
·
)
*
表示为共轭操作算子，e{
·
}表示期望，表示的共轭，j为快
拍数；并取出第m行中所有元素用以定义具有toeplitz结构的矩阵表达式为：
[0022][0023]
其中，和分别表示中的第m行第1列、第m行第l-1列、第m行第l列、第m行第l+1列、和第m行第列元素，为内插ula的势，为虚拟均匀子阵个数；同时，根据进一步构造二值标识矩阵c(m)，其中二值标识矩阵c(m)的第p行第q列的元素[c(m)]
p,q
为：
[0024][0025]
其中，表示中的第p行第q列元素。
[0026]
进一步的，步骤s4包括以下步骤：
[0027]
s41、因为理论上由的零元素补全得到的具有完整元素的toeplitz矩阵f(m)的秩等于相干信源个数，所以在不减少孔径大小的情况下实现解相干的同时具有低秩性；令m＝1，f(1)为在基础上补全缺失元素后所得的矩阵；令是在基础上补全缺失元素后所得的矩阵；利用f(1)和的低秩性，分别考虑下面秩最小化优化问题：
[0028][0029]
和
[0030][0031]
其中，为直积算子，d(1)和表示去噪二值标识矩阵用于去噪操作，分别定义为：
[0032][0033]
和
[0034][0035]
其中，c(1)为借助第1行所构造二值标识矩阵，为借助第行所构造
二值标识矩阵，0
(l-1)
×
l
为(l-1)
×
l全零矩阵，01×
(l-1)
为1
×
(l-1)全零行向量；通过约束等式项和忽略中第一行第一列元素和和第行第列元素对重构的作用；
[0036]
s42、由于公式(6)和(7)所示优化问题是非凸的，因此将其转换为核范数最小化nnm优化问题：
[0037][0038]
和
[0039][0040]
其中，||
·
||
*
为核范数，同时借助半正定规划理论将公式(10)和(11)分别等价于如下半正定规划形式：
[0041][0042]
和
[0043][0044]
其中，w
11
、w
12
、和均为临时变量矩阵，trace(
·
)为求迹运算符，约束条件和表明矩阵满足半正定，(
·
)h表示为共轭转置操作算子；
[0045]
然后通过cvx凸优化工具包求解nnm问题。
[0046]
进一步的，步骤s5具体为：
[0047]
分别对根据cvx工具包所求解得出的f(1)和进行特征值分解(evd)：
[0048][0049]
和
[0050][0051]
其中u1和σ1分别为f(1)下的标准正交阵和对角矩阵，和分别为下的标准正交阵和对角矩阵；令u
1n
＝u1(:,k+1:l)和分别为f(1)和对应的噪声子空间，使得谱峰搜索步长变化，分别通过谱峰搜索得到f(1)对应的噪声子空间的第k个相干信源角度参数和对应的噪声子空间的第k个相干信源角度参数其中k
＝1,
…
,k：
[0052][0053]
和
[0054][0055]
其中，p1和分别为对应的空间谱函数。
[0056]
进一步的，步骤s6中对f(1)对应的噪声子空间的第k个相干信源角度参数和对应的噪声子空间的第k个相干信源角度参数进行平均操作，最终得到第k个相干信源的估计doa参数为：
[0057][0058]
有益效果：与现有技术相比，本发明具有以下优点：
[0059]
(1)与现有的基于核范数最小化(nuclear norm minimization,nnm)和基于原子范数最小化(atom norm minimization,anm)的互质相干信源估计方法相比，本发明通过引入去噪操作，保证低信噪比条件下的高精度相干信源估计性能，进而实现了高精度doa估计。
[0060]
(2)与传统的波束形成方法相比，本发明借助核范数最小化问题实现了相干信源超分辨doa估计。以图2(a)所示m＝2，n＝5的互质阵为例，设置信噪比(snr)等于-5db，快拍数为250，图3(a)的p1和图3(b)的空间谱函数图明显可以估计出角度范围差小于角分辨率等于0.886/(m+n-1)＝7.27
°
的两个相干信源(θ1＝10
°
，θ2＝17
°
)，而传统波束形成方法应用在等物理阵元个数的均匀线阵下其角度分辨率无法突破7.27
°
。
附图说明
[0061]
图1是本发明方法流程图；
[0062]
图2是本发明的不同阵列阵型的阵元位置排布图，其中(a)为m＝2和n＝5配置下的互质阵，(b)为m＝2和n＝5配置下的内插ula；
[0063]
图3是本发明针对m＝3和n＝5配置下的互质阵下所绘的超分辨doa估计谱峰搜索图，其中(a)为p1空间谱函数图二维图，(b)为空间谱函数图；
[0064]
图4是本发明针对m＝3和n＝5配置下的互质阵下所绘的空间谱函数图，其中(a)为p1空间谱函数图二维图，(b)为空间谱函数图；
[0065]
图5是本发明针对m＝3和n＝5配置下的互质阵下所绘的估计误差曲线图，其中(a)为均方根误差(rmse)随snr变化曲线图，(b)为rmse随快拍数变化曲线图。
具体实施方式
[0066]
下面结合附图和具体实施例对本发明进行详细说明。
[0067]
如图1所示，本发明的种基于去噪核范数最小化的互质相干信源估计方法，包括以下步骤：
[0068]
s1、假设现有一对稀疏均匀子阵所构成互质阵列，其中第一个子阵由间距为nd的m个阵元构成，第二个子阵由间距为md的n个阵元构成，d＝λ/2为半波长，λ为波长。因此该互质阵列位置集合标记为同时k个远场相干信源按不同的doa均向该互质阵列发射相同的电磁波信号s1(t)，因此可以得到该互质阵列的接收信号为：
[0069][0070]
其中，为噪声，为阵列流型，和分别表示为第1个和第k个相干信源的导行矢量，θ1和θk分别为第1个和第k个相干信源的真实doa；便于统一标识，此处给出第k(k＝1,
…
,k)个信源的导行矢量以及θk为第k个相干信源的真实doa，l2和分别为第2个和第个阵元所在位置，为该互质阵列的势，α∈ck为囊括了第k(k＝1,
…
,k)个元素为非零复数αk的列向量；(
·
)h表示为共轭转置操作算子。
[0071]
s2、借助内插技巧和互质阵列构造内插均匀线阵(uniform linear array,ula)，其中内插均匀线阵ula的阵元位置集合可标记为l为中间变量。基于此，内插ula的接收信号的第l个元素也就是内插ula阵元位置集合中位于ld位置的阵元所对应的接收信号为：
[0072][0073]
其中，表示内插ula阵元位置集合位于ld位置的阵元所对应的接收信号，表示互质阵列阵元位置集合中位于ld位置的阵元所对应的接收信号，为在上的补集。图2中(a)和(b)给出了m＝2和n＝5配置下的互质阵列和内插均匀线阵ula的阵元位置排布图，结合公式(2)可得：元位置排布图，结合公式(2)可得：
[0074]
s3、根据公式(3)计算内插ula的采样协方差矩阵表达式为：
[0075][0076]
其中，(
·
)
*
表示为共轭操作算子，e{
·
}表示期望，表示的共轭，j为快拍数；并取出第行中所有元素用以定义具有toeplitz结构的矩阵表达式为：
[0077][0078]
其中，和分别表示中的第m行第1列、第m行第l-1列、第m行第l列、第m行第l+1列、和第m行第列元素，为的势，为虚拟均匀子阵个数。同时，根据进一步构造二值标识矩阵c(m)，其中二值标识矩阵c(m)的第p行第q列的元素[c(m)]
p,q
为：
[0079][0080]
其中，表示中的第p行第q列元素。
[0081]
s4、由于借助插值技巧必然导致由构成的toeplitz结构矩阵缺失部分元素，即部分元素值为零，而缺失的元素和背景噪声的干扰是影响精确测角的两大因素。为了克服上述两大因素的影响，首先将中的第一行和最后一行取出，分别得到含有零值元素的和然后利用矩阵补全和低秩理论分别重构出具有完整元素的toeplitz矩阵f(1)和重构具体步骤如下：
[0082]
s41、因为理论上将中零元素补全后得到的具有完整元素的toeplitz矩阵f(m)的秩等于相干信源个数，所以可以在不减少孔径大小的情况下实现解相干的同时具有低秩性。令m＝1，f(1)定位为在基础上补全缺失元素后所得的矩阵。令是在基础上补全缺失元素后所得的矩阵。利用f(1)和具有低秩性这一特性，分别考虑下面秩最小化优化问题：
[0083][0084]
和
[0085][0086]
其中，为直积算子，d(1)和表示去噪二值标识矩阵用于去噪操作，分别定义为：
[0087][0088]
和
[0089][0090]
其中，c(1)为借助第1行所构造二值标识矩阵，为借助第行所构造二值标识矩阵，0
(l-1)
×
l
为(l-1)
×
l全零矩阵，01×
(l-1)
为1
×
(l-1)全零行向量。通过约束等式项和忽略中第一行第一列元素和和第行第列元素对重构的作用(和是带噪观测值)，可以提高对噪声的鲁棒性，进而缓解噪声对测角精度的影响。
[0091]
s42、由于公式(6)和(7)所示优化问题是非凸的，因此将其转换为基于核范数最小化(nnm)优化问题：
[0092][0093]
和
[0094][0095]
其中，||
·
||
*
为核范数，同时本发明借助半正定规划理论将(10)和(11)等价于如下半正定规划形式：
[0096][0097]
和
[0098][0099]
其中，w
11
、w
12
、和均为临时变量矩阵，trace(
·
)为求迹运算符，约束条件和表明矩阵满足半正定，(
·
)h表示为共轭转置操作算子。
[0100]
本发明通过cvx凸优化工具包(m.grant and s.boyd,“cvx:matlab software for disciplined convex programming,version 2.1,”available online at http://cvxr.com/cvx,mar.2014.)以求解nnm问题。
[0101]
s5、分别对根据cvx工具包所求解得出的f(1)和进行特征值分解(evd)：
[0102][0103]
和
[0104][0105]
其中u1和σ1分别为f(1)下的标准正交阵和对角矩阵，和分别为下的标准正交阵和对角矩阵。令u
1n
＝u1(:,k+1:l)和分别为f(1)和对应的噪声子空间，使得谱峰搜索步长变化，分别通过谱峰搜索得到f(1)对应的噪声子空间的第k个相干信源角度参数和对应的噪声子空间的第k个相干信源角度参数
[0106][0107]
和
[0108][0109]
其中，p1和分别为对应的空间谱函数。一般地，假设谱峰搜索步长为0.1
°
，即等间隔0.1
°
变化为了进一步说明所提方法在提高角分辨率方面的有效性，图3(a)的p1和图3(b)的空间谱函数图明显可以估计出角度范围差小于角分辨率等于0.886/(m+n-1)＝7.27
°
的两个相干信源(θ1＝10
°
，θ2＝17
°
)，而传统波束形成方法应用在等物理阵元个数的均匀线阵下其角度分辨率无法突破7.27
°
。图4(a)和图4(b)绘出了snr＝-12db，快拍数等于250，k＝3，三个目标的真实doa定义为θ1、θ2、θ3，其值分别为θ1＝10
°
，θ2＝30
°
，θ3＝60
°
配置下的p1和空间谱函数图，可以明显看出在低信噪比条件下本发明所提方法估计得到的doa离真实doa最近的，即误差最小。
[0110]
s6、对f(1)对应的噪声子空间的第k个相干信源角度参数和对应的噪声子空间的第k个相干信源角度参数进行平均操作，最终得到第k(k＝1,
…
,k)个相干信源的估计doa参数
[0111][0112]
本发明提出的一种去噪核范数最小化的互质相干信源估计方法用以保证低信噪比条件下的高精度相干信源估计性能。仿真结果表明本发明所提方法不仅可以解决了相干信源导致的解相干问题，而且在提高了在低信噪比条件下的相干信源测角精度的同时也实现了超分辨率波达角估计(direction-of-arrival,doa)估计：(1)同物理阵元个数下的互质阵能提供更大的阵列孔径，(2)去噪操作提高对算法对噪声的鲁棒性，(3)采用了核范数最小化实现了超分辨率。图5(a)和(b)分别给出了rmse随snr和快拍数变化曲线图，明显发现本发明所提方法可以实现低信噪比条件下的高精度doa估计，尤其是在-12db至0db。