matlab仿真天线阵代码

matlab仿真天线阵代码

天线阵代码

tail750575.html

一、

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

subplot(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)