run_int_thea.asv

光纤SPR共振峰理论分析

clear
clc
a0=200*10^(-6);
L=10*10^(-3);
NA=0.36;n3=1.38;n4=1.3330;d1=40*10^(-9);d2=0;
a1=0.6961663;a2=0.4079426;a3=0.8774794;b1=0.068403;
b2=0.1162414;b3=9.896161;
t1=cputime;
theta(=linspace((pi/2-asin(NA/1.457))*180/pi,90,50);
theta1(1,:)=theta(1,:)*pi/180;
T1=linspace(400,1000,150);
D_m=linspace(40,40,1);
F2=zeros(length(T1),length(theta1));
R=zeros(2,length(theta1));
r=zeros(1,length(theta1));
for ss2=1:length(D_m)
    d1=D_m(ss2)*10^(-9);
    D=[d1,d2];
for ss1=1:length(T1)
    t=T1(1,ss1)*10^(-3);
    lambda=t*10^(-6);
    n1=sqrt(1+a1*t.^2./(t.^2-b1^2)+a2*t.^2./(t.^2-b2^2)+a3*t.^2./(t.^2-b3^2));
   epsilon2=2.41955+8.66878*t-64.44247*t.^2+7.61924*t.^3+i*(-0.99069+8.09121*t-13.76576*t.^2+9.55063*t.^3);
   n2=sqrt(epsilon2);
   N_index=[n1,n2,n3,n4];
   N_index2=[n1,n2,n3,1.0];
  f1=mat_ME_theta(lambda,N_index,D,theta1,NA,L,a0,R,r);
  f2=mat_ME_theta(lambda,N_index2,D,theta1,NA,L,a0,R,r);
   fm(ss1,:)=f1;
   fair(ss1,:)=f2;
  % t2=cputime-t1
end
F=fm./fair;
end

figure,surf(T1,theta,F')
figure,imagesc(T1,theta,F')