聪明文档网

聪明文档网

最新最全的文档下载
当前位置: 首页> 正在进行安全检测...

正在进行安全检测...

时间:2023-11-21 11:34:26    下载该word文档
第28卷第5期 文章编号:1006—9348(2011)05=0068—04 仿真 011年5月 卫星通信中一种改进的DOA算法 瑜,黄登山,李欢欢 (西北工业大学电子信息学院,陕西西安710072) 摘要:在卫星通信系统的阵列信号处理中,由于阵列孔径的限制,现有DOA解相干算法所能够分辨的信源数以及能够达到 的谱分辨率一直不高,导致不能精确定位信号波达方向,卫星通信无法准确传递信息。因此,提出一种分辨率高、稳定性好 的解相干算法(MSVD),利用人射信号的非圆对称性,对原阵列进行有效的扩展,得到新的虚拟阵列,再用SVD算法对新的 阵列进行谱峰搜索。仿真结果表明,MSVD算法与SVD算法相比,其精度更高、稳定性更好、分辨率更高、且可以估计的信源 数更多。 关键词:波达方向;解相干;高分辨率;高稳定性;非圆对称性 中圈分类号:TP301.6 文献标识码:A A Modifed DoA Algorithm of Satellite Communication ZHANG Yu,HUANG Deng—shan,LI Huan—huan (Dept of Electonics and nfmaon,Nohweser Polechnica Univery,Xi’an Shanxi 710072,China) ABSTRACT:Due to the limitation of the ape ̄ure of the array in array signal processing of satellie communication system,the number ofthe signals,which the existing decorrelate coherent SOUrCes agoritms could hadle,a8 we//as the spectral resolution,is not big,which leads to imprecise location of DOA,then,satelte communication cannot accurately transmi iormaton.So a high resolution and stabiliy algorthm which Ca decorrelate coherent SOUl ̄es is needed.In hi paper,such  new algorhm(MSVD)i proposed,which expands  aperure of he aray ecte- ly when the sials are noncirularly symmetc.Then we cal search spectral peak by SVD algorhm.The computer simulations show that,compared with SVD agorhm,MSVD agoritm has higher accuacy,better stabily ad high- er resoluton.Also,it Ca hadle more si sources than SVD algorthm. KEYWORDS:DOA;ecorrelation;High resoluton;Hi stabiliy;Noncirularly symmetric 引言 在现代卫星通信系统中,由于多经传播或人为干扰,信 号通常高度相关,经典的DOA算法如MUSIC、ESPRIT算法 等由于自相关矩阵降秩,不能精度估计目标信号的来向。针 对这一问题。国内外学者提出了很多去相关算法,其中最常 调制信号的实值特性,Tayem利用阵列接收数据及其共轭信 息,提出了C—ESPRIT算法 J,可处理信号源个数大幅增加, 算法分辨率很高,却仅仅适用于独立信号。本文新提出的 MSVD算法就是利用入射信号的实值特性,首先对原阵列进 行扩展,在信源数不变的情况下,将阵元数由原来的M个扩 展为2M一1个,得到新的虚拟阵列,再利用SVD算法原理对 虚拟阵列进行处理,之后进行谱峰搜索。新算法不仅可以估 计相干信源,且无阵列孔径损失,算法的角度分辨率、稳定 性、估计精度、可估计的的相关信源的数目也有了大幅增加, 用的是空间平滑类算法和矩阵重构类算法。如Shan的空间 平滑算法…,Piai的前后向空间平滑算法 ,王布宏的加权 空间平滑算法M】,高世伟的SVD算法 等等。但这些算法 都是以牺牲天线的有效阵元数来实现相干信源的解相干,虽 然可以估计相干信源,但是由于有效阵元数的减少,造成算 法角度分辨率降低、稳定性变差,可估计的相干信源的个数 下降,阵元的利用率变低,DOA估计性能变差。 后续仿真的结果也表明MSVD算法具有很好的估计性能。 数学模型 考虑在理想情况下,r个远场窄带平面波信号入射到阵 卫星通信系统中广泛应用的BPSK和MASK等调制信 号都具有实值的特点 】,即信号具有非圆对称性。根据这些 元数为M(Ⅳ≤M)的均匀线阵上,且各阵元间距为d;信号源 入射方向为0k(k=1,2,…,Ⅳ),以阵列的第一个阵元为基 准,则阵列示意图如图1所示。 收稿El期:2010—0—22修回日期:2010—05—31 各窄带信号在基准点的复包络分别为s (f), (t),…, 68一 
其中,以是一个m XP的由奇异值组成的矩阵,U是左奇 异矩阵, 是右奇异矩阵。 SVD算法具体步骤如下: 1)通过预处理得到式(5); 2)利用数据矢量按式(6)重构矩阵y; M 3)利用MUSIC算法进行DOA估计。 图1 均匀线阵阵列示意图 4 基于SVD解相干的新算法(MSVD算法) 上节的SVD算法可以分辨相干信号,但是由于是降维处 理,阵列孔径有所下降。因此,该算法最多可以估计M/2个  (t),则第z阵元接收到的信号为 Ⅳ )=∑s  )e斧 +I(t (1) 相干信源,角分辨能力下降、稳定性变差和估计精度降低,即 式中,n (t)表示第l阵元上的高斯白噪声。 将式(1)写成向量形式 X(t)=A(0)S(t)+n(t) (2) 其中 X(t)=[ l(t), 2(t),…, M(t)] A(0)=[o(0 ),o(0 ),…,n(0 )] S(t)=[s,(t), (£),…, (t)] n( )=[1,e— …,e音(…  n(t)=[n。(t),n2(t),…,nM(t)] 式中,上标 表示向量或矩阵转置。 假设信号S(t)和噪声n(t)均为零均值的广义平稳随机 过程,且互不相关,各阵元噪声也互不相关且功率为 :,则阵 列输出数据协方差矩阵和噪声协方差矩阵分别为 R:E[XX ]=ARsA + :, (3) E[n(t)n(t ]=tr21. (4) 式中,R =E[SS ]为信号的自相关矩阵,,为单位矩阵,H表 示向量的复共轭转置。 3 矢量奇异值算法(SVD算法) 由于阵列接收数据的协方差矩阵的最大特征值对应的 特征向量包含所有信号的信息,因此可以直接利用最大特征 矢量重构协方差矩阵实现解相干。矢量奇异值分解算法 (SVD算法)便利用了此原理 ' ,具体原理如下: 对式(3)进行特征值分解得到最大特征值对应的特征 矢量 N eI∑a1(n)n( )=(e …,M) (5) 重构如下矩阵 Y= (6) 式中,P=M—m+1,>N,P>N。 对y进行奇异值分解有 Y=UA (7) DOA估计性能变差。为解决这个问题,本文在SVD算法的基 础上,提出了一种新的MSVD算法。它主要利用卫星通信系 统中信号的非圆对称性,对阵列进行有效扩展得到新的虚拟 阵列 再利用SVD算法的原理对新的阵列进行谱峰搜 索,理论上可以增大阵列孔径,以此达到更高的分辨力和稳 定性,并且可以分辨更多的信源。 根据文献[3],若构造2M一1维矢量y为: (8) 式中,X =[ (t , 一 (£) ,…, (£ ] 。这里的信号 都具有一维的非圆对称性 (如AM、BPSK和MASK信号), 即5 (t :s (t)。 所以Y的解析结构为 Y= S(t)+N(t)  ).㈩+( ) 式中, =diag(exp(j¥ )…ep(j ̄M)),N (t)=(nⅣ (t) …n (t) ),B为 的前(M一1)行数,则A就为(2M一 1)行,即阵元数得到扩展。 所以,接收的数据协方差矩阵R 为 Rrr=E[1 ] (10) 对R 进行特征值分解可得: Ryr=UEUH (11) 其中,最大特征值A.对应的特征矢量 ∑ (n)n( )=( …,M-))(12) 由于最大特征矢量包含所有信号的信息,因此,重构协 方差矩阵Z,[6j: p Zl= lp+ (13)  ● lf2M 式中,m,P为奇数,且m+P=2M,m>N,P>N。 根据MMUSIC 算法得到协方差矩阵z Zz=JZl . (14) 则重构协方差矩阵z 69一 
Z=÷(Z1+Z2) 
对z进行奇异值分解有 
Z= (15) 当信号的入射角度相隔比较远时,SVD和MSVD算法都可以 分辨出信号的方向,但是MSVD算法的角分辨率明显高于 SVD算法;而当信号入射角度相隔比较近时,SVD算法性能 (16) 变差,以至不能正确估计信源的来波方向。因此,MSVD算 法的角分辨能力高于SVD算法。 式中: : o1 、0 0/  diag(Al,A2,…,A,  其中,r=min{m,P},A。 A 2≥…≥A , 为奇异值组 成的m×P矩阵,U为m×m的左奇异矩阵, 为P P的右奇  t   VV#Y- ̄ ̄ 异矩阵。 根据文献[6,0]可知,小特征值对应得左奇异值矩阵 中的矢量组成的空间即是噪声子空间 =U(:,N+1:m) (17) 其对应的导向矢量 ’ 为 0(0)=(o墨 ,…,血2 ,nl, 2,-一,血坐 ) (18) 根据下式进行谱峰搜索可以估计信号源的来波方向 一I口(p) U_ Ⅱ()I 0   ) 在上面讨论的基础上,可以得到MSVD算法的具体实现 步骤: 1)根据式(8)重构数据接收矩阵Y; 2)由式(10)计算阵列接收数据的自相关矩阵R ; 3)利用式(12)、式(13)、式(14)和式(15)重构协方差矩 阵Z; 4)由式(19)对谱峰进行搜索,估计来波方向。 根据以上的算法分析实现过程可知,MSVD算法具有以 下两个亮点: 1)高效的解相干能力。本文利用人射信号的非圆对称 性,将原来的M阵元的阵列扩展为(2M一1)个,增加了阵元 的自由度。同等条件下,MSVD算法的解相干能力较好,具 有高稳定性、高分辨率和高估计精度; 2)可估计相干信源的数目提高。对于M阵元的阵列, SVD算法最多可以估计rM/2]一1个相干信源;而MSVD算 法在M为偶数时最多可以估计(M一2)个相干信源,M为奇 数时可以估计(M—I)个相干信源。因此,MSVD算法可分辨 的相干信源数目比SVD算法多rM/2卜1个(rM/2]表示向 上取整,即取比M/2大的最小整数)。所以,在相同条件下, MSVD算法比SVD算法可估计的相干信源更多。 仿真实验和性能分析 实验条件:若无特殊说明,试验均采用M:8阵元的均匀 线阵,阵元间距d=k/2,窄带远场相干信源N=3,信噪比SNR =10dB,快拍数为K=100,噪声为方差是1的高斯白噪声。 实验1:图2和图3为SVD和MSVD算法独立实验1O次 的结果对比,图1的三个信源的来波方向分别为l0。,3O。, 35。,图2来波方向分别为l0。,20。,30。。由两图可以看出, 70一 呈 . 嚣  :        信号入射角 图2波达方向为10。.3o。。35。   L  —MSD鼻jVD ̄      篡 r…     |、 ~ i— …… ....…… 信号入射角 ● 图3波达方向为10。。20。。30。 实验2:图4为8阵元均匀线阵,6个相干信号源的来波 方向分别为10。,20。,30。,40。,50。,6o。时MSVD算法的仿真 图,此时SVD算法由于解相干时阵元孔径的降低,不能分辨 信源的方向;而MSVD算法由于将阵元有效扩展为M =2M一 1=15个,在一定程度上降低了阵元孔径的损失,因此,可以 正确估计信源的方向。此实验验证了上节的分析:与SVD 算法相比,MSVD算法可以估计更多的相干信源数。 实验3:本实验的2个相干信号源的人射角度分别为 1O。,15。,图5、图6、图7分别为两种算法的成功概率、估计偏 差和估计方差与信噪比的关系比较,两种方法各进行了200 次独立的实验。由三图可以看出,MSVD的成功概率、估计 偏差、估计方差均优于SVD算法,并且可以看出SVD算法适 用于高信噪比的情况,只有信噪比达到一定值时才能分辨信 号的来波方向;而MSVD适用于低信噪比的情况,且具有高 稳定性、高分辨率和高精度。 
0 
2O 
_40 
60 
写 
80 
担 
艇 甜 雷 l0 2O 30 4O 50 60 70 80 90 信号入射角 图4信源数较多时MSVD算法仿真  /’\ MSVr): 浩   姗 r\ 蠼  vD算法  ‘ .r 信噪比/dB 图5 SVD与MSVD成功概率比较  D算  / / 、/  粤  M: yD算法 坦 } / / : 信噪比/B 图6 SVD与MSVD估计偏差比较 结论 本文根据MASK、BPSK和AM等调制信号的实值特性, 并在SVD解相干算法的基础上,提出了一种新的解相干算 法。上一节的仿真结果表明,与SVD解相干算法相比,此算 法由于虚拟阵列的构建,在一定程度上减小了阵列孑L径的损 失,具有高分辨率、高稳定性、高精度及可以估计更多的相干 信源数的优点。它对卫星通信系统的正常运行有重要的现 实意义,为SVD算法在阵列信号出来的高效解相干领域的 应用提供了新的途径。  算j   / 枷 极 坦 l MS 法 / /、 一一 一 信噪比/dB 图7 SVD与MSVD估计方差比较 参考文献: [1]T J Shan,MWax andTKaiath.On spatal smoothingfor diecton of—arrival estmaton of coherent signals.IEEE Trans[J].0n ASSP,1985,33(8):806—811. [2] S UPilai and B H Kwon.Forward/backward spat smoothing te- chiques for coherent sial identcation[J].IEEE Trans.on AS- SP,1989,37(1):8—15. [3]黄蕾,张曙.一种DOA估计新算法及其求根形式[J].系统工 程与电子技术,2007,29(12):2026-2029. [4] 王布宏,王永良,陈辉.相干信源波达方向估计的加权空间平 滑算法[J].通信学报,2003,24(4):31-40. [5] 高世伟,保铮.利用数据矩阵分解实现对空间相干源的超分辨 处理[J].通信学报,1988,9(1):4—13. [6] 汤春荣,黄登山.一种新的高分辨稳定阵列信号估计方法 [J].计算机仿真,2009—10,26(10):115-118. [7]Habt Abeida,Jean—Pierre Delmas.MUSIC-Like Estmaton of Di rection of Arriva f Noncircular Sources[C].IEEE Transactons on Signal Processing.2006,54:2678—2690. [8]许斌,吕元恒.MUSIC与MMUSIC算法对DOA估计性能的比 较[J].火控雷达技术,2008,3(9):52—55. [9]N Tayem,H M Kwon.Conjugat ESPRIT[J].IEEE Trns On Antennas and Propagaton,2004,52(10):2618-2624. [10]王永良,等.空间谱估计理论与算法[M].北京:清华大学出 版社,2004—11. [作者简介] 瑜(1987一),女(汉族),湖北襄樊人,硕士研究 生,主要从事通信信号处理方面研究; 黄登山(1962一),男(汉族),湖南省汨罗市人,教 授,中国电子学会高级会员,主要从事通信系统、信 息处理等方面研究; 李欢欢(1988一),女(汉族),河南洛阳人,硕士研究生,主要从事微 波通信方面的研究。 

免费下载 Word文档免费下载: 正在进行安全检测...

  • 29.8

    ¥45 每天只需1.0元
    1个月 推荐
  • 9.9

    ¥15
    1天
  • 59.8

    ¥90
    3个月

选择支付方式

  • 微信付款
郑重提醒:支付后,系统自动为您完成注册

请使用微信扫码支付(元)

订单号:
支付后,系统自动为您完成注册
遇到问题请联系 在线客服

常用手机号:
用于找回密码
图片验证码:
看不清?点击更换
短信验证码:
新密码:
 
绑定后可用手机号登录
请不要关闭本页面,支付完成后请点击【支付完成】按钮
遇到问题请联系 在线客服