電磁數(shù)值計(jì)算作業(yè)

functionZin=%Linput%ainput%lambdainputU0(4*pi)*10^- E0=(1/(36*pi))*(10^-9); K=2*pi/lambda;w= xyzlable1= xyzlable2= xyzlable3= psi1=ones(N,N);%計(jì)算psi1矩陣psi2ones(N,N);%計(jì)算psi2psi3ones(N,N);%計(jì)算psi3psi4ones(N,N);%計(jì)算psi4psi5ones(N,N);%計(jì)算psi5fork=1:Nform=ifm==psi1(k,m)=1/(2*pi*DeltaL)*log(DeltaL/a)-j*K/(4*pi);psi2(k,m)=1/(2*pi*DeltaL)*log(DeltaL/a)-j*K/(4*pi);psi3(k,m)=1/(2*pi*DeltaL)*log(DeltaL/a)-psi1(k,m)=exp(-psi2(k,m)=exp(-psi3(k,m)=exp(-

fork=form=if(k-m)==-psi4(k,m)=1/(2*pi*DeltaL)*log(DeltaL/a)-psi4(k,m)=exp(-

fork=form=if(k-psi5(k,m)=1/(2*pi*DeltaL)*log(DeltaL/a)-psi5(k,m)=exp(-

Z fork=form=Zin=1/Ii((N+1)/2);%對(duì)稱陣子輸入阻抗%矩量法計(jì)算對(duì)稱陣子參數(shù)雙位法EfunctionZin=%Linput%ainput%lambdainputU0(4*pi)*10^- E0=(1/(36*pi))*(10^-9); K=2*pi/lambda;w=xyzlable1= xyzlable2= xyzlable3= psi1=ones(N,N);%計(jì)算psi1矩陣psi2ones(N,N);%計(jì)算psi2psi3ones(N,N);%計(jì)算psi3psi4ones(N,N);%計(jì)算psi4psi5ones(N,N);%計(jì)算psi5fork=1:Nform=ifm==psi1(k,m)=1/(2*pi*DeltaL)*log(DeltaL/a)-j*K/(4*pi);psi2(k,m)=1/(2*pi*DeltaL)*log(DeltaL/a)-j*K/(4*pi);psi3(k,m)=1/(2*pi*DeltaL)*log(DeltaL/a)-psi1(k,m)=exp(-psi2(k,m)=exp(-psi3(k,m)=exp(-

fork=form=if(k-m)==-psi4(k,m)=1/(2*pi*DeltaL)*log(DeltaL/a)-psi4(k,m)=exp(-

fork=form=if(k-psi5(k,m)=1/(2*pi*DeltaL)*log(DeltaL/a)-psi5(k,m)=exp(-

Z fork=form=Zin=1/Ii((N+1)/2);%對(duì)稱陣子輸入阻抗 theta=-form=1:length(theta)fory=theta1(x(2))-theta1(x(1))F2([1length(theta)])=0; %E面歸一化方向圖（直角坐標(biāo)）holdonplot([theta1(x(1))theta1(x(1))],[0F1(x(1))],'b-plot([theta1(x(2))theta1(x(2))],[0F1(x(2))],'b-.'))grid holdon functionZin=%Linput%ainput%lambdainputU0(4*pi)*10^- E0=(1/(36*pi))*(10^-9); K=2*pi/lambda;w=xyzlable1= xyzlable2= xyzlable3= xyzlable4=zeros(N+2,1);psi1=ones(N,N);%計(jì)算psi1矩陣psi2ones(N,N);%計(jì)算psi2psi3ones(N,N);%計(jì)算psi3psi4ones(N,N);%計(jì)算psi4psi5ones(N,N);%計(jì)算psi5fork=1:Nform=ifm==psi1(k,m)=1/(2*pi*DeltaL)*log(DeltaL/a)-j*K/(4*pi);psi2(k,m)=1/(2*pi*DeltaL)*log(DeltaL/a)-j*K/(4*pi);psi3(k,m)=1/(2*pi*DeltaL)*log(DeltaL/a)-psi1(k,m)=exp(-psi2(k,m)=exp(-psi3(k,m)=exp(-

fork=form=if(k-m)==-psi4(k,m)=1/(2*pi*DeltaL)*log(DeltaL/a)-psi4(k,m)=exp(-

fork=form=if(k-psi5(k,m)=1/(2*pi*DeltaL)*log(DeltaL/a)-

psi5(k,m)=exp(-Z fork=form=%電流分布曲線xyzlable4([1N+2])=[-L/2,L/2];I([1N+2])=0;x=(xyzlable4+L/2)*pi/L;y=sin(x);axis([-gridonf0=1.5e8a=0.000001*[110100100010000forn=1:length(a)Zin(n,m)=Xin(n,m)=forG(1,m)=legend('a/L=0.000001','a/L=0.00001','a/L=0.0001','a/L=0.001','a/L=0.01','a/L=0.1',0)legend('a/L=0.000001','a/L=0.00001','a/L=0.0001','a/L=0.001','a/L=0.01','a/L=0.1',0)text(1,-2,'\fontsize{16}L=1m{f}_{0}=150MHz')gridonforG(1,m)=legend('a/L=0.000001','a/L=0.00001','a/L=0.0001','a/L=0.001','a/L=0.01','a/L=0.1',0)text(0.25,4000,'\fontsize{16}L=1mf}_{0}=150MHz')legend('a/L=0.000001'