FDTD脚本吸收代码

#在脚本编辑区里的吸收
f=getdata("R","f");
T=-transmission("T");
R=transmission("R");
A=1-R-T;
plot(c/f*1e6,R,T,A,"Wavelength","R T");
legend("R","T","A");

#在分析组里的吸收，注意要把监视器T也放到分析组里
make_plots = 0;  
m="E";
f=getdata(m,"f");

T = transmission("T");
R =abs(transmission("R"));
if (make_plots) 
{
  plot(f/1e9,T,R,"f (GHz)","Power");
  legend("T","R");
}

R=getresult("R","T");
R=abs(R.T);

T=getresult("T","T");
T=abs(T.T);

A=1-R-T;
A_matrix=A;
A = matrixdataset("A");
A.addparameter("lambda",c/f,"f",f);
A.addattribute("A",A_matrix);

R_matrix=R;
R = matrixdataset("R");
R.addparameter("lambda",c/f,"f",f);
R.addattribute("R",R_matrix);

T_matrix=T;
T = matrixdataset("T");
T.addparameter("lambda",c/f,"f",f);
T.addattribute("T",T_matrix);