%%%%%%%%%%%%% 10*10平面陣的指向性  %%%%%%%%%%%%%
clear 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=30*pi/180;%入射信號的方位角－45～45，與平面陣的夾角
faii=0*pi/180;%入射信號的俯仰角0～90，與水平夾角
%%% 正弦信號入射
% L=2000;
% Ln=1:L;
% ss=exp(j*(2*pi*f/fs*Ln+2*pi*rand(1,L)));

for fai=-45:45
    for thet=0:90
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+46,thet+1)=Q;
    end
end
 Pvm=(abs(Pm)/max(max(abs(Pm))));
x=-45:1:45;
y=0:1:90;
figure;
mesh(y,x,Pvm);grid on;
%axis([0 90 -45 45 -60 0]);
ylabel('方位角(deg)');
xlabel('俯仰角(deg)');