本发明涉及阵列信号处理,特别是指一种基于双松弛求解的未知相干信源测向方法。
背景技术:
1、空间谱测向技术主要包括music算法、子空间拟合算法等,但是上述算法需要预知或者估计出信源的准确数量才能完成波达角估计。在实际工程应用中,信源数量往往是未知的,此外,传统信源估计技术如信息论法、盖世圆法等,是针对非相干信源条件提出的。因此,当存在相关信源时,往往无法准确的提供信源数量信息,从而对空间谱算法的性能产生严重影响,限制了这些算法的工程应用。
技术实现思路
1、有鉴于此,本发明提出一种基于双松弛求解的未知相干信源测向方法。本发明能够为无源探测领域提供一种针对存在相关信源且信源数量未知条件的高精度测向方法。
2、本发明采用的技术方案为:
3、一种基于双松弛求解的未知相干信源测向方法,包括以下步骤:
4、步骤1:计算阵列信号的协方差矩阵r:
5、
6、上式中,x为m阵元的接收数据矩阵,n为采样点数;
7、步骤2:对r进行特征值分解,得到特征值序列∑n;
8、步骤3:对∑n进行排序,记其中最大特征值和最小特征值分别为λmax和λmin,其中λmax对应的特征矢量记为us,λmin对应的特征矢量记为un;
9、步骤4:构造导向矢量字典矩阵a(θp)及复系数序列c,字典矩阵中包含p个角度,p>>k,k为未知信源数量;
10、步骤5:引入双松弛因子εs和εn,且满足εs>0,εn>0,利用松弛因子,构造不等式凸优化条件,并计算得到系数序列c;
11、步骤6:构造空间谱p(θ):
12、
13、其中,pi∈p,ci∈c,a(θi)为a(θp)的列向量,上标h表示共轭转置,通过搜索空间谱p(θ)中的谱峰即得到信号方位角。
14、进一步地,步骤5中构造的不等式凸优化条件为:
15、minimize||c||2
16、subject to
17、||a(θp)c-us||≤εs
18、||[a(θp)c]h·un||≤εn
19、或
20、minimize||c||1
21、subject to
22、||a(θp)c-us||≤εs
23、||[a(θp)c]h·un||≤εn。
24、本发明与现有技术相比具有如下优点:
25、1、与传统空间谱测向算法方法相比,可实现在未知信源数量前提下对相干信源的准确测向;
26、2、与传统music算法方法相比,具有更高的旁瓣抑制比;
27、3、与传统基于单快拍的稀疏重构测向算法相比,充分利用了阵列的快拍数据,具有更高的算法稳定性和空间分辨率。
1.一种基于双松弛求解的未知相干信源测向方法,其特征在于,包括以下步骤:
2.根据权利要求1所述的一种基于双松弛求解的未知相干信源测向方法,其特征在于,步骤5中构造的不等式凸优化条件为: