m=8;
j=sqrt(-1);
ad=exp(-j*pi*[0:m-1]'*sin(30*pi/180)); % 期望方向陣列響應(yīng)矢量
wop=ad;
thetas=[-90:90];
tm=thetas*pi/180;
am=exp(-j*pi*[0:m-1]'*sin(tm));
A=abs(wop'*am);  %陣列響應(yīng)
A=A/max(A);
figure,polar(tm,A)
A=10*log10(A);  %對數(shù)圖
hold on,title('歸一化陣列響應(yīng)幅值極坐標圖，八陣元')
figure,plot(thetas,A);
hold on,title('八陣元')
hold on,xlabel('入射角/度')
hold on,ylabel('歸一化 A=10*log10(A);')
grid on 
hold on,plot(-45,-160:0.1:0,'r');
hold on,plot(0,-30:0.1:0,'r');
hold on,plot(30,-30:0.1:0,'r');
hold on,plot(60,-30:0.1:0,'r');