2022年2022年雷达信号matlab仿真 .pdf

《2022年2022年雷达信号matlab仿真 .pdf》由会员分享，可在线阅读，更多相关《2022年2022年雷达信号matlab仿真 .pdf（28页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、雷达系统分析大作作者： 陈雪娣学号： 0410420727 1. 最大不模糊距离：,max1252urCRkmf距离分辨率：1502mcRmB2. 天线有效面积：220.07164eGAm半功率波束宽度：36.44odbG3. 模糊函数的一般表示式为22*2;dtetstsfdfjd对于线性调频信号21Rejtppts tcteTT则有：221;ReRepjt TjtdpppttfctcteedtTTTsin1;11dppdpdppfTTfTfTT分别令0,0df可得220;,; 0dfsin0;dpddpf Tff T名师资料总结 - - -精品资料欢迎下载 - - - - - - - -

2、- - - - - - - - - - 名师精心整理 - - - - - - - 第 1 页，共 28 页 - - - - - - - - - sin1;011pppppTTTTT程序代码见附录 1 的 T_3.m, 仿真结果如下 : 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 2 页，共 28 页 - - - - - - - - - 4. 程序代码见附录 1 的 T_4.m, 仿真结果如下 : 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - -

3、 - - - - - - - 名师精心整理 - - - - - - - 第 3 页，共 28 页 - - - - - - - - - 通过比较得知， 加窗后的主副瓣比变大， 副瓣降低到 40db以下，但主瓣的宽度却增加了，约为未加窗时的1.5 倍，主瓣也有一定的损失。5. 由雷达方程221340(4 )tPGTeSNRKT LFR计算可得1196.55 40logSNRR db 作图输出结果如下，程序代码见附录1 的 T_5.m 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 4 页

4、，共 28 页 - - - - - - - - - 在 R=70km时, 计算得单个脉冲的SNR1=2.7497 db, 要达到要求的检测性能则需要12.5dB 的最小检测输入信噪比，而M个相参脉冲积累可以将信噪比提高M倍， 故10) 1(SNRDM=9.4413 因此要达到要求就需要10 个以上的相参脉冲进行积累。可求得可积累脉冲数为：3256dbraNf其中,a为天线的搜索速度等于30o/s.rf是重复频率为 1200hz. 故满足要求 . 6. 设 t 时刻弹舰径向与目标航向的夹角为a （t ）, 目标偏离弹轴方向的夹角为t（ ）, 在t=0 时,31oo（ ）, 1o（0）. 由几何关

5、系知 , sincosooOMRMPR经 t 秒后, cossinsaaM PMPV tV tO MOMV t( )O MtarctanM P( )( )ttsin( )sin( )sin( )aOMV tO MR tttcos ( )cos ( )dasVVtVt又因为cos( )cos( ( )ttcos ( )cos sin( )sin31cos( )sin( )22tttt故31cos ( )sin( )cos( )22dasVVttVt名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - -

6、 - 第 5 页，共 28 页 - - - - - - - - - ( )dft31cos ( )sin( )22dasVVVtt2231cos ( )sin( )22ddasVfVVtt利用以上的关系式即可计算出第i 个重复周期弹目间的距离( )iR t和回波信号的多普勒频率( )dft. 仿真程序代码见附录1 的 T_6.m.实验结果如下 : 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 6 页，共 28 页 - - - - - - - - - 由仿真结果可知, ( )dft的

7、变化不大 , 这表明相对速度的变化不大，同时可求得688/Vdm s. 7. (1)相干积累 : 由于相对速度的变化不大，所以在仿真时取定值688/Vdm s。仿真程序代码见附录1 的 T_7_1.m.实验结果如下 : 相干积累前后的信噪比情况如下图所示：名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 7 页，共 28 页 - - - - - - - - - 由仿真结果知，积累前匹配滤波器输出的信噪比为约12dB。已知 M 个脉冲相参积累可以将信噪比提高M 倍，所以， 64 个脉冲相

8、参积累后的信噪比将提高64 倍（18db） 。相干积累后输出的信噪比约30db，与预期效果相符。(2) 非相干积累：双极点滤波器的时域框图如下：由此可的它的频域响应：22exp2/1dk其中212exp/1cos()ddk22exp2/1dk式中：0.63，2.2dN，N是半功率波束宽度。仿真程序代码见附录1 的 T_7_2.m.实验结果如下 : 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 8 页，共 28 页 - - - - - - - - - 非相干积累前后的信噪比情况如下图

9、所示：名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 9 页，共 28 页 - - - - - - - - - 由仿真结果知， 积累前匹配滤波器的信噪比为约12dB。非相干积累后输出的信噪比约20db。将非相干的结果与相干积累的效果进行比较，可知，相干积累的效果明显优于非相干积累。附录 1 程序代码第 3 题: % T_3.m % clear all clc clf taup=1; %脉冲宽度 100us b=10; %带宽up_down=-1; %up_down=-1正斜率， up

10、_down=1负斜率x=lfm_ambg(taup,b,up_down); %计算模糊函数taux=-1.1*taup:.01:1.1*taup; fdy=-b:.01:b; figure(1) mesh(100*taux,fdy./10,x) %画模糊函数xlabel(Delay - mus) ylabel(Doppler - MHz) zlabel(| chi ( tau,fd) |) title(模糊函数 ) figure(2) contour(100.*taux,fdy./10,x) %画等高线xlabel(Delay - mus) ylabel(Doppler - MHz) titl

11、e(模糊函数等高线 ) grid on N_fd_0=(length(fdy)+1)/2; % fd=0 的位置x_tau=x(N_fd_0,:); % 时间模糊函数figure(3) plot(100*taux,x_tau) axis(-110 110 0 1) xlabel(Delay - mus) ylabel(| chi ( tau,0) |) title( 时间模糊函数 ) grid on N_tau_0=(length(taux)+1)/2; % tau=0 的位置x_fd=x(:,N_tau_0); % 速度模糊函数figure(4) plot(fdy./10,x_fd) 名师资

12、料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 10 页，共 28 页 - - - - - - - - - xlabel(Doppler - MHz) ylabel(| chi ( 0,fd) |) title( 速度模糊函数 ) grid on x_db=20*log10(x+eps); I,J=find(abs(x_db+6)0.09); %取 6db 点的位置I=(I-b/.01)/(1/.01); %Doppler维 坐标变换J=(J-1.1*taup/.01)/(1/.01);

13、%时间维 坐标变换figure(5) %6db 的等高线plot(J*100,I/10,.) axis(-110 110 -1 1) xlabel(Delay - mus) ylabel(Doppler - MHz) title(模糊函数 6db 的等高线 ) grid on %- - - - 模糊函数 - - - function x=lfm_ambg(taup,b,up_down) % taup 脉冲宽度 ; % b 带宽; %up_down=-1 正斜率， up_down=1负斜率eps=0.0000001; i=0; mu=up_down*b/2./taup; for tau=-1.1

14、*taup:.01:1.1*taup i=i+1; j=0; for fd=-b:.01:b j=j+1; val1=1-abs(tau)/taup; val2=pi*taup*(1-abs(tau)/taup); val3=(fd+mu*tau); val=val2*val3+eps; x(j,i)=abs(val1*sin(val)/val); end end % 第 4 题: %T_4.m % 利用频域处理方法进行脉冲压缩 % clear all clc clf eps = 1e-10; Te=100e-6; %脉冲带宽名师资料总结 - - -精品资料欢迎下载 - - - - - - -

15、 - - - - - - - - - - - 名师精心整理 - - - - - - - 第 11 页，共 28 页 - - - - - - - - - Bm=1e6; %调频mu=Bm/Te; %调频斜率Ts=1/(2*Bm); %采样周期Ns=fix(Te/Ts); %采样点数Nf=1024; % fft点数t=0:Ts:Te-Ts; y=exp(j*pi*mu*t.2); %脉冲压缩前的线形调频信号yfft = fft(y,Nf) ; h=zeros(1,Ns); for i=1:Ns h(i)=conj(y(Ns-i+1); end hfft= fft(h,Nf); % 匹配滤波器的频域

16、响应ycomp =abs(ifft(yfft .*hfft); %脉冲压缩maxval = max (ycomp); ycomp = eps + ycomp ./ maxval; % 利用最大值归一化ycomp_db=20*log10(ycomp); % 取对数% 加窗处理 % win = hamming(Ns); h_w=h.*win; % 加窗hfft_w=fft(h_w,Nf); % 加窗的匹配滤波器的频域响应ycomp_w = abs(ifft(yfft .*hfft_w); %脉冲压缩maxval1 = max(ycomp_w); val=ycomp_w ; ycomp_w = ep

17、s + ycomp_w ./ maxval; % 利用 ycomp的最大值归一化ycomp_w1 = eps + val./ maxval1; % 利用 ycomp_w 的最大值归一化ycomp_w_db=20*log10(ycomp_w); % 取对数ycomp_w1_db=20*log10(ycomp_w1); % 取对数% tt =0:Ts:2*Te-Ts; figure(1) plot (tt,ycomp_db(1:2*Ns),g) axis(.2*Te 1.8*Te -60 0 ) xlabel (t - seconds ); ylabel( db) title(没有加窗的脉冲压缩输

18、出 ) grid on figure(2) plot (tt,ycomp_w1_db(1:2*Ns),r) axis(.2*Te 1.8*Te -60 0 ) xlabel (t - seconds ); ylabel( db) title(加窗的脉冲压缩输出 ) grid on 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 12 页，共 28 页 - - - - - - - - - figure(3) plot (tt,ycomp_db(1:2*Ns),g,tt,ycomp_w_

19、db(1:2*Ns),r) axis(.7*Te 1.3*Te -60 0 ) xlabel (t - seconds ); ylabel( db) legend(未加窗 ,加窗); title(脉冲压缩输出对比 ) grid on % 第 5 题: % T_5.m %SNR与距离的关系 % clear all clc eps=1e-10; c = 3.0e+8; % speed of light lambda =0.03; % 波长pt=20; %峰值功率lambda=0.03; % 波长tao=100e-6; %发射脉冲宽度G_db=30; %天线增益 in db sigma=1000;

20、%RCS k=1.38e-23; % Boltzmans constant To=290; %标准室温F_db=2; % 噪声系数 in db L_db=5; % 系统损失 in db R=70e3:-100:0; %距离val=10*log10(pt*tao*lambda2*sigma)/(4*pi)3*k*To)+ 2*G_db-F_db-L_db; SNR=val-40*log10(R); figure(1) plot(R./1e3,SNR) title(SNR与距离的关系 ) xlabel(距离 - km) ylabel(SNR - db) grid on SNR1=val-40*lo

21、g10(70e3)% 计算 R=70km 时的 SNR % 第 6 题: % T_6.m clear all clc lembda=0.03; %波长fr=1200; %重复频率名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 13 页，共 28 页 - - - - - - - - - tra=180/pi; % 度到弧度的转化量alpha=31/tra; % t=0时弹舰径向与目标航向的夹角alpha_p=30/tra;% 导弹运动方向与目标航向的夹角Ro=70e3; % t=0时的

22、弹舰距离Vs=10; % 舰船速度Va=680; % 导弹速度OM=Ro*sin(alpha); % t=0时 弹舰垂直距离MP=Ro*cos(alpha); % t=0时 弹舰垂直距离% t=0:1/fr:10; OM_t=OM-0.5*Va.*t; % t时刻 弹舰垂直距离MP_t=MP-Vs.*t-sqrt(3)*Va.*t/2;% t时刻 弹舰垂直距离alpha_t=atan(OM_t./MP_t); %t时刻弹舰径向与目标航向的夹角R_t=OM_t./sin(alpha_t); %t时刻弹舰距离% t 时刻弹舰径向速度vd_t=(sqrt(3)/2*Va+Vs).*cos(alpha

23、_t)+0.5*Va.*sin(alpha_t); fd_t=2*vd_t/lembda; %t时刻多普勒频移% figure(1) plot(t,R_t) title (t时刻弹舰距离 ) xlabel(时间 - s) ylabel(弹舰距离 - m) figure(2) plot(t,fd_t) title(t时刻多普勒频移 ) xlabel(时间 - s) ylabel(多普勒频移 - hz) % 第 7 题: %T_7_1.m % 相干积累clear all clc clf c=3e8; % speed of light Te=100e-6; % 发射脉冲宽度Be=1e6; %带宽mu

24、=Be/Te; %调频斜率Ts=1/(2*Be); %采样频率Ro=70e3; % 起始距离fo=c/0.03; % 中心频率Vr=688; %径向速度名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 14 页，共 28 页 - - - - - - - - - t=0:Ts:Te-Ts; W=exp(j*pi*mu*t.2); Wf=fft(W,1024); % nnn=fix(Ro-30e3)/75);%采样的起始位置 , 从 30km开始采样R=0:75:15e3-75; %在 3

25、0km和 45km之间采样，采样间隔75m for i=1:200 for k=1 :64 Ri(k,i)=R(i)-Vr*Ts*(k-1); end end taoi=2*Ri/c; echo=100.225*0.707*(randn(64,1024)+j*randn(64,1024); j=sqrt(-1); for i=1:64 %回波信号 echo(i,nnn:nnn+199)=echo(i,nnn:nnn+199)+. exp(-j*2*pi*fo*taoi(i,:)+j*pi*mu*taoi(i,:).2); end for i=1:64 %脉冲压缩 sp2(i,:)=ifft(f

26、ft(echo(i,:),1024).*conj(Wf),1024); end for k=1:1024 % 相干积累 sct(:,k)=abs(fftshift(fft(sp2(:,k),256); end sct=sct./max(max(sct);%归一化 sp=sp2./max(max(sp2);%归一化 %积累前后信噪比输出 figure(1) plot(20*log10(abs(sp) ylabel(-db) title(相干积累前 ) axis(1 1024 -30 0) figure(2) plot(20*log10(sct) ylabel( - db) title(相干积累输

27、出 ) axis(1 1024 -30 0) %积累结果输出r=(1:1024)*75+30e3)./1e3; dp=(-128:127)*(Be/128)/1e3; figure(1) mesh(r,dp,sct) xlabel(距离 km) ylabel(Doppler - kHz) 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 15 页，共 28 页 - - - - - - - - - title(相干积累输出结果 ) figure(2) contour(r,dp,sct)

28、axis(30 100 -200 200) xlabel(距离 km) ylabel(Doppler - kHz) title(R-fd 等高线 ) grid on dp=(-32:31)*(Be/32)/1e3; figure(3) mesh(r,dp,abs(echo)/max(max(abs(echo) xlabel(距离 km) ylabel(Doppler - kHz) title(相干积累前的结果 ) % T_7_2.m % 非相干积累clc clear all c=3e8; % speed of light Te=100e-6; % 发射脉冲宽度Be=1e6; %带宽mu=Be/

29、Te; %调频斜率Ts=1/(2*Be); %采样频率Ro=70e3; % 起始距离fo=c/0.03; % 中心频率Vr=688; %径向速度fr=1200; %重复频率t=0:Ts:Te-Ts; W=exp(j*pi*mu*t.2); Wf=fft(W,1024); %- 双极点滤波器 -% sheta_3_db=6.4; % 半功率波束宽度v=30; % 天线的搜索速度N=sheta_3_db*fr/v; wd_tao=2.2/N; xi=0.63; k1=2*exp(-xi*wd_tao/sqrt(1-xi2)*cos(wd_tao); k2=exp(-2*xi*wd_tao/sqrt

30、(1-xi2); NN=64; w=-pi:pi/NN:pi-pi/NN; j=sqrt(-1); H=exp(-j.*w)./(1-k1*exp(-j*w)+k2*exp(-2*j.*w); h=ifft(H,64); 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 16 页，共 28 页 - - - - - - - - - %- - - - 信号处理 - - - % nnn=fix(Ro-30e3)/75);%采样的起始位置 , 从 30km开始采样R=0:75:15e3-75;

31、 %在 30km和 45km之间采样，采样间隔75m for i=1:200 for k=1 :64 Ri(k,i)=R(i)+Vr*Ts*(k-1); end end taoi=2*Ri/c; echo=100.275*0.707*randn(64,1024)+j*randn(64,1024); j=sqrt(-1); for i=1:64 %回波信号，加随机相位模拟非相干信号 echo(i,nnn:nnn+199)=echo(i,nnn:nnn+199). +exp(-j*2*pi*fo*taoi(i,:)+j*pi*mu*taoi(i,:).2 . +j*2*pi*rand*ones(1

32、,200); end for i=1:64 %脉冲压缩 sp2(i,:)=ifft(fft(echo(i,:),1024).*conj(Wf),1024); end for i=1:1024 %用双极点滤波器进行非相干积累 isct(:,i)=conv(sp2(:,i),h); end for k=1:1024 sct(:,k)=abs(fftshift(fft(isct(:,k),256); end sct=sct./max(max(sct); %归一化 sp2=sp2./max(max(abs(sp2); %归一化 %积累前后信噪比输出 figure(1) plot(20*log10(ab

33、s(sp2) ylabel(-db) title(非相干积累前 ) axis(1 1024 -30 0) figure(2) plot(20*log10(sct) ylabel( - db) title(非相干积累输出 ) axis(1 1024 -30 0) %积累结果输出 r=(1:1024)*75+30e3)./1e3; dp=(-128:127)*(Be/128)./1e3; figure(3) mesh(r,dp,sct) xlabel(距离 km) 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - -

34、 - - - - - 第 17 页，共 28 页 - - - - - - - - - ylabel(Doppler - kHz) title(非相干积累输出结果 ) figure(4) contour(r,dp,sct) axis(30 100 -200 200) xlabel(距离 km) ylabel(Doppler - kHz) title(R-fd 等高线 ) grid on dp=(-32:31)*(Be/32)/1e3; figure(5) mesh(r,dp,abs(echo)/max(max(abs(echo) xlabel(距离 km) ylabel(Doppler - kH

35、z) title(非相干积累前的结果 ) % % 二 翻译11.2比幅单脉冲比幅单脉冲跟踪类似于对于圆形区域而言需要四个斜的波束来测量目标的角度位置。区别就是这四个波束是同时产生而不是相继产生的。因为这个目的，一种特殊的天线馈电正是利用一个单独的脉冲产生四个波束原理，因此它的名字叫做单脉冲。并且，单脉冲跟踪更加精确并且不容易受圆柱形的近点角的影响，比如调幅干扰和获得的倒置的电子对抗措施。最后，连续并且圆锥形变化的雷达回波降低了跟踪的精度；尽管如此，自从一个单独的脉冲被用来产生错误信号后者对单脉冲技术已经不是一个问题。单脉冲跟踪雷达既可以用天线反射器又可以用相控阵天线。图 11.7 显示了一种典

36、型的单脉冲天线的模型。A,B,C,D 四个波束描绘出四个圆锥形扫描波束的位置。这四个方向，大体上呈喇叭状，用来产生单脉冲天线的模式。振幅单脉冲处理器要求这四个信号有相同的相位但却有不同的名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 18 页，共 28 页 - - - - - - - - - 幅度。图 11.7 单脉冲天线模型一个解释振幅单脉冲技术的好方法就是用天线轴线的圆形中心来表示目标回波信号，正如图11.8a所示，在图中四个象限表示四个波束。在这种情况下，这四个喇叭接受相等的能

37、量，这就显示目标位于天线跟踪的轴线上。然而，当目标不在轴线上时（图11.8b-d） ，各喇叭中的能量就会不平衡。这个不平衡的能量用来产生差信号来驱动伺服控制系统。单脉冲处理包括计算天线模型的一个加两个不同的差（方位角和仰角）然后除去差信道的电压核信道的电压，信号的角度就可以确定了。图 11.8用图形解释单脉冲的概念(a) 目标在轴线上（b-d）目标不在轴线上雷达通过不断比较返回波束的幅度和相位来判断目标对轴线的位移。四个信号的相位相位在传输和接受模式都是连续的，这是关键性的。因此，任何名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - -

38、- 名师精心整理 - - - - - - - 第 19 页，共 28 页 - - - - - - - - - 数字电路或是微波比较电路都可以利用。图11.9 显示了一个典型的微波比较仪的方框图，图中三个接收通道分别被称为加通道，仰角差通道和方位角差通道。图 11.9 单脉冲比较器为了产生仰角差波束，我们可以用波束差（A-D ）或（B-C） 。然而，通过先形成和的形式 (A+B) 和(D+C)然后计算 (A+B) -(D+C)的差，我们获得一个更强的仰角差信号，el。同样，通过首先形成和的形式(A+B) 和(D+C) 然后计算(A+B) -(D+C) 的差，一个更强的方位角差信号，az，也产生了

39、。一个简单的单脉冲雷达的结构方框图如图11.10.所示。加通道被用作发射和接收。在接收机中加通道为其他两个不同的通道提供相位参考。加通道还可以用来测距。为了说阐明加通道和差通道的天线模型是如何形成的，我们将会假设一个sin/为天线的单独的一部分和斜角0。归一化的加信号（方位角或仰角）如下所示名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 20 页，共 28 页 - - - - - - - - - 0000sin()sin()( )（11.7）图 11.10 简单比幅单脉冲雷达框图而归

40、一化的差信号为0000sin()sin()( )（11.8）MATLAB 函数“ mono_pulse.m ”函数“mono_pulse.m ” 的功能是实现等式 (11.7)和(11.8)。它的输出包括和模式和差模式天线方向图以及差和比。函数形式为：mono_pulse(phi0) 其中phi0 为斜角，单位是弧度。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 21 页，共 28 页 - - - - - - - - - MATALB 函数 mono_pulse.m程序如下：fun

41、ction mono_pulse(phi0) % phi0 为跟踪轴偏离波束轴的偏角eps = 0.0000001; angle = -pi:0.01:pi; y1 = sinc(angle + phi0); y2 = sinc(angle - phi0); ysum = y1 + y2; % 计算和模式天线方向图ydif = -y1 + y2; % 计算差模式天线方向图figure (1) plot (angle,y1,b,angle,y2,r); grid; xlabel (方位（弧度） ) ylabel (斜视模式 ) figure (2) plot(angle,ysum,k); gri

42、d; xlabel (方位（弧度） ) ylabel (和模式 ) figure (3) plot (angle,ydif,k); grid; 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 22 页，共 28 页 - - - - - - - - - xlabel (方位（弧度） ) ylabel (差模式 ) angle = -pi/4:0.01:pi/4; y1 = sinc(angle + phi0); y2 = sinc(angle - phi0); ydif = -y1 +

43、 y2; ysum = y1 + y2; dovrs = ydif ./ ysum; % 计算差和比figure(4) plot (angle,dovrs,k); grid; xlabel (方位（弧度） ) ylabel (电压增益 ) 图 11.11(a-c)是对应于斜角00.15rad的和模式和差模式天线方向图。图 11.12(a-c)和图 11.11类似，只是对应的斜角00.75rad。很显然，和差模式非常依赖于斜角。用一个相对小的斜角比用一个大的角度能产生一个更好的和模式。并且，小斜角情况下差模式的斜率更陡峭。名师资料总结 - - -精品资料欢迎下载 - - - - - - - -

44、- - - - - - - - - - 名师精心整理 - - - - - - - 第 23 页，共 28 页 - - - - - - - - - 图 11.11a 两种斜模式斜角00.15rad图 11.11b 和图 11.11a对应的和模式名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 24 页，共 28 页 - - - - - - - - - 图 11.11c 和图 11.11a对应的差模式图 11.12a 两种斜模式斜角00.75rad名师资料总结 - - -精品资料欢迎下载

45、- - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 25 页，共 28 页 - - - - - - - - - 图 11.12b 和图 11.12a对应的和模式图 11.12c 和图 11.12a对应的差模式差通道给我们指明了目标是否在跟踪轴上。然而，信号的振幅不仅依靠于目标的角度位置，它也跟目标的距离和有效散射面积有关。由于这个原因，名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 26 页，共 28 页 - - -

46、- - - - - - 比值/被用来精确的估计误差角度，这个误差角度仅仅依靠目标的角度位置。让我们来看一看怎样计算这个误差角度，首先来看一下方位角误差角度。我们定义信号1S和2S为1SAD（11.9）2SBC（11.10）和信号是12SS。方位角差信号是12azSS。如果12SS，那么这两个通道有相同的相位o0(因为和通道被用来做相位参考)。相反的，如果12SS，那么这两个通道相位相差o180。可对仰角通道做同样的分析，信号相应的为1SAB和2SDC。因此，误差信号输出为| |cos| |（11.11）这里代表和通道和差通道之间的相位角，它等于o0或o180。如果0，那么目标在跟踪轴上； 否则

47、它就不在跟踪轴上。 图 11.13(a,b)给出了在图 11.11和 11.12所示的单脉冲雷达和通道和差通道的比值/图名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 27 页，共 28 页 - - - - - - - - - 图 11.13a 对应于图 11.11a的差和比图 11.13b 对应于图 11.12a的差和比名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 28 页，共 28 页 - - - - - - - - -