仿真题目一 单脉冲和差测角仿真.docx

仿真题目一单脉冲和差测角仿真题目要求：采用高斯型天线方向图绘制单脉冲和差测角的和、差波束及-波形, 并将-数据制表，以便找出偏离等信号轴的角度，给定-的值即可给出偏离角 度1.单脉冲和差测角原理雷达测角的物理基础是电波在均匀介质中传播的直线性和雷达天线的方向 性，分为振幅法和相位法两大类,其中振幅法测角又分为最大信号法和等信号法, 等信号测角采用两个相同且彼此部分重叠的波束，其方向图如下图1所示，若目 标处在两波束的交叠轴OA方向，则两波束收到的信号强度相等，否则一个波束 收到的信号强度高于另一个，故常称OS为等信号轴当两个波束收到的回波信 号相等时，等信号轴所指的方向即为目标方向若目标处在OB方向，波束2的 回波比波束1的强，处在OC方向时，贝V与之相反，因此比较两个波束回波的强 弱就可以判断目标偏离等信号轴的方向，并可用查表的方法估计出偏离等信号轴 的大小M方向 则方向Cb)Ca)图1等信号测角(a)波束(b)显示器画面设天线电压方向性函数为F(9 )，等信号轴OA的指向为0，则波束1、2的方向 0性函数可分别写为F (0)二 F(0 )二 F(0+0 -0 )1 1 k oF (0)二 F(0 )二 F(0-0 -0 )2 2 k o0为0与波束最大值方向的偏角。

k 0用等信号法测量时，波束1接收到的回波信号u二KF (0)二KF(0 -0 )，波束21 1 k t收到的回波电压值u二KF (0)二KF(-0 -0 )二KF(0 +0 )，式中0为目标偏离等2 2 k t k t t信号轴0的角度，这里对u和u信号进行和差法处理，可以获得目标信号0的信0 1 2 t息由u及u可以求得其差值A(6)及和值2(9)，即1 2A(6)二 u (6) -u (6)二 K[F(6 -6 ) - F(6 +6 )]1 2 kt k t2(6)二 u (6) + u (6)二 K[F(6 -6 ) + F(6 +6 )]1 2 k t k t在等信号轴附近差信号及和信号分别可近似表示为A(6 ) q 26 dF(6)| kt t d6 6=6o2 (6 ) q 2 F (6 ) kt o即可求得其和差波束2(6)及A(6)，如图2所示归一化的和差值为A 6 dF (6), = t | ，2 F (6 ) d6 6=6oo由于△正比于目标偏离6的角度6，故可用它来判读6的大小及方向2F*崗 IvF图2和差测角法(a)两波束的方向图(b)差波束响应(c)和波束响应2.单脉冲和差测角仿真高斯型天线方向图函数为f (0)二e」3863(e”)2, 0为半功率波束宽度，在8mm波段，rk = = 730mm,d = 160mm，天线波束宽度为0 = 1.2 ,波束倾角0 =0 /3。

入 -3db d k -3dbMatlab仿真程序及结果如下：k=0.730;%参数设定d=0.160;labda=2*pi/k; %波长 theta_3db=1.2*labda/d; % 天线波束宽度 theta_k=theta_3db/3; %相对等场强方向的波束倾斜角 theta=-2*theta_3db:0.2:2*theta_3db;f1=exp(-1.3863*(theta-theta_k).入2/theta_3db入2); f2=exp(-1.3863*(theta+theta_k).入2/theta_3db入2);sigma=f1+f2;delta=f1-f2; figure,subplot(221),plot(theta,f1,'r-'),grid on hold on,plot(theta,f2),xlabel('角度\theta' ),ylabel('两个响应'); subplot(222),plot(theta,sigma),xlabel('角度 \theta' ),ylabel('和波束 \Sigma'),grid onsubplot(223),plot(theta,delta),xlabel('角度 \theta' ),ylabel('差波束 \Delta'),grid onsubplot(224),plot(theta,(delta./sigma)),grid on xlabel('角度\theta' ),ylabel( '\Delta/\Sigma')角度B1 ； ；—角度B归一化和差倩1 ; : ;ftv度角3. A数据制表，E由于最大单值测角范围为有限，因此只考虑0在[-30° ,30° ]范围内的数 据，本文中为简化内容，只讨论[-15°，15°]范围内的数值。

利用一次回归曲 线拟合，得到对应的一次曲线Matlab 程序及结果如下：k=0.730;d=0.160;%参数设定labda=2*pi/k;% 波长theta_3db=1.2*labda/d;% 天线波束宽度 theta_k=theta_3db/3; % 相对等场强方向的波束倾斜角 theta=-15:1:15;f1=exp(-1.3863*(theta-theta_k).入2/theta_3db入2);f2=exp(-1.3863*(theta+theta_k).入2/theta_3db入2);sigma=f1+f2;% 和波束函数delta=f1-f2;%差波束函数t=[-15:1:15]';f=[-0.2115 -0.1978 -0.1840 -0.1701 -0.1562 -0.1422 -0.1281 -0.1140 -0.0999-0.0857-0.0715-0.0572 -0.0429 -0.0286 -0.014030.01430.02860.04290.05720.07150.08570.09990.11400.12810.14220.1562 0.1701 0.1840 0.1978 0.2115]';x=[ones(size(t)) t];a=x\ff3=[ones(size(t)) t]*a;figure，plot(theta,(delta./sigma),'g-')xlabel('角度\theta' ),ylabel( '\Delta/\Sigma');grid onhold onplot(t,f3,'r-')a =0.00000.0142得到直线方程为△辽=0.01420从程序运行结果可以看出，利用二次回归曲线拟合得到的圆点直线和原始数据绘 制出来的直线基本上重合，可以看出，结果的误差是很小的。

得到下面的表格:9-15-14-13-12-11-10-9-8-7-6-5-4-3-2-1算得-0.-0.-0.-0.-0.-0.-0.-0.-0.-0.-0.-0.-0.-0.-0.A211918417015614212811409908070504020115780122109571572298643拟合-0.-0.-0.-0.-0.-0.-0.-0.-0.-0.-0.-0.-0.-0.-0.A 得—211918417015614112711309908070504020127853108642510967258442接上表012345678910111213141500.00.00.00.00.00.00.00.10.10.10.10.10.10.10.214328642957271585799914028425670849711122108500.00.00.00.00.00.00.00.10.10.10.10.10.10.10.2142284425567709851992134274156708498126801357由上表也可以看出计算得到的和曲线拟合出来的结果误差很小。