表面等离子体共振的Matlab模拟，可以计算出共振峰的准确位置.rar

【项目代码】表面等离子体共振的Matlab模拟，可以计算出共振峰的准确位置.rar

%多层结构最终反射系数曲线计算

n0=1.65;%棱镜折射率
e0=n0^2;%棱镜介电常数

n1=1.56;%波片 
e1=n1^2;%波片介电常数

e2=-13.4+1.4i;%金属膜介电常数

n3=1.35;%1.40;1.45;1.50;%被测物折射率
e3=n3^2;

n4=1.33;%缓冲液体
e4=n4^2;

lamta=632.8;%波长 

theta_degree=linspace（4070100）;%入射光度数离散化
theta_deg=pi/180*theta_degree;  %化为弧度

for number=1:100;
theta=theta_deg（number）;

kz0=2*pi/lamta*sqrt（e0-e0*sin（theta）^2）;
kz1=2*pi/lamta*sqrt（e1-e0*sin（theta）^2）;
kz2=2*pi/lamta*sqrt（e2-e0*sin（theta）^2）;
kz3=2*pi/lamta*sqrt（e3-e0*sin（theta）^2）;
kz4=2*pi/lamta*sqrt（e4-e0*sin（theta）^2）;

zeta0=e0/kz0;
zeta1=e1/kz1;
zeta2=e2/kz2;
zeta3=e3/kz3;
zeta4=e4/kz4;


r34=（zeta4-zeta3）/（zeta4+zeta3）;
r23=（zeta3-zeta2）/（zeta3+zeta2）;
r12=（zeta2-zeta1）/（zeta2+zeta1）;
r01=（zeta1-zeta0）/（zeta1+zeta0）;


d3=0;  %被测物厚度
d2=60; %金属膜厚度
d1=40; %波片厚度


r24=（r23+r34*exp（2*i*d3*kz3））/（1+r23*r34*exp（2*i*d3*kz3））;
r14=（r12+r24*exp（2*i*d2*kz2））/（1+r12*r24*exp（2*i*d2*kz2））;
r04=（r01+r14*exp（2*i*d1*kz1））/（1+r01*r14*exp（2*i*d1*kz1））;

r=abs（r04）^2;

R（number）=r;%最终反射系数

end

plot（theta_degreeR‘r‘）

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1153  2006-05-24 08:59  RSPR.m

     文件     256894  2006-09-21 17:15  基于SPR光谱分析的液体折射率测量研究.pdf

----------- ---------  ---------- -----  ----

               258265                    3