coner.m

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

%************************
%
%线性插值法求顶点处的场强
%************************

function Lx=coner(flag,F) 
global I J c dx dy dt
switch flag
    %中心点为(I+1)/2,(J+1)/2
case 1
%(1,1)
    L=(1-2*c*dt/sqrt(dx^2+dy^2))*F(1,1,1)+2*c*dt/sqrt(dx^2+dy^2)*F(2,2,1);
case 2
%(I,1)
    L=(1-2*c*dt/sqrt(dx^2+dy^2))*F(I,1,1)+2*c*dt/sqrt(dx^2+dy^2)*F(I-1,2,1); 
case 3
%(I,J)
    L=(1-2*c*dt/sqrt(dx^2+dy^2))*F(I,J,1)+2*c*dt/sqrt(dx^2+dy^2)*F(I-1,J-1,1); 
case 4
%(1,J)
    L=(1-2*c*dt/sqrt(dx^2+dy^2))*F(1,J,1)+2*c*dt/sqrt(dx^2+dy^2)*F(2,J-1,1); 
end
d0=sqrt((I-1)^2+(J-1)^2)/2-1;
a0=sqrt(d0/(d0+1));%  衰减常数
Lx=a0*L;