相控阵雷达的发射方向图.m

相控阵雷达的重要性就不用多说了。本人仔细编写了相控阵雷达的仿真波束图

%%%%%%%%%%%%%%%%%% 相控阵的发射方向图
clear,clc;
fi=1:90;
fi0=0;
sita=1:180;
sita0=90;
PI=pi/180;
M=16;            %%%%%%%%%%%%%%% M行
N=16;            %%%%%%%%%%%%%%% N列

%%%%%%%%%%%%%%%% 行阵发射方向图
I_m=chebwin(M,20); %%%%%%%%%%%%% 20dB切比雪夫加权
f_fi=zeros(M,90);
for i=1:M
    f_fi(i,:)=I_m(i,1)*exp(j*pi*(i-1)*(sind(fi)-sind(fi0)));
end
for m=1:90
%     F_fi(1,j)=10*log10(sum(f_fi(:,j)));
    F_fi(1,m)=sum(f_fi(:,m));
end
subplot(2,1,1);
plot(fi,F_fi);

%%%%%%%%%%%%%%%% 列阵发射方向图
I_n=chebwin(N,20); %%%%%%%%%%%%% 20dB切比雪夫加权
f_sita=zeros(N,180);
for i=1:N
    f_sita(i,:)=I_n(i,1)*exp(j*pi*(i-1)*(cosd(sita)-cosd(sita0)));
end
for m=1:180
%     F_fi(1,j)=10*log10(sum(f_fi(:,j)));
    F_sita(1,m)=sum(f_sita(:,m));
end
subplot(2,1,2);
plot(sita,F_sita);

%%%%%%%%%%%%%%%% 整个阵面发射方向图
F_fi_sita=F_fi.'*F_sita;
F_fi_sita=abs(F_fi_sita);
F_fi_sita=F_fi_sita/max(max(F_fi_sita));
F_fi_sita=20*log10(F_fi_sita);
for i1=1:90,
    for i2=1:180,
        if F_fi_sita(i1,i2)<-60,
            F_fi_sita(i1,i2)=-60;
        end
    end
end        
figure(2);
% surf(fi,sita,F_fi_sita);
% [X,Y]=meshgrid(fi,sita);
meshz(F_fi_sita);
axis([0 180 0 90 -60 0]);