m=8;
j=sqrt(-1);
ad=exp(-j*pi*[0:m-1]'*sin(30*pi/180)); % array response vector in the direction of desired signal.期望方向陣列響應矢量
wop=ad;
thetas=[-90:90];
tm=thetas*pi/180;
am=exp(-j*pi*[0:m-1]'*sin(tm));
A=abs(wop'*am);  %array response陣列響應
A=A/max(A);
figure,polar(tm,A)
A=10*log10(A);  %log figure對數圖
hold on,title('歸一化陣列響應幅值極坐標圖，八陣元')
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,-160:0.1:0,'r');
hold on,plot(30,-160:0.1:0,'r');
hold on,plot(60,-160:0.1:0,'r');