1、浅析应用空间谱估计技术测向摘要：空间谱估计测向技术是一种不同于传统的振幅测向法和相位测向法的全新测向方法，它是近三十年在经典谱估计理论基础上发展起来的，是一种以多元天线阵结合现代数字信号处理技术为基础的新型测向技术。本文介绍了空间谱估计的原理，阐述了空间谱估计的五种典型算法并通过分析针对其中三种基础算法的性能进行比较，提出了关于空间谱测向领域发展的见解。1、引言空间谱估计测向技术是近三十年来发展起来的一门新兴的测向处理技术，这种测向技术因为采用了先进的数字信号处理方法，具有传统测向体制无可比拟的技术优势，展现出良好的应用前景，成为国际无线电侦测领域的研究热点。1979 年美国人 R.O.Sch
2、midt 提出著名的 MUSIC（Multiple Signal Classification 多信号分类）算法，标志着空间谱估计测向进入了繁荣发展的阶段，经过三十年的发展，其理论已经比较成熟，但是到目前为止见诸报道的实用空间谱估计测向系统并不多。近年来，随着无线电通信技术的进步，空间谱估计测向技术的研究不断深入，推动了无线电测向领域的长足发展。2、空间谱估计的原理对电磁信号方向的测量也就是要测量辐射源入射电磁波的同相位波前，因此与传统测向系统类似，整个空间谱估计设备系统由三部分组成，即空间辐射源、空间接收阵列及算法处理器。空间谱估计测向并不是从各天线阵元所接收到信号幅度或相位的简明的数学公式
3、上直接求出 DOA，而是根据各阵元的输出信号来估计空间频率，进而求出 DOA 等参数，充分利用了各阵元信号所含综合信息，通过统计处理方法来估计 DOA。 由于电磁波不仅是时间的函数，还是空间位置的函数，如果对各天线单元的输出同时进行采样，其样本既是时域样本，也是对信号进行空间采样的样本，即对各天线单元的输出同时进行的一次采样具有空间域特性。对位于不同位置的各天线单元输出同时采样的样本进行处理，利用信号的时域函数和功率谱密度函数的特性，就能得到它的空间谱。采用空间谱估计方法计算信号的空间频率，进而确定空间入射波相对于接收天线的传播方向。 空间谱估计的一般方法是构造一个以信号方向为参数的“谱函数”
4、 ，并使得在信号的到达方向上具有尖锐的峰值，这样在进行谱分析时，其峰值就指示了信号的 DOA。3、空间谱估计的算法3.1 阵列信号的数学模型1)关于辐射源的假设辐射源是位于远场且各向同性的点源，因其传输距离远大于接收阵列的尺寸，所以可认为接收阵列处感应的是平面波。辐射源是窄带信号，即相对于信号的载频而言，信号包络的带宽很窄(包络慢变)，因此，该类信号对各阵元的影响仅有由其到达各阵元的波程差而引起的相位差异。2)关于接收天线阵的假设接收阵列由位于空间已知坐标处的阵元按一定的形式排列而成。假设阵元的接收特性仅与其位置有关而与其尺寸无关，增益均相等，相互之间的互耦忽略不计。3)关于噪声的假设假设为加
5、性高斯白噪声，各阵元上的噪声相互统计独立，且噪声与信号是统计独立的。在以上假设的基础上，建立描述阵列信号的数学模型。设 D 个辐射源入射到由 M(MD)个阵元组成的天线阵列，记第 i 个入射信号为 ，其入射角度)(tsi为 ( 包括方位角 和仰角 )，i=1，2，D。记第 k 个阵元相对于参考i点的相位差为 ( )，其输出为 (t)，噪声为 (t)，k=1，2，M ，则：kkxkn(t)= + (t)kxDits1)(ikje以 T 表示矩阵的转置，并记：(t)= (t) (t)X)(1tx2MxT(t)= (t) (t)N1n2MnT(t)= (t) (t)S)(1ts2DsT( )= (
6、) ( )A1a2aD( )= )1je2)MjeT则阵列输出的数学模型可表示为： )()(tNSAtX式中 与阵列的形状和所有信号的入射角度有关。若阵列阵形已知，则只包含信号的方向信息，因此称其为方向矩阵，相应地称 为方向矢量。A )(a3.2 MUSIC 算法设阵列输出信号的空间自相关矩阵为 R，其计算式如下IAPtXERHH2)(式中 H 表示矩阵的共轭转置，E 表示计算统计平均，P 是信号的自相关矩阵， 是噪声的自相关矩阵， 是噪声功率，I 是 MM 单位阵。I22R 是非奇异的正定 Hermitian 矩阵，可以利用酉变换实现对角化，其相似对角阵由 M 个不同的正实数组成，与之对应的
7、 M 个特征矢量是线性独立的。将R 按照特征分解的结果写成 HnHSHUU式中 是由按从大到小顺序排列的特征值组成的对角阵，U 是与特征值相对应的特征矢量矩阵。 是 中 D 个较大特征值构成的对角阵， 则由 中S nM-D 个小特征值构成， 和 是 U 中对应于 和 的部分。由于 的列矢nSnU量张成 A 的正交补空间， 则张成 A 的值域空间 RA(矩阵的值域指由矩阵S各列的线性组合构成的矢量集合，当各列为最大线性无关组时就张成了矩阵的值域空间) ，所以称 为信号子空间, 为噪声子空间。进行 DOA 估计时，信SUnU号个数 D 一般是未知的，需要先行作出正确的判断，常用的方法是信息准则(A
8、IC)。MUSIC 算法进行 DOA 估计需要在所有的 取值中寻找在 空间有),(SU最大投影的 D 个 或在其正交空间 投影的倒数最大的 D 个 ，就可以)(anU)(a得到信号的来向。在 和 上的投影分别计算为：nUS2)(STUads2)(1nTUanMUSIC 算法步骤清晰，性能稳定，需要进行一次特征分解和二维的谱峰搜索，计算复杂度和运算量适中，适合实际应用。3.3 其他谱估计算法基于特征结构的谱估计算法很多，如空间平滑、旋转子空间不变(ESPRIT)、零点预处理和加权子空间拟和(WSF)等，这些算法在实际应用中各有优缺点，具体算法如下：1)空间平滑算法将具有相同阵型、但是在阵列所处平
9、面存在平移量的若干子阵列的输出信号的自相关矩阵相加之后，对得出的矩阵进行特征分解、谱峰搜索得到DOA 估计，原理的数学描述如下：设 表示第 i 个子阵输出信号的自相关矩阵，iR则： IBPfXERHiHii 2)(式中 是第 i 个子阵的输出，B 是该子阵的方向矩阵， 是iNSBf if该子阵相对于参考阵列的旋转因子矩阵，)()(21Diiijji eedagf将所有的 相加得到： ，式中iRMiHIBQR12是满秩的。MiHPfQ1经过平滑处理后 R 具有和 MUSIC 算法中讨论的自相关矩阵一样的特性，可以按照 MUSIC、ESPRIT 等算法来进行 DOA 估计。设空间两相干来波的方位是
10、( )、( )，采用双排直线阵，空20,13,6间平滑处理后的 MUSIC 谱图如图 1 所示：图 1 相干信号谱估计图可见平滑前 MUSIC 算法完全失效，而空间平滑算法就是利用阵元的冗余性，以降低阵列有效孔径的手段达到去相干的效果。2)ESPRIT 算法ESPRIT 算法利用一对对偶阵列输出信号的自相关矩阵得到对信号空间的两个估计，通过求解这两个空间的旋转矩阵 计算信号的来向，即： ，ffVUS式中 和 是由一对对偶阵列的输出求得的对信号空间的两个估计。在最小二SUV乘的约束下得到： ,f 是对角阵，矩阵元素的定义与空间平滑中SHSUVf1)(的元素一样，据此可以计算出信号的来向。考虑到
11、和 并非精确已知，而if SVU是存在估计误差，因此实际应用中使用总体最小二乘法计算 f。ESPRIT 算法可以通过简单的三角公式计算信号的来向，不需要谱峰搜索， MUSIC 算法运算效率高，但是 ESPRIT 算法要求阵列中存在对偶阵限制了它的实际应用，并且当多个信号同时到达阵列时，ESPRIT 算法需要将方位角和仰角进行正确配对才能得到正确的 DOA 估计。改进的 ESPRIT 算法可以解决多个信号方位角和仰角的正确配对问题。将阵列划分为三个子阵：)()(321tNGSAtXF式中 F，G 为旋转矩阵，利用 F，G 可以由子阵 1 的阵列输出数据矩阵得到子阵 2、3 的阵列输出数据矩阵，并
12、且位于 F 和 G 相同位置的 和 包含了同kfg一个来波的方向信息 。采用矩阵束的方法，由矩阵 可以计算k )3,21(itX出矩阵 F，G 的对角元素，从而得到来波的 DOA 估计。3)零点预处理零点预处理算法对阵型及信号是否相干不敏感，它通过加权矩阵 W 使天线在 D 个来波信号中的某 D-1 个来向上形成零陷，阵列输出数据 X 经 W 加权即 之后对 Y 进行 DOA 估计就简化成一个信号的估计问题，可以借助XWYHMUSIC 等算法实现，因此算法关键在于求得 W。一般的解法是求解欠定的线性方程组 ，其中 为： ，我们采用特征0AHA)()(121Daa分解的方法也能得到 W。零点预处
13、理算法的具体步骤如图 2 所示。图 2 表示零点预处理算法需要知道 D 个来波方位的初始值，并经过多次迭代才能得到准确结果，因此算法成败的关键在于初始值和收敛条件的选取。初值的获取需要来波方位的先验知识或者结合多个阵列实现(例如由可以处理相干信号的阵列给出 DOA 估计作初值，再用其它阵列进行迭代达到准确结果)，不仅增加了算法的复杂性，也限制了算法的适用范围。4)WSF 算法加权子空间拟和(Weighted Subspace Fitting)算法具有比 MUSIC 算法更强的同频多信号的分辨能力(信号相干与否不影响算法的性能)和更高的精度，但实现困难。WSF 算法归结为求解下式所示的非线性最优
14、化问题：221T,minarg， FSAWU即寻找参数 使得 与 在最小平方规则下最接近。式中 W 是T，SURA一正定加权对角阵，T 是一哑变量。当 A 已知时， ，那么 ，式中211-HA)(WUTSHSUWtraceminAtrace 表示求矩阵的迹， 是 A 的正交补空间的投影矩阵。选择 ,则上式的估计具有最小渐近方差。12SSIW综上所述，以上各算法的性能比较如下表：算法 统计性能 实 现MUSIC 好 阵型不限，需特征分解和谱峰搜索，简单ESPRIT 好 阵型受限，需特征分解和三角函数计算，简单WSF 很好 阵型不限，需特征分解和迭代，复杂4、空间谱测向与传统测向方法比，具有以下突
15、出优点（1）多信号测向能力，既可以对不相关或部分相关的多个同频来波信号进行同时测向，也可以通过预处理对几个相干信号同时测向（抗多径测向能力）；（2） 测向分辨力高，突破了瑞利极限，能分辨出落入阵列同一波束内的多个信号（超分辨测向能力）；（3） 测向精度高，采用阵列信号处理方法，可以充分利用复杂的数学工具，获得更高的测向精度；（4） 测向灵敏度高，在短数据低信噪比条件下亦能获得良好测向性能。5、结论当前，测向技术领域的发展特点是：在短波大基础测向中依然是Wullenveber 体制风光依旧，而在高速宽带测向中则是相关干涉仪和 Watson-Watt 体制各有千秋，尤其是 FFT 多信道测向技术应
16、用最广泛。虽然空间谱估计测向技术拥有一些技术优越性，但是要真正在实际工作中发挥出来，还需要在多个方面下功夫。如何提高实用效果和发挥技术优势，是空间谱估计测向技术面临的最大挑战。对此，笔者认为，空间谱估计测向技术要进一步解决好以下几个问题：（1） 如何将空间谱估计测向应用到宽带侦测中去，即发展宽带多信号的空间谱估计测向技术，只有这样才能适应目前日益复杂的信号环境，这一条至关重要；（2） 空间谱估计测向必须结合 DBF 对多个来波进行空域滤波，加强其信号识别及对感兴趣信号的抗干扰接收能力。用户需要的不是一堆示向度，真正需要的是到底各个信号是从哪个方向来的，有可能还关心信号的具体信息是什么，所以一套
17、完整的空间谱估计测向系统必须具有 DBF 能力，才能充分发挥其效力；（3） 研究更稳健的、普适性更高的信号个数估计算法；（4） 进一步提高其对通道、天线失配的宽容性。浅 议 空 间 谱 估 计 测 向 技 术 的 实 用 化 问 题摘要：本文简单介绍了空间谱估计测向技术的发展历程和基本原理，详细探讨了其在实际应用中遇到的技术难题，并给出一种实用性较强的短波空间谱估计测向系统的设计思路。 关键词：空间谱估计测向MUSIC 算法 短波测向 引言 空间谱估计测向技术是近三十年来发展起来的一门新兴的测向处理技术，这种测向技术因为采用了先进的数字信号处理方法，具有传统测向体制无可比拟的技术优势，展现出良
18、好的应用前景，成为国际无线电侦测领域的研究热点。1979年美国人 R.O.Schmidt 提出著名的 MUSIC（Multiple Signal Classification 多信号分类）算法，标志着空间谱估计测向进入了繁荣发展的阶段，经过三十年的发展，可以说其理论已经比较成熟，但是到目前为止见诸报道的在实用空间谱估计测向系统并不多，这不能不说是一个遗憾。尤其是近年来，随着无线电通信技术的不断发展，无线电测向技术有了长足进步，各生产厂家纷纷推出了各种新型无线电测向设备，但是从中很难见到空间谱估计测向设备的身影。 下表是目前见诸公开报道的国际上主流公司具有代表性的一些侦测产品，可见只有美国的 DRS 和以色列的 Rafae 两家公司采用了空间谱估计测向技术
