matlab对菲涅尔公式画图.docVIP

  • 30
  • 0
  • 约6.84千字
  • 约 8页
  • 2017-02-02 发布于江苏
  • 举报
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)

文档评论(0)

1亿VIP精品文档

相关文档