matlab-对菲涅尔公式画图.doc

2010～2011学年第一学期 MATLAB课程 用MATLAB 语言，绘制SPR理论中 rp、rs、R对于入射角和折率比值n=n1/n2的二维和三维图形,其中n的范围为（0.6，1.5），并结合每一个图形分析其特点及物理意义。 (其中rp为反射光中平行分量的反射系数；rs为反射光中垂直分量的反射系数； R-反射率；n1、n2分别表示两种不同介质的折射率) 答案： 程序： （1）%反射系数r，和反射率R与折射率之比n的关系； clear; clc; clf; n=0.6:0.05:1.5; zeta1=pi/10; %入射角 zeta2=real(asin(n.*sin(zeta1))); %折射角 rpz=-n.*cos(zeta2)+cos(zeta1); %平行分量反射部分分子 rpm=n.*cos(zeta2)+cos(zeta1); %平行分量反射部分分母 rp=rpz./rpm; %平行分量反射系数 rsz=n.*cos(zeta1)-cos(zeta2); %垂直分量反射部分分子 rsm=n.*cos(zeta1)+cos(zeta2); %垂直分量反射部分分母 rs=rsz./rsm; %垂直分量反射系数 Rp=rp.^2; %平行分量反射率 Rs=rs.^2; %垂直分量反射率 t=0.6:0.005:1.5; y=0; y1=-0.4:0.005:0; t1=1; plot(n,rp,c,n,rs,*:,n,Rp,k-,n,Rs,g,t,y,-,t1,y1,-) text(1.3,-0.1,rp); text(1.2,0.09,rs); text(0.65,0.023,Rp); text(0.65,0.065,Rs); text(0.9,-0.2,(n1)); text(1.2,-0.2,(n1)); title(relationships between rp,rs,Rp,Rs and n=n1/n2); xlabel(n=n1/n2); ylabel(r or R); （1）分析其物理意义：（取入射角pi/10，小于Brewster角，小于临界角） 当n1（光疏到光密介质）时，反射光中垂直分量（s）的反射系数rs0，说明反射光中s分量与入射光中的s分量相位相反；而p分量，即平行分量，其反射系数rp0，说明反射光中p分量与入射光中的p分量相位相同。 当n=1时，相当于光束在同一种介质中传输，因此无反射，反射系数为零； 当n1（光密到光疏介质）时，反射光中rs0，说明反射光中s分量与入射光中的s分量相位相同；而p分量rp0，说明p分量在两种介质中相位相反。 %反射系数r，和反射率R与入射角x的关系（n1）；clear; clc; clf; n=0.6:0.05:1.5; zeta1=pi/10; %入射角 zeta2=real(asin(n.*sin(zeta1))); %折射角rpz=-n.*cos(zeta2)+cos(zeta1); %平行分量反射部分分子 rpm=n.*cos(zeta2)+cos(zeta1); %平行分量反射部分分母 rp=rpz./rpm; %平行分量反射系数 rsz=n.*cos(zeta1)-cos(zeta2); %垂直分量反射部分分子 rsm=n.*cos(zeta1)+cos(zeta2); %垂直分量反射部分分母 rs=rsz./rsm; %垂直分量反射系数 Rp=rp.^2; %平行分量反射率 Rs=rs.^2; %垂直分量反射率 Brewster=acot(n)