带特定形状挡板的平行板传输线.docVIP

带特定形状挡板的平行板传输线 对于开槽形状如图8的平行传输线，在开槽处添加上边界条件Ex=0就可以对其进行仿真。参数选取与上面均匀传输线仿真相同，并在同一观测面下记录总电压，用公式，求解反射电压和透射电压。 程序和图如下： d=0.18;c=3e8;L=6;mu=4*pi*1e-7;ep=1/(4*pi*c^2)*1e7;T=0.5e-9;t0=3*T;delt_z=0.015; delt=delt_z/(2*c); A=(c*delt-delt_z)/(c*delt+delt_z); Nz=L/delt_z; Nx=d/delt_z; B=zeros(1,Nx); C=zeros(1,Nx); F=zeros(1,120); Ex=zeros(Nx,Nz+1); Ez=zeros(Nx+1,Nz); Hy=zeros(Nx,Nz); Vref1total=zeros(1,2000); Vref2total=zeros(1,2000); Vref1re=zeros(1,2000); Vref2tr=zeros(1,2000);%参数的设置 for n=1:120 F(n)=exp(-(n*delt-t0)^2/T^2); end for n=1:2000 if n=120 for i=1:Nx Ex(i,1)=F(n);%加激励源 end end for i=1:Nx for k=1:Nz Hy(i,k)=Hy(i,k)+(delt/(mu*delt_z))*(Ez(i+1,k)-Ez(i,k)+Ex(i,k)-Ex(i,k+1)); %蛙跳公式 end end for i=1:Nx B=Ex(i,2); C=Ex(i,Nz);%n时刻变量的保存 end for i=1:Nx for k=2:Nz Ex(i,k)=Ex(i,k)+(delt/(ep*delt_z))*(Hy(i,k-1)-Hy(i,k));%蛙跳公式 end end if n120 for i=1:Nx Ex(i,1)=B+A*(Ex(i,2)-Ex(i,1));%入口边界条件的切换 end end for i=2:Nx for k=1:Nz Ez(i,k)=Ez(i,k)+(delt/(ep*delt_z))*(Hy(i,k)-Hy(i-1,k));% end end for i=1:Nx/3%挡板条件的设置 Ex(i,Nz/2+1)=0; Ex(i,Nz/2+5)=0; end for i=Nx/3:2*Nx/3 Ex(i,Nz-2)=0; end for k=Nz/2+1:Nz/2+5 Ez(Nx/3+1,k)=0; end for k=Nz-2:Nz+1 Ez(Nx/3,k)=0; Ez(2*Nx/3,k)=0; end for k=1:Nz Ez(Nx+1,k)=0; Ez(1,k)=0; end for i=1:Nx Ex(i,Nz+1)=0;%挡板条件的设置 end if k==101 for i=1:Nx Vref1total(n)=Vref1total(n)+Ex(i,101)*delt_z;%参考电压 end end if k==301 for i=1:Nx Vref2total(n)=Vref2total(n)+Ex(i,301)*delt_z;%参考电压 end end if mod(n,50)==0 Vref1re(n)= Vref1total(n)-Vref1in(n);%反射电压 Vref2tr(n)=Vref2total(n);%透射电压 %plot(Vref1re); %plot(Vref2tr); subplot(3,1,1) mesh(Ex) subplot(3,1,2) mesh(Ez) subplot(3,1,3) mesh(Hy) pause(0.3) end end