一种近场信源多参数联合估计的降维music方法

文档序号:9749666阅读:789来源:国知局
一种近场信源多参数联合估计的降维music方法
【技术领域】
[0001] 本发明涉及阵列信号处理领域,尤其是指一种近场信源多参数联合估计的降维方 法。
【背景技术】
[0002] 电磁矢量传感器阵元大多假设由2至6个共点放置相互正交的电偶极子或磁偶极 子构成,它可以接收入射电磁波全部的电场分量和磁场分量,获得了更多的入射信号信息, 从而相较于标量传感器具有更强的抗干扰能力和更高的分辨率等诸多优势,具有重要的军 事、民事应用价值以及广阔的应用前景。
[0003] 由电磁矢量传感器组成的阵列也称为极化敏感阵列,它不仅可以接收入射信号的 空域信息,同时可以获得入射电磁波固有的极化信息。而极化敏感阵列下的信源定位多参 量联合估计是目前阵列信号处理领域的重要研究内容之一,现有的空域-极化域参数联合 估计方法主要基于子空间类方法,如极化MUSIC(Multiple Signal Classification)方法、 极化ESPRIT (Estimation of Signal Parameters via Rotational Invariance Techniques)类方法等,上述方法多考虑目标信源位于阵列的远场区域,并假设部分参量已 知,旨在减少待估计参量个数,避免多维谱峰搜索带了的复杂运算量。若入射信源在阵列的 近场区域即菲涅尔区域内,不仅需要考虑信源方位角信息同时还要考虑信源的距离信息, 加上极化幅角与极化相位差两个极化参数,面临着需要进行四维谱峰搜索,需要付出巨大 的运算量,难以实现。

【发明内容】

[0004] 本发明提供一种近场信源多参数联合估计的降维MUSIC方法,以解决基于MUSIC方 法近场信源方位角、距离、极化辅角、极化相位差联合估计需要进行四维谱峰搜索,计算复 杂度很高,难以工程化应用的问题。
[0005] 本发明采取的技术方案是,包括下列步骤:
[0006] 步骤一:选取电偶极子与磁偶极子组成的电磁矢量传感器COLD阵列接收入射信 号,获得电偶极子与磁偶极子子阵列入射信号的协方差矩阵1^ 8]与1?[11];
[0007] 步骤二:将电偶极子和磁偶极子的接收数据协方差矩阵相加求和,获得仅包含信 源方位角与距离参数的新协方差矩阵;
[0008] 步骤三:建立谱函数,通过二维谱峰搜索获得方位角与距离估计值成;
[0009] 步骤四:通过两次一维谱峰搜索获得极化辅角与极化相位差的估计值六々
[0010] 本发明所述步骤一获得电偶极子与磁偶极子子阵列入射信号的协方差矩阵R[gg] 与1?[11]的步骤是:
[0011] K个信源入射到Μ对COLD阵元组成的线性阵列,当入射信源俯仰角Φ =90°,即入射 信源投影在固定在y轴上,由y〇z平面入射到传感器阵列,d为阵元间距且等长,设位于坐标 原点处阵元为参考阵元,第m个阵元在某t个采样时刻的接收数据为:
[0012]
[0013]
[0014] 式中,sk(t)代表第k个入射信号,和4?(f):为第m个传感器的噪声,cok = -2Ji dsin9k/A,(J)k = jid2cos9k/Ark,A代表信号波长,0k,rk, yk,nk分别为信源方位角,距离,极化 辅角和极化相位差;
[0015] 将电偶极子与磁偶极子子阵列的接收数据〃 [f](小〃以?)写成向量形式:
[0016] u[g](t)=As[g](t)+n[g](t)
[0017] u[1](t)=As[1](t)+n[1](t)
[0018] 其中,电偶极子与磁偶极子子阵列接收的K个信号分别为s[g](t)、s[1](t)
[0019]
[0020] s[1](t)=-[si(t)cos( γ i)---sk(t)cos( yk)-,,SK(t)cos( γκ)]τ
[0021] n[g](t)和n[1](t)分别是电偶极子和磁偶极子传感器噪声的向量形式,AzUPh
[· ]τ是矩阵的转置运算;
[0022] 根据电偶极子阵列接收数据和磁偶极子阵列接收数据丨i*;),计算出各自的 自协方差矩阵1^8]与1? [11]:
[0023]
[0024] Κ~?
[0025] 其中5(0) = :1(6^)3^((?.),Pk为信号功率,〇2为噪声功率,Ε{ ·}表不对矩阵求均值运 算,I表示是ΜΧΜ的单位矩阵,(·)Η表示矩阵的共辄转置运算。
[0026]本发明所述步骤二包括:
[0027] 将协方差矩阵1?[%]与1?[11]相加求和得到新的协方差矩阵R:
[0028]
[0029]容易看出,矩阵R仅包含信源方位角信息01{,实现多参量分步求解以降低算法计算 复杂度目的,对矩阵R进行特征值分解获得信号和噪声子空间,由于假设有Κ个信号,故特征 值分解结果有Κ个大特征值对应的特征矢量构成信号子空间,记为??;,:Μ-Κ个小特征值对应 的特征矢量构成噪声子空间,记为?.Β。
[0030]本发明所述步骤三包括:
[0031 ]由于阵列信号矢量张成子空间和信号子空间相同,并且与噪声子空间正交,表不 为:
[0032]
[0033]
[0034] Smnn表示向量的张成方向。构津信源方位角与距离参量的谱函数,表示为:
[0035]
[0036] 在成f平面内进行二维谱峰搜索得到信源的角度和距离参数。
[0037]本发明所述步骤四包括:
[0038] 对电偶极子接收数据自协方差矩阵1?[#]进行特征值分解,将M-K个小特征值张成 的噪声子空间记为t TBl,此谱函数中包含信源方位角、距离与极化幅角三个参量,代入步骤 三所得估计值么?,有新的谱函数
[0039]
[0040] 通过一维谱峰搜索获得极化幅角的估计值^同理,对电偶极子与磁偶极子接收数 据互协方差矩阵1^1]进行特征值分解,并代入估计值么卩,^,有谱函数
[0041]
[0042] 其中U"2为矩阵R[gl]特征值分解后得到的噪声子空间,通过一维谱峰搜索获得极 化相位差的估计值f。
[0043] 本发明将四维谱峰搜索问题分解为一个二维谱峰搜索和两个一维谱峰搜索,降低 了算法的运算量,优点是:
[0044] (1)本发明利用电偶极子与磁偶极子组成的电磁矢量传感器接收数据自协方差矩 阵的特性,将方位角与距离参数从多参量中提取出来,通过二维谱峰搜索率先求解。
[0045] (2)本发明通过分步求解过程,使得MUSIC方法在多参量联合估计下计算量级从0 {η 4}降为〇{n2+2n}(n为搜索范围内点数),大大提高了系统的运算效率。
【附图说明】
[0046] 图1是本发明的流程图;
[0047] 图2是本发明采用的电磁矢量传感器阵列的示意图。
【具体实施方式】 [0048] 包括下列步骤:
[0049] 步骤一:选取电偶极子与磁偶极子组成的电磁矢量传感器(concentered orthogonal loop and dipole,简称COLD)阵列接收入射信号,获得电偶极子与磁偶极子子 阵列入射信号的协方差矩阵1^8]与1?[11]:
[0050] K个信源入射到Μ对COLD阵元组成的线性阵列,当入射信源俯仰角Φ =90°,即入射 信源投影在固定在y轴上,由y〇z平面入射到传感器阵列,d为阵元间距且等长,设位于坐标 原点处阵元为参考阵元,第m个阵元在某t个采样时刻的接收数据为:
[0053] 式中,sk(t)代表第k个入射信号,n,[,f](i)和旧(0为第m个传感器的噪声,《 k = -23i dsin9k/A,(J)k = jid2cos9k/Ark,A代表信号波长,0k,rk, yk,nk分别为信源方位角,距离,极化 辅角和极化相位差。[0054] 将电偶极子与磁偶极子子阵列的接收数据"[,外)写成向量形式:[0055] u[g](t)=As[g](t)+n[g](t)[0056] u[1](t)=As[1](t)+n[1](t)[0057] 其中,电偶极子与磁偶极
当前第1页1 2 
网友询问留言 已有0条留言
  • 还没有人留言评论。精彩留言会获得点赞!
1