空天混合双基地雷达杂波建模.docx

1、空天混合双基地雷达杂波建模ISSN 1000-0054CN 11-2223/N 清华大学学报(自然科学版J T singh ua Un iv (Sci &Tech ,2007年第47卷第10期2007,V o l.47,N o.1013/401606-1609空天混合双基地雷达杂波建模李华,汤俊,彭应宁(清华大学电子工程系,北京100084收稿日期:2006-09-26作者简介:李华(1981,女(汉,四川,博士研究生。通讯联系人:彭应宁,教授,E -mail :peng yn hsp .ee .ts inghua .edu .cn摘要:为了研究空天混合双基地雷达中的杂波特性和杂波抑制算法,提

2、高对慢速目标的检测性能,该文分析了空天混合双基地雷达的几何关系,获得杂波等距离环的解析解并建立了杂波模型。同时基于该模型,对几种特定空间几何配置下的杂波等距离环和杂波角度Do ppler 线进行了仿真分析。结果表明混合模式的杂波特性对卫星平台速度和位置更加敏感,合理的双基地几何配置可以提高空时二维信号处理的性能,达到更好的最小可检测速度。关键词:双基地雷达;空天混合双基地雷达;杂波建模中图分类号:T N 957文献标识码:A文章编号:1000-0054(200710-1606-04Clutter modeling for Hybrid bistaticspace -based radarLI

3、Hua ,TANG Jun ,PE NG Yingning(Department of Electronic Engineering ,Tsinghua University ,Beij ing 100084,China Abstract :Hybrid bistatic space-bas ed rad ar (HB-SBRuses atran smitter on a s atellite w ith th e receiver on th e airplane.Cluttermodeling is introduced to r educe clu tter and improve sl

4、ow -moving target detection in HB-S BR.Simulations of the clutter isorange contours and clutter angle-doppler traces for varies g eometr ies sh ow s that the clutter char acteristics of th e h ybrid bistatic mode are very sen sitive to th e s atellite ch aracteristics and a better m inimum detection

5、 velocity can be achieved w ith the pr op er bis tatic configuration.Key words :bis tatic radar;hybr id bis tatic space-bas ed radar;cluttermodel天基雷达系统由于具有全天时、全天候、全球大范围区域的信息获取能力,近年来正逐渐成为一个新的研究热点。单基雷达系统与双/多基地雷达系统相比有着系统成本低,抗电子干扰、抗摧毁等很多优点。实际的天基双/多基地雷达系统可以采用空天(飞机/卫星混合平台、天基(卫星/卫星平台等多种形式1,2。由于天基雷达系统采用下视工作

6、的方式,同时由于卫星平台具有极高的运动速度,使得天基雷达系统面临着复杂的杂波背景。双基地雷达系统的杂波特性比单基地系统的杂波特性更为复杂,而系统杂波的建模研究工作对于后续的杂波抑制算法研究具有至关重要的意义。近年来,国内外已经开展了一些机载双基地雷达杂波的建模和仿真工作并分析杂波特性3-6。本文提出了针对空天混合双基地雷达的杂波建模方法,该模型的接收机采用无人机平台,发射机放置在卫星平台上。并在该模型基础上,对几种特定空间几何配置下的空天混合双基地系统的杂波等距离环和杂波角度Doppler 关系进行了仿真和分析。1空天混合双基地雷达几何关系分析在单基地雷达系统中,接收机在一个距离门内接收到的杂

7、波回波来自于一个等距离圆环。在双基地雷达系统下,对应于一个距离门的杂波等距离环为以接收机和发射机为焦点的椭球与地球球面的交线。由于两个曲面交线的求解比较复杂,本文将利用合理的近似来求出该等距离环的解析解。由于地球是球体,接收机对地面的最大覆盖范围是一个球冠。但考虑到接收机是无人机,它的飞行高度低,对地面的覆盖范围很小,因此可以认为接收机对地球的覆盖范围为一个圆平面,定义为接收圆。来自接收圆外的回波信号将无法被接收机接收到。如图1所示,我们以接收机与地心O 的连线为Z 轴,与地面的交点O 为原点,与Z 轴垂直的平面为X OY 平面。发射机和地心连线与X OY 平面交于P 点,O P 的方向即为X

8、 轴,收发机到地心连线的夹角为 。椭球与X OY 平面的交线为杂波等距离环,它是一个椭圆。要注意的是,当 较大时,在大多数情况下,只有部分椭圆环在接收圆内,因此杂波环的积分范围发生变化,不再是一个完整椭圆,使得杂波回波只存在部分方位角内。这与机载双基地雷达不同,在机载情况下,发射机和接收机在地面的投影一般认为在同一平面内,杂波等距离环则为该平面内的完整椭圆,即杂波信号是从各个方向来的。从后面的分析还将看出, 不仅影响等距离环在接收圆内的形状,还将影响杂波谱随距离的变化率,是 一个很重要的参数。图1空天混合双基地几何关系2杂波等距离环的建立杂波等距离环的位置和形状与发射机和接收机的坐标及双基地距

9、离L 有关。接收机和发射机离地面的高度分别为H r 和H t ,它们到地心连线的夹角为 。则它们的坐标为:R :R x =0,R y =0,R z =H r ;T :T x =(H t +R e sin ,T y =0,T z =(H t +R e cos -R e .其中R e 为地球半径。只需先求出椭圆在OX Y 平面内的4个顶点A BCD 的位置,就可以计算出焦点,从而写出椭圆方程。由分析可知,该椭圆关于X 轴对称 ,长轴顶点和焦点都在X 轴上。先求解长轴顶点AB 。如图2所示,写出椭球在X OZ 平面内的椭圆方程,焦点为R 和T 。然后令Z 为零,就得到A 和B 的坐标:A x ,B

10、x =-( 2a 2S z sin -2b 2S x co s 2(b 2cos 2 +a 2sin 2.(1其中:a =L /2,b =a 2- D RT /2 2,D RT 为接收机和发射机之间的距离。=arctanT z -R zT x -R x,=(2a 2S z sin -2b 2S x cos 2-4(b 2cos 2+a 2sin 2(b 2S x +a 2S z -a 2b 2,S x =(R x +T x /(2cos +(R z +T z /(2sin ,S z =-(R x +T x /(2sin +(R z +T z /(2cos .再利用椭圆的对称性,得到短轴顶点C

11、和D 的坐标为A x +B x2,C y ,仅有y 坐标未知,代入D RC+D TC =2L 可解。D RC 、D TC 分别为C 点到R 和T 的距离。图2双基地椭球解出C 、D 的坐标为C y ,D y =E 2-A x +B x2-T x2-T 2z .(2其中E =4L 2+(A x +B x (R x -T x +T 2x +T 2z -R 2x -R 2z4L.3空天混合双基地雷达地杂波模型双基地雷达与单基地雷达最大的区别是杂波单元相对接收机的Doppler 频率f t 与发射机和接收机的相对关系和运动方向都有关系,杂波等距离环上任意一个散射点M 的Doppler 频移为f t =

12、V r !cos r +V t !cos t ,(3其中:!为波长;V r 、V t 分别为接收机R 和发射机T 的速度矢量,其方向一般认为平行于X OY 平面。假设卫星轨道为圆轨道,由向心力可知,V t 只能在垂直O R 的平面内。#r 、#t 分别为V r 、V t 在X OY 平面内的投影V r 、V t 与X 轴的夹角。r 、t 分别为RM 、TM 与V r 、V t 所成的锥角,其中RM 、TM 分别为接收机和发射机到杂波单元M 的向量。锥角采用内积的方法得到co s r =V r RM V r RM ,cos t =V t TM V t TM .设发射机和接收机都采用面阵,阵元之间

13、等间1607李华,等:空天混合双基地雷达杂波建模距,且阵面都采用正侧式放置。发射机和接收机采用的阵元数目分别Q t N t 和Q r N r ,Q 表示行数,N 表示列数。将杂波等距离环分成N c 个杂波单元,则第i 个杂波单元相对于接收机的空间频率为f s i =d !cos r i ,(4归一化的时间频率为f t i =V r !f r cos r +V t !f rcos t ,(5其中f r 为脉冲重频。定义a i 和b i 为空域Fo urier 导引矢量和时域Fo urier 导引矢量:a i =1,ej2 f s i,ej2 fs i,ej2(N r-1 fs iT,b i =1

14、,ej2 ft i,ej2 f t i,ej2(K -1 ft iT.则二维Four ier 导引矢量S i 即定义为矢量a i 和b i的Kronecker 积,即S i =b i a i ,S i 为N K 1维,K 为脉冲数。由上面建立的导引矢量可以将雷达接收到的杂波单元回波的一组空时快拍数据表示为x i =%i S i ,%i 为该杂波单元信号回波的幅度,其平均功率为 i ,根据雷达方程,i =P t G t G r !2&i d A(43D 2t i D 2r i.(6式中:P t 为发射功率,G t 、G r 为在该杂波单元上的发射、接收天线增益;!为波长;&i 为该杂波单元的散射

15、系数,它与发射接收角度有关,此时为简便,假设相同的&i 。D t i 、D r i 是发射和接收机到杂波单元的距离;d A 为该微小杂波单元的面积。雷达接收到的该距离环上的杂波信号为N c 个杂波单元回波的和,即X c =Nci =1%iS i.(7X c 称为该距离单元上的一个时空快拍。杂波协方差阵为R c =E X c X H c=Nci =1i S iSH i.(84空天混合双基地雷达杂波特性分析STAP 算法通常需要很多训练数据来估计杂波协方差矩阵,这些数据来自与待检测距离单元相邻的距离单元。如果该数据随距离变化是不均匀的,则将增加杂波协方差矩阵的估计误差,降低ST AP 算法的杂波抑

16、制性能。杂波随距离的均匀性可以由杂波角度-Do ppler 线来描述,该线反映的是杂波谱背脊线的空间频率与时间频率的关系5,不同的距离门可以画出相应的角度-Doppler 谱线,这些谱线的重合程度反映了杂波信号随距离变化的特性。在单基地模式下,如果天线采用正侧视放置,杂波角度-Doppler 谱呈线性且在不同的距离单元都是重合的。在双基地模式下,杂波角度-Doppler 谱随距离变化不再重合,加大了利用空时二维信号处理的难度。对于机载双基地雷达,发射机与接收机的速度大小差别不大,由式(5可以看出一个杂波单元的时间频率大小将同时受到发射机和接收机运动的影响。然而在空天混合双基地模式下,发射机放置

17、在卫星上,它的运动速度可以达到每秒几千米。相比之下,接收机的速度非常小,几乎可以忽略。因此,对于空天混合这种模式,杂波时间频率将主要受发射机径向速度的影响。双基地的几何关系随发射机和接收机的高度、飞行方向和夹角 而变化。下面着重研究由 变化引起的杂波特性的改变。这里假设#r 和#t 都为0,发射机和接收机的高度分别为850km 和10km ,速度分别为7.408km /s 和100m /s 。f r 为3kHz ,!为0.3m 。待测目标P 点位于接收机的正侧视方向且距离接收机200km 。为了比较杂波随距离变化的平稳性还选取了另外2个距离环进行比较,它们与P 点所在的待检测距离环相距50km

18、 。当 =2时,杂波等距离环的形状如图3a 所示,随着距离的增加,等距离环逐渐超出接收机圆的范围,其相应的杂波角度Do ppler 线如图3b 所示。在天基情况下,由于杂波谱宽远远超过了f r ,将发生严重的Doppler 模糊。原本随着距离的增加,杂波谱的范围增加,Do ppler 模糊应变得更加严重,但是由于当距离增加到一定值后,杂波环上只有部分回波存在,因此杂波谱的范围实际减小,Doppler 模糊缓和了。同时可以看出,对于不同的距离,杂波角度Doppler 线不重合,表明杂波数据的不均匀性。当 =15时,杂波等距离环的形状如图4a 所示,距离环都只有部分杂波信号存在。杂波角度Doppl

19、er 线如图4b 所示,对于不同的距离,杂波角度Doppler 线接近重合,表明杂波的不均匀性大大减弱。同时杂波谱宽也减小了,几乎没有Doppler 模糊。这是因为当发射机和接收机距离很远后,接收圆内的杂波单元相对于发射机的角度变化小,而杂波的Doppler 多由发射机飞行参数决定。这与机载双基地情况不同。机载时,杂波Doppler 由发射机和接收机共同决定,发射机和接收机离的越近,越倾向于单基地,杂波谱线间距越小,杂波数据越均匀。1608清华大学学报(自然科学版2007,47(10 图3 =2 杂波时间空间频率特性分析图4 =15杂波时间空间频率特性分析图5比较了2种 角度下的信干噪比损失。

20、当=15时的凹口更窄,最小可检测速度明显好于 =2。这是由于 增加,训练数据越均匀,使得杂波协方差阵的估计误差更小,ST AP 性能越好。因此,在进行慢速目标检测时,可以利用这一特点获得更好的检测性能。图5信干噪比损失5结论本文分析了当发射机放在卫星平台上,接收机采用无人机这种混合空天双基地工作模式下的杂波等距离环的特点,并提出了新的杂波模型的建立方法。在该模型基础上,本文还对几种特定空间几何配置下的空天混合双基地系统的杂波等距离环和角度Doppler 线进行了仿真,并分析了空天混合模式下杂波呈现的新特点。参考文献(References 1Hartnett M P.Ground and air

21、borne tar get detection w ith bistatic adap tive s pace based radar J .A er osp ace and E le ctronic S ystems M ag az ine ,1999,14(10:39-44.2Hartnett M P,Davis M E.Bis tatic su rveillance concept of operations C/Proceedings IEEE 2001Radar conference.U SA,2001:75-80.3王炎,吴曼青,陈学.空载双极地雷达地杂波模拟J .系统工程与电子技

22、术,2000,22(3:23-26.W ANG Yan ,WUM anqing,C HEN Xuequan.Clutter s imulationofair borneb istaticradarJ.S ystemsE ngineering and E le ctronic s ,2000,22(3:23-26.(in Ch ines e4Klemm R.Comp aris on b etw een m onostatic and bis tatic anten na configurations for ST AP J.IE EE Tr ans A erosp E le ctron Sy st ,2000,36(2:596-608.5ZHANG Yuhong,Him ed B.E ffects of geometry on clutter characteris ticsofbistatic radarsC/IEEERadarCon ference.US A,2003:417-424.6Him ed B.Effects of bistatic clutter dis pers ion on S TAP s ystem s J.I E E Pr oc Rad ar S onar N av ig ,2003,150(1:28-32.1609李华,等:空天混合双基地雷达杂波建模