%%%%%%%%%%%%% 10*10平面陣的指向性  %%%%%%%%%%%%%
clear all;
close all;
clc;
f=60*10^3;
fs=5*f;
w=2*pi*f;
c=1500;
lmd=c/f;
d=lmd/2;
%d=0.015;
theti=90*pi/180;%入射信號的方位角－45～45，與平面陣的夾角
faii=0*pi/180;%入射信號的俯仰角0～90，與水平夾角
m=5;
for fai=0:90
    for thet=0:180
for i=1:10
    for k=1:10
        gx=exp(-j*(i-1)*2*pi*d/lmd*(cos(thet*pi/180)*cos(fai*pi/180)-cos(theti)*cos(faii)));
        gz=exp(-j*(k-1)*2*pi*d/lmd*(sin(fai*pi/180)-sin(faii)));
G(i,k)=gx*gz;
    end
end
P=sum(G);
Q=sum(P);
Pm(fai+1,thet+1)=Q;
    end
end
 Pvm=20*log10(abs(Pm)/max(max(abs(Pm))));
for fai=0:90
    for thet=0:180
 if Pvm(fai+1,thet+1)<(-40)
     Pvm(fai+1,thet+1)=-40;
 end
    end
end
     
x=0:180;
y=0:1:90;
figure;
%plot(x,t);
mesh(x,y,Pvm);grid on;
%axis([0 180 0 90 -60 0]);
xlabel('方位角(deg)');
ylabel('俯仰角(deg)');