TFSF边界条件

待续

 %1D FDTD simulation with a simple absorbing boundary condition

 % and a TFSF boundary between hy[] and ez[].

 SIZE = ;

 ez = zeros(,SIZE);

 hy = zeros(,SIZE);

 imp0=377.0;

 maxTime = ;

 for qTime =  :  : maxTime

      hy(SIZE)=hy(SIZE-);%

     for mm =  :  : SIZE -

         hy(mm) = hy(mm) + (ez(mm+)- ez(mm))/imp0;

     end

 % correction for Hy adjacent to TFSF boundary */

  hy()=hy()-exp(-(qTime - )*(qTime - ) / ) / imp0;

     figure(); t_hy = plot(hy);

     title('Magnetic Field');

     ez()=ez();%

     for mm =  :  : SIZE

          ez(mm) = ez(mm) + (hy(mm)- hy(mm-))*imp0;

     end

     figure(); t_ez = plot(ez);

     title('Electric Field')

     ez() = ez() + exp(-(qTime+0.5-(-0.5)-30.0)*(qTime+0.5-(-0.5)-30.0)/100.0);

 %     pause(0.01);

 % %     refreshdata(t_hy);

 %     refreshdata(t_ez);

 %     drawnow;

 end