pp=reshape(pp,N,N-2); mat=[pp,q',s'];%构造矩阵 I=mat\\t';
II=[0;I(1:N-2);0];%加上两端零电流 Current=abs(II);
x=linspace(-L/2,L/2,N); figure(1);
string=['b','g','r','y','c','k','m','r'];
string1=['ko','bo','yo','co','mo','ro','go','bo']; plot(x,Current,string(ii),'linewidth',1.3); xlabel('L/\\lambda'),ylabel('电流分布'); grid on hold on
%legend('L=0.1\\lambda','L=0.2\\lambda','L=0.3\\lambda','L=0.4\\lambda','L=0.5\\lambda','L=0.6\\lambda','L=0.7\\lambda','L=0.8\\lambda','L=0.9\\lambda','L=1\\lambda') legend('L=0.1\\lambda','L=0.3\\lambda','L=0.5\\lambda','L=0.7\\lambda','L=0.9\\lambda','L=1.1\\lambda','L=1.3\\lambda','L=1.5\\lambda') Zmn=1/I((N+1)/2);%%%%%%V=1v theta=linspace(0,2*pi,360);
for m=1:360 for n=1:N
F1(m,n)=II(n).*exp(j*k*z(n)*cos(m*pi/180))*len*sin(m*pi/180); end end
F2=-sum(F1');
F=F2/max(F2);%%%归一化 figure(2);
polar(theta,abs(F),string(ii)); title('E面归一化方向图') view(90,-90)
%legend('L=h\\lambda','L=0.3\\lambda','L=0.3\\lambda','L=0.4\\lambda','L=0.5\\lambda','L=0.6\\lambda','L=0.7\\lambda','L=0.8\\lambda','L=0.9\\lambda','L=1\\lambda') legend('L=0.1\\lambda','L=0.3\\lambda','L=0.5\\lambda','L=0.7\\lambda','L=0.9\\lambda','
L=1.1\\lambda','L=1.3\\lambda','L=1.5\\lambda') hold on figure(3) kk=1;
for phi=0:pi/180:2*pi for n=1:N
FF(n)=II(n)*len*exp(i*k*len*n*cos(pi/2))*sin(pi/2); end;
FFF(kk)=sum(FF); kk=kk+1; end;
phi=0:pi/180:2*pi;
polar(phi,FFF/max(abs(FFF)),string(ii));title('pattern,F({\\theta},{\\phi}),\\theta=90');
legend('L=0.1\\lambda','L=0.3\\lambda','L=0.5\\lambda','L=0.7\\lambda','L=0.9\\lambda','L=1.1\\lambda','L=1.3\\lambda','L=1.5\\lambda') hold on figure(4)
polar(phi,FFF/max((FFF)),string(ii));title('pattern,F({\\theta},{\\phi}),\\theta=90'); hold on figure(5) mm=1;
for theta=0:0.01*pi:pi; for n=1:N
E(1,n)=2*pi*c*u0*len/(4*pi*1)*(exp(-i*k*1)*exp(i*k*len*n*cos(theta))*sin(theta)); end
EE=E*II;
G(mm)=(4*pi*1^2)/ata/abs(II((N-1)/2+1))^2/(-real(Zmn))*abs(EE)^2; mm=mm+1; end
归
一
化
H-plane
不
同
L/\\lambda
H-plane