矩量法matlab程序設(shè)計(jì)實(shí)例

pp=reshape(pp,N,N-2);mat=[pp,q',s'];%構(gòu)造矩陣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('電流分布');gridonholdon%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=1vtheta=linspace(0,2*pi,360);form=1:360forn=1:NF1(m,n)=II(n).*exp(j*k*z(n)*cos(m*pi/180))*len*sin(m*pi/180);endendF2=-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')holdonfigure(3)kk=1;forphi=0:pi/180:2*piforn=1:NFF(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('不同L/\lambdaH-planepattern,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')holdonfigure(4)polar(phi,FFF/max((FFF)),string(ii));title('歸一化H-planepattern,F({\theta},{\phi}),\theta=90');holdonfigure(5)mm=1;fortheta=0:0.01*pi:pi;forn=1:NE(1,n)=2*pi*c*u0*len/(4*pi*1)*(exp(-i*k*1)*exp(i*k*len*n*cos(theta))*sin(theta));en