雷达大作业-振幅和差角度测量及仿真

1、雷达原理大作业单脉冲自动测角的原理及应用学院：电子工程学院作者：2016年5月21日单脉冲自动测角的原理及应用 单脉冲自动测角的原理及应用一摘要单脉冲测角法是属于振幅法测角中的等信号法中的一种，其测角精度高，抗干扰能力强，在现实中得到了广泛的应用。而其中对于接收支路要求不太严格的双平面振幅和差式单脉冲雷达，更是备受青睐。本文首先讲述了单平面振幅和差式单脉冲雷达自动测角的原理，再简述了双平面振幅和差式单脉冲雷达自动测角的结构框图，接着简述了本文仿真所用的一些原理和公式推导，包括天线方向图函数及其导数的推导，最后做了基于高斯形天线方向图函数的单脉冲自动测角，基于辛克函数形天线方向图函数的单脉冲自动

2、测角，和基于高斯形天线方向图函数的双平面单脉冲自动测角。源代码在附录里。二重要的符号说明符号含义误差角00.5天线方向图函数半功率波瓣宽度，也就是归一化幅度值为0.707时对应的角度值00天线方向图函数零功率波瓣宽度5单平面内，两波束相对天线轴线的偏角0P双平面内，子波束在方位面相对天线轴线的偏角0P双平面内，子波束在俯仰面相对天线轴线的偏角三单平面振幅和差式单脉冲自动测角原理单脉冲测角法是属于振幅法测角中的等信号法中的一种。在单平面内，两个相同的波束部分重叠，交叠方向即为等信号轴的方向。将这两个波束接收到的回波信号进行比较就可以在一定范围内，一定精度要求下测到目标的所在角度。因为两个波束同时

3、接到回波，故单脉冲测角获得目标角误差信息的时间可以很短，理论上只要分析一个回波脉冲即可，所以称之为“单脉冲”。因取出角误差的具体方式不同，单脉冲雷达种类很多，其中应用最广的是振幅和差式单脉冲雷达，其基本原理说明如下：角误差信号两馈源形成的波束(b)和波束(c)差波束振幅和差式单脉冲雷达取得角误差信号基本方法是将这两个波束同时收到的信号进行和差处理，分别得到和信号和差信号。其中差信号即为该角平面内角误差信号。若目标处在天线轴方向(等信号轴)，误差角=0，则两波束收到的回波信号振幅相同，差信号等于0。目标偏离等信号轴而有一个误差角&时，差信号输出振幅与&成正比而其符号则由偏离方向决定。和差比较器这

4、里主要使用双T插头，示意图如下图2(a)所示。它有四个端口：和端，差端和1,2端。假定四个端都是匹配的，则从和端输出信号时，1,2端输出等幅同相的信号，差端无输出；从1,2端输入同相信号时，和端输出两信号之和，差端输出两信号之差。图2双T接头和差比较器示意图双T接头(b)和差比较器示意图在发射信号时，从发射机来的信号加在和端，故1,2端输出等幅同相的信号，两波束在空间各点产生的场强同相相加，形成发射和波束的天线方向性函数为F)。接收时，回波脉冲同时加到1,2端，此时在和端，输出两个回波信号同相相加之和，记为Ee；在差端，输出两信号反相相加之和，记为Ea。假设两个波束方向性函数完全相同，记为F(

5、0)，两波束衰减倍数为k,两波束相对天线轴线的偏角为8，则对于0方向的目标来说：和信号振幅为：E=kF(0)F(8-0)+kF(0)F(8+0)=kF2(0)EA差信号振幅为：E=kF(0)F(8-0)-kF(0)F(8+0)=kF(0)F(0)其中：F(0)二F(8-0)+F(8+0)，F(0)二F(8-0)F(8+0)。EA实际情况下，0是很小的，可以对F(8-0)和F(8+0)在8附近做一阶泰勒展开:F(6-0)二F(6)-F)0+o(02)沁F(6)-F)0F(6+0)二F(6)+F)0+o(02)二F(6)+F)0则.E=kF2(0)沁kF(6)+F(6)0+F(6)-F(6)02=

6、4kF2(6)工工E=kF(0)F(0)沁kF(6)+F(6)0+F(6)F(6)0(F(6)F(6)0)-(F(6)+F(6)0)A2A=4kF(6)F(6)0E则：E2据此,-4kF(6)F(6)04kF2(6)F(6)0n0硕0占其中耳是个常数，耳=F(6)F6)1EAnE2可以求得0=，这就是单平面振幅和差式单脉冲雷达测角的公式。四双平面振幅和差式单脉冲自动测角简述为了对空中目标进行自动方向跟踪，必须在方位角和俯仰角两个角平面上进行角跟踪,因而必须获得方位角和俯仰角的误差信号。为此，需要用4个馈源照射一个反射体，以形成四个对称的相互部分重叠的波束。双平面上天线四个子波束在三维空间中的辐

7、射图图3雷达天线空间子波束分布图假设其中一个子波束的三维方向图函数为f(0,e)，其中0代表方位角，代表俯仰角，0和分别为子波束在方位面和俯仰面相对中心轴的偏离角。pp则图3所示的4个子波束的数学表达式为：子波束A：F(0,)二F(0-0,-)App子波束B：F(0,)二F(0+0,-)Bpp子波束C：F(0,)二F(0+0,+)Cpp子波束D：F(0,)二F(0-0,+)Dpp*A+B+C+Dyf(A+M-伯+方宜角误差暗号双平面振幅和差式单脉冲自动测角的原理框图如下图4所示：图4.双平面振幅和差式单脉冲自动测角的原理框图(AB)-(C+I俯忡曲燥建信号AB1Dc其中和信号为F(0,)=F(

8、0,)+F(0,)+F(0,)+F(0,)SABCD俯仰角误差信号为：Fa(0，)=Fa(0，)+Fb(0，)-FC(0，)-Fd(0，)方位角误差信号为：5(0，)=Fa(0，)+Fd(0，)-Fb(0，)-FC(0，)再根据之前的单平面振幅和差式单脉冲雷达测角的公式有：1EaES其中，耳=af(0,)f(0,)app0=0,=pp五仿真所用原理和公式推导1.单平面单脉冲雷达方向图函数A.高斯形方向图函数单脉冲自动测角的原理及应用92 其方向图函数为：F(9)=e-1.4話2.8乂-1.4仝其导函数为F(8)=盯8e9-oB.辛克函数形方向图函数其方向图函数为：0.5F(9)二Sa(2兀9)

9、,9罟8cos(|8)-sin(誥8)其导函数为F(8)=02k82902.双平面脉冲雷达方向图函数此处对三维天线使用简化模型，将三维天线看做由两个二维平面(方位面和俯仰面)方向图相乘的结果，即F(9,)二F(9)F(帕。TOC o 1-5 h z HYPERLINK l bookmark78 o Current Document 9由于仿真中只以高斯形方向图函数为例做了仿真，所以此处只讨论三维的高斯形方向图函数，设9和分别为子波束在方位面和俯仰面相对中心轴的偏离角：pp9221.41.4F(9,)=e9蟲*e蟲其偏导函数为：QF(9,)399=9,=pp2.8921.4e92920.51.4

10、e02.50.53F(9,)9=9,=pp2.820.51.4le025*p1.4e92T0.5关于两波束相对天线轴线的偏角的选取：两波束相对天线轴线的偏角8一般选择为单脉冲自动测角的原理及应用 v是半功率波瓣宽度，也就是归一化幅度为0.707时候对应的0值20.50.5在单平面上，直接取;20.5在双平面上，可以取0=0=0。pp20.5六仿真结果1.基于高斯形天线方向图函数的单脉冲自动测角下图是两馈源形成的波束，和波束和差波束波形图0.60.&0.40.2-323-101角度值（rad0角搜值（rad0F（悶）尸阳）波形下图是单平面上目标所在角度和测得角度的比较曲线：诣斯力向函数的白动测加

11、黠刃杠所在方时角嵐捎比较TJBMB.腿fcEK足四整吸邛更经二2.基于辛克函数形天线方向图函数的单脉冲自动测角下图是两馈源形成的波束，和波束和差波束波形图：F（斟E+FW波形20.86；：E。斗0.20-0.2-3-1Q1ffiWL(rad)角度ri(rad)fl下图是单平面上目标所在角度和测得角度的比较曲线：辛克方向雷数的口动测仙结呆和H标所在方向仙度值比较(P皿)切昵嘅欺sss矍嘅卑主堕匸AO.-0.3-0.2-0.100.1O.20.3角度值（旧3.基于高斯形天线方向图函数的双平面单脉冲自动测角下图是双平面上天线四个子波束的方向图函数图像：230-2四个子波束的天线方向图函数俯仰ft(r

12、ad)_2方位ft(rad)四个子披束的天线方向图函数下图是三维和函数的图像：和團数的-维图像下图是方位角差函数图像：2200-20.5方位角签團数的二维图像-05、谢仰角（rad揷-2方位ft(rad)i?下图是俯仰角差函数图像：2200-20.5俯仰角并團数的二维图似-05、谢仰角（rad揷-2方位ft(rad)i?下图是目标方位角与俯仰角之积的三维曲面和测得方位角与俯仰角之积的三维曲面：o4)20.60.20.60.400.2-020-0.1汀标方位需切俯仰角之积耳测得的方仗角门傭仰角之觀的辨團傻俯仰角(rad艸方讥角(rad)0-0.4-0.64总下图是目标方位角与测得的方位角的比较：

13、目标的方位角和测得的方位角比较-目标的方位角(rad)fl测得的方位角(前)禺(pg1心虺翼慝嘅0股呂炬匸-0.3-02-0.10010.203方位角(怡日下图是目标俯仰角与测得的俯仰角的比较：冃标的俯仰角和测得的俯仰角比较-目标的俯仰Hj(rad)测得的嵋仰角“胡诂-0,3-0200.10,203俯抑角(rad)七附录源码ld.m(单平面的仿真代码)：clear;closeall;clc%高斯函数的方向图函数绘图和自动测角%parameters%theta_r=0.4;%theta_p=0.2;%k=1;%operations%theta=(-pi):(2*pi/1000):pi;%f1=e

14、xp(-1.4*(theta-theta_p)42)/theta_rA2);%f2=exp(-1.4*(theta+theta_p).A2)/theta_rA2);%sigma=f1+f2;delta=f1-f2;%figure(1);%subplot(3,1,1);plot(theta,f1,theta,f2,-,Linewidth,2);axistight;%xlabel(角度值(rad)theta);ylabel(幅度);title(F(theta+delta)+F(thetaAdelta)波形);%holdon;plot(theta,0.7*ones(size(theta),k);te

15、xt(1,0.8,归一化幅度=0.707);%subplot(3,1,2);plot(theta,sigma,Linewidth,2);axistight;%xlabel(角度值(rad)theta);ylabel(幅度);title(F(theta-delta)+F(theta+delta)波形)；%subplot(3,1,3);plot(theta,delta,Linewidth,2);axistight;%xlabel(角度值(rad)theta);ylabel(幅度);title(F(theta-delta)-F(theta+delta)波形);%target=(-pi/8):(2*p

16、i/1000):(pi/8);%result=zeros(size(target);%Fsum=exp(-1.4*(target-theta_p*ones(size(target).A2)/theta_rA2)+exp(-1.4*(target+theta_p*ones(size(target).A2)/theta_rA2);%Fdelta=exp(-1.4*(target-theta_p*ones(size(target).A2)/theta_rA2exp(-1.4*(target+theta_p*ones(size(target).A2)/theta_rA2);%Esum=k*Fsum.A

17、2;%Edelta=k.*Fsum.*Fdelta;%F_ori=exp(-1.4*(theta_pA2)/(theta_rA2);%F_dao=-2.8*theta_p*exp(-1.4*(theta_p.A2)/theta_rA2)/(theta_rA2);%res=-(Edelta*F_ori)./(Esum*F_dao);%figure;plot(target,target,target,res,-r,Linewidth,2);axistight;legend(目标所在角度,测角所得角度);%xlabel(角度值(rad);ylabel(目标所在角度值和测角所得角度值(rad);%ti

18、tle(高斯方向函数的自动测角结果和目标所在方向角度值比较)；%辛克函数的方向图函数绘图和自动测角%parameterstheta_r=0.8;theta_p=0.18;k=1;%operationstheta=(-pi):(2*pi/1000):pi;f1=sin(2*pi*(theta-theta_p)/theta_r)./(2*pi*(theta-theta_p)/theta_r);f2=sin(2*pi*(theta+theta_p)/theta_r)./(2*pi*(theta+theta_p)/theta_r);sigma=f1+f2;delta=f1-f2;figure(1);s

19、ubplot(3,1,1);plot(theta,f1,theta,f2,-,Linewidth,2);axistight;holdon;plot(theta,0.7*ones(size(theta),k);text(1,0.8，归一化幅度=0.707);xlabel(角度值(rad)theta);ylabel(幅度);title(F(theta+delta)+F(theta-delta)波形);subplot(3,1,2);plot(theta,sigma,Linewidth,2);axistight;xlabel(角度值(rad)theta);ylabel(幅度);title(F(thet

20、a-delta)+F(theta+delta)波形)；subplot(3,1,3);plot(theta,delta,Linewidth,2);axistight;xlabel(角度值(rad)theta);ylabel(幅度);title(F(theta-delta)-F(theta+delta)波形);target=(-pi/8):(2*pi/1000):(pi/8);result=zeros(size(target);Fsum=sin(2*pi*(target-theta_p*ones(size(target)/theta_r)./(2*pi*(target-theta_p*ones(s

21、ize(target)/theta_r)+sin(2*pi*(target+theta_p*ones(size(target)/theta_r)./(2*pi*(target+theta_p*ones(size(target)/theta_r);Fdelta=sin(2*pi*(target-theta_p*ones(size(target)/theta_r)./(2*pi*(target-theta_p*ones(size(target)/theta_r)-sin(2*pi*(target+theta_p*ones(size(target)/theta_r)./(2*pi*(target+t

22、heta_p*ones(size(target)/theta_r);Esum=k*Fsum42;Edelta=k.*Fsum.*Fdelta;F_ori=sin(2*pi*theta_p/theta_r)/(2*pi*theta_p/theta_r);F_dao=(2*pi*theta_p*cos(2*pi*theta_p/theta_r)/theta_r-sin(2*pi*theta_p/theta_r)/(2*pi*(theta_pA2)/theta_r);res=-(Edelta*F_ori)./(Esum*F_dao);figure;plot(target,target,target,

23、res,-r,Linewidth,2);axistight;legend(目标所在角度,测角所得角度);xlabel(角度值(rad);ylabel(目标所在角度值和测角所得角度值(rad);title(辛克方向函数的自动测角结果和目标所在方向角度值比较)；ld2.m(双平面的仿真代码)：clear;closeall;%parameterstr=0.4;tp=0.2;k=1;%高斯方向图函数三维图像绘制%setax=(-pi):(2*pi/1000):pi;%setay=(-pi):(2*pi/1000):(pi);%f1=(exp(-1.4*(setax-tp*ones(size(setax

24、).A2)/trA2)*(exp(-1.4*(setay-tp*ones(size(setay).A2)/trA2);%A%f2=(exp(-1.4*(setax+tp*ones(size(setax).A2)/trA2)*(exp(-1.4*(setay-tp*ones(size(setay).A2)/trA2);%B%f3=(exp(-1.4*(setax+tp*ones(size(setax).A2)/trA2)*(exp(-1.4*(setay+tp*ones(size(setay).A2)/trA2);%C%f4=(exp(-1.4*(setax-tp*ones(size(setax

25、).A2)/trA2)*(exp(-1.4*(setay+tp*ones(size(setay).A2)/trA2);%D%figure;%meshc(setax,setay,f1);%holdon;meshc(setax,setay,f2);%holdon;meshc(setax,setay,f3);%holdon;meshc(setax,setay,f4);%xlabel(方位角(rad)theta);ylabel(俯仰角(rad)phi);zlabel(幅度);title(四个子波束的天线方向图函数);%Fsum=f1+f2+f3+f4;%figure;meshc(setax,setay

26、,Fsum);xlabel(方位角(rad)theta);ylabel(俯仰角(rad)phi);%zlabel(幅度);title(和函数的三维图像);%Fdelst=f1+f4-f2-f3;%figure;meshc(setax,setay,Fdelst);xlabel(方位角(rad)theta);ylabel(俯仰角(rad)phi);%zlabel(幅度);title(方位角差函数的三维图像)；%Fdelph=f1+f2-f3-f4;%figure;meshc(setax,setay,Fdelph);xlabel(方位角(rad)theta);ylabel(俯仰角(rad)phi);

27、%zlabel(幅度);title(俯仰角差函数的三维图像)；%目标自动测角targetx=(-pi/8):(pi/1000):(pi/8);targety=(-pi/8):(pi/1000):(pi/8);f1=(exp(-1.4*(targetx-tp*ones(size(targetx).A2)/trA2)*(exp(-1.4*(targety-tp*ones(size(targety).A2)/trA2);%Af2=(exp(-1.4*(targetx+tp*ones(size(targetx).A2)/trA2)*(exp(-1.4*(targety-tp*ones(size(targety).A2)/trA2);%Bf3=(exp(-1.4*(targetx+tp*ones(size(targetx).A2)/trA2)*(exp(-1.4*(targety+tp*ones(size(targety).A2)/trA2);%Cf4=(exp(-1.4*(