%%%%%%%%%%%%%%%%%% 相控陣的發射方向圖
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]);