基于PD-ALM算法的稀疏阵列DOA估计方法技术

本发明专利技术公开了一种基于PD‑ALM算法的稀疏阵列DOA估计方法，包括以下步骤：将稀疏阵列每次快拍获得的采样数据重排成托普利兹矩阵x

【技术实现步骤摘要】
基于PD-ALM算法的稀疏阵列DOA估计方法
本专利技术属于雷达信号处理技术，具体涉及一种基于PD-ALM算法的稀疏阵列DOA估计方法。
技术介绍
近年来，随着压缩感知和稀疏重构理论在雷达信号处理领域中的不断发展，基于压缩感知和稀疏重构的DOA估计方法得到了广泛的研究。矩阵填充理论由压缩感知理论衍生出来，它能够对缺失数据的矩阵进行填充从而得到完整的矩阵。矩阵填充主要研究当矩阵中仅观测到部分数据或者矩阵中存在部分数据缺失时，利用已知矩阵元素的相关性，来对缺失的数据进行填充。矩阵填充在数学形式上可以描述为一个仿射秩最小化问题，由于秩函数具有非光滑性和非凸性的特点，所以这类问题通常是NP难问题。一些研究人员利用核范数最小化来替代秩最小化问题，得到了很多有效的算法，包括内点法、奇异值阈值法(SVT)、近似奇异值分解的不动点延拓法、类直线加速策略(APGL)的近似梯度加速法、低秩矩阵拟合方法等。但这些算法都必须满足严格的恢复条件，包括限制等距性质(RIP)，零空间性质(NSP)，s-good性质等。所以，研究人员也开始将目光投向于直接求解秩最小化问题的有效算法。专利CN201910464969.7公开了一种基于矩阵填充的嵌套阵阵元失效下的波达方向估计方法及装置，该方法保留了嵌套阵自身的优势，并利用矩阵填充算法填充了更多的阵元来进行DOA估计。该方法虽有效提高了稀疏阵列的DOA估计精度，但只适用于嵌套阵，并不适用于任意的稀疏阵。专利CN201810110343.1公开了一种矩阵补全方法，该方法首先利用逼近函数近似计算矩阵的秩，然后建立一种逼近矩阵秩的低秩重建模型，接着提出基于非凸函数的低秩矩阵重建模型的求解算法。该方法重建的矩阵精度高，易于操作，可以从少量数据中恢复出完整信号。但该方法适用于秩较小的大型矩阵，当应用在小型稀疏天线阵列中并不能得到预期的效果。
技术实现思路
本专利技术的目的在于提供一种基于PD-ALM算法的稀疏阵列DOA估计方法。实现本专利技术目的的技术解决方案为：一种基于PD-ALM算法的稀疏阵列DOA估计方法，包括以下步骤：步骤1、设置稀疏阵列的快拍总数为imax，i＝1；步骤2、将稀疏阵列第i次快拍获得的采样数据x(t)重排构造成托普利兹矩阵xT，该矩阵即为低秩稀疏矩阵；步骤3、运用PD-ALM算法对低秩稀疏矩阵xT进行基于秩最小化的优化矩阵填充，得到满阵x'T；步骤4、获取矩阵x'T中的第一行数据，该数据即为补全后的第i次快拍数据，然后将该数据作为稀疏阵列接收数据矩阵X'的第i列数据；步骤5、i＝i+1，重复步骤2～步骤4直到i＝imax，这时所有快拍采样数据补全完毕，得到稀疏阵列接收数据矩阵X'；步骤6、运用DOA估计算法对数据矩阵X'进行DOA估计。本专利技术中基于秩最小化的矩阵填充算法与现有的用核范数逼近秩项的方法不同，采用罚分解法直接求解秩最小化问题，其显著优点为：1)在小规模稀疏矩阵中，用该算法可以更快更好的恢复出满阵；2)在矩阵秩相对较大时，该算法依然具有很好的性能；3)将该算法应用在小型稀疏阵列DOA估计中，可以减少阵元间的互耦效应，更加准确地估计来波方向。下面结合附图对本专利技术作进一步详细描述。附图说明图1为本专利技术基于秩最小化的矩阵填充方法——PD-ALM算法流程图。图2为本专利技术矩阵填充方法在稀疏阵中的应用方法流程图。图3为稀疏线阵模型图。图4是本专利技术实施例1的均方根误差对比图。图5是本专利技术实施例2的均方根误差对比图。图6是本专利技术实施例3的均方根误差对比图。具体实施方式结合图1和图2，本专利技术为一种基于秩最小化的矩阵填充方法——PD-ALM算法及其在稀疏阵列DOA估计中的应用，其中稀疏阵列模型如图3所示，DOA估计方法包括以下步骤：步骤1、设置稀疏阵列的快拍总数为imax，i＝1；步骤2、将稀疏阵列第i次快拍获得的采样数据x(t)重排构造成托普利兹矩阵xT，其中，M为阵列的阵元数目，xm(t)为第m个阵元在t时刻的信号接收数据，m＝1,2,...M，矩阵xT为低秩稀疏矩阵。步骤3、运用PD-ALM算法对低秩稀疏矩阵xT进行基于秩最小化的优化矩阵填充，得到满阵x'T。引入一个辅助变量Y，该矩阵填充算法的优化模型为其中，E为等效缺失阵列，x'T为满阵。PD-ALM算法采用罚分解法和增广拉格朗日乘子法求解上述优化问题，设罚参数为ρ，μ为一个给定的正数，Z为拉格朗日乘子，优化问题可描述为：当ρ→∞时，上式的解收敛于优化问题模型的解。该算法通过外循环和内循环实现，外循环不断增大ρ，内循环更新x'T,Y。具体步骤为：步骤3-1、输入：观测矩阵xT，最大迭代次数kmax、lmax，μ和ρ的递增步长t1＞1、t2＞1，内、外循环终止条件中的参数ε1＜＜1、ε2＜＜1步骤3-2、迭代过程：(1)外循环初始化：ρ0＝ρ,l＝1(2)内循环初始化：μ0＝μ,E0＝0,Z0＝0,k＝1(3)计算(4)计算(5)计算(6)计算μk+1＝t1μk(7)计算其中，(8)计算若且k≤kmax，则进入(9)；若或k＞kmax，则进入(10)(9)k＝k+1，重复(3)-(8)(10)计算若且l≤lmax，则进入(11)；若或l＞lmax，则跳出步骤2-2，进入步骤2-3(11)计算ρl+1＝t2ρl(12)取Yk+1的第一行，重排成托普利兹矩阵得到新的Y0(13)l＝l+1，重复(2)-(12)步骤3-3、输出：步骤4、获取矩阵x'T中的第一行数据，该数据即为补全后的第i次快拍数据，然后将该数据作为稀疏阵列接收数据矩阵X'的第i列数据；步骤5、i＝i+1，重复步骤2～步骤4直到i＝imax，这时所有快拍采样数据补全完毕，得到稀疏阵列接收数据矩阵X'；步骤6、运用DOA估计算法对数据矩阵X'进行DOA估计，下面以root-MUSIC算法为例，具体为：步骤6-1、对矩阵X'进行协方差计算：R＝E{[X'][X']H}步骤6-2、对矩阵R进行特征分解，确定信号子空间US和噪声子空间UN；由于特征子空间具有方向矩阵A与噪声子空间UN正交的性质，故有：aH(θ)UN＝0步骤6-3、设多项式f(θ)为：步骤6-4、将步骤6-2中求得的噪声子空间UN代入到步骤6-3的多项式f(θ)中，求该多项式的根，所得结果即为要求的来波方向。下面结合三个实施例对本专利技术作进一步详细描述。实施例1设置均匀线阵的阵元个数分别为16、25、36、49、64，随机关闭阵元的个数为总个数的0.3。同时设置快拍数为1本文档来自技高网...

【技术保护点】
1.一种基于PD-ALM算法的稀疏阵列DOA估计方法，其特征在于，包括以下步骤：/n步骤1、设置稀疏阵列的快拍总数为i

【技术特征摘要】
1.一种基于PD-ALM算法的稀疏阵列DOA估计方法，其特征在于，包括以下步骤：
步骤1、设置稀疏阵列的快拍总数为imax，i＝1；
步骤2、将稀疏阵列第i次快拍获得的采样数据x(t)重排构造成托普利兹矩阵xT，该矩阵即为低秩稀疏矩阵；
步骤3、运用PD-ALM算法对低秩稀疏矩阵xT进行基于秩最小化的优化矩阵填充，得到满阵x'T；
步骤4、获取矩阵x'T中的第一行数据，该数据即为补全后的第i次快拍数据，然后将该数据作为稀疏阵列接收数据矩阵X'的第i列数据；
步骤5、i＝i+1，重复步骤2～步骤4直到i＝imax，这时所有快拍采样数据补全完毕，得到稀疏阵列接收数据矩阵X'；
步骤6、运用DOA估计算法对数据矩阵X'进行DOA估计。


2.根据权利要求1所述的基于PD-ALM算法的稀疏阵列DOA估计方法，其特征在于，步骤2中构造的托普利兹矩阵xT为：



其中，M为阵列的阵元数目，xm(t)为第m个阵元在t时刻的信号接收数据，m＝1,2,...M，矩阵xT为低秩稀疏矩阵。


3.根据权利要求1所述的基于PD-ALM算法的稀疏阵列DOA估计方法，其特征在于，步骤3中引入一个辅助变量Y，则基于秩最小化的矩阵填充算法的优化模型为：



其中，E为等效缺失阵列，x'T为满阵；
PD-ALM算法采用罚分解法和增广拉格朗日乘子法求解上述优化问题，设罚参数为ρ，μ为一个给定的正数，Z为拉格朗日乘子，优化问题可描述为：



当ρ＞∞时，上式的解收敛于优化问题模型的解；
该算法通过外循环和内循环实现，外循环不断增大ρ，内循环更新x'T,Y，具体步骤为：
步骤3-1、输入：观测矩阵xT，最大迭代次数kmax、lmax，μ和ρ的递增步长t1＞1、t2＞1，内、外循环终止条件中的参数ε1＜＜1、ε2＜＜1；
步骤3-2、迭代过程：
(1)...

【专利技术属性】
技术研发人员：芮义斌，高进盈，谢仁宏，李鹏，高媛，李雨航，杨恺文，季宇豪，
申请(专利权)人：南京理工大学，
类型：发明
国别省市：江苏;32