%%%%%%%%%%%%% 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)');
评论0