FDTD脚本吸收代码

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);