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

雷达原理大作业评分单脉冲自动测角的原理及应用学院：电子工程学院作者：2016年5月21日单脉冲自动测角的原理及应用1单脉冲自动测角的原理及应用一．摘要单脉冲测角法是属于振幅法测角中的等信号法中的一种，其测角精度高，抗干扰能力强，在现实中得到了广泛的应用。而其中对于接收支路要求不太严格的双平面振幅和差式单脉冲雷达，更是备受青睐。本文首先讲述了单平面振幅和差式单脉冲雷达自动测角的原理，再简述了双平面振幅和差式单脉冲雷达自动测角的结构框图，接着简述了本文仿真所用的一些原理和公式推导，包括天线方向图函数及其导数的推导，最后做了基于高斯形天线方向图函数的单脉冲自动测角，基于辛克函数形天线方向图函数的单脉冲自动测角，和基于高斯形天线方向图函数的双平面单脉冲自动测角。源代码在附录里。二．重要的符号说明符号含义误差角0.5天线方向图函数半功率波瓣宽度，也就是归一化幅度值为0.707时对应的角度值0天线方向图函数零功率波瓣宽度单平面内，两波束相对天线轴线的偏角p双平面内，子波束在方位面相对天线轴线的偏角p双平面内，子波束在俯仰面相对天线轴线的偏角三．单平面振幅和差式单脉冲自动测角原理单脉冲测角法是属于振幅法测角中的等信号法中的一种。在单平面内，两个相同的波束部分重叠，交叠方向即为等信号轴的方向。将这两个波束接收到的回波信号进行比较就可以在一定范围内，一定精度要求下测到目标的所在角度。因为两个波束同时接到回波，故单脉冲测角获得目标角误差信息的时间可以很短，理论上只要分析一个回波脉冲即可，所以称之为“单脉冲”。因取出角误差的具体方式不同，单脉冲雷达种类很多，其中应用最广的是振幅和差式单脉冲雷达，其基本原理说明如下：1.角误差信号雷达天线在一个平面内有两个重叠的部分，如下图1所示：图1.振幅和差式单脉冲雷达波束图单脉冲自动测角的原理及应用2(a)两馈源形成的波束(b)和波束(c)差波束振幅和差式单脉冲雷达取得角误差信号基本方法是将这两个波束同时收到的信号进行和差处理，分别得到和信号和差信号。其中差信号即为该角平面内角误差信号。若目标处在天线轴方向（等信号轴），误差角0，则两波束收到的回波信号振幅相同，差信号等于0。目标偏离等信号轴而有一个误差角时，差信号输出振幅与成正比而其符号则由偏离方向决定。2.和差比较器这里主要使用双T插头，示意图如下图2（a）所示。它有四个端口：和端，差端和1，2端。假定四个端都是匹配的，则从和端输出信号时，1,2端输出等幅同相的信号，差端无输出；从1,2端输入同相信号时，和端输出两信号之和，差端输出两信号之差。图2.双T接头和差比较器示意图(a)双T接头(b)和差比较器示意图在发射信号时，从发射机来的信号加在和端，故1,2端输出等幅同相的信号，两波束在空间各点产生的场强同相相加，形成发射和波束的天线方向性函数为()F。接收时，回波脉冲同时加到1,2端，此时在和端，输出两个回波信号同相相加之和，记为E；在差端，输出两信号反相相加之和，记为E。假设两个波束方向性函数完全相同，记为()F，两波束衰减倍数为k，两波束相对天线轴线的偏角为，则对于方向的目标来说：和信号振幅为：2()()()()()EkFFkFFkF差信号振幅为：()()()()()()EkFFkFFkFF其中：()()()FFF，()()()FFF。实际情况下，是很小的，可以对()F和()F在附近做一阶泰勒展开：单脉冲自动测角的原理及应用322()()()()()()()()()()()()FFFoFFFFFoFF则：222()[()()()()]4()EkFkFFFFkF()()[()()()()][(()())(()())]EkFFkFFFFFFFF4()()kFF则：24()()()4()()EkFFFEkFF，其中是个常数，()()FF。据此，可以求得1EE，这就是单平面振幅和差式单脉冲雷达测角的公式。四．双平面振幅和差式单脉冲自动测角简述为了对空中目标进行自动方向跟踪，必须在方位角和俯仰角两个角平面上进行角跟踪，因而必须获得方位角和俯仰角的误差信号。为此，需要用4个馈源照射一个反射体，以形成四个对称的相互部分重叠的波束。双平面上天线四个子波束在三维空间中的辐射图图3.雷达天线空间子波束分布图假设其中一个子波束的三维方向图函数为(,)F，其中代表方位角，代表俯仰角，p和p分别为子波束在方位面和俯仰面相对中心轴的偏离角。则图3所示的4个子波束的数学表达式为：子波束A：(,)(,)AppFF子波束B：(,)(,)BppFF单脉冲自动测角的原理及应用4子波束C：(,)(,)CppFF子波束D：(,)(,)DppFF双平面振幅和差式单脉冲自动测角的原理框图如下图4所示：图4.双平面振幅和差式单脉冲自动测角的原理框图其中和信号为(,)(,)(,)(,)(,)ABCDFFFFF俯仰角误差信号为：(,)(,)(,)(,)(,)ABCDFFFFF方位角误差信号为：(,)(,)(,)(,)(,)ADBCFFFFF再根据之前的单平面振幅和差式单脉冲雷达测角的公式有：1EE1EE其中，,1(,)(,)ppppFF，,1(,)(,)ppppFF。五．仿真所用原理和公式推导1.单平面单脉冲雷达方向图函数A．高斯形方向图函数单脉冲自动测角的原理及应用5其方向图函数为：220.51.4()Fe其导函数为220.51.420.52.8()Fe。B．辛克函数形方向图函数其方向图函数为：0()(2)FSa，其导函数为00020222cos()sin()()2F2.双平面脉冲雷达方向图函数此处对三维天线使用简化模型，将三维天线看做由两个二维平面（方位面和俯仰面）方向图相乘的结果，即(,)()()FFF。由于仿真中只以高斯形方向图函数为例做了仿真，所以此处只讨论三维的高斯形方向图函数，设p和p分别为子波束在方位面和俯仰面相对中心轴的偏离角：22220.50.51.41.4(,)*Fee其偏导函数为：22220.50.51.41.42,0.5(,)2.8*pppppFee，22220.50.51.41.420.5,(,)2.8*pppppFee。3.关于两波束相对天线轴线的偏角的选取：两波束相对天线轴线的偏角一般选择为：单脉冲自动测角的原理及应用60.5120.5是半功率波瓣宽度，也就是归一化幅度为0.707时候对应的值在单平面上，直接取0.512；在双平面上，可以取0.512pp。六．仿真结果1.基于高斯形天线方向图函数的单脉冲自动测角下图是两馈源形成的波束，和波束和差波束波形图下图是单平面上目标所在角度和测得角度的比较曲线：单脉冲自动测角的原理及应用72.基于辛克函数形天线方向图函数的单脉冲自动测角下图是两馈源形成的波束，和波束和差波束波形图：单脉冲自动测角的原理及应用8下图是单平面上目标所在角度和测得角度的比较曲线：3.基于高斯形天线方向图函数的双平面单脉冲自动测角下图是双平面上天线四个子波束的方向图函数图像：单脉冲自动测角的原理及应用9下图是三维和函数的图像：下图是方位角差函数图像：单脉冲自动测角的原理及应用10下图是俯仰角差函数图像：下图是目标方位角与俯仰角之积的三维曲面和测得方位角与俯仰角之积的三维曲面：单脉冲自动测角的原理及应用11下图是目标方位角与测得的方位角的比较：下图是目标俯仰角与测得的俯仰角的比较：单脉冲自动测角的原理及应用12七．附录源码ld.m（单平面的仿真代码）：===========================================================================================clear;closeall;clc%%高斯函数的方向图函数绘图和自动测角%%parameters%theta_r=0.4;%theta_p=0.2;%k=1;%%operations%theta=(-pi):(2*pi/1000):pi;%f1=exp(-1.4*((theta-theta_p).^2)/theta_r^2);%f2=exp(-1.4*((theta+theta_p).^2)/theta_r^2);%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(\theta-\delta)波形');%holdon;plot(theta,0.7*ones(size(theta)),'k');text(1,0.8,'归一化幅度=0.707');%subplot(3,1,2);plot(theta,sigma,'Linewidth',2);axistight;单脉冲自动测角的原理及应用13%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*pi/1000):(pi/8);%result=zeros(size(target));%Fsum=exp(-1.4*((target-theta_p*ones(size(target))).^2)/theta_r^2)+exp(-1.4*((target+theta_p*ones(size(target))).^2)/theta_r^2);%Fdelta=exp(-1.4*((target-theta_p*ones(size(target))).^2)/theta_r^2)-exp(-1.4*((target+theta_p*ones(size(target))).^2)/theta_r^2);%Esum=k*Fsum.^2;%Edelta=k.*Fsum.*Fdelta;%F_ori=exp(-1.4*(theta_p^2)/(theta_r^2));%F_dao=-2.8*theta_p*exp(-1.4*(theta_p.^2)/theta_r^2)/(theta_r^2);%res=-(Edelta*F_ori)./(Esum*F_dao);%figure;plot(target,target,target,res,'--r','Linewidth',2);axistight;legend('目标所在角度','测角所得角度');%xlabel('角度值(rad)');ylabel('目标所在角度值和