基于MIMO雷达降维酉求根MUSIC算法角度估计.pdf
《基于MIMO雷达降维酉求根MUSIC算法角度估计.pdf》由会员分享,可在线阅读,更多相关《基于MIMO雷达降维酉求根MUSIC算法角度估计.pdf(6页珍藏版)》请在咨信网上搜索。
1、2023 年第 7 期计算机与数字工程收稿日期:2023年1月21日,修回日期:2023年2月17日基金项目:国家自然科学基金项目(编号:61801170,61801435)资助。作者简介:尹大雨,男,硕士研究生,研究方向:MIMO雷达目标定位。章飞,男,博士,副教授,研究方向:水下自主导航、多传感器信息融合、随机集理论及应用、电力电子技术与装置等。姬传堂,男,硕士研究生,研究方向:现代测控技术。1引言多输入多输出(multiple-input-multiple-outputMIMO)雷达16采用多个天线发射和接收相互正交的雷达信号,因此能够产生虚拟阵元,扩大阵列孔径,提高测量精度,从而引起了
2、国内外学者的广泛关注。近些年来,MIMO雷达波达方向711(Direction of Arrival,DOA)估计问题一直是个热点研究方向。文献 12 研究了基于均匀线阵与非均匀线阵MUSIC 的相关算法,但需要谱峰搜索,计算量较大。文献 13 提出了基于双基地MIMO雷达的求根MUSIC方法,通过传统求根MUSIC算法分别估计出DOD和DOA,虽避免了谱峰搜索,但当信噪比较低时,精度不高。文献 14 采用了降维求根MUSIC算法,该算法虽然降低了计算维度,但还需要进行复值协方差矩阵的特征值分解和计算。文献总第 405期2023 年第 7 期计算机与数字工程Computer&Digital E
3、ngineeringVol.51No.7基于MIMO雷达降维酉求根MUSIC算法角度估计尹大雨章飞姬传堂(江苏科技大学电子信息学院镇江212100)摘要针对均匀线阵的单基地多输入多输出(MIMO)雷达波达方向(DOA)问题,提出一种基于波束空间降维酉变换的实值求根MUSIC算法。该算法在传统MUSIC算法的基础上通过降维变换与波束空间方法对接收数据进行降维处理,接着酉变换构建实值协方差矩阵,然后进行实值特征值分解得到噪声子空间。通过构建多项式,用多项式求根的方法来代替传统MUSIC算法中的谱峰搜索,多项式所得的根即是目标的DOA估计。该方法不需要进行谱峰搜索,计算量大大减少,仿真结果验证了算法
4、的有效性。关键词MIMO雷达;酉变换;降维变换;波束空间;求根MUSIC算法中图分类号TN953+.5DOI:10.3969/j.issn.1672-9722.2023.07.002Angle Estimation of Music Algorithm for Dimension Reduction andUnitary Rooting Based on MIMO RadarYIN DayuZHANG FeiJI Chuantang(School of Electronic Information,Jiangsu University of Science and Technology,Zhe
5、njiang212100)AbstractTo solve the Direction of Array(DOA)problem of single-base multi-input multi-output(MIMO)radar in uniformlinear array,a real value MUSIC algorithm based on unitary transformation of dimensionality reduction in beam space is proposed.Based on the traditional MUSIC algorithm,the r
6、eceived data are processed by dimensionality reduction transformation and beamspace method,and then the real value covariance matrix is constructed by unitary transformation,and then the noise subspace is obtained by real value eigenvalue decomposition.By constructing polynomial,the method of findin
7、g the root of polynomial is used to replace the spectral peak search in the traditional MUSIC algorithm.The root obtained by polynomial is the DOA estimate of the target.The proposed method does not need to search the spectral peak,and the computation is greatly reduced.The simulation results verify
8、 the effectiveness of the proposed method.Key WordsMIMO radar,unitary transformation,dimensionality reduction transformation,beam space,root MUSIC algorithmClass NumberTN953+.51457第 51 卷15 通过酉求根 MUSIC算法,将协方差矩阵转化为实值协方差矩阵,在文献 13 的基础上降低了运算量,但当阵元数较多时,计算量依然比较可观。文献 16 提出了一种单基地MIMO雷达降维酉ESPRIT的方法,降低了传统 ESPR
9、IT算法的计算量,但需要较多快拍数来提高角度估计精度。文献17 提出了波束空间算法,以减轻角度估计算法的计算负担,且有着较好的效果。文献 18 提出一种降维变换,用于通过合并MIMO雷达虚拟阵列中的重复元素来降低 MIMO 雷达角度估计的复杂性。文献 19 提出了一种基于单基地MIMO雷达波束空间降维 MUSIC空间谱的方法,该方法取得了较高的角度分辨率,但仍然需要进行谱峰搜索,计算量大。针对均匀线阵MIMO雷达MUSIC算法存在谱峰搜索计算量大,本文提出了一种基于单基地MIMO 雷达均匀线阵的波束空间降维酉求根 MUSIC算法。该算法在传统MUSIC算法的基础上通过降维变换与波束空间方法对接
10、收数据进行降维处理,接着酉变换构建实值协方差矩阵,然后进行实值特征值分解得到噪声子空间。通过构建多项式,用多项式求根的方法来代替传统MUSIC算法中的谱峰搜索,多项式所得的根即是目标的DOA估计。该方法不需要进行谱峰搜索,计算量大大减少。2单基地均匀线阵MIMO雷达数据模型考虑包含一个均匀线阵的单基地MIMO雷达系统,如图1所示。发射阵列和接收阵列均是为间距半波长的均匀线阵,其中发射阵列包含 M个阵元,接收阵列包含N个阵元。假设存在K个不相关的信号,则接收阵列匹配滤波器的输出信号为x(t)=ar()1at()1ar()kat()ks(t)+n(t)=As(t)+n(t)(1)式中的k为第K个目
11、标的角度估计 A=ar()1at()1ar()kat()k,ar()k和at()k分别是对应于k的接收和发射导向矢量。ar()kat()k表示Kronecker积;s(t)=s1(t)s2(t)sk(t)TCk1,sk(t)=kej2fkt,fk是多普勒频率,k表示幅度;n()t表示零均值方差为2IMN的MN1维的复高斯白噪声向量。由L个快拍数构成的数据矩阵为X=AS+N(2)式中的矩阵X为X=x()t1x()t2x()tL,S=s()t1s()t2s()tL,N=n()t1n()t2n()tL。目标发射阵列接收阵列图1单基地MIMO雷达示意图3波束空间降维变换处理接收-发射联合导向矢量ar(
12、)kat()k可以表示为ar()kat()k=Gb()k(3)其中:b(k)=1exp(-jsink)exp(-j(M+N-2)ink)T则矩阵A可表示为A=GB。B=b(1)b(2)b(k)C()M+N-1 K(4)根据式(4),定义W=GHG,同时可以表示为W=diag(12 min(MN)min(MN)|M-N+121)(5)对于接收信号x()t使用降维变换W12GH得y(t)=W(-1/2)GHx(t)=W(-1/2)GHAs(t)+W(-1/2)GHn(t)=W(-1/2)Wb(1)b(2)b(k)s(t)+W(-1/2)GHn(t)=W(1/2)Bs(t)+W(-1/2)GHn(t
13、)(6)G=1000001000000100M010000100000010M0010000010000001MCMN()M+N-1(7)因为降维变换矩阵是稀疏矩阵且降维变换没有引起色噪声,所以变换带来的计算复杂度很小。波束空间变换20也是降低计算复杂度的一种方法,将之用于式(6)。定义矩阵Q:尹大雨等:基于MIMO雷达降维酉求根MUSIC算法角度估计14582023 年第 7 期计算机与数字工程Q=leftrightbbH()d()(8)当leftrigt,假 设 波 束 形 成 矩 阵TC()M+N1 F由对应于矩阵Q的较大特征值F(FM+N-1)的特征向量构成14。当F为波束空间的波束数
14、的时候,根据式(6)中的波束空间处理,得到z(t)=THy(t)=THW(-1/2)GHAs(t)+THW(-1/2)GHn(t)=THW(1/2)Bs(t)+V(9)其中V=W(-1/2)GHn(t),根据文献 17,式(9)中波束空间处理的波束形成增益定义为g()=(THW-12GHA()HTHW-12GHA()AH()A()(10)式中A()=ar()at(),Z(t)CF1,V=THV,由于THT=IF,V是一个F1的高斯白向量,均值为0,方差为2IF。对于信号模型(9)协方差矩阵RZ可由L个快排得到即RZ=1Ll=1Lz(tl)zH(tl)(11)信号子空间是由RZQS的K个最大特征
15、值所对应的特征向量EkKk=1组成。在无噪声的情况下,信号子空间可以表示为QS=THW12BDS,其中DS为满秩矩阵。可见信号子空间是THW12B的线性组合。噪声子空间是由RZQn的 F-K个较小特征值所对应的特征向量EkFk=K+1所组成的,输入信号所对应的转向矢量THW12b()位于信号子空间QS中,因此与噪声子空间QS正交。由此可得波束空间低复杂度MUSIC空间谱如下:fmusic()=1 THW12b()HQnQHnTHW12b()(12)4酉变换实值协方差矩阵显然当RZ为中心厄密特矩阵时,具有RZ=JMNR*ZJMN这一性质。一般情况下,我们采用前后平均式的方法来改善角度估计可以得到
16、:RQH=12RZ+JMNRZ*JMN(13)其中JMN表示MN维副对角线元素均为1,其余元素均为0的交换矩阵。对RQH进行酉变换处理RU=OHMNRQHOMN(14)OMN为稀疏酉矩阵,其偶数维和奇数维分别为O2K=12IKjIKJK-jJK(15)O2K+1=12IK0jIK0T20TJK0-jJK(16)由式(14)推导可知:RU=12OHMNRZOMN+JMN()OHMNRZOMN*JMN=1/2OHMNRZOMN+()O*MNHR*ZO*MN=Re()OHMNRZOMN(17)5单基地均匀线阵 MIMO 酉变换ROOT-MUSIC 算法对实值协方差矩阵RU进行特征值分解:RU=ESU
17、SEHS+EnUnEHn(18)式中US为PP的对角矩阵,其对角元素包含P个较大的特征值;Un为对角元素较小的MNP个特征值的对角矩阵;ES由P个较大特征值对应的特征向量组成;En由余下的特征向量组成。由传统ROOT-MUSIC 算法的原理可以求到单基地 MIMO雷达的空间波束降维酉求根多项式:f(z)=W12b(z)HTOMN()VVH-ESEHS(TOMN)HW12b(z)(19)令f()z=0,即只需求出式(19)的N个接近于单位圆上的根即可,也就是对于等距均匀线阵有i=arcsin2dargzi i=1N(20)由式(20)可以估计出第 i个目标的 DOA。结合以上理论分析,本文的方法
18、步骤可以总结如下:Step 1:根据单基地MIMO雷达搭建数据模型。Step 2:对接收信号进行降维变换,即y(t)=W(-1/2)GHx(t)。Step 3:基于式(8)构建波束空间变换矩阵T,利用波束空间处理获得z(t)。Step 4:对接收矩阵z(t)构造协方差矩阵RZ。Step 5:对协方差矩阵进行酉变换处理。Step 6:对实值协方差矩阵进行Root-MUSIC算法运算,估计出目标的DOA。6仿真实验本论文采用Monte Carlo实验来衡量算法的角1459第 51 卷度估计性能。定义均方根误差(Root Mean SquareError,RMSE)为RMSE=1k1k1JJ=11(
19、)KJ-k2(21)式中J=300表示Monte Carlo试验次数,KJ表示第J次实验中k的DOA估计值,k为角度真实值。本文采用的是均匀线阵的单基地MIMO雷达,阵元间距为半波长,发射阵元数M=6,接收阵源数N=6,快拍数=300,存在3个信号目标其角度分别为1=10,2=30,3=45。实验一:在信噪比SNR=5的情况下对本文角度估计算法进行100次仿真实验,如图2显示,可以看出本文算法能够同时对多个独立目标进行准确角度估计,且适用于低信噪比情况,无需进行谱峰搜索。图2本文算法在SNR=5的角度估计实验二:3个相互独立的目标,信噪比由-5dB到20dB,间隔为5dB。图3中显示的是本文算
20、法与MIMO MUSIC算法,ESPRIT算法,UESPRIT算法随信噪比变动RMSE变化趋势。表1则给出具体对比数据,可以看出,随着信噪比的不断增大,四种算法的均方根误差都逐渐变小,且本文算法精度更高,更稳定。图3角度估计精度与信噪比的变化趋势实验三:3个相互独立目标信噪比均为 5dB,收发阵元数均为6,进行蒙特卡洛实验次数为300次,图4显示了在各个算法下的角度估计均方根误差与快拍数之间的关系,表 2 给出了具体对比数据,可以看出各个算法随着快拍数的增加均方根误差逐渐减小。本文算法精度最高,且在快拍数低的时候有比MUSIC算法,ESPRIT算法和UESPRIT算法更高的精度。表1随信噪比变
21、化RMSE对比SNR-505101520Our method0.23090.11020.05620.02640.01530.0101MUSIC0.45290.23950.12930.07110.04170.0229ESPRIT0.6450.35170.18650.10950.06860.0321UESPRIT0.25020.1280.07080.04030.02110.012图4角度估计精度与快拍数的变化关系表2随快拍数变化RMSE对比Snapshots50100150200250300Our method0.03350.02720.02670.02510.02450.0237MUSIC0.0
22、6730.04980.04380.03840.03680.034ESPRIT0.08140.06510.04960.0440.03860.0335UESPRIT0.04840.04520.04490.04140.03680.0323实验四:快拍数=50,信噪比=10dB,在阵元数分别为10、20、30、40的情况下进行300次蒙特卡洛实验。图5显示了在不同阵元数的情形下角度估计精度的变化趋势,表3给出了具体数据对比,可以看出,随着阵元数的增加,各算法精度都在逐渐提高。在低阵元数的情况下本文算法有着比MUSIC算法,ESPRIT算法和UESPRIT算法更高的精确度。图5角度估计精度与阵元数的变化
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于 MIMO 雷达 降维酉 求根 MUSIC 算法 角度 估计
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。