untitled2.m

光纤SPR共振峰理论分析

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