matlab仿真天线阵代码

1、天线阵代码一、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、./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等幅等矩阵列

3、的归一化阵因子 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.*

4、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改变的天线

5、方向图 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

6、/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;

7、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).*(

8、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); subplot(132);polar(t,k8);xlabel(f=3ghz,n=20,d=1/2波长,a=/2); subplot(133);polar(t,k9);xlabel(f=3ghz,n=20,d=0.7波长,a=/2); %- %-3d-天线方向图 n_tehta = 130; %- 采样

9、视角点的仰角 n_phi = 130; %-采样点的方向角 tehta,phi=meshgrid(eps:pi./(n_tehta-1):pi,. %meshgrid函数为矩形区域的设定范围是epftehta 0phi2 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

10、 =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=radi

11、o_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，采用伪彩色表示曲面的高度 camli

12、ght 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); figure(7) surf(x3,y3,z3) camlight right light colorbar axis image rotate3d on title(f=3ghz,n=12,d=1/2波长,a=/2)

13、; %- %-只有阵列单元方向角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);

14、%- b_3d=a2+(beta.*d1.*cos(t3); y_1_3d=(n4/2).*(b_3d); y_2_3d=(1/2).*(b_3d); b_1=sin(y_1_3d)./(n4.*sin(y_2_3d); %非归一化的阵因子k5 k_5=abs(b_1); radio_5 =k_5; x5=radio_5.*sin(tehta).*cos(phi); y5=radio_5.*sin(tehta).*sin(phi); z5=radio_5.*cos(tehta); %- c_3d=a3+(beta.*d1.*cos(t3); v_1=(n4/2).*(c_3d); v_2=(1

15、/2).*(c_3d); c_1=sin(v_1)./(n4.*sin(v_2); %非归一化的阵因子k6 k_6=abs(c_1); radio_6 =k_6; x6=radio_6.*sin(tehta).*cos(phi); y6=radio_6.*sin(tehta).*sin(phi); z6=radio_6.*cos(tehta); %- figure(8) surf(x4,y4,z4); %三维绘图函数surf，采用伪彩色表示曲面的高度 camlight right light colorbar axis image rotate3d on title(f=3ghz,n=10,d

16、=1/4波长,a=0); figure(9) surf(x5,y5,z5); camlight right light colorbar axis image rotate3d on title(f=3ghz,n=10,d=1/4波长,a=/2); figure(10) surf(x6,y6,z6) camlight right light colorbar axis image rotate3d on title(f=3ghz,n=10,d=1/4波长,a=/2+/10); %- %-只有阵列单元间隔d改变的天线方向3d图 n5=20;d2=lambda/4;d3=lambda/2;d4=0

17、.7*lambda;a4=pi/2; d_3d=a4+(beta.*d2.*cos(t3); p_1=(n5/2).*(d_3d); p_2=(1/2).*(d_3d); d_1=sin(p_1)./(n5.*sin(p_2); %非归一化的阵因子k7 k_7=abs(d_1); radio_7 =k_7; x7=radio_7.*sin(tehta).*cos(phi); y7=radio_7.*sin(tehta).*sin(phi); z7=radio_7.*cos(tehta); %- e_3d=a4+(beta.*d3.*cos(t3); q_1=(n5/2).*(e_3d); q_

18、2=(1/2).*(e_3d); e_1=sin(q_1)./(n5.*sin(q_2); %非归一化的阵因子k8 k_8=abs(e_1); radio_8 =k_8; x8=radio_8.*sin(tehta).*cos(phi); y8=radio_8.*sin(tehta).*sin(phi); z8=radio_8.*cos(tehta); %- f_3d=a4+(beta.*d4.*cos(t3); r_1_3d=(n5/2).*(f_3d); r_2_3d=(1/2).*(f_3d); f_1=sin(r_1_3d)./(n5.*sin(r_2_3d); %非归一化的阵因子k9

19、 k_9=abs(f_1); radio_9 =k_9; x9=radio_9.*sin(tehta).*cos(phi); y9=radio_9.*sin(tehta).*sin(phi); z9=radio_9.*cos(tehta); %- figure(11) surf(x7,y7,z7); %三维绘图函数surf，采用伪彩色表示曲面的高度 camlight right light colorbar axis image rotate3d on title(f=3ghz,n=20,d=1/4波长,a=/2); figure(12) surf(x8,y8,z8); camlight ri

20、ght light colorbar axis image rotate3d on title(f=3ghz,n=20,d=1/2波长,a=/2); figure(13) surf(x9,y9,z9) camlight right light colorbar axis image rotate3d on title(f=3ghz,n=20,d=0.7波长,a=/2);二、%-均匀直线阵列天线的应用之一：边射阵 clc clear all f=3e10; %30ghz,厘米波 i=1; %天线电流值 lambda=(3e8)/f; %lambda=c/f 波长 beta=2.*pi/lambd

21、a; n=15; t=0:0.01:2*pi; d1=lambda/4; %没有栅瓣效应的边射阵，即间隔d波长 w2=beta.*d2.*cos(t); %定义kdcos(方向角) z3=(n/2).*w2; z4=(1/2).*w2; f2=sin(z3)./(n.*sin(z4); k2=abs(f2); figure(1) subplot(121);polar(t,k1);title(边射阵 f=30ghz,n=15,d=1/4波长); subplot(122);polar(t,k2);title(边射阵（有栅瓣） f=30ghz,n=15,d=1.5倍波长); 三、%-均匀直线阵列天线的应用之二：普通端射阵 clc clear all f=3e10; %30ghz,厘米波 i=1; %天线电流值 lambda=(3e8)/f; %lambda=c/f 波长 beta=2.*pi/lambda; n=15; t=0:0.01:2*pi; d1=lambda/4; %没有栅瓣效应的普通端射阵，即间隔d1/2波长 w2=beta.*d2.*cos(t)