s6赛季-

2021年1月21日发(作者：钱朝鼎)

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.*(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(W,y1) grid on; %
绘出
N=4
等幅等矩阵列的归一化

阵因子

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

subplot(312);plot(W,y2) grid on; %
绘出
N=8
等幅等矩阵列的归一化

阵因子

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

subplot(313);plot(W,y3) grid on; %
绘出
N=12
等幅等矩阵列的归一

化阵因子

xlabel('f=3GHz,N=12,d=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(t,K1);xlab el('f=3GHz,N=4,d=1/2
波长,a=π/2');

subpl ot(132);polar(t,K2);xlabel('f=3GHz,N=8,d=1/2
波长 ,a=π/2');

,a=π/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
波长
,a=0');
subplot( 132);polar(t,K5);xlabel('f=3GHz,N=10,d=1/4
波长,a =π/2');



subplot(133);polar(t,K6 );xlabel('f=3GHz,N=10,d=1/4
波
长,a=π/2+π/10') ;

%-------------------------------------- -----------------------------
-
--
%
只有阵列单元间隔
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=π/2');
s ubplot(132);polar(t,K8);xlabel('f=3GHz,N=20,d=1/2< br>波长
,
a=π/2');
subplot(133);polar(t,K9 );xlabel('f=3GHz,N=20,d=0.7
波长,a=π/2');
%----------------------------------------------- --------------------
-
------
%------- --------------------------------3D-
天线方向图

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

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

[tehta,p hi]=meshgrid(eps:pi./(n_tehta-1):pi,... %meshgrid
函数为矩形
区域的设定范围是
epf
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);
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(5)
surf(X1,Y1,Z1); %
三维绘图函数
surf
，采用伪彩色表示曲面的高度
camlight
right
light
colorbar
axis image
rotate3D on
title('f=3GHz,N=4,d=1/2
波长,a=π/2'); figure(6)

surf(X2,Y2,Z2);


camlight right
light
colorbar
axis image
rotate3D on
title('f=3GHz,N=8,d=1/2
波长,a=π/2'); fi
gure(7)
surf(X3,Y3,Z3)
camlight right
light
colorbar
axis image
rotate3D on
title('f=3GHz,N=12,d=1/2
波长,a=π/2');
< br>%----------------------------------------------
%--------------------
只有阵列单元方向角
a
改变的 天线方向图

N4=10;d1=lambda/4;a1=0;a2=pi/2;a3=pi/2+pi/10;
A_3d=a1+(beta.*d1.*cos(t3));
x_1=(N4/2).*(A_3d);
x_2=(1/2).*(A_3d);
A_1=sin(x_1)./(N4.*sin(x_2)); %
非归一化的阵因子
K4
K_4=abs(A_1);
radio_4 =K_4;
X4=radio_4.*sin(tehta).*cos(phi); Y4=radio_4.*sin(tehta).*sin(phi);
Z4=radio_4.*cos(tehta);

s6赛季-

s6赛季-

s6赛季-

s6赛季-

s6赛季-

s6赛季-

s6赛季-

s6赛季-