電磁數值計算作業(yè)男11 2

ifn==N

z(m,n-1)=cos(k*zm);elseifn==N+1

z(m,n-1)=sin(k*zm);

end

end

else

z(m,n-1)=p(m,n);end

z1=inv(z);t=zeros(N,1);

form=1:N

zm=-0.5*L+m*l;t(m)=1/(j*2*377)*sin(k*abs(zm));

endzz=z1*t;form=1:N

ifm==1

I(m)=0;

elseifm==N

I(m)=0;

end

elseI(m)=zz(m-1);end

zin(pp)=1/I((N+1)/2);end

x=zeros(1,N);y1=zeros(1,N);form=1:N

x(1,m)=w0/w(m);

end

forq=1:N

y1(1,q)=imag(zin(q));

endplot(x,y1,'k');holdon;

endend

、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、、

E面圖源程序：

functionhallenZin2

N=201;L=0.5;l=L/N;t=8.85*10^(-12);u=4*3.14*10^(-7);r=0.00001*L;

w=2/3*L;

f=3*10^8/w;k=2*3.14/w;form=1:N

forn=1:N

ifm==n

p(m,n)=1/(2*3.14)*log(l/r)-j*k*l/4/3.14;

end

end

else

end

R=((m-n)^2*l^2+r^2)^0.5;p(m,n)=l*exp(-j*k*R)/(4*3.14*R);

form=1:N

zm=-0.5*L+m*l;forn=2:N+1

zn=-0.5*L+n*l;ifn==N

z(m,n-1)=cos(k*zm);elseifn==N+1

z(m,n-1)=sin(k*zm);

end

end

else

z(m,n-1)=p(m,n);end

z1=inv(z);t=zeros(N,1);

I=zeros(N,1);

form=1:N

zm=-0.5*L+m*l;t(m)=1/(j*2*377)*sin(k*abs(zm));

endzz=z1*t;form=1:N

ifm==1

I(m)=0;

elseifm==N

I(m)=0;

elseI(m)=zz(m-1);

end

endEE=zeros(N,1);t=zeros(N,1);

E=zeros(N,1);t=linspace(0,2*pi,N);forq=1:N

forp=1:N

E(p)=l*sin(t(q))*exp(j*k*p*l*cos(t(q)));

end

endEE(q)=I'*E;

forp=1:N

F(p)=abs(EE(p)/max(EE));

endF

polar(t,F,'b--');view(90,-90);holdon;

end

end

endEE(q)=I'*E;

forp=1:N

F(p)=