- 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)
您可能关注的文档
最近下载
- 河北广谋新材料科技有限公司年产600万平方米复合革及隔音布、400万平方米新型复合网项目环境影响报告表.docx VIP
- 部编版语文五年级下册第七单元教材解读大单元集体备课.pptx VIP
- 部编版语文五年级下册第六单元教材解读大单元集体备课.pptx VIP
- 部编版语文五年级下册第五单元教材解读大单元集体备课.pptx VIP
- 部编版语文五年级下册第四单元教材解读大单元集体备课.pptx VIP
- 部编版语文五年级下册第三单元教材解读大单元集体备课.pptx VIP
- 部编版语文五年级下册第二单元教材解读大单元集体备课.pptx VIP
- 部编版语文五年级下册第一单元教材解读大单元集体备课.pptx VIP
- 车工 (数控车床)理论知识考核要素细目表四级.pdf VIP
- 本科课程《基础护理学》教案,第十二章给药.doc VIP
原创力文档

文档评论(0)