run_intens.m

光纤SPR共振峰理论分析

clear
clc
%load('F_air.mat')
n4=1.3480;d1=40*10^(-9);d2=0;
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);
%lambda1=1.7614*10^(-5); lambda2=1.4541*10^(-7);
lambda1=0.89342*10^(-5); lambda2=1.6826*10^(-7);%lambda1_c,lambda2_p
t1=cputime;
theta=56;
T1=linspace(400,1000,601);
D_m=linspace(40,40,1);
F2=zeros(length(D_m),length(T1));
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));
    n3=sqrt(1+aa1*t.^2./(t.^2-bb1)+aa2*t.^2./(t.^2-bb2)+aa3*t.^2./(t.^2-bb3));
   %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);
   epsilon2=1-(lambda)^2*lambda1/((lambda2)^2*(lambda1+i*lambda));
   n2=sqrt(epsilon2);
   N_index=[n1,n2,n3,n4];
   N_index2=[n1,n2,n3,1.0];
  f1=mat_ME(lambda,N_index,D,theta1,NA,L,a0,R,r);
  f2=mat_ME(lambda,N_index2,D,theta1,NA,L,a0,R,r);
   fm(1,ss1)=f1;
   fair(1,ss1)=f2;
  % t2=cputime-t1
end
F(ss2,:)=fm./fair;
%F(ss2,:)=fm;
F_min(1,ss2)=min(F(ss2,:));
[x2,y2]=min(F_min(1,ss2));
T1(1,y2);
half_x=(max(F(ss2,:))+min(F(ss2,:)))/2;
F2=F(ss2,:)-half_x;
F3=F2(1:end-1).*F2(2:end);
[x1,y1]=find(F3<=0);
F_half(1,ss2)=T1(1,y1(2))-T1(1,y1(1));
 t2=cputime-t1
end
%figure,plot(D_m,F_min)
%figure,plot(D_m,F_half)
%figure,surf(D_m,T1,F')
%figure,imagesc(D_m,T1,F')
%figure,plot(T1,fm./fair)