本发明属于信号处理领域,尤其涉及一种稀疏表示模型下优化正则参数的目标角度估计方法。
背景技术:
目标角度/波达方向估计directionofarrival,doa是阵列信号处理中的一项重要研究内容,典型的doa估计算法有capon提出的最小方差谱估计法minimumvariancedistortionlessresponse,mvdr和以schimidt提出的多重信号分类multiplesignalclassification,music为代表的子空间算法,然而这些算法在低信噪比、低快拍的情况下,性能会下降很多,特别是在信号源相关的情况下,算法的估计精度会变得更差。
近年来由candes和donoho等人建立起来的压缩感知compressivesensing,cs理论引起了研究者们的极大的兴趣,稀疏表示问题成为研究的热点。在波达方向估计方面,目标信源在空域具有天然的稀疏性,可以通过稀疏表示对目标角度进行估计,不需要大量的时间采样,并且对信源的相关性也没有特别的限定,可以获得高分辨率。基于稀疏表示的目标角度估计,首先建立稀疏表示模型下的目标函数,对目标函数求解,得到与过完备角度字典一一对应的稀疏向量,通过此稀疏向量的峰值,可实现对目标角度的估计。但在对稀疏表示模型下的目标函数进行求解的过程中,涉及到正则参数的选取问题,参数选取合适与否,对稀疏表示算法的性能至关重要。
文献“malioutovd,cetinm,willskyas.asparsesignalreconstructionperspectiveforsourcelocalizationwithsensorarrays[j].ieeetransactionsonsignalprocessing,2005,53(8):3010-3022”通过偏差原理discrepancyprinciple对目标函数中的正则参数进行选取,偏差原理要求噪声的统计特性先验已知,这限制了基于稀疏表示的目标角度估计方法的实际应用。
申请号为cn201310331365.8的专利公开了一种基于软稀疏表示的doa估计方法,基于稀疏性的前提下,应用迭代加权最小方差法求解软稀疏解,来估计目标源方位。首先选取好初始值和正则化参数,确定迭代结束条件;然后将选择好的初始值和参数代入到软稀疏表示迭代公式中进行迭代;最后满足迭代结束条件,退出迭代,得到软稀疏解,确定来波方向,也即实现了信号的doa估计。但该方法在正则化参数的选取中,依赖于噪声功率,若未知噪声功率,则该方法不能施行。
技术实现要素:
本发明的目的在于公开准确度高、不需要实际噪声统计特性的一种稀疏表示模型下优化正则参数的目标角度估计方法。
实现本发明的技术思路是:首先构造一个稀疏表示模型下的目标函数;然后给定一个初始正则参数ui;对稀疏表示下的目标函数求解,得到稀疏向量
一种稀疏表示模型下优化正则参数的目标角度估计方法,包含如下步骤:
步骤(1):从m个传感器阵列获取t次采样数据,第t次采样得到一个m×1维的信号向量x(t),t次采样得到一个m×t维的信号矩阵x,x=[x(t1),x(t2),…,x(tt)];
步骤(2):设定目标所在角度空域范围[θd,θu],将目标所在角度空域范围等间隔划分为n个角度;设
步骤(2.1):设定目标所在角度空域范围[θd,θu],将目标所在角度空域范围等间隔划分为n个角度,设
步骤(2.2):将目标所在角度空域范围[θd,θu]等间隔划分为n个角度后,得到过完备的角度字典:
步骤(2.3):获得过完备阵列流形矩阵
上式中,
步骤(3):构造目标函数:
上式中,||·||f表示frobenius范数,
步骤(4):根据集中度对目标函数中的正则参数μ进行选取,得到优化的正则参数μopt:
步骤(4.1):正则参数μ初始化:对迭代变量初始化i=0,对正则参数μ设定一个初始值
步骤(4.2):用
步骤(4.3):利用凸优化工具箱求解代替后的目标函数,得到稀疏表示矩阵
步骤(4.4):定义集中度efi为:
上式中,spi(l)表示
步骤(4.5):取一个正实数ε,判断:如果1-efi>ε,使正则参数按照
步骤(5):根据优化的正则参数μopt,对目标角度进行估计,得到目标角度的估计值。
根据优化的正则参数μopt,通过凸优化工具箱求解目标函数:
得到优化的稀疏表示矩阵
本发明的有益效果为:
本发明提出稀疏表示向量的集中度概念,利用稀疏表示向量的集中度优化正则参数,保证目标角度估计的正确性;本发明不需要已知噪声的统计特性,在实际噪声统计特性未知的情况下,仍能实现正确的目标测向。
附图说明
图1是一种稀疏表示模型下优化正则参数的目标角度估计方法流程图;
图2是本发明方法选取的优化的正则参数与随机选取的正则参数下得到的空间谱的对比;
图3是信噪比为-5db时本发明方法和偏差原理方法在噪声方差等于1和噪声方差等于0.77时得到的空间谱;
图4是信噪比为5db时本发明方法和偏差原理方法在噪声方差等于1和噪声方差等于0.77时得到的空间谱;
图5是本发明方法和偏差原理方法在噪声方差等于1和噪声方差等于0.77时得到的信源角度估计的均方根误差曲线。
具体实施方式
下面结合附图来进一步描述本发明:
实现本发明的技术思路是:首先构造一个稀疏表示模型下的目标函数;然后给定一个初始正则参数ui;对稀疏表示下的目标函数求解,得到稀疏向量
如图1,一种稀疏表示模型下优化正则参数的目标角度估计方法,包含如下步骤:
步骤(1)m个传感器阵列获取t次采样数据,其中第t次采样得到一个m×1维的信号向量x(t);t次采样得到一个m×t维的信号矩阵为x,x=[x(t1),x(t2),…,x(tt)];
步骤(2)设定目标所在角度空域范围[θd,θu],将该角度空域范围划分为n个角度,其中第n个角度值为
对角度空域范围[θd,θu]等间隔划分为n个角度后,得到过完备的角度字典:
过完备阵列流形矩阵
其中,
步骤(3)构造目标函数
步骤(4)根据集中度对目标函数中的正则参数μ进行选取,得到优化的正则参数μopt,具体步骤为:
步骤(4.1)正则参数μ初始化:即对迭代变量初始化i=0,正则参数μ设定一个小的初始值
步骤(4.2)用
利用凸优化工具箱求解此目标函数,得到稀疏表示矩阵
其中,spi(l)表示
步骤(4.3)如果1-efi>ε时(ε为一个很小的正实数),使得正则参数按照式:
步骤(4.4)如果1-efi≤ε时,停止迭代,
步骤(5)根据优化的正则参数μopt,对目标角度进行估计。
具体为:通过凸优化工具箱求解目标函数:
得到优化的稀疏表示矩阵
实施例1:
以下结合仿真实验,对本发明作进一步详细说明:
(一)仿真条件
本发明仿真实验中软件仿真平台为matlabr2016a,在仿真实验中阵列是均匀直线阵,阵元间距为半波长;信源是非相关信源;噪声是独立同分布的高斯白噪声,信号与噪声不相关;在空间区域的[-90°,90°]范围等间隔进行网格划分,角度间隔为1°,建立过完备角度字典。在仿真实验中设定小的正实数ε=10-3。
(二)仿真内容及结果
仿真1,本实验用于比较本方明方法选取的正则参数和随机选取的正则参数下得到空间谱。假设空间远场有三个窄带信号源入射到阵列,入射角度分别为10°、30°和50°,三个信号源信噪比都为10db,阵元数为8,采样数是50。本发明发方选取的正则参数μ=26.8,随机选取的正则参数μ=60,这两个正则参数下的空间谱图如图2所示。
参照图2,本发明方法选取的正则参数得到的空间谱与随机选取的正则参数得到的空间谱相比较,其主瓣更加尖锐,旁瓣更低,且能够对目标信源精确定向;而随机选取的参数得到的空间谱,主瓣较宽,旁瓣很高;可以看出,本发明方法对正则参数的选取是可行的。
仿真2,本实验用于观察噪声方差的改变对本方明方法和偏差原理方法的空间谱的影响。假设空间远场有两个窄带信号源入射到阵列中,入射的角度分别为-25°和65°,阵元数为8,采样数是50,两个信号具有相同的信噪比;定义噪声方差为
参照图3和图4,当
仿真3,我们把利用偏差原理进行正则参数选取的稀疏表示方法简称为偏差原理方法;本实验用于观察噪声的方差变化对本发明方法和偏差原理方法的目标角度估计精度的影响。偏差原理中的正则参数β按照
参照图5,当
从图2~图5,可以看出,本发明对正则参数的选取是可行的,且角度估计精度较高,能够在噪声的统计特性未知情况下对正则参数进行合适选取,从而实现目标角度的正确估计。
以上具体实施例子,是为了说明本发明在均匀直线阵是能够实现的,并不是仅仅局限在这种仿真条件下;对上述具体实施方式进行修改或者对技术特征进行替换,其技术方案的本质仍然没能够脱离本发明技术方案的所属范围。
以上所述并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
综上,本发明公开了一种稀疏表示模型下优化正则参数的目标角度估计方法,其属于阵信号处理领域。本发明的优点是在噪声统计特性未知情况下,仍能得到合适的正则参数,实现目标角度的正确估计。本发明的主要步骤为:首先构造一个稀疏表示模型下的目标函数;然后给定一个初始正则参数ui;对稀疏表示模型下的目标函数求解,得到稀疏表示向量