clear?all;
close?all;
clc
N=3;
M=3;%陣元個數
f=300000000;%信號頻率
c=300000000;
lambda=c/f;
d=0.5*lambda;
theta=linspace（0360361）;?
phi=linspace（0360361）;
phi0=30;???%俯仰角
theta0=180;?%預先設定的方向，方位角
k?=?1;
p?=?zeros（3M*N）;%3*9的矩陣
for?iii?=?1?:?N
????for?jjj?=?1:?M
????????p（:k）?=?[iii*djjj*d0]‘;
????????k?=?k+1;
????end
end???%陣元的位置
%p=[?0.03????0.03????0.03????0.06????0.06????0.06????0.09????0.09????0.09???
%????0.03????0.06????0.09????0.03????0.06????0.09????0.03????0.06????0.09
%????0???????0???????0???????0???????0???????0???????0???????0???????0????]

figure（1）;
plot3（p（1:）p（2:）p（3:）‘ko‘）;
hold?on;
xlabel（‘\it?x‘）;
ylabel（‘\it?y‘）;
zlabel（‘\it?z‘）;


k0=-2*pi/lambda*[sin（phi0*pi/180）*cos（theta0*pi/180）sin（phi0*pi/180）*sin（theta0*pi/180）cos（phi0*pi/180）].‘;
w0=exp（1i*k0.‘*p）;
for?ii=1:length（theta）
????for?jj=1:length（phi）
????????k=-2*pi/lambda*[sin（phi（jj）*pi/180）*cos（theta（ii）*pi/180）sin（phi（jj）*pi/180）*sin（theta（ii）*pi/180）cos（phi（jj）

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????1458??2018-05-29?15:38??方向圖\mianzhen.m

?????文件???????1279??2018-05-31?07:38??方向圖\UCL.m

?????文件???????1217??2018-05-08?15:37??方向圖\ULA.m

?????目錄??????????0??2018-07-02?13:05??方向圖

-----------?---------??----------?-----??----

?????????????????3954????????????????????4