matlab计算传统表面等离子体共振角度谱的反射率曲线.zip

%注释：lamda 波长； clear; c=3.0*10^8; for p=0:1 d=(32)*10^-9; %金膜厚度 for o=1 %折射率步数 for m=1:7000 %角度步数 for n=0 %波长步数 lamda=365*10^-9+p*10^-9;%波长 theta=(10+m*0.01)*pi/180;%角度 e0=1.53627^2; %棱镜介电常数 e1=-2.6+0.43i; % e1=LD(lamda,'Au','LD'); % e1=-12.98+1.41i; %金膜介电常数 % e1=-26.0779 + 1.8907i; % e1=-12.598+1.15117i; % e1=(0.185+5.11i)^2; e2=(1.00026)^2; %被测介质 Xp0=e0/sqrt(e0-e0*(sin(theta))^2); Xp1=e1/sqrt(e1-e0*(sin(theta))^2); Xp2=e2/sqrt(e2-e0*(sin(theta))^2); Xs0=2*pi*sqrt(e0-e0*(sin(theta))^2)/lamda; Xs1=2*pi*sqrt(e1-e0*(sin(theta))^2)/lamda; Xs2=2*pi*sqrt(e2-e0*(sin(theta))^2)/lamda; rp01=(Xp0-Xp1)/(Xp0+Xp1); rp12=(Xp1-Xp2)/(Xp1+Xp2); rs01=(Xs0-Xs1)/(Xs0+Xs1); rs12=(Xs1-Xs2)/(Xs1+Xs2); kz1=2*pi*sqrt(e1-e0*(sin(theta))^2)/lamda; kz2=2*pi*sqrt(e2-e0*(sin(theta))^2)/lamda; rp02=(rp01+rp12*exp(2i*kz1*d))/(1+rp01*rp12*exp(2i*kz1*d)); rs02=(rs01+rs12*exp(2i*kz1*d))/(1+rs01*rs12*exp(2i*kz1*d)); Rp(m,p+1)=abs(rp02); Rp2(m,p+1)=abs(rp02)^2; Rs(m,p+1)=abs(rs02); Php(m,p+1)=angle(rp02); Phs(m,p+1)=angle(rs02); end end end end plot(Rp2);