帶特定形狀擋板的平行板傳輸線.doc

帶特定形狀擋板的平行板傳輸線對(duì)于開槽形狀如圖8的平行傳輸線，在開槽處添加上邊界條件Ex=0就可以對(duì)其進(jìn)行仿真。參數(shù)選取與上面均勻傳輸線仿真相同，并在同一觀測(cè)面下記錄總電壓，用公式，求解反射電壓和透射電壓。程序和圖如下：d=0.18;c=3e8;L=6;mu=4*pi*1e-7;ep=1/(4*pi*c2)*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);%參數(shù)的設(shè)置for n=1:120 F(n)=exp(-(n*delt-t0)2/T2);endfor n=1:2000 if n120 for i=1:Nx Ex(i,1)=B+A*(Ex(i,2)-Ex(i,1);%入口邊界條件的切換 endend 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);% endend for i=1:Nx/3%擋板條件的設(shè)置 Ex(i,Nz/2+1)=0; Ex(i,Nz/2+5)=0;endfor 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; endfor k=1:Nz Ez(Nx+1,k)=0; Ez(1,k)=0; endfor i=1:Nx Ex(i,Nz+1)=0;%擋板條件的設(shè)置endif k=101for i=1:Nx Vref1total(n)=Vref1total(n)+Ex(i,101)*delt_z;%參考電壓endend if k=301for i=1:Nx Vref2total(n)=Vref2total(n)+Ex(i,301)*delt_z;%參考電壓 endendif mod(n,50)=0 Vref1re(n)= Vref1total(n)-Vref1in(n);%反射電壓 Vref2tr(n)=Vref2total(n);%透射電壓 %plot(Vref1re); %plot(Vref2tr); subplot(3,1,1)mesh(E