《西安邮电大学光学实验matlab仿真结果分析与程序.doc》由会员分享，可在线阅读，更多相关《西安邮电大学光学实验matlab仿真结果分析与程序.doc（44页珍藏版）》请在金锄头文库上搜索。
1、西安邮电大学光学实验matlab仿真结果分析与程序 光学实验 实验报告 课程名称： 姓 名： 学 院： 系 部： 专 业： 年 级： 学 号： 指导教师： 光学实验 伍金霄 电子工程学院 光电子技术系 电子科学与技术 科技1201 05122012 刘娟 2014年 12 月 24 日 光波在介质中界面上的反射及透射特性 一实验目的： 1掌握反射系数及透射系数的概念； 2掌握反射光与透射光振幅和相位的变化规律； 3掌握布儒斯特角和全反射临界角的概念。 二实验原理： 1 反射定律和折射定律 光由一种介质入射到另一种介质时，在界面上将产生反射和折射。现假设二介质为均匀、透明、各向同性介质，分界面为
2、无穷大的平面，入射、反射和折射光均为平面光波，其电场表示式为 El?E0le?i(?lt?kl?r) l=i, r, t 式中，脚标i、r、t分别代表入射光、反射光和折射光；r是界面上任意点的矢径，在图2-1所示的坐标情况下，有 r=ix+jy 图2-1 平面光波在界面上的反射和折射 图2-2 ki、kr、kt三波矢关系 根据电磁场的边界条件，可以得到如下关系 ?i?r?t (ki?kr)?r?0 (ki?kt)?r?0 这些关系表明： 入射光、反射光和折射光具有相同的频率； 入射光、反射光和折射光均在入射面内，ki、kr和kt波矢关系如图2-2所示。 进一步可得 kisin?i?krsin?
3、r kisin?i?ktsin?t 或 nisin?i?nrsin?rnisin?i?ntsin?t 即介质界面上的反射定律和折射定律，它们给出了反射光、折射光的方向。折射定律又称为斯涅耳(Snell)定律。 2 菲涅耳公式 s分量和p分量 通常把垂直于入射面振动的分量称做s分量，把平行于入射面振动的分量称做p分量。为讨论方便起见，规定s分量和p分量的正方向如图2-3所示。 图2-3 s分量和p分量的正方向 反射系数和透射系数 假设介质中的电场矢量为 El?E0le?i(?t?kl?r) l=i, r, t 其s分量和p分量表示式为 Elm?E0lme?i(?t?kl?r) m=s,p 则定义
4、s分量、p分量的反射系数、透射系数分别为 rm? tm?E0rmE0imE0tm E0tm 菲涅耳公式 假设界面上的入射光、反射光和折射光同相位，根据电磁场的边界条件及s分量、p分量的正方向规定，可得 Eis?Ers?Ets 和 Hipcos?1?Hrpcos?1?Htpcos?2 利用?H?E，上式变为 (Eis?Ers)n1cos?1?Etsn2cos?2 再利用折射定律，消去Ets，经整理可得 Erssin(?2?1)? Eissin(?2?1) 根据反射系数定义，得到 rs?sin(?1?2) sin(?1?2) ts?2n1cos?1 n1cos?1?n2cos?2 将所得到的表示式
5、写成一个方程组，就是著名的菲涅耳公式： rs?E0rssin(?1?2)n1cos?1?n2cos?2tan?1?tan?2? E0issin(?1?2)n1cos?1?n2cos?2tan?1?tan?2 E0rp E0ip?tan(?1?2)n2cos?1?n1cos?2sin2?1?sin2?2 ?tan(?1?2)n2cos?1?n1cos?2sin2?1?sin2?2rp? ts? tp?E0ts2cos?1sin?22n1cos?1?E0issin(?1?2)n1cos?1?n2cos?2E0tp E0ip?2cos?1sin?22n1cos?1?sin(?1?1)cos(?1?2
6、)n2cos?1?n1cos?2 这些系数首先是由菲涅耳用弹性波理论得到的，所以又叫做菲涅耳系数。 于是，如果已知界面两侧的折射率n1、n2和入射角1，就可由折射定律确定折射角2，进而可由上面的菲涅耳公式求出反射系数和透射系数。图2-4绘出了在n1n2(光由光疏介质射向光密介质)和n1n2(光由光密介质射向光疏介质)两种情况下，反射系数、透射系数随入射角1的变化曲线。 图2-4 rs、rp、ts、tp随入射角1变化曲线 反射系数与透射系数不仅反映了反射光和透射光相对于入射光的振幅改变，它还反映了反射光和透射光相对于入射光的相移。图2-6给出了反射光随入射角产生的相位改变。 图2-5 ?rs、?
7、rp随入射角1变化。(a)（b）为光疏到光密的情况；(c)（d）为光密到光疏的情 况 3 反射率和透射率 菲涅耳公式给出了入射光、反射光和折射光之间的场振幅和相位关系。不计吸收、散射等能量损耗，入射光能量在反射光和折射光中重新分配，而总能量保持不变。 图2-6 光束截面积在反射和折射时的变化(在分界面上光束截面积为1) 如图2-6所示，若有一个平面光波以入射角1斜入射到介质分界面，平面光波的强度为Ii，则每秒入射到界面上单位面积的能量为 Wi=Ii cos1 由此可以得到反射率、透射率的表达式分别为 R?Wr?r2 Wi T?Wtn2cos?22?t Win1cos?1 将菲涅耳公式代入，即可
8、得到入射光中s分量和p分量的反射率和透射率的表示式分别为 sin2(?1?2) Rs?r?2sin(?1?2)2 s tan2(?1?2) Rp?r?2tan(?1?2)2 p Ts?n2cos?22sin2?1sin2?2ts? 2n1cos?1sin(?1?2) Tp? 显然有 n2cos?22sin2?1sin2?2tp?2 2n1cos?1sin(?1?2)cos(?1?2) Rs?Ts?1 Rp?Tp?1 综上所述，光在介质界面上的反射、透射特性由三个因素决定：入射光的偏振态，、入射角、界面两侧介质的折射率。 反射率随入射角的变化关系见图2-7. 图2-7 R随入射角1变化曲线 三实
9、验流程 2、实验程序 光密到光疏 clc clear all n01=1.52; n02=1; thta1=0:pi/90:pi/2; %以每两度取点 thta2=asin(sin(thta1)*1.52/1); %a为透射角，b为入射角 %c=b*pi/180; n=46; for i=1:n if thta1(i)>(41.8*pi/180); %rs1(i)=(1.52*cos(thta1(i)-1*cos(thta2)./(1.52*cos(thta2(i)+1*cos(thta2); rs=abs(1.52*cos(thta1)-1*cos(thta2)./(1.52*cos(
10、thta1)+1*cos(thta2); rp=abs(1*cos(thta1)-1.52*cos(thta2)./(1*cos(thta1)+1.52*cos(thta2); ts=abs(2*1.52*cos(thta1)./(1.52*cos(thta1)+1*cos(thta2); tp=abs(2*1.52*cos(thta1)./(1*cos(thta1)+1.52*cos(thta2); else rs=(1.52*cos(thta1)-1*cos(thta2)./(1.52*cos(thta1)+1*cos(thta2); rp=(1*cos(thta1)-1.52*cos(t
11、hta2)./(1*cos(thta1)+1.52*cos(thta2); ts=(2*1.52*cos(thta1)./(1.52*cos(thta1)+1*cos(thta2); tp=(2*1.52*cos(thta1)./(1*cos(thta1)+1.52*cos(thta2); end end subplot(2,2,1) plot(thta1,rs,g-*) hold on plot(thta1,rp,g:*) hold on plot(thta1,ts,b-*) hold on plot(thta1,tp,b:*) hold on xlabel(入射角); ylabel(反射系
12、数/透射系数); legend(rs,rp,ts,tp) title(反射光与透射光振幅的变化,fontname,宋体,color,blue,fontsize,16); grid on %第二个图rs %m=43; %布儒斯特角大概为30度 for i=1:n rs=(1.52*cos(thta1)-1*cos(thta2)./(1.52*cos(thta2)+1*cos(thta2); if thta1(i)<=(41.8*pi/180); %小于全反射角时 ph(i)=0; %elseif thta1(i)>30*pi/180&&thta1(i)>41.8*pi/180; % ph=0; else % thta1(i)>41.8*pi/180; ph(i)=angle(rs(i); end end subplot(2,2,2) plot(thta1,ph,b-+) set(gca,YDir,reverse) legend(ph(rs) title(反射光rp的相位变化,fontname,宋体,color,blue,fontsize,16); %第三个图 %m=32; %布儒斯特角大概为30度 for i=1:n rp=(n02*cos(thta1)-n01*cos(thta2)./(n02*cos(thta1)+n
