海面多路径效应对被动微波测角的影响分析

海面多路径效应对被动微波测角的影响分析作者：樊鹏飞 欧阳中辉 来源：现代电子技术2014年第 24期摘 ;要： 在拦截超低空目标时，海面多路径效应会严重影响被动制导型舰空导弹的测角精 度。首先阐述海面多路径效应产生机理，然后结合旋转式相位干涉仪的测角原理，通过建立海 面多路径镜面反射几何模型，给出了双天线直达信号和镜面反射信号表达式。在此基础上，选 取典型超低空弹道，根据仿真结果代入多路径计算模型，分析超低空弹道多路径效应对测角精 度的影响。仿真结果表明，海面多路径效应对被动微波体制导引头测角影响显著，为下一步研 究抑制多路径影响措施提供了重要参考。关键词：海面多路径效应; 被动雷达导引头;镜面反射;超低空弹道中图分类号：TN97- 34 ;文献标识码：A ;文章编号： 1004- 373X (2014) 24- 0056- 05Simulation analysis on influence of sea surface multipath effect on PRS angle measurementFAN Peng- fei, OUYANG Zhong- hui( Department of Ordnance Science and Technology， Naval Aeronautical and Astronautical University, Yantai 264001 , China)Abstract： As attacking the ultralow- altitude targets, the sea surface multipath effect has great influence on passive radar seeker ( PRS) angle measuring precision of ship- to- air missile. The occurring principle of sea surface multipath effect is described. The specular reflection geometric model of the sea surface multipath effect is built in combination with the angle measuring principle of rotating phase interferometer. By this model, the expressions of direct signal and specular reflection signal of dual- atenna are given. On this basis, typical low- altitude ballistic trajectories are selected to analyze the influence of the multipath effect on the angle measuring accuracy according to the results of the simulation calculation in multipath effect model. The simulation results show that sea surface multipath effect affects PRS angle measuring accuracy significantly. The achievement provided an important reference for the next research on restraining multipath effect.Keywords： sea surface multipath effect; passive radar seeker; specular reflection; ultralow- altitude ballistic trajectory0 ;引 ;言电磁波在传播过程中经地面、海面或其他物体反射后，往往存在包括直线传播以外的多条 路径。在雷达接收机处，来自同一辐射源的电磁波经不同路径传播后形成的多径信号，会给雷 达对目标正常的截获和跟踪性能造成影响，这种现象被称为多路径效应1。末端舰空导弹的拦截目标多为超低空反舰导弹，跟踪雷达的工作状态往往是低仰角姿态， 通常会受到海面反射信号的干扰，因此被动雷达天线极易受多路径效应影响2。当反射表面 比较光滑时，这种误差更为强烈，严重时会使导弹跟踪到镜像目标，造成拦截失败。为了提高 舰空导弹的拦截效能，必须详细分析海面多路径效应对被动微波测角的影响，从而为进一步研 究抑制多路径效应措施提供参考。近年来，国内外已对多路径效应做了许多工作，研究集中于反射模型建立、半实物仿真、 三维地形建模、抑制措施等方面。文献3采用面元KA法及微扰法（SPA）建立了粗糙海面多 路径电磁散射模型;文献4提出了基于半实物仿真的多路径仿真系统工程实现的方法;文献5采 用随机分形插值算法对起伏地面和海面进行了三维地形建模;文献6将传统的多目标分辨算法 （C2 算法）与偏差补偿技术相结合应用于低角多路径环境下跟踪目标俯仰角的测量。总体上 看，现有海面多路径效应研究的不足之处主要包括：对舰载雷达、主动雷达研究较多，缺乏对 被动雷达导引头影响的研究;对被动雷达影响的研究大多未考虑相位干涉仪双天线的旋转;与导 弹飞行弹道的联系不够紧密等。基于此，本文结合旋转式相位干涉仪的测角原理，建立了海面多路径镜面反射模型，选取 典型超低空弹道，详细分析了超低空下海面多路径效应对测角精度的影响。1 ;海面多路径效应分析雷达侦察装备在复杂海面背景条件下接收到的多路径信号也是非常复杂的。按照反射表面 起伏不平的程度，可分为平滑表面和粗糙表面。在光滑平坦的表面，表面反射主要为镜面反射; 在粗糙不平的表面上，还会产生漫反射。根据之前的研究结果，海情较小时镜面反射下多路径 效应影响会较为严重;而海情较大时多路径反射比较杂乱，起伏较快，可通过取平均值的方法 减小偏差3。因此，本文着重研究镜面反射情况下海面多路径效应对被动微波测角的影响， 镜面反射多路径示意图如图 1所示。图1 镜面反射多路径示意图由于目标和镜像目标分别为两个矢量点源，两个矢量点源通过干涉叠加后进行矢量合成， 形成一个新的矢量点源，也就是两个目标的合成中心。由于镜像目标的影响，导引头测到的角 度既不是目标角度也不是镜像目标角度，而是这个新的合成矢量的角度，这个合成中心位置随 着弹目相对运动在目标与镜像的连线上改变，这就是目标的角闪烁现象，给导引头测角带来严 重影响。2 ;镜面反射模型2.1 ;旋转式相位干涉仪测角原理在被动雷达测角中相位干涉仪是较常用的一种角度测量方法，其基本原理是通过鉴别不同 天线接收到的平面电磁波信号的相位差，然后经过角度变换计算出目标辐射源的视线角7。被动微波子系统采用相位干涉仪测角面临的基本问题是难以同时满足不模糊测角范围和测 角精度的要求。旋转式双天线跟踪测角系统能够较好地解决这个问题，在保证较高测角精度、 较大不模糊测角范围的基础上，利用一对测向天线实现对目标的空间定向。旋转式相位干涉仪 实质上是采用时延跟踪环路形成电子角度跟踪系统，使角度信息转换为弹体旋转频率的交流幅 度与相位信息，其测角双天线运动模型如图2所示。有关旋转式相位干涉仪测角模型详见文献 8，在此不再赘述。图2 旋转式双天线运动模型2.2 ;镜面反射几何模型当海面起伏高度差满足瑞利判据9时，认为反射信号来自反射点附近的第一菲涅尔区， 此时主要产生镜面反射信号。镜面反射几何模型如图3所示，辐射源位于位置T,坐标为 (0，0，hT)，ala2为基线长为d啲天线，以旋转频率阴绕平行于y轴的水平线旋转， 旋转面在xOz面上的投影如图4所示，天线中心在yOz面内，高度为hM,且到y轴的投影 点坐标为(0，L，0)。以天线在xOz面上的投影与z轴重合时(al在上a2在下)为初始 时刻，经过时间t后，al， a2的坐标为(xal, yal， zal)和(xa2, ya2， za2)。假设在 整个过程中不考虑微波辐射目标的俯仰、方位和滚动，也不考虑接收导引头的俯仰变化。图 3 镜面反射几何模型图4旋转天线在xOz面上的投影以天线al为例，此时辐射源T与天线al的距离为R，a为直达信号与辐射源主波束的 夹角，I为镜面反射点，反射点到T，al的距离分别为Rl、R2, al、a2分别为反射信号 与T和al主波束的夹角，为反射信号与海平面的夹角。各个相应的参数求解如下：由旋转 天线在xOz面上投影的几何关系，天线al的坐标(xal, yal, zal)的求解公式为：xal=-d2sin (rt) yal=Lzal=hM+d2cos (rt) ; ; ; ; ; ; (l)辐射源T到天线al距离：R=xal2+yal2+(zal-hT) 2 ; ; ; ; ;(2)直达路径与辐射源主波束的夹角：a=arccosyalR; （3）镜面入射波与水平面的夹角即掠射角：0=arctanhT+zalxal2+yal2 ; ; ; ; ; （4）进一步，反射点I到辐射源T、天线的距离Rl、R2分别为：Rl=hMsinO ; ; ; ; ; ; ;（5）R2=za1sin0 ; ; ; ; ; ; ; （6）镜面反射点I的坐标（xI, yI，0）为：xI=Rlcos0 （xa12+ya12） xalyI=Rlcos0 （xa12+ya12） yal ; ; ; ; ; （7） 镜面反射波与辐射源主波束的夹角：al=arccosyIRl ; ; ; ; ; ; ;（8）镜面入射波与天线a1 主波束的夹角：a2=arccosyal-yIR2 ; ; ; ; ; （9）直达波的延时td、反射波的延时tr为：td=Rc ; ; ; ; ; ; ; ;（10）tr=R1c+R2c ; ; ; ; ; ;（11）式中c为光速。3 ; 天线接收信号模型3.1 ;辐射源发射信号模型雷达发射信号为：St （t） =AtGtft（0）exp （jct）v（t） ; （ 12）式中：At为发射信号的幅度;Gt为发射天线的电压增益;ft (0)为发射天线方向图函 数;c为载频;v(t)为调制函数，是Np个宽度为Tp的矩形脉冲构成的脉冲串。若不考虑 脉间捷变频和线性调频，则：St (t) =AtGtft(0)exp (ject) k=ONp-lrectt-kTrTp ; (13)3.2 ;天线接收信号模型参数定义如图1所示，Sd为直达波，Si为入射波，Sr为反射波，镜面反射满足入射角等于反射 角，即如=屮，结合电磁理论易得导引头处总的接收信号为10：S=Sd+Sr=Atf (0t) +ArpDf (0r) ; ; (14)式中：At为直达信号的幅度;Ar为反射信号的幅度;f (0t)为直达波方向图;f (0r)为 反射波方向图；D为扩散因子，考虑到地球曲率的影响，反射波照到凸起的地球表面会引起扩 散，使得电磁波能量密度衰减；p 为表面反射系数，与反射表面的散射特性、入射波的入射角 屮小雷达的工作波长和极化方式有关，下面对p的求解进行讨论。对于理想的光滑平坦表面，反射系数为菲涅尔反射系数p0,即p=pO。可以利用极化形 式、入射余角n 和雷达工作频率通过菲涅尔方程计算得到，不同极化形式下的菲涅尔反射系 数p0 表达式如下所示。对于垂直极化有：pO=sin屮i-cos2isin屮i+-cos2屮i ; ; ; (15)对于水平极化有：pO=sin屮i-cos2屮isin屮i+-cos2屮i ; ; ; ; (16)对于圆极化有：p0 是垂直极化与水平极化时的中值。式(15)和式(16)中，是复介电常数，由下式确定：=K0-jGW0=，-j，-j60XG ； (17)式中：0是自由空间的介电常数；K是相容率，K0是反射面的相对介电常数;g是电导 率。文献9给出了一些典型海面、地表的回、和叫的数值。对于有一定粗糙度的反射面，认为其属于相对平坦面，即在第一菲涅尔反射区内，表面高 度变化Ah满足瑞利判据，这时p=pOps,其中ps为镜面散射因子，是表征反射面的粗糙度对 镜面反射幅度衰减影响的参数。镜面散射因子通常用其均方根值(Root Mean Square，RMS) 表示，它与反射面粗糙度因子r 的关系为：ps=exp-2 (2nD 2，0r0.1 ; ; ; (18)r=ohsin屮认;(19)式中：oh为表面起伏高度的标准差。因此，对于海面有一定粗糙度的情况，应用镜面反射模型，天线总的接收信号为：S=Sd+Sr=Atf (Ot)T+ArpOpsf (Or) ; ; (20)3.3 ;天线接收信号模型求解过程对于起伏较小的海面，考虑只有镜面反射时，天线接收的信号包括直达信号、镜面反射信 号。假设辐射和接收天线方向图ft (Ot)，fr (Or)均为高斯函数，将在镜面反射几何模型 和接收信号模型中求得的参数代入公式，即得t时刻天线a1 接收到的直达信号、反射信号。(1) t时刻天线a1 接收到的直达信号Sd (t)为：将所求得的a, td代入式(12)，即得直达波信Sd (t)： Sd (t)=AtGtft (a)exp (jec (t-td)v (t-td)fr (a)Gr ; ; ; =AtGtft (a)exp (jec (t-td) ; ; ; ; ;k=0Np-1rect t- kTr-tdTpfr (a)Gr式中：At, wc和旦(t)的定义与前面表述一致，ft (a)， fr (a)为直达波在辐射 源T和接收天线a1处的方向图。(2) t时刻天线a1 接收到的海面反射信号Sr (t)为：将所求得的a1, a2，O和tr代入式(15)和式(13)，即得反射波信号Sr (t)：Sr(t) =AtGtft(a1) exp(jec(t-tr) v(t-tr) p0(Otr) fr(a2) Gr ; ; ; =AtGtft(a1) exp(jec(t-tr) ; ; ; ; ;k=0Np-1rectt-kTr-trTpp0(O) fr(a2) Gr式中：ft (a1) ，fr (a2)为反射波在T和a1处的方向图，pO (O)为掠射角为O时 的海面的反射系数；td, tr分别为直达波和反射波到达天线的延时;Gr为接收天线的电压增 益。(3) t时刻天线al总的接收信号为：Sa1(t) =Sd1(t) +Sr1(t) ; ; ; ; ;(23)式中：Sdl (t)为直达信号;Sr1 (t)为反射信号。同理可以求出t时刻天线a2收到的直达波信号为Sd2 (t)、反射波信号为Sr2 (t), 天线a2总的接收信号为Sa2 (t)。将所求得的Sal, Sa2经过比相、解模糊后，可得目标 角度的仿真结果。4 ; 舰空导弹运动模型导弹运动的动态过程通常是由导弹的运动方程组积分得到的，为方便起见，在讨论导弹的 运动学时可将导弹当成质点来考虑，重量不变，在航迹坐标系中给出运动学模型。为建立相应 的运动学方程组，可将导弹质心的速度投影到地面坐标系中。本文的建模和仿真中，均采用大 地坐标系，并以舰空导弹发射时的初始位置作为坐标原点。反舰导弹可以采用匀速飞行也可以 采用变速飞行模式，这通常与导弹的机动策略有关。为了简化分析，假设反舰导弹作等速直线 运动。设反舰导弹当前坐标为( xt， k， yt， k， zt， k) ，舰空导弹当前坐标为( xm， k， ym, k，zm, k)，飞行速度为Vm,弹道角为申k,速度矢量倾角为Ok,则舰空导弹在下 一个时刻 k+l 时的位置为：xm, k+1=xm, k+vm cos0k- cospk- Atym, k+1=ym, k+vm sin0k- Atzm, k+1=zm, k+vm cos0k- sing At ; ; ; (24)式中At为仿真时间增量。以 xOy 平面为例,仿真时刻 k+1 时的弹目视线角为qxy( k+1) =arctanyt, k+1-ym, k+1xt, k+1-xm, k+1 ; ; ; ;( 25)仿真时刻k+1的舰空导弹速度矢量转动角速度为Nxy (k+1)。本文以比例导引法为例说 明导引规律模型。比例导引法是指导弹在攻击目标的导引过程中,导弹速度矢量的旋转角速度 与目标线的旋转角速度成比例的一种导引方法11,其导引关系方程为：屮k=Kqk ; ; ; ; ; ; (26)式中K为比例系数，一般取26。仿真过程中，屮k+1, qk+1可由下式计算：qk+1=qk+1-qkAt屮k+l=k+l-屮kAt ; ; ; ; ; (27)式中：qk,屮k分别为k时刻的弹目视线角、导弹速度与基准线夹角;qk+1,屮k+1分 别为下一时刻k+1时的弹目视线角、导弹速度与基准线夹角，At为跟踪仿真间隔。5 ; 海面多路径效应仿真5.1 ;仿真参数设置对典型海面超低空目标多路径效应进行仿真,仿真初始参数设置情况如下：（1） 舰空导弹。舰空导弹仿真初始时刻的坐标为（0, 0，0）,发射倾角为gm=15。,航 向角为屮m=15,飞行速度vm=600 ;m/s,被动雷达制导时的比例导引系数为4。（2） 目标特性。反舰导弹初始坐标为（10 000, 10, 200），飞行速度vt=300 ;m/s,在 水平面上作等速直线飞行，航向角屮t=180,雷达辐射波长0.03 m,辐射功率30 W,垂直极 化，发射天线增益20 dB,反舰导弹主动雷达天线在指向舰空导弹方向上的幅度与在反射点方 向上的幅度相等。（3）环境设定。在末端舰空导弹拦截反舰导弹过程中,辐射源与接收天线的距离比较 短，可以不考虑扩散因子的影响，D-1。海水介电常数k=65-652.1i。如果海面起伏均方差满足瑞利条件,则认为是光滑的。与雷达波长相比,带有小的毛细波 的水面都被认为是光滑的,因此起伏较小海面的反射系数即为菲涅尔反射系数。对于采用垂直极化的辐射源,综合式（15）和式（18）得到光滑海面的反射系数为：p=p0ps=exp-22nohsin呱2p0（0） ; （24）式中：0为入射余角（掠射角）为海表的复介电常数;对于起伏较小的海面oh可近似 为0,因此，海面的反射系数仍为p=p0 （0）。5.2 ;仿真结果及分析选取典型超低空仿真弹道,根据仿真计算结果代入多路径计算模型,分析超低空弹道的多 路径效应对测角精度的影响,仿真结果如图5,图6所示,其中图5为舰空导弹拦截目标三维 弹道仿真,图6为舰空导弹被动雷达受海面多路径影响出现的测角误差。图5舰空导弹拦截目标三维弹道仿真图6多路径效应对被动雷达测角影响从仿真计算结果可以得出以下结论：(1) 海面多路径效应对被动微波接收体制导引头有影响，导引头受多路径效应影响测角 出现了明显偏差;(2) 对于超低空目标，在弹道中段由多路径引起的误差很小，误差都在1以内;在弹道 末段1 km左右区域由多路径引起的误差增长至2左右，被动雷达测角精度受到较大影响，导 弹脱靶距离较大。(3) 多路径对测角精度的影响主要与弹目距离和弹目相对海面高度有关，目标和导引头 相对高度如果变化，目标/镜像合成中心将围绕着真实目标做上下运动变化;(4) 多路径效应对测角的影响与目标飞行航路无明显关系。6 ;结 ; 语被动微波导引头的制导精度取决于测角精度，但由于海面多路径效应的存在使得测角精度 大大降低。本文从舰空导弹攻防对抗机理出发，系统阐述了海面多路径的产生原理，建立了结 合旋转相位干涉仪测角原理的海面多路径镜面反射几何模型，最后选取典型超低空弹道进行仿 真，详细分析了多路径对被动微波导引头测角的影响。为下一步研究抑制多路径影响措施、多 模复合交班、被动微波导引头滤波器设计打下基础。参考文献1 任子西，多路径效应对反辐射导弹被动雷达导引头性能的影响分析J.战术导弹技术，2009 (3)： 1- 5.2 陈鑫，王浩丞，唐勇，等多径环境中被动雷达导引头测向性能分析J.电子信息对抗技 术， 2011， 26(4)： 1- 4.3 王晓冰，童广德，梁子长.海面多路径下比相测角误差仿真分析J.制导与引信，2009， 30(2)： 11- 15.4 袁兴鹏，张金全，杨立永海面雷达信号环境多路径效应建模与仿真技术研究J.舰船电 子对抗， 2013， 36(6)： 68- 72.5 段世忠，周荫清，张孟，等.主动雷达导引头多路径效应的数字仿真J.北京航空航天大 学学报， 2002， 28( 4)： 447- 450.6 吕韶昱，万建伟，辛勤海面多径环境下雷达目标俯仰角测量提取的研究与应用J国防 科技大学学报， 2007， 29( 5)： 48- 53.7 吴宝东，陈舒基于相位干涉仪测向系统的相位误差分析J.舰船电子对抗，2008, 31 (3)： 74- 76.8 沈康一种旋转式相位干涉仪测角系统研究J.现代电子技术，2011，34 (15)： 26- 28.9 M W朗陆地和海面的雷达波散射特性M.薛德镛，译北京：科学出版社，1981.10 刘亚平多径传播对反反舰导弹相位干涉仪测角影响的建模与仿真D长沙：国防科学 技术大学， 2007.11 钱杏芳，林瑞雄，赵亚男导弹飞行力学M.北京：北京理工大学出版社，2006.