基于MATLAB的直线天线阵的仿真

基于MATLAB的直线天线阵的仿真，（附带仿真图），时间为2018年4月21日。

%************  天线阵代码   ****************

%%  程序界面清空
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