www.pudn.com > mainamplitudecomparison.rar > mainamplitudecomparison.m, change:2017-04-13,size:1087b


 
clc 
clear all 
close all 
 
phi0 = 10*pi/180; %指向间隔 
angle = -pi:0.01:pi;  %成像方位角范围 
y1 = sinc((angle+phi0));  %阵列一sinc方向图 
y2 = sinc((angle-phi0));  %阵列二方向图 
ydif = y2 - y1;  %差通道 
ysum = y1 + y2;  %和通道 
dovrs = ydif./ysum;  %和差比幅曲线 
angle_xx = angle/pi*180;  %成像方位角范围 
 
% figure 
% plot(angle_xx,y1,'k',angle_xx,y2,'k'); 
% grid; 
% xlabel('Angle-radians') 
% ylabel('Squintedpattern') 
% figure 
% plot(angle_xx,ysum,'k'); 
% grid; 
% xlabel('Angle-radians') 
% ylabel('Sumpatterns') 
% figure 
% plot(angle_xx,ydif,'k'); 
% grid; 
% xlabel('Angle-radians') 
% ylabel('Differencepattern') 
 
 
figure 
plot(angle_xx,dovrs,'k') 
grid; 
xlabel('Angle-radians') 
axis([-40 40 -0.7 0.7]) 
ylabel('voltagegain') 
 
phi0 = 10*pi/180; %指向间隔 
theta = 10*pi/180; 
angle = -pi:0.01:pi;  %成像方位角范围 
y1 = sinc((angle+theta+phi0));  %阵列一sinc方向图 
y2 = sinc((angle+theta-phi0));  %阵列二方向图 
ydif = y2 - y1;  %差通道 
ysum = y1 + y2;  %和通道 
k = ydif./ysum; 
est_theta = k/ 0.0112; 
figure 
plot(angle_xx,est_theta);