MATLAB仿真天线阵代码

clc clear all f=3e9;

N1=4;N2=8;N3=12; a=pi/2; % 馈电相位差 i=1; % 天线电流值

lambda=(3e8)/f; %lambda=c/f d=lambda/2;

beta=2 、 *pi/lambda;

W=-2*pi:0 、 001:2*pi; y1=sin((N1 、 *W 、 /2)) 、

/(N1 、 y1=abs(y1);

r1=max(y1);

y2=sin((N2 、 *W 、 /2)) 、

/(N2 、 y2=abs(y2);

r2=max(y2);

y3=sin((N3 、 *W 、 /2)) 、

/(N3 、

y3=abs(y3);

r3=max(y3);

%归一化阵因子绘图程序

figure(1) subplot(311);plot(W,y1) ; grid on; % 阵因子

xlabel('f=3GHz,N=4,d=1/2 波长,a= n /2') subplot(312);plot(W,y2) ; grid on; % 阵因子

xlabel('f=3GHz,N=8,d=1/2 波长,a= n /2')

subplot(313);plot(W,y3) ; grid on; % 化阵因子 xlabel('f=3GHz,N=12,d=1/2 波长,a= n /2')

% --------------- %只有参数N 改变的天线方向

图 t=0:0 、0 1 :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);

% ---------------

天线阵代码

波长

*(sin(W 、/2))); % *(sin(W 、/2))); % *(sin(W 、/2))); % 归一化阵因子 归一化阵因子 归一化阵因子 绘出N=4等幅等矩阵列的归一化 绘出N=8等幅等矩阵列的归一化 绘出N=12等幅等矩阵列的归一

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(t,K1);xlabel('f=3GHz,N=4,d=1/2 subplot(132);polar(t,K2);xlabel('f=3GHz,N=8,d=1/2 subplot(133);polar(t,K3);xlabel('f=3GHz,N=12,d=1/2 % %只有阵列单元方向角 a 改变的天线方向图

N4=10;d1=lambda/4;a1=0;a2=pi/2;a3=pi/2+pi/10;

A=a1+(beta 、*d1 、*cos(t)) ; x1=(N4/2) 、*(A); x2=(1/2) 、*(A);

A1=sin(x1) 、/(N4 、*sin(x2)); % 非归一化的阵因子 K4 K4=abs(A1);

% -------------------

B=a2+(beta 、*d1、*cos(t));

y_1=(N4/2) 、*(B);

y_2=(1/2) 、*(B);

B1=sin(y_1) 、/(N4 、 *sin(y_2)) ; % 非归一化的阵因子 K5 K5=abs(B1); % -------------------

C=a3+(beta 、*d1、*cos(t));

v1=(N4/2) 、*(C);

v2=(1/2) 、*(C);

C1=sin(v1) 、/(N4 、*sin(v2)); % 非归一化的阵因子 K6 K6=abs(C1);

% ------------------- 绘图函数

figure(3)

subplot(131);polar(t,K4);xlabel('f=3GHz,N=10,d=1/4 subplot(132);polar(t,K5);xlabel('f=3GHz,N=10,d=1/4 subplot(133);polar(t,K6);xlabel('f=3GHz,N=10,d=1/4 长,a= n /2+ n /10');

% ------------------------------------------------

波长,a= n /2'); 波长,a= n /2'); 波长,a= n 波长,a=0'); 波长,a= n /2');

波

%只有阵列单元间隔 d 改变的天线方向图

N5=20;d2=lambda/4;d3=lambda/2;d4=0 、7*lambda;a4=pi/2;

D=a4+(beta 、*d2 、*cos(t));

p1=(N5/2) 、*(D);

p2=(1/2) 、*(D);

D1=sin(p1) 、/(N5 、*sin(p2)); % 非归一化的阵因子 K7

K7=abs(D1);

% ----------------------

E=a4+(beta 、*d3、*cos(t));

q1=(N5/2) 、*(E);

q2=(1/2) 、*(E);

E1=sin(q1) 、/(N5 、*sin(q2)); % 非归一化的阵因子 K8 K8=abs(E1);

% ----------------------

F=a4+(beta 、*d4、*cos(t));

r_1=(N5/2) 、*(F);

r_2=(1/2) 、*(F);

F1=sin(r_1) 、/(N5 、 *sin(r_2)); % 非归一化的阵因子 K9 K9=abs(F1); % ---------------- 绘图函数

figure(4)

subplot(131);polar(t,K7);xlabel('f=3GHz,N=20,d=1/4 波长,a= n /2'); subplot(132);polar(t,K8);xlabel('f=3GHz,N=20,d=1/2

波长,a= n /2'); subplot(133);polar(t,K9);xlabel('f=3GHz,N=20,d=0

% ------------------------------------------------ % ---------------------------- 3D- 天线方向图

n_tehta = 130; % ------------- 采样视角点的仰角

n_phi = 130; % --------------- 采样点的方向角

[tehta,phi]=meshgrid(eps:pi 、/(n_tehta-1):pi, 、、、 %meshgrid

为矩形区域的设定范围就是epfvtehtav n 0

0: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);

7 波长,a= n /2');

函数