1、2008 年 5 月 Journal on Communications May 2008第 29 卷第 5 期 通 信 学 报 Vol.29 No.5共轭传播算子测向算法刘剑 1, 王延伟 1, 黄知涛 2, 周一宇 2(1. 空军工程大学 电讯工程学院，陕西 西安 710077；2. 国防科技大学 电子科学与工程学院，湖南 长沙 410073)摘 要： 针对一维信号提出了共轭传播算子测向算法(COP)，新算法利用了阵列输出数据的共轭，通过对阵列孔径的扩展，可对多于阵元数的信号进行测向，其分辨力和测角精度优于 OPM(正交传播算子测向算法)和 MUSIC算法。分析了新算法的均方误差性能和计算
2、复杂度，得到了均方误差的解析表达式。仿真实验验证了 COP 算法的优良性能，均方误差的理论值与仿真值相符。关键词：阵列信号处理；测向；阵列扩展；传播算子方法中图分类号：TN911.7 文献标识码：A 文章编号：1000-436X(2008)05-0013-06Conjugate propagator method for direction findingLIU Jian1, WANG Yan-wei1, HUANG Zhi-tao2, ZHOU Yi-yu2(1. Telecommunication Engineering Institute, Airforce Engineering Un
3、iv., Xian 710077, China;2. School of Electronic Science and Engineering, National Univ. of Defense Technology, Changsha 410073, China)Abstract: A new direction finding algorithm called conjugate propagator method (COP) for one-dimension signals was proposed. The new algorithm uses the array output a
4、nd its conjugate form to extend the array aperture so that it can detect more signals than sensors and has better performance (in terms of resolution and angular precision) than OPM (orthonormal-PM) and MUSIC algorithms. The mean square error (MSE) performance of COP was analyzed and the analytical
5、expression of MSE was given. The computational complexity of COP was analyzed also. Simulation results show the better performance of COP than OPM and MUSIC algorithms and the accordance between the theoretical results and the simulation results. Key words: array signal processing; direction finding
6、; aperture extension; propagator method1 引言以 MUSIC 算法 1为代表的子空间类高分辨测向算法可对多信号同时测向，测向精度也远高于传统测向方法，但这些算法需要对协方差矩阵进行特征分解(ED)或对数据矩阵进行奇异值分解(SVD)。传播算子方法(PM) 2不需要 ED 或 SVD，信号数远小于阵元数时计算量明显小于 MUSIC 算法。正交传播算子方法(OPM)在中高信噪比下的性能接近于 MUSIC 算法，低信噪比下对信号个数适当过估计其性能也接近于 MUSIC 算法 2。但对OPM，可测向信号数小于阵元数、低信噪比小阵列情况下，若信号个数仅略小于阵元数
7、，则 OPM算法难以最佳过估 计 信 号 数 ， 从 而 性 能 下 降 。 利 用高 阶 累 积 量 对 阵 列 进 行 扩 展 3可 以 解 决 这 2 个 问 题 ，但 估 计 高 阶 累 积 量 矩 阵 的 计 算 量 很 大 。 为 增 加 可 测向 信 号 数 并 提 高 测 角 精 度 和 分 辨 力 ， 同 时 保 持 较 小的 计 算 量 ， 本 文 针 对 一 维 信 号 提 出 了 共 轭 PM 算 法收稿日期：2006-08-21；修回日期：2008-03-01基金项目：国家自然科学基金资助项目 (60502040)Foundation Item: The Nation
8、al Natural Science Foundation of China (60502040)14 通 信 学 报 第 29 卷(COP)。2 数据模型及传播算子方法的测向原理2.1 数据模型假设有 D 个不相关窄带信号入射到 M 元阵列上，信号入射方向与阵列共平面，阵列位于远场。文中认为信号数 D 已知或已估计得到了。上标“*”表示共轭， “T”表示转置， “H” 表示共轭转置。阵列输出矩阵为(1)T011,MXxASN其中 为信号矩阵， DSs 0,n为噪声矩阵，T11,Mn为方向矩阵，01(),()Aaa为信号矢量，is为噪声矢量, ,0,0,1,122j(cosin)j(cossi
9、n)T()e,exy xMyddda为方向矢量， 为波长，( )为阵元 i 的坐标xiy(假设阵元 0 为参考阵元，位于原点，即)。假设信号为平稳的零均值随机过,0,xyd程，信号 i( )的功率为 ，噪声的实1,D 2,si部与虚部不相关，每个阵元上的噪声功率均为 ，2n信号与噪声不相关。2.2 传播算子的定义假设导向矩阵 A 是列满秩的，则 A 中有 D 行是线性独立的，其他行可由这 D 行线性表示。以下分析中，总是假设其前 D 行是线性独立的。将导向矩阵分块为(2)ABM其中， 和 分别为 DD 维和(M D)D 维矩阵。A假设 非奇异，传播算子定义为由 MD 维复空间 到 D 维复空间
10、 的惟一线性算子MCCP，P 满足(3)HABP或(4)(),MDMDIQO式中， 表示 的单位矩阵。()可以证明 2，Q 的列张成的空间就是噪声子空间。2.3 传播算子方法的测向原理式(4)表明，导向矢量 与 Q 的列向量正交，)(a则 DOA 为式（5）的解(5)H()1MDO由此可得到类似 MUSIC 算法的空间谱函数(6)HPM()()fa搜索式(6)空间谱的负峰就可以得到信号的波达方向(DOA) 。如果噪声是白的，为提高算法性能，用矩阵 Q 的标准正交化形式代替 Q，则有(7)H1/20()正交传播算子方法(OPM)的空间谱函数为 2(8)OPM()fa0()a其中， 为 的估计值，
11、估计方法如下。0Q2.4 传播算子的估计引入数据矩阵如下形式的分块(9)ABDX其中， 和 分别为 DL 维和( MD)L 维矩阵，AL 为快拍数。传播算子的估计值 由最小化下面的代价函P数得到(10)()J2HBAX其中， 表示 Frobenius 范数。代价函数 为 J的二次函数，最小化 的最佳解为P(11)H1()AABP亦可由阵列输出的协方差矩阵估计得到，方法与以上过程类似，参见文献2。3 共轭传播算子测向算法3.1 阵列扩展测向中阵列扩展的概念首先是由 Dogan 等针对四阶测向算法提出的 3。Dogan 等利用虚拟互相关概念得到了四阶累积量与互相关的关系，将四阶累积量矩阵映射到协方
12、差矩阵，相当于阵列孔径扩展。假设信号是非高斯的而噪声是高斯的，则协方差矩阵和四阶累积量矩阵可分别写为 4 第 5 期 刘剑等：共轭传播算子测向算法 15(12)H21(,)()DxSiinMiRaI(13)*H1(,()xSiiiii Qa其中， 和 分别为信号矢量的协方差矩阵和四阶累积量矩阵， ， 是由四阶累积HESRSQ量 组成的矩阵。相对于式(12)，式*,iklmcus(13)的导向矢量得到了扩展，而导向矢量与阵元位置是对应的，所以利用四阶累积量可以虚拟地增加阵元、扩展阵列孔径，从而提高算法的测向精度和分辨力 4,5。阵列扩展是通过导向矢量的扩展实现的，这一思想也可用于传播算子测向算法
13、。传播算子测向算法中，要求导向矩阵是列满秩的，所以信号数 D 必然小于阵元数 M。由式(2)的矩阵分块，要想增大可测向的信号个数 D，必须增加 A 的不相关的行数，一个简单的方法就是将 A 的共轭与 A同时使用，构成扩展的导向矩阵。参考阵元在 A中对应的行与其共轭是相同的，将这一行去掉。定义扩展导向矢量为 ，其中，)(bTT1(,)a)(1a T，其中，,(1),(1)2jcossine,xMydd ,1,12j(cosin)exyd， ， 。这就ii, iyi, M相当于虚拟出了坐标为的阵元，阵列扩),()( 11(,)1, xyx展的示意图见图 1。图 1 中实点表示实际阵元，空心点表示虚
14、拟出的阵元。扩展导向矩阵为。对 B 如下分块011(),()DBbb(14)2ABM假设 是非奇异的，于是得到扩展传播算子 ，A EP满足EP(15)HEABP即(16)HE21E(21),MDMDIQO图 1 阵列扩展示意图利用 QE 与扩展导向矢量 的正交性可进行)(b测向，如何由阵列输出得到 QE 请见 3.2 节。3.2 共轭传播算子测向算法(COP)对于 BPSK 和 AM 等一维信号， 6。*S记 ，V1011(),()Daa，符号“”表示无噪声量。将*VXSA式(14)乘以 可得(17)V21ABMDXYSY其中， 和 分别为 DL 维和 (2MD1)L 维矩A阵。此时有(18)
15、HEABPY存在噪声时，扩展传播算子的估计值 由式EP(19)给出(19)H1E()AB其中， 为存21VDMYXYN在噪声的数据。在此注意到，扩展后噪声 的统EN计特性满足(20)H2E,kpnLkpI(21)*,MN其中， 为 的第 i 行。类似于式(16)，可得到Ei(22)EE21MDPQI对其进行标准正交化16 通 信 学 报 第 29 卷(23)0EQHE(1/2)则共轭传播算子测向算法(COP)的空间谱函数为(24)COP0E()fb()b由 PM 的原理及 3.1 节阵列扩展的概念可知，只要扩展导向矩阵的前 D 行是线性独立的， COP算法就适用。也就是说，当信号入射角度各不相
16、同时，COP 算法适用于均匀线阵以及所有满足导向矩阵前 D 行线性独立的其它阵列。3.3 COP 算法测向的均方误差性能记无噪声时 COP 算法的零谱函数为(25)H0E()()fbQ其中， ， ，1/20EEQE21MDPI。为书写方便，将存在噪声时H1()ABPCOP 算法的零谱函数式(24) 记为 ，则 对)(f)(f的一阶导数为()fH0E()dbH0EQ(26)Q0E2Re0()d其中， 为 对 的一阶导数。由式 (16)，()，则 的二阶导数为Hb)(fHH0E0Ee()()f db(27)2d记 为 的估计值，则 ，由一阶泰勒f展开式得 7 0()ff()f(28)式(28)中最
17、后一步忽略的是 的高阶小量。所以(29)()fAHRe()()bQdd0E式(28)中， Hb0EHTN0E2。所以，()Qd0E(30)HEe()()ddQ0式(30)中， ， 。记1(2)AMDLYT0ABS，则标量 u 满足HE()ubNdH(31)*1ReeR2uu记 和 HE121,MTNgH012,EQh第 5 期 刘剑等：共轭传播算子测向算法 17，则 ，所以21Mh21H()()Mkkubghd21HHT*1Re() ()()kpkpMkpkk bghddh(32)其中，。所以THT2TE()kpEkpngN,2kMp212T*HT*2Re()()()()Mn kku bdhd
18、(33)类似于式(32)，有(34)*21HH21Re()()()eMkppkppkkpkpu bghdbhg其中， ，HHEE()kpgTN2nH,kpT又 ，所以0QQ00*2 0Re()()()n EubdQd(35)显然， ，H1(21)H(2)ADMMDLBT0,1TT(21)(2)(21)AD 0其中， 。利用H,0,EdiagssS式(33)、式(34)和式(31) 及式(30)，可得信号 i 的测向均方差为(36)HH12 21T1TT*HT*2H()(E()Re()()()()()AiAiii MAiAAAiikMikiiLL SNSNSSbBRbdQdb dhd0E0E其中
19、， 为 的前 D 个元素组成的列向量，()Aii，011dag,SNSNSN为信号 i 的信噪比。2,isinR3.4 COP 算法的计算复杂度传播算子类测向算法最大的优点是不需要 ED或 SVD，从而降低计算量。MUSIC 算法中，对协方差矩阵进行 ED 所需的计算量约为 次复)(3MO数乘法，由数据矩阵估计协方差矩阵约需 次L2复数乘法，直接对数据矩阵进行 SVD 的计算量为次复数乘法。通常 LM，所以 MUSIC)(2LMO算法的主要计算量为 次复数乘法。COP)(2算法中，由式(19) 计算扩展传播算子的计算量为次复数乘法，对 进)1(3D EQ行 Gramm-Schmidt 正交化需
20、要 (1)2O次复数乘法。当 时，COP 与2)DLMSUIC 算法的计算量之比约为 。这说/明，阵元较多而信号较少时，COP 的计算量小于MUSIC 算法。4 仿真结果仿真中采用均布 4 元线阵，阵元间距半波长，信号形式 BPSK，信号间不相关，高斯白噪声，快拍数 50，主要验证 COP 算法的阵列扩展能力、分辨力与测角精度、信号数过估计的影响。实验 1 阵列扩展。假设有 6 个 BPSK 信号入射到阵列，波达方向为40，60 ，80 ，100，120，140 ，信噪比均为 10dB。COP 算法用了阵列流形的共轭，等效于阵列孔径扩展，故可对多于阵元数的信号进行测向。在均布线阵条件下，COP
21、 算法可将 M 元阵列扩展到 (2M1)元阵列，可测向信号数为(2M2)个。图 2 为 COP 算法 20 次空间谱的叠加。第 5 期 刘剑等：共轭传播算子测向算法 17图 2. COP 算法的空间谱实验 2 分辨力与测角精度。2 个信号从85，95入射，等信噪比 220dB，分辨力和测角精度由 300 次测向结果的统计得到。仿真中，空间谱只有一个负峰时认为分辨不成功，并令两个角的估计值都等于这个负峰对应的角度，第 i ( i=1,2)个入射角的均方根误差定义为，总的测角均方根误差为3021()/30iij j两个信号均方根误差的平均，分辨概率定义为成功分辨两个源的次数与总实验次数的比值。由图
22、3 可以看出 COP 算法的均方根误差明显小于 OPM和 MUSIC 算法，而较低信噪比下分辨概率明显高于 OPM 和 MUSIC 算法，并且由式(36)得到的COP 算法测角均方误差的理论值与实际仿真值很接近。实 验 3 信 号 数 过 估 计 对 COP 算 法 分 辨 力 和 测角 均 方 根 误 差 的 影 响 。 信 噪 比 为 2dB， 其 他 仿 真 条件 同 实 验 2， 图 4 给 出 信 号 数 估 计 值 为 26 时 COP算 法 的 性 能 。 文 献 2已 指 出 ， 低 信 噪 比 下 对 信 号 数适 当 过 估 计 会 提 高 传 播 算 子 类 算 法 的
23、测 向 性 能 。 由图 4 可 知 信 号 数 估 计 为 5 时 COP 算 法 的 分 辨 力 和 测角 均 方 根 误 差 性 能 均 达 到 最 佳 。 文 献 2提 出 ， 信 噪比 较 低 时 信 号 数 适 当 过 估 计 可 使 OPM 的 性 能 接 近于 MUSIC 算 法 ， 而 图 3 中 OPM 性 能 明 显 比MUSIC 差 ， 这 是 因 为 阵 列 孔 径 较 小 ， 可 以 对 信 号数 过 估 计 的 余 地 很 小 。 而 COP 算 法 虚 拟 出 了 更 多 的阵 元 ， 可 由 过 估 计 信 号 数 进 一 步 提 高 算 法 性 能 ， 如图
24、 4 所 示 。(a) 分辨力(b) 均方根误差图 3 分辨力与测角均方根误差(a) 分辨力(b) 均方根误差图 4 信号数过估计对 COP 算法性能的影响18 通 信 学 报 第 29 卷5 结束语提 出 了 适 用 于 一 维 信 号 的 共 轭 传 播 算 子 测 向 算法 (COP)， 新 算 法 利 用 了 阵 列 输 出 数 据 的 共 轭 ， 扩展 了 阵 列 孔 径 ， 可 对 多 于 阵 元 数 的 信 号 进 行 测 向 ，其 分 辨 力 和 测 角 精 度 优 于 OPM 和 MUSIC 算 法 ， 低信 噪 比 下 对 信 号 数 适 当 过 估 计 可 进 一 步 提
25、 高 算 法 性能 。 通 过 对 COP 算 法 零 谱 函 数 的 一 阶 泰 勒 展 开 分 析了 其 均 方 误 差 性 能 ， 得 到 了 均 方 误 差 的 解 析 表 达 式 ；分 析 了 COP 算 法 的 计 算 复 杂 度 。 仿 真 实 验 验 证 了COP 算 法 的 优 良 性 能 ， 均 方 误 差 的 理 论 值 与 仿 真 值相 符 。参考文献：1 SCHMIDT R O. Multiple emitter location and signal parameter estimation J. IEEE Trans，AP-34, 1986, (3): 276-2
26、80.2 MARCOS S, MARSAL A, BENIDIR M. The propagator method for source bearing estimationJ. Signal Processing-42, 1995, (2): 121-138.3 DOGAN M C, MENDEL J M. Application of cumulants to array processing-part I: aperture extension and array calibrationJ. IEEE Trans SP-43, 1995, (5): 1200-1216.4 CHEVALI
27、ER P, FERREOL A. On the virtual array concept for the fourth-order direction finding problemJ. IEEE Trans SP-47, 1999, (9): 2592-2595. 5 PORAT B, FRIEDLANDER B. Direction finding algorithms based on high-order statisticsJ. IEEE Trans SP-39, 1991, (9): 2016-2023.6 TAYEM N, KWON H M. Conjugate ESPRIT(
28、C-SPRIT)J. IEEE Trans AP-52, 2004, (10): 2618-2624.7 STOICA P, NEHORAI A. MUSIC, maximum likelihood, and cramer-rao boundJ. IEEE Trans ASSP-37, 1989, (5): 720-741.作者简介：黄知涛（1976-），男，湖北荆州人，博士，国防科技大学副教授，主要研究方向为循环平稳信号处理、综合电子战。周一宇（1948-），男，上海人，博士，国防科技大学教授，主要研究方向为综合电子战。刘剑（1978-），男，山东潍坊人，博士，空军工程大学电讯工程学院讲师，主要研究方向为阵列信号处理、综合电子战、卫星通信。王延伟（1978-），男，河南平顶山人，硕士，空军工程大学电讯工程学院讲师，主要研究方向为无线通信理论。
