run_intens.asv

光纤SPR共振峰理论分析

clear
clc
%load('F_air.mat')
a0=200*10^(-6);
L=10*10^(-3);
NA=0.36;n3=1.38;n4=1.3731;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;
theta1=linspace(pi/2-asin(NA/1.457),89*pi/180,50);
T1=linspace(0.4,1,150);
D_m=linspace(40,60,20);
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);
    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(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_min(1,ss2)=min(F(ss2,:));
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)