浅析应用空间谱估计技术进行测向1

1、浅析应用空间谱估计技术测向摘要：空间谱估计测向技术是一种不同于传统的振幅测向法和相位测向法的 全新测向方法，它是近三十年在经典谱估计理论基础上发展起来的， 是一种以多 元天线阵结合现代数字信号处理技术为基础的新型测向技术。本文介绍了空间谱估计的原理，阐述了空间谱估计的五种典型算法并通过分析针对其中三种基础算 法的性能进行比较，提出了关于空间谱测向领域发展的见解。1、引言空间谱估计测向技术是近三十年来发展起来的一门新兴的测向处理技术，这种测向技术因为采用了先进的数字信号处理方法，具有传统测向体制无可比拟的技术优势，展现出良好的应用前景，成为国际无线电侦测领域的研究热点。For pers onal

2、 use only in study and research; not for commercial use1979年美国人Signal Classification多信号分类）算法，标志着空间谱估计测向进入了繁荣发展的阶段，经过三十年的发展，其理论已经比较成熟，但是到 目前为止见诸报道的实用空间谱估计测向系统并不多。近年来，随着无线电通信技术的进步，空间谱估计测向技术的研究不断深入，推动了无线电测向领域的长足发展。2、空间谱估计的原理对电磁信号方向的测量也就是要测量辐射源入射电磁波的同相位波前，因此与传统测向系统类似，整个空间谱估计设备系统由三部分组成，即空间辐射源、 空间接收阵列及算法处

3、理器。For pers onal use only in study and research; not for commercial use空间谱估计测向并不是从各天线阵元所接收到信号幅度或相位的简明的数 学公式上直接求出 DOA ，而是根据各阵元的输出信号来估计空间频率，进而求 出 DOA 等参数，充分利用了各阵元信号所含综合信息，通过统计处理方法来估 计 DOA 。由于电磁波不仅是时间的函数， 还是空间位置的函数， 如果对各天线单元的 输出同时进行采样， 其样本既是时域样本， 也是对信号进行空间采样的样本， 即 对各天线单元的输出同时进行的一次采样具有空间域特性。 对位于不同位置的各 天

4、线单元输出同时采样的样本进行处理， 利用信号的时域函数和功率谱密度函数 的特性，就能得到它的空间谱。 采用空间谱估计方法计算信号的空间频率， 进而 确定空间入射波相对于接收天线的传播方向。空间谱估计的一般方法是构造一个以信号方向为参数的 “谱函数 ”，并使得在 信号的到达方向上具有尖锐的峰值， 这样在进行谱分析时， 其峰值就指示了信号 的 DOA 。For personal use only in study and research; not for commercial use3、空间谱估计的算法3.1 阵列信号的数学模型1）关于辐射源的假设辐射源是位于远场且各向同性的点源，因其传输距离远

5、大于接收阵列的尺 寸，所以可认为接收阵列处感应的是平面波。 辐射源是窄带信号， 即相对于信号 的载频而言，信号包络的带宽很窄 （包络慢变 ），因此，该类信号对各阵元的影响 仅有由其到达各阵元的波程差而引起的相位差异。2）关于接收天线阵的假设接收阵列由位于空间已知坐标处的阵元按一定的形式排列而成。 假设阵元的 接收特性仅与其位置有关而与其尺寸无关， 增益均相等，相互之间的互耦忽略不 计。3）关于噪声的假设假设为加性高斯白噪声，各阵元上的噪声相互统计独立，且噪声与信号是统 计独立的。在以上假设的基础上，建立描述阵列信号的数学模型。设 D个辐射源入射 到由M（M>D）个阵元组成的天线阵列，记第

6、i个入射信号为Si（t），其入射角度为 6（°包括方位角，和仰角），i=1，2,，D。记第k个阵元相对于参考点的 相位差为*（0），其输出为Xk（t），噪声为氐，k=1，2，M，贝Dx三 Si（t）ej 丞（®）nxk（t）=i 二 e+nk（t）以T表示矩阵的转置，并记：X(t)=Xi(t) X2(t)N (t) = ni(t )n2 (t)nM (t)S(t)=Si(t)S2(t) s D(t)A(汩a(0i)a（心 2）(°D)a(0)=ej 1(° ej 2(° ejT则阵列输出的数学模型可表示为:X(t ) = A(0)S(t )N(

7、t )式中A°)与阵列的形状和所有信号的入射角度有关。若阵列阵形已知，则A3只包含信号的方向信息，因此称其为方向矩阵，相应地称为方向矢量。3.2 MUSIC 算法设阵列输出信号的空间自相关矩阵为 R,其计算式如下R = EX(t)XH(t)l - APY Fl式中H表示矩阵的共轭转置，E表示计算统计平均，P是信号的自相关矩阵，2 2二1是噪声的自相关矩阵， '二是噪声功率，I是 M M单位阵。R是非奇异的正定Hermitian矩阵，可以利用酉变换实现对角化，其相似对 角阵由M个不同的正实数组成，与之对应的 M个特征矢量是线性独立的。将 R 按照特征分解的结果写成RUUH 二U

8、s 上 sU? Un 上 nU：式中上是由按从大到小顺序排列的特征值组成的对角阵，U是与特征值相对 应的特征矢量矩阵。上s是上中D个较大特征值构成的对角阵，上n则由上中M-D 个小特征值构成，Us和Un是U中对应于上S和上n的部分。由于Un的列矢量张成 A的正交补空间，Us则张成A的值域空间RA(矩阵的值域指由矩阵各列的线 性组合构成的矢量集合，当各列为最大线性无关组时就张成了矩阵的值域空间 )， 所以称Us为信号子空间，Un为噪声子空间。进行DOA估计时，信号个数D 一 般是未知的，需要先行作出正确的判断，常用的方法是信息准则(AIC)。MUSIC算法进行DOA估计需要在所有的 巩二一)取值

9、中寻找在Us空间有 最大投影的D个或在其正交空间Un投影的倒数最大的D个，就可以 得到信号的来向。在Un和Us上的投影分别计算为：ds = aT(p)US dn =aT( o)UnMUSIC算法步骤清晰，性能稳定，需要进行一次特征分解和二维的谱峰搜 索，计算复杂度和运算量适中，适合实际应用。3.3其他谱估计算法基于特征结构的谱估计算法很多，如空间平滑、旋转子空间不变 （ESPRIT）、 零点预处理和加权子空间拟和（WSF）等，这些算法在实际应用中各有优缺点， 具 体算法如下：1）空间平滑算法将具有相同阵型、但是在阵列所处平面存在平移量的若干子阵列的输出信号的自相关矩阵相加之后，对得出的矩阵进行

10、特征分解、谱峰搜索得到DOA估计，原理的数学描述如下：设 R表示第i个子阵输出信号的自相关矩阵，贝R 二旦XiXiH）二 BfiPfiHBH 二21式中人=Bfi S + N是第：个子阵的输出，b是该子阵的方向矩阵，fi是该 子阵相对于参考阵列的旋转因子矩阵，f =di aej 2ej 2 ej "门1 M将所有的R相加得到:r =丄' r = bqB1二 2iMy，式中是满秩的经过平滑处理后R具有和MUSIC算法中讨论的自相关矩阵一样的特性， 可以按照MUSIC、ESPRIT等算法来进行DOA估计。设空间两相干来波的方位是（10°，20 ）、（16°，3

11、0 ），采用双排直线阵，空 间平滑处理后的MUSIC谱图如图1所示：-1O100 ralmiliicinuazimuth图1相干信号谱估计图可见平滑前MUSIC算法完全失效，而空间平滑算法就是利用阵元的冗余性, 以降低阵列有效孔径的手段达到去相干的效果。2）ESPRIT 算法ESPRIT算法利用一对对偶阵列输出信号的自相关矩阵得到对信号空间的两 个估计，通过求解这两个空间的旋转矩阵f计算信号的来向，即：L?S =VSf， 式中吃和必是由一对对偶阵列的输出求得的对信号空间的两个估计。在最小二乘的约束下得到：f = (VSVs)JVHUS,f是对角阵，矩阵元素的定义与空间平滑中fi 的元素一样，据

12、此可以计算出信号的来向。考虑到 VS和US并非精确已知，而是 存在估计误差，因此实际应用中使用总体最小二乘法计算f。ESPRIT算法可以通过简单的三角公式计算信号的来向，不需要谱峰搜索， MUSIC算法运算效率高，但是ESPRIT算法要求阵列中存在对偶阵限制了它的 实际应用，并且当多个信号同时到达阵列时，ESPRIT算法需要将方位角和仰角进行正确配对才能得到正确的 DOA估计。改进的ESPRIT算法可以解决多个信号方位角和仰角的正确配对问题。将阵列划分为三个子阵：X1(t ) = A G)S(t ) + N(t )X/t ) = A( G，)FS：t ) N(t )X3(t ) = A 0)G

13、St ) + W )式中F，G为旋转矩阵，利用F，G可以由子阵1的阵列输出数据矩阵得到 子阵2、3的阵列输出数据矩阵，并且位于 F和G相同位置的f k和gk包含了同 一个来波的方向信息°k。采用矩阵束的方法，由矩阵Xi(t)(i =伍，3)可以计算 出矩阵F，G的对角元素，从而得到来波的 DOA估计。3) 零点预处理零点预处理算法对阵型及信号是否相干不敏感，它通过加权矩阵W使天线 在D个来波信号中的某D-1个来向上形成零陷，阵列输出数据 X经W加权即 丫二时X之后对丫进行DOA估计就简化成一个信号的估计问题，可以借助 MUSIC等算法实现，因此算法关键在于求得 W。一般的解法是求解欠

14、定的线性 方程组Wab二0，其中A为：伦=也曲2)，我们采用特征分解 的方法也能得到 w。零点预处理算法的具体步骤如图2所示。图R零点预处理算法流程图图2表示零点预处理算法需要知道 D个来波方位的初始值，并经过多次迭 代才能得到准确结果，因此算法成败的关键在于初始值和收敛条件的选取。初值 的获取需要来波方位的先验知识或者结合多个阵列实现（例如由可以处理相干信 号的阵列给出DOA估计作初值，再用其它阵列进行迭代达到准确结果），不仅增 加了算法的复杂性，也限制了算法的适用范围。4）WSF算法加权子空间拟和（Weighted Subspace Fitting算法具有比MUSIC算法更强的 同频多信号

15、的分辨能力（信号相干与否不影响算法的性能）和更高的精度，但实现 困难。WSF算法归结为求解下式所示的非线性最优化问题:1USW2-A即寻找参数'?，T?使得rus '与RA在最小平方规则下最接近。式中w是 正定加权对角阵，T是一哑变量。当A已知时，T=（aha）-1ahUSW2，那么舀二min咲皿昵，式中trace表示求矩阵的迹，二斤是A的正交补空间的投影矩阵选择W= 一?s一?S,则上式的估计具有最小渐近方差。综上所述，以上各算法的性能比较如下表：算法统计性能实现MUSIC好阵型不限，需特征分解和谱峰搜索，简单ESPRIT好阵型受限，需特征分解和三角函数计算，简单WSF很好阵

16、型不限，需特征分解和迭代，复杂4、空间谱测向与传统测向方法比，具有以下突出优点（1）多信号测向能力，既可以对不相关或部分相关的多个同频来波信号进行 同时测向，也可以通过预处理对几个相干信号同时测向（抗多径测向能力）；（2）测向分辨力高，突破了瑞利极限，能分辨出落入阵列同一波束内的多 个信号（超分辨测向能力）；（3）测向精度高，采用阵列信号处理方法，可以充分利用复杂的数学工具， 获得更高的测向精度；（4）测向灵敏度高，在短数据低信噪比条件下亦能获得良好测向性能。5、结论当前，测向技术领域的发展特点是：在短波大基础测向中依然是 Wullenveber 体制风光依旧，而在高速宽带测向中则是相关干涉仪

17、和 Watson-Watt体制各有千 秋，尤其是FFT多信道测向技术应用最广泛。虽然空间谱估计测向技术拥有一些技术优越性，但是要真正在实际工作中发挥 出来，还需要在多个方面下功夫。如何提高实用效果和发挥技术优势， 是空间谱 估计测向技术面临的最大挑战。对此，笔者认为，空间谱估计测向技术要进一步 解决好以下几个问题：（1）如何将空间谱估计测向应用到宽带侦测中去，即发展宽带多信号的空间谱估计测向技术，只有这样才能适应目前日益复杂的信号环境，这一条至关重要；（2）空间谱估计测向必须结合DBF对多个来波进行空域滤波，加强其信号 识别及对感兴趣信号的抗干扰接收能力。 用户需要的不是一堆示向度，真正需要

18、的是到底各个信号是从哪个方向来的， 有可能还关心信号的具体信息是什么， 所 以一套完整的空间谱估计测向系统必须具有 DBF能力，才能充分发挥其效力；（3）研究更稳健的、普适性更高的信号个数估计算法；（4）进一步提高其对通道、天线失配的宽容性。浅议空间谱估计测向技术的实用化问题摘要：本文简单介绍了空间谱估计测向技术的发展历程和基本原 理，详细探讨了其在实际应用中遇到的技术难题，并给出一种实用性较强的短波空间谱估计测向系统的设计思路。关键词：空间谱估计测向MUSIC算法短波测向引言空间谱估计测向技术是近三十年来发展起来的一门新兴的测向处理技术，这种测向技术因为采用了先进的数字信号处理方法， 具有传统测向体制无可比拟的技术优势，展现出良好的应用前 景，成为国际无线电侦测领域的研究热点。1979 年美国人Signal Classification多信号分类）算法，标志着空间谱估