资源描述：
《仿真题目一-单脉冲和差测角仿真.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、仿真题目一单脉冲和差测角仿真题目要求：采用高斯型天线方向图绘制单脉冲和差测角的和、差波束及波形，并将数据制表，以便找出偏离等信号轴的角度，给定的值即可给出偏离角度。1.单脉冲和差测角原理雷达测角的物理基础是电波在均匀介质中传播的直线性和雷达天线的方向性，分为振幅法和相位法两大类，其中振幅法测角又分为最大信号法和等信号法，等信号测角采用两个相同且彼此部分重叠的波束，其方向图如下图1所示，若目标处在两波束的交叠轴OA方向，则两波束收到的信号强度相等，否则一个波束收到的信号强度高于另一个，故常称OS为等信号轴。当两个波束收到的回波信号相等时，等信号轴所指的方向即
2、为目标方向。若目标处在OB方向，波束2的回波比波束1的强，处在OC方向时，则与之相反，因此比较两个波束回波的强弱就可以判断目标偏离等信号轴的方向，并可用查表的方法估计出偏离等信号轴的大小。图1等信号测角（a）波束（b）显示器画面设天线电压方向性函数为，等信号轴OA的指向为，则波束1、2的方向性函数可分别写为为与波束最大值方向的偏角。用等信号法测量时，波束1接收到的回波信号，波束2收到的回波电压值，式中为目标偏离等信号轴的角度，这里对和信号进行和差法处理，可以获得目标信号的信息。由及可以求得其差值及和值，即在等信号轴附近差信号及和信号分别可近似表示为即可求得
3、其和差波束及，如图2所示。归一化的和差值为，由于正比于目标偏离的角度，故可用它来判读的大小及方向。图2和差测角法（a）两波束的方向图（b）差波束响应（c）和波束响应2.单脉冲和差测角仿真在8mm波段，。Matlab仿真程序及结果如下：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
4、).^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-'),gridonholdon,plot(theta,f2),xlabel('角度theta'),ylabel('两个响应');subplot(222),plot(theta,sigma),xlabel('角度theta'),ylabel('和波束Sigma'),gridonsubplot(223),p
5、lot(theta,delta),xlabel('角度theta'),ylabel('差波束Delta'),gridonsubplot(224),plot(theta,(delta./sigma)),gridonxlabel('角度theta'),ylabel('Delta/Sigma')3.数据制表，由于最大单值测角范围为有限，因此只考虑在[-30°,30°]范围内的数据，本文中为简化内容，只讨论[-15°，15°]范围内的数值。利用一次回归曲线拟合，得到对应的一次曲线。Matlab程序及结果如下：k=0.730;d=0.160;%参数设定la
6、bda=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
7、-0.1281-0.1140-0.0999-0.0857-0.0715-0.0572-0.0429-0.0286-0.014300.01430.02860.04290.05720.07150.08570.09990.11400.12810.14220.15620.17010.18400.19780.2115]';x=[ones(size(t))t];a=xff3=[ones(size(t))t]*a;figure，plot(theta,(delta./sigma),'g-')xlabel('角度theta'),ylabel('Delta/Sigma
8、');gridonholdonplot(t,f3,'r-')a=0
