一种基于矩阵补全的二阶统计量重构DOA估计方法与流程

本发明属于信号处理领域，具体涉及一种基于矩阵补全的二阶统计量重构doa估计方法。

背景技术：

在雷达、移动通信、无线传感器网络和射电天文学等领域，波达方向角估计(directionofarrival,doa)是一类重要的问题，通常被称为定向(directionfinding,df)估计或doa估计。基于信号和噪声子空间的传统doa估计算法如多重信号分类(multiplesignalclassification,music)算法可显著改善doa估计精度和分辨力。以music方法为代表的子空间类算法通常假设噪声是零均值、方差为1的复高斯白噪声。然而当附加在阵列传感器上的噪声为非均匀高斯噪声时，对接收信号协方差进行特征空间分解会引起信号子空间泄漏，诸如子空间类算法的性能将会急剧下降甚至失效。为进一步提高doa的估计精度和分辨力，基于压缩感知理论，liangg等人利用目标信号的空域稀疏性，提出一种稀疏重构l1-svd算法。在信源数量已知条件下，该算法将得到doa高精度估计。然而，在非均匀高斯噪声或未知信源先验信息条件下，该算法将无法分辨两个邻近角度，因而空间分辨力较差。

近年来，随着信号处理技术的快速发展，研究人员相继提出多种非均匀高斯噪声下的doa估计算法。pesaventom等人提出一种非均匀高斯噪声下最大似然(maximumlikelihood,ml)估计算法，该算法通过逐步迭代求解信号和噪声的对数似然函数以实现doa参数估计。然而，ml算法具有的较强初始值依赖性及较大运算复杂度限制该算法的广泛应用。liaob等人提出一种非均匀高斯噪声下基于矩阵补全的doa参数估计算法(matrixcompletionbasedmusic,mc-music)。该算法基于矩阵补全思想重构出无噪声信号协方差，而后采用传统类music算法实现doa参数估计，不仅降低了非均匀噪声的影响，也进一步避免doa估计算法的迭代求解。然而，该算法没有考虑协方差矩阵元素之间相关性，有可能导致算法求解数值不稳定，从而使得算法稳定性较差。

技术实现要素：

本发明所要解决的技术问题是解决传统doa估计算法在非均匀噪声下角度估计精度差和分辨率低的问题，提出了一种基于矩阵补全的二阶统计量重构doa估计方法。

本发明为解决其问题所采用的的技术方案为：提出一种基于矩阵补全的二阶统计量重构doa估计方法(mc-wlosrss)。首先建立阵列接收信号模型；其次基于矩阵补全方法，引入弹性正则化因子将接收信号协方差矩阵重构为无噪声协方差矩阵；而后在二阶统计量域下通过矩阵求和平均将无噪声协方差矩阵多矢量问题转化为单矢量问题；最后利用稀疏重构加权l1范数实现doa估计。该方法的具体步骤为：

步骤1：建立接收信号模型

假设有q个远场窄带信号入射至具有m个阵元的均匀线性阵列，则阵列接收信号模型可表示为：

其中，x(t)为m×1信号接收矢量，sq(t)和a(θq)＝[1e^-jα…e-^j(m-1)α]^t分别为第q个信号源的幅度和阵列导向矢量，α＝2πsin(θq)/λ表示阵元之间的第q个信号到达此阵元时的相移，{θ1,θ2,…,θq,…θq}为q个窄带信号的发射角度，d为阵元间距，λ为载波信号波长，通常d≤λ2，n(t)＝[n1(t),n2(t),…,nm(t)]^t为互不相关的0均值非均匀高斯噪声，即n(t)～cn(0,w)，w为噪声n(t)功率协方差矩阵，信号sq(t)互不相关。

为便于推导，式(1)接收信号模型可进一步改写为：

x(t)＝as(t)+n(t)(2)

其中，阵列导向流型矩阵且m＞＞q，即假设阵元数量远大于信号数量，信号波形矢量

在多快拍数条件下，式(2)接收信号模型可进一步表示为：

x^m×l＝a^m×qs^q×l+n^m×l(3)

其中，x＝[x(1),x(2),…,x(l)]、s＝[s(1),s(2),…,s(l)]、n＝[n(1),n(2),…,n(l)]分别为l个快拍下的接收信号矩阵、信号幅度矩阵、非均匀高斯噪声矩阵，l表示快拍数。

基于式(3)，上述接收信号的协方差可表示为：

其中，r为接收信号协方差矩阵，为信号功率协方差矩阵，pq为单信号功率，r0＝apa^h为不包含噪声的信号协方差矩阵，非均匀高斯噪声功率协方差矩阵且为2单非均匀高斯噪声功率，信号与噪声互不相关。

步骤2：基于矩阵补全的信号协方差重构算法

为避免在非均匀高斯噪声下，信号协方差矩阵特征空间分解会引起信号子空间泄漏，诸如子空间类算法的性能将会急剧下降甚至失效的情况，本发明引入矩阵补全重构算法以消除非均匀高斯噪声对信号协方差r的影响。

(1)矩阵补全重构算法

假设阵列信号模型中的阵元数远大于信号数(m＞＞q)，则无噪声信号协方差的秩rank(r0)＝q＜m，即无噪声信号协方差矩阵r0是一个低秩矩阵，且协方差矩阵r0非主对角线元素值等价于接收信号对应元素值，故可进一步通过矩阵补全理论对r0进行主对角线元素重构，以此消除非均匀高斯噪声对信号协方差r的影响。

对给定矩阵其在相应的子集ω内投影，即采样可表示为：

其中子集表示采样元素索引集合，m,n分别为采样最大索引行数和列数，xij,[pω(x)]ij分别为矩阵x和pω(x)第ij个索引元素，pω(·)为正交投影算子。

鉴于子集ω是随机抽取的，使得采样后的信号空域数据矩阵满足矩阵低秩特性，故而可通过矩阵补全方法重构出其中的未知元素，即重构模型可表示为：

其中为待重构变量，x为其对应的已知变量。

由于秩函数具有非凸性，故将其等价松弛为凸核范数则式(6)重构模型可等价转化为凸优化模型，即：

其中核范数等价于矩阵所有奇异值之和。

(2)引入弹性正则化因子

为避免利用矩阵补全重构算法对强相关性数据进行求解时出现解不稳定的情况，本发明引入弹性正则化项增加矩阵补全重构算法求解稳定性，则式(7)凸优化模型进一步表示为：

其中，τ为和之间的平衡度正则化因子。

基于式(4)和式(8)，矩阵补全凸优化求解模型可改写为：

由无噪声协方差r0是一个半正定矩阵，可得：

τ||r0||*＝τtr(r0)＝tr(τr0)(10)

其中tr(·)表示矩阵的迹，即矩阵对角线元素之和。

又则式(9)最小化项可重写为：

同时，式(9)约束条件pω(r0)＝pω(r)可等价表示为：

jvec(r0-r)＝0(12)

其中j代表一个m(m-1)×m²的选择矩阵。

(3)正则化矩阵补全重构算法

基于式(11)和式(12)，式(9)凸优化模型可进一步改写为：

式(13)凸优化模型是一个线性规划问题，约束条件jvec(r0-r)＝0在实际工程中要求较为苛刻，故可将其等价松弛化为：

||jvec(r0-r)||2≤ξ(14)

其中，ξ表示为与接收信号协方差r相关的一误差常数。

则式(14)凸优化模型可等价表示为：

为便于后续矩阵补全重构算法求解，式(15)凸优化求解问题可进一步表示为：

其中，t为辅助优化变量。

优化问题(16)中约束条件可进一步改写为：

根据schur补定理及tr(x^hx)＝[vec(x)]^hvec(x)可知，式(17)凸优化求解问题可转化为半定规划问题(semidefiniteprogramming,sdp)，即：

上述sdp问题可运用cvx工具包实现高效求解。

步骤3：二阶统计量域下的加权l1稀疏重构算法

本发明基于矩阵补全理论(mc理论)重构无噪声协方差r0，以此消除非均匀高斯噪声影响，而后可采用传统子空间类算法实现doa参数估计。然而，基于无噪声协方差r0的传统类doa估计算法没有考虑到协方差矩阵r0元素之间的的相关性，在实现doa参数估计的同时，会导致算法计算量较大、低信噪比条件下估计稳定性差和分辨率低等问题。

为进一步提高低信噪比条件下来波信号角度估计精度和分辨力，降低稀疏重构算法计算复杂度，本发明基于mc理论，提出一种二阶统计量域下的加权l1范数稀疏重构算法。

(1)协方差矩阵求和平均

无噪声协方差r0中的元素，即第i行和第j列的阵列输出互相关系数可表示为：

对信号协方差r0按行展开，可表示为：

其中，当i-j＝u-v时，互相关系数rij＝ruv，i,j,u,v∈[1,m]，即它们含有相同的系数特征，可通过二阶统计量域下的求和平均将协方差多矢量问题进一步转化为单矢量问题。

定义(2m-1)×1的向量υ，通过对信号协方差r0求和平均，其第个元素可表示为：

(2)接收信号单矢量稀疏表示

υ的向量形式可表示为：

υ＝b(θ)p(22)

其中，p＝[p1,p2,…,pq]^t，为虚拟阵列流型矩阵，其第q列可表示为：

上述单矢量问题可转变为稀疏信号重构问题，将空间波达方向角θ分割为n(n＞＞m)个网格，即则稀疏表示下的接收信号矢量可表示为：

其中，为过完备基矩阵，为k稀疏的信号功率向量。

(3)基于加权l1范数的doa优化问题求解

上述稀疏重构问题可通过l1范数约束最优化算法进行求解，即：

其中，η为估计误差参数因子，||·||f和||·||1表示2范数和1范数。

由式(25)可知，l1范数约束最优化算法对目标信号大系数的约束比对小系数的约束更为严格，导致其事实上是一种有偏的估计。

基于此，本发明采用加权l1范数稀疏重构算法实现非均匀噪声下的信源doa估计，即：

式(26)还可进一步改写为：

式(27)凸优化问题可通过二阶锥规划(second-orderconeprogramming,socp)重构实现doa参数估计。

基于上述讨论，本专利提出的mc-wlosrss算法可表述如下：

ⅰ)初始化：

(1).υ＝avr(r0)，maxiter，ε，η；

(2).

ⅱ)迭代：

(3).权值

(4).解如下优化问题以更新系数；

ⅲ)终止：

(5).重复步骤3、4直至满足如下条件或者m≥maxiter。

其中，avr代表式(21)矩阵求和平均运算；maxiter表示算法最大迭代次数；ε为算法迭代终止参数；为权值w平衡约束因子，防止算法迭代求解时权值分母趋于无穷大；η为估计误差参数因子；m为算法进行第m次迭代；为算法迭代初始化值；为算法进行第m次迭代获得的角度谱峰估计值。

本发明基于矩阵补全理论，引入弹性正则化因子重构无噪声信号协方差矩阵，从而避免非均匀高斯噪声的影响，且提高了mc求解的数值稳定性；此外所提算法在二阶统计量域下通过矩阵求和平均将稀疏多矢量问题转化为单矢量问题以改善低信噪比条件下doa估计性能，一定程度上降低了算法计算复杂度，且在低信噪比条件下具有较高的角度估计精度和分辨力。

附图说明

以下结合附图和具体实施方式对本发明作进一步详细说明。

图1为本发明实现的流程图。

图2为在信噪比snr＝0db和-5db条件下的非相干信号空域谱对比图。

图3为在邻近非相干信号空域谱对比图。

图4为在信源数变化对music、mc-wlosrss算法的影响。

图5为在非相干信号时doa估计均方根误差随信噪比的变化曲线图。

图6为doa估计均方根误差随快拍数的变化曲线图。

具体实施方式

为了使本发明的目的、技术方案和优点更加清楚，下面结合附图和具体的实施方式对本发明进行详细描述。

本发明提出一种基于矩阵补全的二阶统计量重构doa估计方法，该方法的具体步骤为：

步骤1：建立接收信号模型

假设有q个远场窄带信号入射至具有m个阵元的均匀线性阵列，则阵列接收信号模型可表示为：

其中，x(t)为m×1信号接收矢量，sq(t)和a(θq)＝[1e^-jα…e^-j(m-1)α]^t分别为第q个信号源的幅度和阵列导向矢量，α＝2πsin(θq)/λ表示阵元之间的第q个信号到达此阵元时的相移，{θ1,θ2,…,θq,…θq}为q个窄带信号的发射角度，d为阵元间距，λ为载波信号波长，通常d≤λ2，n(t)＝[n1(t),n2(t),…,nm(t)]^t为互不相关的0均值非均匀高斯噪声，即n(t)～cn(0,w)，w为噪声n(t)功率协方差矩阵，信号sq(t)互不相关。

为便于推导，式(1)接收信号模型可进一步重写为：

x(t)＝as(t)+n(t)(29)

其中，阵列导向流型矩阵且m＞＞q，即假设阵元数量远大于信号数量，信号波形矢量

在多快拍数条件下，式(2)接收信号模型可进一步表示为：

x^m×l＝a^m×qs^q×l+n^m×l(30)

其中，x＝[x(1),x(2),…,x(l)]、s＝[s(1),s(2),…,s(l)]、n＝[n(1),n(2),…,n(l)]分别为l个快拍下的接收信号矩阵、信号幅度矩阵、非均匀高斯噪声矩阵，l表示快拍数。

基于式(3)，上述接收信号的协方差可表示为：

其中，r为接收信号协方差矩阵，为信号功率协方差矩阵，pq为单信号功率，r0＝apa^h为不包含噪声的信号协方差矩阵，非均匀高斯噪声功率协方差矩阵且为单非均匀高斯噪声功率，信号与噪声互不相关。

步骤2：基于矩阵补全的信号协方差重构算法

为避免在非均匀高斯噪声下，信号协方差矩阵特征空间分解会引起信号子空间泄漏，诸如子空间类算法的性能将会急剧下降甚至失效的情况，本发明引入矩阵补全重构算法以消除非均匀高斯噪声对信号协方差r的影响。

(1)矩阵补全重构算法

假设阵列信号模型中的阵元数远大于信号数(m＞＞q)，则无噪声信号协方差的秩rank(r0)＝q＜m，即无噪声信号协方差矩阵r0是一个低秩矩阵，且协方差矩阵r0非主对角线元素值等价于接收信号对应元素值，故可进一步通过矩阵补全理论对r0进行主对角线元素重构，以此消除非均匀高斯噪声对信号协方差r的影响。

对给定矩阵其在相应的子集ω内投影，即采样可表示为：

其中子集表示采样元素索引集合，m,n分别为采样最大索引行数和列数，xij,[pω(x)]ij分别为矩阵x和pω(x)第ij个索引元素，pω(·)为正交投影算子。

鉴于子集ω是随机抽取的，使得采样后的信号空域数据矩阵满足矩阵低秩特性，故而可通过矩阵补全方法重构出其中的未知元素，即重构模型可表示为：

其中为待重构变量，x为其对应的已知变量。

由于秩函数具有非凸性，故将其等价松弛为凸核范数则式(6)重构模型可等价转化为凸优化模型，即：

其中核范数等价于矩阵所有奇异值之和。

(2)引入弹性正则化因子

为避免利用矩阵补全重构算法对强相关性数据进行求解时出现解不稳定的情况，本发明引入弹性正则化项增加矩阵补全重构算法求解稳定性，则式(7)凸优化模型进一步表示为：

其中，τ为和之间的平衡度正则化因子。

基于式(4)和式(8)，矩阵补全凸优化求解模型可改写为：

由无噪声协方差r0是一个半正定矩阵，可得：

τ||r0||*＝τtr(r0)＝tr(τr0)(37)

其中tr(·)表示矩阵的迹，即矩阵对角线元素之和。

又则式(9)最小化项可重写为：

同时，式(9)约束条件pω(r0)＝pω(r)可等价表示为：

jvec(r0-r)＝0(39)

其中j代表一个m(m-1)×m²的选择矩阵。

(3)正则化矩阵补全重构算法

基于式(11)和式(12)，式(9)凸优化模型可进一步改写为：

式(13)凸优化模型是一个线性规划问题，约束条件jvec(r0-r)＝0在实际工程中要求较为苛刻，故可将其等价松弛化为：

||jvec(r0-r)||2≤ξ(41)

其中，ξ表示为与接收信号协方差r相关的一误差常数。

则式(14)凸优化模型可等价表示为：

为便于后续矩阵补全重构算法求解，式(15)凸优化求解问题可进一步表示为：

其中，t为辅助优化变量。

优化问题(16)中约束条件可进一步改写为：

根据schur补定理及tr(x^hx)＝[vec(x)]^hvec(x)可知，式(17)凸优化求解问题可转化为半定规划问题(semidefiniteprogramming,sdp)，即：

上述sdp问题可运用cvx工具包实现高效求解。

步骤3：二阶统计量域下的加权l1稀疏重构算法

本发明基于mc理论重构无噪声协方差r0，以此消除非均匀高斯噪声影响，而后可采用传统子空间类算法实现doa参数估计。然而，基于无噪声协方差r0的传统类doa估计算法没有考虑到协方差矩阵r0元素之间的的相关性，在实现doa参数估计的同时，会导致算法计算量较大、低信噪比条件下估计稳定性差和分辨率低等问题。

为进一步提高低信噪比条件下来波信号角度估计精度和分辨力，降低稀疏重构算法计算复杂度，本发明基于mc理论，提出一种二阶统计量域下的加权l1范数稀疏重构算法。

(1)协方差矩阵求和平均

无噪声协方差r0中的元素，即第i行和第j列的阵列输出互相关系数可表示为：

对信号协方差r0按行展开，可表示为：

其中，当i-j＝u-v时，互相关系数rij＝ruv，i,j,u,v∈[1,m]，即它们含有相同的系数特征，可通过二阶统计量域下的求和平均将协方差多矢量问题进一步转化为单矢量问题。

定义(2m-1)×1的向量υ，通过对信号协方差r0求和平均，其第个元素可表示为：

(2)接收信号单矢量稀疏表示

υ的向量形式可表示为：

υ＝b(θ)p(49)

其中，p＝[p1,p2,…,pq]^t，为虚拟阵列流型矩阵，其第q列可表示为：

上述单矢量问题可转变为稀疏信号重构问题，将空间波达方向角θ分割为n(n＞＞m)个网格，即则稀疏表示下的接收信号矢量可表示为：

其中，为过完备基矩阵，为k稀疏的信号功率向量。

(3)基于加权l1范数的doa优化问题求解

上述稀疏重构问题可通过l1范数约束最优化算法进行求解，即：

其中，η为估计误差参数因子，||·||f和||·||1表示2范数和1范数。

由式(25)可知，l1范数约束最优化算法对目标信号大系数的约束比对小系数的约束更为严格，导致其事实上是一种有偏的估计。

基于此，本发明采用加权l1范数稀疏重构算法实现非均匀噪声下的信源doa估计，即：

式(26)还可进一步改写为：

式(27)凸优化问题可通过二阶锥规划(second-orderconeprogramming,socp)重构实现doa参数估计。

基于上述讨论，本专利提出的mc-wlosrss算法可表述如下：

ⅰ)初始化：

(1).υ＝avr(r0)，maxiter，ε，η；

(2).

ⅱ)迭代：

(3).权值

(4).求解如下优化问题以更新系数；

ⅲ)终止：

(5).重复步骤3、4直至满足如下条件或者m≥maxiter。

其中，avr代表式(21)矩阵求和平均运算；maxiter表示算法最大迭代次数；ε为算法迭代终止参数；为权值w平衡约束因子，防止算法迭代求解时权值分母趋于无穷大；η为估计误差参数因子；m为算法进行第m次迭代；为算法迭代初始化值；为算法进行第m次迭代获得的角度谱峰估计值。

本发明的效果可通过仿真进一步说明：

仿真条件：均匀线性阵列空间阵元个数m＝12，信号快拍数l＝500，非相干信号功率为误差常数ξ＝5，τ＝5m，ε＝0.01，网格数量n＝181。信噪比定义为其中为单噪声功率。均方根误差定义为：

其中，k为蒙特卡洛实验重复次数。

非均匀高斯噪声功率协方差定义为：

w＝diag{2.0,10,2.5,5.0,0.5,1.5,3.0,5.0}(56)

其中，diag{·}表示对角化算子。

仿真内容：

仿真1：考虑信噪比snr＝0db和-5db两种情况。图2为在信噪比snr＝0db和-5db条件下的非相干信号空域谱对比图。

从图2(a)可看出，传统music、wl1和l1-svd算法在非均匀噪声条件下不能有效分辨10°和16°两个目标，而mc-music算法和本发明所提mc-wlosrss算法可以有效分辨10°和16°两个目标。图2(b)为三个入射角度分别为-3°、10°和16°、snr＝-5db的信号空域谱输出对比图。从图2(b)可看出，在低信噪比情况下，只有本发明所提mc-wlosrss算法可以有效分辨三个角度，mc-music等其它传统类算法均不能有效分辨10°和16°两个目标。图2表明，与mc-music等其它传统类算法相比，所提算法在非均匀高斯噪声和低信噪比条件下仍具有较好的doa估计性能。另外，从图2还可以看出，本发明所提mc-wlosrss算法具有更窄的主瓣谱峰和更低的旁瓣，进一步验证本发明所提mc-wlosrss算法在非均匀高斯噪声和低信噪比条件下具有较高的角度估计精度和分辨力。

仿真2：考虑三个入射角度分别为-3°、10°和13°的非相干信号，信噪比snr＝5db。图3为在邻近非相干信号空域谱对比图。

由图3可以看出，由于非均匀高斯噪声的影响，传统music算法和wl1算法不能有效分辨10°和13°两个邻近角度，mc-music算法虽然重构出无噪声协方差，但是受限于算法角度分辨率的影响，不能对邻近角度实现有效估计。本发明所提mc-wlosrss算法和理想高斯白噪声下的l1-svd(ideal-l1-svd)算法可以对邻近角度实现有效分辨。然而，需要注意的是，ideal-l1-svd算法虽然正确估计出两个邻近角度，但是相对于本发明所提mc-wlosrss算法而言，ideal-l1-svd算法在两个邻近角度之间的谱峰估计精度较差，而本发明所提mc-wlosrss算法在两个邻近角度之间的谱峰估计精度较高，且具有更低的旁瓣，表明本发明所提mc-wlosrss算法对邻近角度具有较高角度分辨力。

仿真3：考虑三个入射角度分别为-3°、10°和16°的非相干信号，信噪比snr＝5db，信源数k＝1，2，3时，music、mc-wlosrss算法的doa估计有效性。图4为在信源数变化对music、mc-wlosrss算法的影响。

图4(a)为信源数变化对music算法性能影响。从4(a)可以看出，当信源数k分别为1，2时，传统music算法不能有效分辨三个角度。只有当信源数等于待估计角度数量，即k＝3时，music才能有效分辨三个角度。图4(b)为信源数变化对mc-wlosrss算法性能影响。从4(b)可以看出，当信源数k分别为1，2，3时，本发明所提mc-wlosrss算法均能有效分辨三个角度，信源数的变化对本发明所提mc-wlosrss算法几乎没有影响。从图4可以看出，传统子空间类doa估计算法由于受制于信号或信号子空间数目，即信源数目的先验信息，导致子空间类算法的doa估计性能对信源数的变化较为敏感；本发明所提mc-wlosrss算法利用目标空域稀疏性，采用稀疏重构对待估计信号进行稀疏逼近，而无需已知信源先验信息，进一步体现所提mc-wlosrss算法优越性。

仿真4：考虑两个入射角度为-3°和5°的非相干信号，快拍数l＝500，信噪比snr＝[-8:2:12]，进行200次蒙特卡洛独立重复实验。图5为在非相干信号时doa估计均方根误差随信噪比的变化曲线图。

从图5可以看出，传统music算法在低信噪比条件下的doa估计rmse相对较高。mc-music算法采用mc理论重构出无噪声协方差，以消除分均匀噪声的影响，故而mc-music算法在低信噪比条件下doa估计rmse低于传统music算法。在低信噪比条件下，wl1算法和l1-svd算法有一个相对较低的rmse。然而，从图5还可以看出，在给定仿真条件下，本发明所提mc-wlosrss算法估计性能明显优于wl1算法和l1-svd算法，具有较低rmse，特别是在低信噪比条件下，这种优势更加明显，表明所提mc-wlosrss算法具有较好的doa估计性能。

仿真5：考虑两个入射角度为-3°和5°的非相干信号，信噪比snr＝0db，快拍数l＝[100:1200]，进行200次蒙特卡洛独立重复实验。图6为doa估计均方根误差随快拍数的变化曲线图。

从图6可以看出，随着快拍数的增加，本发明所提mc-wlosrss算法及music、mc-music、wl1和l1-svd算法的rmse逐渐降低。然而，需要注意的是，所提mc-wlosrss算法rmse均低于music、mc-music、wl1和l1-svd算法，表明所提mc-wlosrss算法在非均匀高斯噪声下具有更优的doa参数估计性能。

仿真结果表明，与传统的music、mc-music、wl1和l1-svd算法相比，所提算法在非均匀高斯噪声和低信噪比条件下，具有较好的doa参数估计性能。由此，本发明所提算法可以为工程应用中阵列信号处理领域的doa估计性能研究提供坚实的理论与实现依据。

发明基于矩阵补全理论，引入弹性正则化因子重构无噪声信号协方差矩阵，从而避免非均匀高斯噪声的影响，且提高了mc求解的数值稳定性；此外所提算法在二阶统计量域下通过矩阵求和平均将稀疏多矢量问题转化为单矢量问题以改善低信噪比条件下doa估计性能，一定程度上降低了算法计算复杂度，且在低信噪比条件下具有较高的角度估计精度和分辨力。