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

1、带特定形状挡板的平行板传输线对于开槽形状如图8的平行传输线，在开槽处添加上边界条件Ex=O就可以对其进行仿真。参数选取与上面均匀传输线仿真相同，并在同一观测面下记录总 电压，用公式v；ei二v；o;ai(t)亠驚们，求解反射电压和透射电 压。程序和图如下：d=0.18;c=3e8;L=6;mu=4*pi*1e-7;ep=1/(4*pi*cA2)*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=z

2、eros(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:120F(n )=exp(-( n*delt-t0f2/TA2);endfor n=1:2000 if n<=120 for i=1:NxEx(i,1)=F( n);%加激励源endendfor i=1:Nx for k=

3、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);式endendfor i=1:NxB=Ex(i,2);C=Ex(i,Nz);%用寸刻变量的保存endfor i=1:Nx%蛙跳公式for k=2:NzEx(i,k)=Ex(i,k)+(delt/(ep*delt_z)*(Hy(i,k-1)-Hy(i,k);endendif n>120for i=1:NxEx(i,1)=B+A*(Ex(i,2)-Ex(i,1);%入口边界条件的切换endendfor i=2:Nxfor k=1:NzE

4、z(i,k)=Ez(i,k)+(delt/(ep*delt_z)*(Hy(i,k)-Hy(i-1,k);endendfor i=1:Nx/3%挡板条件的设置Ex(i,Nz/2+1)=0;Ex(i,Nz/2+5)=0;endfor i=Nx/3:2*Nx/3Ex(i,Nz-2)=0;endfor k=Nz/2+1:Nz/2+5Ez(Nx/3+1,k)=0;endfor k=Nz-2:Nz+1Ez(Nx/3,k)=0;Ez(2*Nx/3,k)=0;endfor k=1:NzEz(Nx+1,k)=0;Ez(1,k)=0;endfor i=1:NxEx(i,Nz+1)=0;%挡板条件的设置endif

5、k=101for i=1:Nx%参考电压%参考电压Vref1total(n)=Vref1total(n)+Ex(i,101)*delt_z;endendif k=301 for i=1:NxVref2total(n)=Vref2total(n)+Ex(i,301)*delt_z;endendif mod(n,50)=0Vref1re(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 end45040035Q20015050106820.235030003&#