%************??天線陣代碼???****************

%%??程序界面清空
clc??
clear?all??
close?all
%%??
f?=?3e9;???????%通信頻率為3*10^9HZ
N1?=?4;???????%對應4元均勻線陣
N2?=?8;???????%對應8元均勻線陣
N3?=?12;?????%對應12元均勻線陣
a?=?pi/2;??????%饋電相位差??
i?=?1;???????????%天線電流值?
c?=?3e8;??????%光速
lambda?=?c/f;????????%波長??
d?=?lambda/2;???????%相鄰單元之間的距離
beta?=?2.*pi/lambda;??
W?=?-2*pi:0.001:2*pi;??
y1?=?sin（（N1.*W./2））./（N1.*（sin（W./2）））???%歸一化陣因子??
y1?=?abs（y1）;??
r1?=?max（y1）;??
?
y2?=?sin（（N2.*W./2））./（N2.*（sin（W./2）））;????%歸一化陣因子??
y2?=?abs（y2）;??
?
r2?=?max（y2）;??
y3?=?sin（（N3.*W./2））./（N3.*（sin（W./2）））;???%歸一化陣因子??
y3?=?abs（y3）;??
r3?=?max（y3）;??
?
?%歸一化陣因子繪圖
figure（1）??
subplot（311）;??plot（Wy1）;??grid?on;?????%繪出N=4等幅等矩陣列的歸一化陣因子??
xlabel（‘f=3GHzN=4d=1/2波長a=π/2‘）??
subplot（312）;??plot（Wy2）;??grid?on;?????%繪出N=8等幅等矩陣列的歸一化陣因子??
xlabel（‘f=3GHzN=8d=1/2波長a=π/2‘）??
subplot（313）;??plot（Wy3）;??grid?on;????%繪出N=12等幅等矩陣列的歸一化陣因子??
xlabel（‘f=3GHzN=12d=1/2波長a=π/2‘）??
?
%只有參數N改變的天線方向圖??
t?=?0:0.01:2*pi;???
W?=?a+（beta.*d.*cos（t））;??
z1?=?（N1/2）.*（W）;??
z2?=?（1/2）.*（W）;??
W1?=?sin（z1）./（N1.*sin（z2））;??%非歸一化的陣因子K1?
K1?=?abs（W1）;??
%------------------------------??
W?=?a+（beta.*d.*cos（t））;??
z3?=?（N2/2）.*（W）;??
z4?=?（1/2）.*（W）;??
W2?=?sin（z3）./（N2.*sin（z4））;??%非歸一化的陣因子K2??
K2?=?abs（W2）;??
%------------------------------??
W?=?a+（beta.*d.*cos（t））;??
z5?=?（N3/2）.*（W）;??
z6?=?（1/2）.*（W）;??
W3?=?sin（z5）./（N3.*sin（z6））;??%非歸一化的陣因子K3??
K3?=?abs（W3）;??

%非歸一化陣因子繪圖??
figure（2）??
subplot（131）;polar（tK1）;xlabel（‘f=3GHzN=4d=1/2波長a=π/2‘）;??
subplot（132）;polar（tK2）;xlabel（‘f=3GHzN=8d=1/2波長a=π/2‘）;??
subplot（133）;polar（tK3）;xlabel（‘f=3GHzN=12d=1/2波長a=π/2‘）;

n_tehta?=?130;???????????%采樣視角點的仰角??
n_phi?=?130;??????????????%采樣點的方向角????

%meshgrid函數為矩形區域的設定范圍是epf[tehtaphi]?=?meshgrid（eps:pi/（n_tehta-1）:pi0:2*pi/（n_phi-1）:2*pi）;?????
t3=tehta;??
%----?只有參數N改變的天線方向3D圖?---------??
M?=?a+（beta.*d.*cos（t3））;????
%---------------------------
N1?=?4;
N2?=?8;
N3?=12;??
?z_1?=?（N1/2）.*（M）;??
z_2?=?（1/2）.*（M）;??
M1?=?sin（z_1）./（N1.*sin（z_2））;??%非歸一化的陣因子K1??
K_1?=?abs（M1）;??
radio_1?=?K_1;???
X1?=?radio_1.*sin（tehta）.*cos（phi）;???
Y1?=?radio_1.*sin（tehta）.*sin（phi）;???
Z1?=?radio_1.*cos（tehta）;??
%-----------------------------------??
M?=?a+（beta.*d.*cos（t3））;??
?z_3?=?（N2/2）.*（M）;??
z_4?=?（1/2）.*（M）;??
M2?=?sin（z_3）./（N2.*sin（z_4））;???%非歸一化的陣因子K2??
K_2?=?abs（M2）;??
radio_2?=?K_2;???
X2?=?radio_2.*sin（tehta）.*cos（phi）;???
Y2?=?radio_2.*sin（tehta）.*sin（phi）;???
Z2?=?radio_2.*cos（tehta）;?
%-------------------------------------??
M?=?a+（beta.*d.*cos（t3））;??
?z_5?=?（N3/2）.*（M）;??
z_6?=?（1/2）.*（M）;??
M3?=?sin（z_5）./（N3.*sin（z_6））;??%非歸一化的陣因子K3?
K_3?=?abs（M3）;??
radio_3??=?K_3;???
X3?=?radio_3.*sin（tehta）.*cos（phi）;???
Y3?=?radio_3.*sin（tehta）.*sin（phi）;???
Z3?=?radio_3.*cos（tehta）;??

%3D繪圖函數??
figure（3）??
surf（X1Y1Z1）;???%三維繪圖函數surf采用偽彩色表示曲面的高度???
camlight?right??
light??

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2018-04-22?04:11??代碼\
?????文件???????43518??2017-04-17?22:32??代碼\2D_1.png
?????文件???????55822??2017-04-17?22:33??代碼\2D_2.png
?????文件??????150478??2017-04-17?22:27??代碼\3D_1.png
?????文件??????180840??2017-04-17?22:31??代碼\3D_11.png
?????文件??????158454??2017-04-17?22:27??代碼\3D_2.png
?????文件??????183956??2017-04-17?22:31??代碼\3D_22.png
?????文件??????201374??2017-04-17?22:27??代碼\3D_3.png
?????文件??????216234??2017-04-17?22:31??代碼\3D_33.png
?????文件????????3618??2017-04-17?22:37??代碼\Linear_array_antenna_logodo.m