MATLAB仿真天线阵代码.pdf

天线阵代码 一、 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) ; gridon; %绘出N=4 等幅等矩阵列的归一化 阵因子 xlabel(f=3GHz,N=4,d=1/2 波长,a=π/2) subplot(312);plot(W,y2) ; gridon; %绘出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); subplot(133);polar(t,K3);xlabel(f=3GHz,N=12,d=1/2 波长,a=π/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); % -- %只有