xcdt.m

来自「二维光子晶体弯曲波导的fdtd程序,matlab语言 ,TE波, PML边界条件」· M 代码 · 共 7 行

function L=xcdt(i0,im,F)% x=0/xm 的边界条件
global J c dt dx dy
jb=2:J-1;
L=-F(im,jb,1)...
   +(c*dt-dx)/(c*dt+dx)*(F(im,jb,3)+F(i0,jb,1))...
   +2*dx/(c*dt+dx)*(F(i0,jb,2)+F(im,jb,2))...
   +(c*dt)^2*dx/(2*dy^2*(c*dt+dx))*(F(i0,jb+1,2)-2*F(i0,jb,2)+F(i0,jb-1,2)+F(im,jb+1,2)-2*F(im,jb,2)+F(im,jb-1,2));