对称阵子天线matlab

方向图clearlamda=10;%自由空间的波长n=1.125;l=n*lamda;k=2*pi/lamda;%自由空间的相移常数theta0=[0.0001:0.1:360];theta=theta0*pi/180;fori=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));辐射阻抗clearlamda=10;%自由空间的波长n=0.1:0.01:0.9;l=n*lamda;k=2*pi/lamda;%自由空间的相移常数gama=0.5772156;fori=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)阻抗特性clcclearlamda=100;%自由空间的波长n=0.1:0.01:0.9;l=n*lamda;a=[30,300,3000];%l和a的比值k=2*pi/lamda;%自由空间的相移常数gama=0.5772156;forj=1:length(a)fori=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;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)*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('辐射电阻随l/λ的变化曲线');xlabel('l/λ');ylabel('Rr(Ω)');figureplot(n,Xr)%辐射电抗title('l/a不同时辐射电抗随l/λ的变化曲线');xlabel('l/λ');ylabel('Xr(Ω)');figureplot(n,Rin)%输入电阻title('不同平均特性阻抗时输入电阻随l/λ的变化曲线');xlabel('l/λ');ylabel('Rin(Ω)');figureplot(n,Xin)%输入电抗title('不同平均特性阻抗时电抗随l/λ的变化曲线');xlabel('l/λ');ylabel('Xin(Ω)');对称阵子clcclearlamda=10;%自由空间的波长n=0.25;l=n*lamda;k=2*pi/lamda;%自由空间的相移常数theta0=[0.0001:0.1:360];theta=theta0*pi/180;fori=1:length(theta0)fe(i)=abs((cos(k*l*cos(theta(i)))-cos(k*l))/sin(theta(i)));fh(i)=1-cos(k*l);endfigurepolar(theta,fe/max(fe));figurepolar(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);forj=1:nnz=(x./sqrt(x.^2+y.^2)).^2.*(...-sin(sqrt(x.^2+y.^2)-j*2*pi/nn)...-cos(sqrt(x.^2+y.^2)-j*2*pi/nn)./sqrt(x.^2+y.^2)...);v=[-1-0.9-0.8-0.75-0.5-0.25-0.1-0.0500.050.10.250.5...0.750.80.91];%v=[-0.75-0.5-0.2500.250.50.75];contour(z,v);%contour(z,30);%colorbaraxisequal%gridonM(:,j)=getframe;endmovie(M,100,10)clearallN=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',[100100640480])forn=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),'elements,'...,num2str(d),'\lambdaapart'],'fontsize',18)xlabel(['\alpha=',num2str(alpha/framemax*(n-1))],'Color','k','fontsize',18)holdonplot(N.*cos(phi),N.*sin(phi),'b','LineWidth',2);plot(dip-N/2-0.5,dip*0,'o','linewidth',3)holdoffaxisequalM(:,n)=getframe(gcf);endclfresetset(gcf,'Position',[100100800600])axisoffmovie(M,1,6)clearallN=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',[100100640480])forn=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),'elements,'...,num2str(d),'\lambdaapart'],'fontsize',18)xlabel(['\alpha=',num2str(alpha/framemax*(n-1))],'Color','k','fontsize',18)holdonplot(N.*cos(phi),N.*sin(phi),'b','LineWidth',2);plot(dip-N/2-0.5,dip*0,'o','linewidth',3)holdoffaxisequalM(:,n)=getframe(gcf);endclfresetset(gcf,'Position',[100100800600])axisoffmovie(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*(i_point-1);z=linspace(-l,l,100);%积分离散fori_point=1:point;%求解A向量r=((matrix_wavelength(i_point)-z).^2+a^2).^(1/2);%场源距离离散g=exp(-j*k*r)./r;%格林函数离散A(i_point)=trapz(z,sin(k*(1/4-abs(z))).*g);%A元素确定endfori_point=1:point;%求解B向量r=((matrix_wavelength(i_point)-z).^2+a^2).^(1/2);%场源距离离散g=exp(-j*k*r)./r;%格林函数离散B(i_point)=trapz(z,sin(2*k*(1/4-abs(z))).*g);%B元素确定endfori_point=1:point;%求解C向量r=((matrix_wavelength(i_point)-z).^2+a^2).^(1/2);%场源距离离散g=exp(-j*k*r)./r;%格林函数离散C(i_point)=trapz(z,sin(3*k*(1/4-abs(z))).*g);%B元素确定endfori_point=1:point;%求解D向量r=((matrix_wavelength(i_point)-z).^2+a^2).^(1/2);%场源距离离散g=exp(-j*k*r)./r;%格林函数离散D(i_point)=trapz(z,sin(4*k*(1/4-abs(z))).*g);%B元素确定endfori_point=1:point;%求解E向量r=((matrix_wavelength(i_point)-z).^2+a^2).^(1/2);%场源距离离散g=exp(-j*k*r)./r;%格林函数离散E(i_point)=trapz(z,sin(5*k*(1/4-abs(z))).*g);%B元素确定endfori_point=1:point;%求解F向量F(i_point)=cos(k*matrix_wavelength(i_point));%B元素确定endimpedance_matrix=[A.',B.',C.',D.',E.',F.'];%阻抗矩阵确定fori_point=1:point;%电压矩阵确定voltage_matrix(i_point)=(-j*v0/60)*sin(k*abs(matrix_wavelength(i_poin