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

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

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

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

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

5、波束衰减倍数为 k,两波束相对天假设两个波束方向性函数完全相同，记为线轴线的偏角为 ，则对于方向的目标来说和信号振幅为：EkF ( )F()差信号振幅为：EkF ( )F()其中：F （）F()F()，实际情况下，是很小的，可以对 F（kF ()F()kF2()kF ()F()kF ( )F ()F ()F()F( )。)和F(）在附近做一阶泰勒展开为E ;在差端，输出两信号反相相加之和，记为 EF( ) F( ) F ( ) o( 2) F( ) F ()F( ) F( ) F ( ) o( 2) F( ) F ()4. 双平面振幅和差式单脉冲自动测角简述为了对空中目标进行自动方向跟踪，必须

6、在方位角和俯仰角两个角平面上进行角跟踪，因而必须获得方位角和俯仰角的误差信号。 为此，需要用4个馈源照射一个反射体，以形成双平面上天线四个子波束在三维空间中的辐射图四个对称的相互部分重叠的波束。图3.雷达天线空间子波束分布图假设其中一个子波束的三维方向图函数为 F(,)，其中 代表方位角， 代表俯仰角，p和p分别为子波束在方位面和俯仰面相对中心轴的偏离角。则图3所示的4个子波束的数学表达式为：子波束 A: Fa( , ) F( p , p )子波束 B： Fb( , ) F( p子波束 C: Fc( , ) F( p , p )子波束 D: Fd( , ) F( p , p )双平面振幅和差式

7、单脉冲自动测角的原理框图如下图 4所示:其中和信号为F (,)Fa(,)Fb( , ) Fc(,)Fd(,)俯仰角误差信-号为：F (,)Fa( , ) Fb(,)Fc(,)Fd(,)方位角误差信-号为：F (,)Fa(，) Fd(,)Fb(,)Fc(,)再根据之前的单平面振幅和差式单脉冲雷达测角的公式有:1E1EE|E|其中,1 F (,)F ( p, p)p p1 F (,)F ( p, p)p p5. 仿真所用原理和公式推导1单平面单脉冲雷达方向图函数A.高斯形方向图函数其方向图函数为：F ()21.4 _20.50.52.82_0.5其导函数为F ()B.辛克函数形方向图函数其方向图函

8、数为：F()Sa(2 ),02 2 cos( ) sin( )其导函数为F ()0 022将三维天线看做由两个二维平面 (方位面和俯仰面)方2双平面脉冲雷达方向图函数 此处对三维天线使用简化模型,向图相乘的结果，即 F( , ) F)F ()。由于仿真中只以高斯形方向图函数为例做了仿真,所以此处只讨论三维的高斯形方向图函数，设p和p分别为子波束在方位面和俯仰面相对中心轴的偏离角:21.420.51 4 20.5其偏导函数为:3关于两波束相对天线轴线的偏角的选取： 两波束相对天线轴线的偏角 一般选择为:0.50.5是半功率波瓣宽度，也就是归一化幅度为0.707时候对应的 值在单平面上，直接取6.

9、 仿真结果1.基于高斯形天线方向图函数的单脉冲自动测角 下图是两馈源形成的波束，和波束和差波束波形图F图是单平面上目标所在角度和测得角度的比较曲线:2基于辛克函数形天线方向图函数的单脉冲自动测角 下图是两馈源形成的波束，和波束和差波束波形图：F(fl+A+F(M)波矗10.501.30.60.20-0.2 一化啣叮=0.7070和度fiLgdWF偲d)*F佩*刖波形0旬度值(填F(4F锻*股形辛克力向菌数的口动测角结果和訂标所乔方向角度ffi比较F图是单平面上目标所在角度和测得角度的比较曲线：0.40.30.20,10-0.1-0.2-0,3-0.4-0.3 -0.2 -0.1 0 0.1 0

10、.2角度B(rad)3.基于高斯形天线方向图函数的双平面单脉冲自动测角 下图是双平面上天线四个子波束的方向图函数图像：H牛波柏山电V川胸函敢3四个子波束的天线方向图函F图是三维和函数的图像:和函数的二樂国悵F图是方位角差函数图像:方位角签團数的三堆图像F图是俯仰角差函数图像:俯仰角赛團数的F图是目标方位角与俯仰角之积的三维曲面和测得方位角与俯仰角之积的三维曲面:目标方位角与俯仲角2积与刮得的方位坤耳帕何坤之段的怨国俊D.60.30.20.10-0.1-0.2-0.3F图是目标方位角与测得的方位角的比较：口标的方位毎和测得的方位角比较 H标的厅检角(rad)0 测得的方悅爲(用日)心0.2 0.

11、3-0.3 -0.2 -0.1 0 0.12/1 ft (rad)F图是目标俯仰角与测得的俯仰角的比较：H标的俯仰角和测得的K仰角比较0.30.20J0-0.1-0.2-0.3-0.3 02 -0.1 0 0.1 0.2 0.3俯仰ft(rad)7. 附录源码Id.m （单平面的仿真代码）clear; close all; clc% %高斯函数的方向图函数绘图和自动测角% % parameters% theta_r=0.4;% theta_p=0.2;% k=1;% % operati ons% theta=(-pi):(2*pi/1000):pi;% f1=exp(-1.4*(theta -

12、theta_p).A2)/theta_rA2);% f2=exp(-1.4*(theta+theta_p).A2)/theta_A2);% sigma=f1+f2; delta=f1 -f2;% figure(1);% subplot(3,1,1); plot(theta,f1,theta,f2, -,Linewidth,2); axis tight;% xlabel(角度值(rad)theta); ylabel(幅度);title(F(theta+delta)+F(theta -delta)波形);% hold on; plot(theta,0.7*ones(size(theta),k);

13、text(1,0.8,归一化幅度=0.707);% subplot(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(F(theta -delta) -F(theta+d

14、elta) 波形 );% target=(-pi/8):(2*pi/1000):(pi/8);% result=zeros(size(target);%Fsum=exp(-1.4*(target -theta_p*ones(size(target).A2)/theta_A2)+exp( -1.4*(target+theta_p*o n es(size(target).A2)/theta_rA2);%Fdelta=exp(-1.4*(target -theta_p*ones(size(target).A2)/theta_rA2) -exp(-1 .4*(target+theta_p*o nes(size(target).A2)/theta