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

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

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

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

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

5、方向性函数完全相同，记为，两波束衰减倍数为k，两波束相对天线轴线的偏角为，则对于方向的目标来说：和信号振幅为：差信号振幅为：其中：，。实际情况下，是很小的，可以对和在附近做一阶泰勒展开：则：则：， 其中是个常数，。据此，可以求得，这就是单平面振幅和差式单脉冲雷达测角的公式。四双平面振幅和差式单脉冲自动测角简述为了对空中目标进行自动方向跟踪，必须在方位角和俯仰角两个角平面上进行角跟踪，因而必须获得方位角和俯仰角的误差信号。为此，需要用4个馈源照射一个反射体，以形成四个对称的相互部分重叠的波束。双平面上天线四个子波束在三维空间中的辐射图 图3.雷达天线空间子波束分布图假设其中一个子波束的三维方向图

6、函数为，其中代表方位角，代表俯仰角，和分别为子波束在方位面和俯仰面相对中心轴的偏离角。则图3所示的4个子波束的数学表达式为：子波束A：子波束B：子波束C：子波束D：双平面振幅和差式单脉冲自动测角的原理框图如下图4所示：图4. 双平面振幅和差式单脉冲自动测角的原理框图其中和信号为俯仰角误差信号为：方位角误差信号为：再根据之前的单平面振幅和差式单脉冲雷达测角的公式有：其中，。五仿真所用原理和公式推导1.单平面单脉冲雷达方向图函数A高斯形方向图函数其方向图函数为：其导函数为。B辛克函数形方向图函数其方向图函数为：， 其导函数为2.双平面脉冲雷达方向图函数此处对三维天线使用简化模型，将三维天线看做由两

7、个二维平面（方位面和俯仰面）方向图相乘的结果，即。由于仿真中只以高斯形方向图函数为例做了仿真，所以此处只讨论三维的高斯形方向图函数，设和分别为子波束在方位面和俯仰面相对中心轴的偏离角：其偏导函数为：，。3.关于两波束相对天线轴线的偏角的选取：两波束相对天线轴线的偏角一般选择为： 是半功率波瓣宽度，也就是归一化幅度为0.707时候对应的值在单平面上，直接取；在双平面上，可以取。六仿真结果1. 基于高斯形天线方向图函数的单脉冲自动测角下图是两馈源形成的波束，和波束和差波束波形图 下图是单平面上目标所在角度和测得角度的比较曲线：2.基于辛克函数形天线方向图函数的单脉冲自动测角下图是两馈源形成的波束，

8、和波束和差波束波形图：下图是单平面上目标所在角度和测得角度的比较曲线：3. 基于高斯形天线方向图函数的双平面单脉冲自动测角下图是双平面上天线四个子波束的方向图函数图像：下图是三维和函数的图像：下图是方位角差函数图像：下图是俯仰角差函数图像：下图是目标方位角与俯仰角之积的三维曲面和测得方位角与俯仰角之积的三维曲面：下图是目标方位角与测得的方位角的比较：下图是目标俯仰角与测得的俯仰角的比较：七附录源码ld.m（单平面的仿真代码）：=clear; close all; clc% % 高斯函数的方向图函数绘图和自动测角% % parameters% theta_r=0.4; % theta_p=0.2

9、; % k=1;% % operations% theta=(-pi):(2*pi/1000):pi;% f1=exp(-1.4*(theta-theta_p).2)/theta_r2);% f2=exp(-1.4*(theta+theta_p).2)/theta_r2);% sigma=f1+f2; delta=f1-f2;% figure(1);% subplot(3,1,1); plot(theta,f1,theta,f2,'-','Linewidth',2); axis tight;% xlabel('角度值(rad)theta'); y

10、label('幅度'); title('F(theta+delta)+F(theta-delta)波形');% hold on; plot(theta,0.7*ones(size(theta),'k'); text(1,0.8,'归一化幅度=0.707');% subplot(3,1,2); plot(theta,sigma,'Linewidth',2); axis tight;% xlabel('角度值(rad)theta'); ylabel('幅度'); title('F

11、(theta-delta)+F(theta+delta)波形');% subplot(3,1,3); plot(theta,delta,'Linewidth',2); axis tight;% xlabel('角度值(rad)theta'); ylabel('幅度'); title('F(theta-delta)-F(theta+delta)波形');% target=(-pi/8):(2*pi/1000):(pi/8);% result=zeros(size(target);% Fsum=exp(-1.4*(target

12、-theta_p*ones(size(target).2)/theta_r2)+exp(-1.4*(target+theta_p*ones(size(target).2)/theta_r2);% Fdelta=exp(-1.4*(target-theta_p*ones(size(target).2)/theta_r2)-exp(-1.4*(target+theta_p*ones(size(target).2)/theta_r2);% Esum=k*Fsum.2;% Edelta=k.*Fsum.*Fdelta;% F_ori=exp(-1.4*(theta_p2)/(theta_r2);% F

13、_dao=-2.8*theta_p*exp(-1.4*(theta_p.2)/theta_r2)/(theta_r2);% res=-(Edelta*F_ori)./(Esum*F_dao);% figure; plot(target,target,target,res,'-r','Linewidth',2); axis tight; legend('目标所在角度','测角所得角度'); % xlabel('角度值(rad)'); ylabel('目标所在角度值和测角所得角度值(rad)');% t

14、itle('高斯方向函数的自动测角结果和目标所在方向角度值比较');% 辛克函数的方向图函数绘图和自动测角% 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=

15、f1-f2;figure(1);subplot(3,1,1); plot(theta,f1,theta,f2,'-','Linewidth',2); axis tight;hold on; 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)波形');subpl

16、ot(3,1,2); plot(theta,sigma,'Linewidth',2); axis tight;xlabel('角度值(rad)theta'); ylabel('幅度'); title('F(theta-delta)+F(theta+delta)波形');subplot(3,1,3); plot(theta,delta,'Linewidth',2); axis tight;xlabel('角度值(rad)theta'); ylabel('幅度'); title('

17、;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(size(target)/theta_r)+sin(2*pi*(target+theta_p*ones(size(target)/theta_r)./(2*pi*(target+theta_p*ones(size(target)/t

18、heta_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+theta_p*ones(size(target)/theta_r);Esum=k*Fsum.2; Edelta=k.*Fsum.*Fdelta;F_ori=sin(2*pi*theta_p/theta_r)/(2*pi*theta

19、_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_p2)/theta_r);res=-(Edelta*F_ori)./(Esum*F_dao);figure; plot(target,target,target,res,'-r','Linewidth',2); axis tight; legend('目标所在角度','测角所得角度'); xlabel('角度值(ra

20、d)'); ylabel('目标所在角度值和测角所得角度值(rad)');title('辛克方向函数的自动测角结果和目标所在方向角度值比较');=ld2.m（双平面的仿真代码）：=clear; close all;% 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) ).2 )/tr2 )'*(

21、exp( -1.4*( ( setay-tp*ones(size(setay) ).2 )/tr2 ); %A% f2=(exp( -1.4*( ( setax+tp*ones(size(setax) ).2 )/tr2 )'*(exp( -1.4*( ( setay-tp*ones(size(setay) ).2 )/tr2 ); %B% f3=(exp( -1.4*( ( setax+tp*ones(size(setax) ).2 )/tr2 )'*(exp( -1.4*( ( setay+tp*ones(size(setay) ).2 )/tr2 ); %C% f4=(e

22、xp( -1.4*( ( setax-tp*ones(size(setax) ).2 )/tr2 )'*(exp( -1.4*( ( setay+tp*ones(size(setay) ).2 )/tr2 ); %D% figure;% meshc(setax,setay,f1); % hold on; meshc(setax,setay,f2);% hold on; meshc(setax,setay,f3);% hold on; meshc(setax,setay,f4);% xlabel('方位角(rad)theta'); ylabel('俯仰角(rad)

23、phi'); zlabel('幅度'); title('四个子波束的天线方向图函数');% Fsum=f1+f2+f3+f4;% figure; meshc(setax,setay,Fsum); xlabel('方位角(rad)theta'); ylabel('俯仰角(rad)phi'); % zlabel('幅度'); title('和函数的三维图像');% Fdelst=f1+f4-f2-f3;% figure; meshc(setax,setay,Fdelst); xlabel(

24、9;方位角(rad)theta'); ylabel('俯仰角(rad)phi'); % zlabel('幅度'); title('方位角差函数的三维图像');% Fdelph=f1+f2-f3-f4;% figure; meshc(setax,setay,Fdelph); xlabel('方位角(rad)theta'); ylabel('俯仰角(rad)phi'); % zlabel('幅度'); title('俯仰角差函数的三维图像');% 目标自动测角targetx=(-

25、pi/8):(pi/1000):(pi/8); targety=(-pi/8):(pi/1000):(pi/8); f1=(exp( -1.4*( ( targetx-tp*ones(size(targetx) ).2 )/tr2 )'*(exp( -1.4*( ( targety-tp*ones(size(targety) ).2 )/tr2 ); %Af2=(exp( -1.4*( ( targetx+tp*ones(size(targetx) ).2 )/tr2 )'*(exp( -1.4*( ( targety-tp*ones(size(targety) ).2 )/t

26、r2 ); %Bf3=(exp( -1.4*( ( targetx+tp*ones(size(targetx) ).2 )/tr2 )'*(exp( -1.4*( ( targety+tp*ones(size(targety) ).2 )/tr2 ); %Cf4=(exp( -1.4*( ( targetx-tp*ones(size(targetx) ).2 )/tr2 )'*(exp( -1.4*( ( targety+tp*ones(size(targety) ).2 )/tr2 ); %DFsum=f1+f2+f3+f4; Fdelst=f1+f4-f2-f3; Fdelph=f1+f2-f3-f4;Esum=k*Fsum.2; Edelst=k.*Fsum.*Fdelst; Edelph=k.*Fsum.*Fdelph;F_ori=(exp(-1.4*(tp2)/(tr2)2;F_dao=-2.8*tp*exp(-1.4*