天线阵电流激励幅度相等相位成规律的MATLAB程序

%Ö±ÏßÕó£¬ÕóÔªµÄµçÁ÷·ù¶ÈÏàµÈ£¬Ïàλ³ÉijÖÖ¹æÂÉ·Ö²¼Ê±ÌìÏߵķ½Ïòͼ%ÕóÔª¼ä¾àΪ0.5±¶²¨³¤clcclearN=5;c=3e8;f=10e9;lamda=c/f;k=2*pi/lamda;theta=linspace(0,pi,400);phi=linspace(0,2*pi,400);d=0.5*lamda;z=linspace(0,(N-1)*d,N);%*****************ÏàλÏßÐÔ·Ö²¼************************k1=-240/360*2*pi;alpha1=k1*z./((N-1)*d);I=1;f=0;forn=0:N-1temp=I*exp(1i*k*n*d.*cos(theta)+1i*alpha1(n+1));f=f+temp;endfa=abs(f);x1=(fa.*sin(theta))'*cos(phi);y1=(fa.*sin(theta))'*sin(phi);z1=(fa.*cos(theta))'*ones(size(phi));figure(1)subplot(2,1,1)mesh(x1,y1,z1)title('µçÁ÷ÏàλÏßÐÔ·Ö²¼')subplot(2,1,2)polar(theta,fa)%***************Ïàλ¶þ´Î¹æÂÉ·Ö²¼*********************k2=48/360*2*pi;alpha2=k2*((2*z-(N-1)*d)/((N-1)*d)).^2;I=1;f=0;forn=0:N-1temp=I.*exp(1i*k*n*d.*cos(theta)+1i*alpha2(n+1));f=f+temp;endfa=abs(f);x2=(fa.*sin(theta))'*cos(phi);y2=(fa.*sin(theta))'*sin(phi);z2=(fa.*cos(theta))'*ones(size(phi));figure(2)subplot(2,1,1)mesh(x2,y2,z2)title('µçÁ÷¶þ´Î¹æÂÉ·Ö²¼')subplot(2,1,2)polar(theta,fa)%%***************ÏàλÈý´Î¹æÂÉ·Ö²¼*********************k3=-19.2/360*2*pi;alpha3=k3*((2*z-(N-1)*d)/((N-1)*d)).^3;I=1;f=0;forn=0:N-1temp=I.*exp(1i*k*n*d.*cos(theta)+1i*alpha3(n+1));f=f+temp;endfa=abs(f);x3=(fa.*sin(theta))'*cos(phi);y3=(fa.*sin(theta))'*sin(phi);z3=(fa.*cos(theta))'*ones(size(phi));figure(3)subplot(2,1,1)mesh(x3,y3,z3)title('µçÁ÷Èý´Î¹æÂÉ·Ö²¼')subplot(2,1,2)polar(theta,fa)