untitled2.m

来自「光纤SPR共振峰理论分析」· M 代码 · 共 12 行

t=0.6328
a1=0.6961663;a2=0.4079426;a3=0.8974794;b1=4.679148*10^(-3);
b2=1.31206*10^(-2);b3=93.934;
aa1=0.48755108;aa2=0.39875031;aa3=2.31203530;
bb1=1.8821784*10^(-3);bb2=8.95188847*10^(-3);bb3=5.66135591*10^(2);
 n1=sqrt(1+a1*t.^2./(t.^2-b1)+a2*t.^2./(t.^2-b2)+a3*t.^2./(t.^2-b3))
 n3=sqrt(1+aa1*t.^2./(t.^2-bb1)+aa2*t.^2./(t.^2-bb2)+aa3*t.^2./(t.^2-bb3))
 lambda=t*10^(-6)
 lambda1=1.7614*10^(-5)
 lambda2=1.4541*10^(-7)
 epsilon2=1-(lambda)^2*lambda1/((lambda2)^2*(lambda1+i*lambda))