getcon.m

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

%----------------------------------------------
%   输入参数in为波长（单位：m）
%   输出out为对应的电导率
%----------------------------------------------
function out=getcon(in)
nn=in*1e6;
coon=(-6.21198035452752*nn.^5+ 36.9470450716619*nn.^4-83.1702487377279*nn.^3+...
   94.5079035255191*nn.^2-48.61825472792*nn+9.64235361750187);
eee=(-108.909148353719*nn.^5+ 671.044293885632*nn.^4-1547.0900731642*nn.^3+ ...
   1674.14670244375*nn.^2-773.494216929895*nn+130.284211701206);

e0=(1e-9)/(36*pi);
w=2*pi*3e8./(in);     % 光频率

out=(coon+i*eee)*e0.*w;