《对称阵子天线matlab》由会员分享，可在线阅读，更多相关《对称阵子天线matlab（6页珍藏版）》请在金锄头文库上搜索。
1、方向图clearlamda=10;%自由空间的波长n=1.125;l=n*lamda;k=2*pi/lamda;%自由空间的相移常数theta0=0.0001:0.1:360; theta=theta0*pi/180;for i=1:length(theta0) fe(i)=abs(cos(k*l*cos(theta(i)-cos(k*l)/sin(theta(i);%fh(i)=1-cos(k*l);end% figurepolar(theta,fe/max(fe);xlabel(l=1.125 入)；% figure% polar(theta,fh/max(fh);辐射阻抗clearlam
2、da=10;%自由空间的波长n=0.1:0.01:0.9;l=n*lamda;k=2*pi/lamda;%自由空间的相移常数gama=0.5772156;for i=1:length(n)Rr1(i)=sin(2*k*l(i)*(Si(4*k*l(i)-2*Si(2*k*l(i); Rr2(i)=cos(2*k*l(i)*(gama+log(k*l(i)+Ci(4*k*l(i)-2*Ci(2*k*l(i);Rr3(i)=2*(gama+log(2*k*l(i)-Ci(2*k); Rr(i)=(Rr1(i)+Rr2(i)+Rr3(i)*30;endplot(n,Rr)阻抗特性clcclearla
3、mda=100;%自由空间的波长n=0.1:0.01:0.9;l=n*lamda;a=30,300,3000;%l 和 a 的比值k=2*pi/lamda;%自由空间的相移常数gama=0.5772156;for j=1:length(a) for i=1:length(n)Rr1(i)=sin(2*k*l(i)*(Si(4*k*l(i)-2*Si(2*k*l(i);Rr2(i)=cos(2*k*l(i)*(gama+log(k*l(i)+Ci(4*k*l(i)-2*Ci(2*k*l(i);Rr3(i)=2*(gama+log(2*k*l(i)-Ci(2*k);Rr(i)=(Rr1(i)+Rr
4、2(i)+Rr3(i)*30;Xr1(i)=2*Si(2*k*l(i);Xr2(i)=sin(2*k*l(i)*(gama+log(k*l(i)+Ci(4*k*l(i)-2*Ci(2*k*l(i)-2*log(a(j); Xr3(i)=cos(2*k*l(i)*(2*Si(2*k*l(i)-Si(4*k*l(i);Xr(i,j)=(Xr1(i)+Xr2(i)+Xr3(i)*30;Zca(i)=120*(log(2*(a(j)-1); afa(i)=Rr(i)/(Zca(i)*l(i)*(1-sin(2*k*l(i)/(2*k*l(i);Rin(i,j)=Zca(i)*(sinh(2*afa(i
5、)*l(i)-afa(i)/k*sin(2*k*l(i)/(cosh(2*afa(i)*l(i)-cos(2* k*l(i);Xin(i,j)=-Zca(i)*(afa(i)/k*sinh(2*afa(i)*l(i)+sin(2*k*l(i)/(cosh(2*afa(i)*l(i)-cos(2 *k*l(i);endendfigureplot(n,Rr)%辐射电阻title (辐射电阻随1/入的变化曲线)；xlabel(l/ 入);y1abe1(Rr (Q);figurep1ot(n,Xr)%辐射电抗tit1e(1/a不同时辐射电抗随1/入的变化曲线)；x1abe1(1/ 入);y1abe1(
6、Xr (Q);figurep1ot(n,Rin)%输入电阻tit1e(不同平均特性阻抗时输入电阻随1/入的变化曲线)；x1abe1(1/ 入);y1abe1(Rin (Q);figurep1ot(n,Xin)%输入电抗tit1e(不同平均特性阻抗时电抗随1/入的变化曲线);x1abe1(1/ 入);y1abe1(Xin (Q);对称阵子c1cc1ear1amda=10;%自由空间的波长n=0.25;1=n*1amda;k=2*pi/1amda;%自由空间的相移常数theta0=0.0001:0.1:360; theta=theta0*pi/180;for i=1:length(theta0)
7、fe(i)=abs(cos(k*l*cos(theta(i)-cos(k*l)/sin(theta(i); fh(i)=1-cos(k*l);endfigure polar(theta,fe/max(fe);figure polar(theta,fh);用 MATLAB 语言编程计算半波对称振子场图，电流分布和输入阻抗 n=2*pi;d=n/50+0.001;x,y = meshgrid(-n:d:n,-n:d:n);nn=30j=0;M = moviein(nn); for j=1:nnz=(x./sqrt(x.A2+y42).人2.*(. -sin(sqrt(x42+y.A2)-j*2*p
8、i/nn).-cos(sqrt(x.A2+y.A2)-j*2*pi/nn)./sqrt(x.A2+y.A2).);v=-1 -0.9 -0.8 -0.75 -0.5 -0.25 -0.1 -0.05 0 0.05 0.1 0.25 0.5 .0.75 0.8 0.9 1;%v= -0.75 -0.5 -0.25 0 0.25 0.5 0.75 ;contour(z,v);%contour(z,30);%colorbaraxis equal%grid onM(:,j) = getframe; endmovie(M,100,10) clear allN= 10;d= 0.5;alpha = 2*p
9、i;dip=1:N; kd=2*pi*d;phi=0:2*pi/800.1:2*pi; framemax = 48;M = moviein(framemax); set(gcf,Position,100 100 640 480) for n=1:framemaxu=kd*cos(phi)-alpha/framemax*(n-1); F=abs(sin(N.*u./2+N*1e-15)./sin(u./2+1e-15);plot(F.*cos(phi),F.*sin(phi),r,LineWidth,3); %polar(phi,abs(F);title( num2str(N), element
10、s, .,num2str(d), lambda apart,fontsize, 18)xlabel(alpha= ,num2str(alpha/framemax*(n-1),Color,k,fontsize, 18) hold onplot(N.*cos(phi),N.*sin(phi),b,LineWidth,2); plot(dip-N/2-0.5,dip*0,o,linewidth,3) hold offaxis equalM(:,n) = getframe(gcf); end clf reset set(gcf,Position,100 100 800 600) axis off mo
11、vie(M,1,6)clear allN= 10;d= 0.5;alpha = 2*pi;dip=1:N; kd=2*pi*d;phi=0:2*pi/800.1:2*pi; framemax = 48;M = moviein(framemax); set(gcf,Position,100 100 640 480)for n=1:framemax u=kd*cos(phi)-alpha/framemax*(n-1); F=abs(sin(N.*u./2+N*1e-15)./sin(u./2+1e-15); plot(F.*cos(phi),F.*sin(phi),r,LineWidth,3);
12、%polar(phi,abs(F);title( num2str(N), elements, .,num2str(d), lambda apart,fontsize, 18) xlabel(alpha= ,num2str(alpha/framemax*(n-1),Color,k,fontsize, 18) hold onplot(N.*cos(phi),N.*sin(phi),b,LineWidth,2); plot(dip-N/2-0.5,dip*0,o,linewidth,3) hold offaxis equalM(:,n) = getframe(gcf);end clf reset s
13、et(gcf,Position,100 100 800 600) axis off movie(M,1,6) clear;clc; range=1/4;point=6; wave_length_value=1;measurement=7.022*(10)人(-3);v0=1;k=(2*pi)/wave_length_value; a=measurement*wave_length_value;l=range*wave_length_value; step=l/(point-1);i_point=1:point; % 测试点赋值 matrix_wavelength(i_point)=step*(
14、i_point-1);z=linspace(-l,l,100); % 积分离散for i_point=1:point; % 求解 A 向量 r=(matrix_wavelength(i_point)-z).人2 + aA2).A(1/2); % 场源距离离散 g=exp(-j*k*r)./r; % 格林函数离散A(i_point)=trapz(z,sin(k*(1/4-abs(z).*g); % A 元素确定endfor i_point=1:point; % 求解 B 向量r=(matrix_wavelength(i_point)-z).人2 + aA2).A(1/2); % 场源距离离散 g=exp(-j*k*r)./r; % 格林函数离散B(i_point)=trapz(z,sin(2*k*(1/4-abs(z).*g); % B 元素确定endfor i_point=1:point; % 求解 C 向量 r=(matrix_wavelength(i_point)-z).A2 + aA2).A(1/2); % 场源距离离散 g=exp(-j*k*r)./r; % 格林函数离散C(i_point)=trapz(z,sin(3*k*(1/4-abs(z).*g); % B 元素确定endfor i_point=1:point; % 求解 D
