setcon.m

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

C0=1/(sqrt(u0/e0)*(g0/2+e0/dt));
D0=sqrt(u0/e0)*(g0/2-e0/dt);

C1=1/(sqrt(u0/e0)*(g1/2+e1/dt));
D1=sqrt(u0/e0)*(g1/2-e1/dt);

C2=1/(sqrt(u0/e0)*(g2/2+e2/dt));
D2=sqrt(u0/e0)*(g2/2-e2/dt);

%////////////////////////////////////
%  C10,C20,D10,D20代表空气的介质参数
%  C11,C21,D11,D21代表玻璃的介质参数
%  C12,C22,D12,D22代表银膜的介质参数
%////////////////////////////////////
%  空气介质参数
C(1:I,1:J)=C0;
D(1:I,1:J)=D0;