矩形阵列指向性MATLAB程序

clear

M=64;c=340;%ÕóÔª×ÜÊý

theta=-2:0.02:2;%ÈëÉäÑö½Ç

N=length(theta);

phi=-2:0.02:2;%ÈëÉä·½Ïò½Ç

omiga=1000;%ÈëÉä½ÇƵÂÊ

d=0.25

for k=1:N

for m=1:N

for n=0:7

for l=0:7

t(n+1,l+1)=omiga*(n*d*(cos(phi(m))*sin(theta(k))-1/2)+l*d*(sin(phi(m))*sin(theta(k))-1/2))/c;%µÚiºÅÕóÔªÏà¶ÔÕóÁÐÖÐÐĵÄʱ¼äÑÓ³Ù

A(n+1,l+1)=exp(-1i*t(n+1,l+1));

end

end

B(k,m)=1/M*sum((sum(A)));

end

end

[xx,yy]=meshgrid(theta,phi);

surf(xx,yy,abs(B));

title('1000Hz入射方向图');

xlabel('入射仰角 x/m');

ylabel('入射方向角 y/m');

figure

contour(B,8);

title('1000Hz 等高线图');

xlabel('入射仰角 x/m');

ylabel('入射方向角 y/m');