ycdt.m

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

function L=ycdt(j0,jm,F)% y=0/ym 的边界条件
global I c dt dx dy
ib=2:I-1;
L=-F(ib,jm,1)...
   +(c*dt-dy)/(c*dt+dy)*(F(ib,jm,3)+F(ib,j0,1))...
   +2*dy/(c*dt+dy)*(F(ib,j0,2)+F(ib,jm,2))...
   +(c*dt)^2*dy/(2*dx^2*(c*dt+dy))*(F(ib+1,j0,2)-2*F(ib,j0,2)+F(ib-1,j0,2)+F(ib+1,jm,2)-2*F(ib,jm,2)+F(ib-1,jm,2));