四．实验结果和分析
由光疏到光密的仿真结果
图1我们可以看出rs和rp及ts和tp的变化，也可以清楚的看出在布儒斯特角处rp分量为零，这正好验证了在布儒斯特角处无p分量，因为反射率是反射系数的平方，反射系数为零，所以反射率就为零。
图2我们可以看出rs的相位是一直没有变的，因为在菲涅耳公式中rs为两个复振幅之比，我们可以根据第一个图看出来，它们之比的符号一直未变，而一个确定的复振幅可以写成余弦函数的形式，+π或-π时它们都会变，所以综上我们可以推出rs的相位没有变化，而图中也证明了我们的推断是正确的。
图3，反应了在布儒斯特角出它的相位发生了π的跃变，而根据一个确定的波的表达式来看它是由余弦函数的的变化来确定的，而rp在菲涅耳表达式中是两个确定的余弦函数之比，所以rp由正变为负的时候，其中有一个余弦函数肯定相位发生了变化（奇变偶不变，符号看象限），正好可以说明此问题。
由光密到光疏的仿真结果
11 / 42
图1我们可以看出rs和rp及ts和tp的变化，也可以看到当r分量变为1时，p分量变为0，但不是直接变为0，用衰逝波的存在也可以解释该现象，也可以验证透射率与反射率这和为1，的定义。
图2我们可以看出rs的相位在全反射角处发生了变化，而且是慢慢变到-π的，此时入射角已经逐渐大于全反射临界角了，它满足一个全反射相位变化公式，随着入射角的增大，相位是一个逐渐变化的过程。
图3，反应了在布儒斯特角处它的相位发生了π的跃变，而根据一个确定的波的表达式来看它是由余弦函数的的变化来确定的，而rp在菲涅耳表达式中是两个确定的余弦函数之比，所以rp由正变为负的时候，其中有一个余弦函数肯定相位发生了变化（奇变偶不变，符号看象限），且在布儒斯特角处，而在全反射角处也会发生变化，而且是逐渐变化的，这是因为当入射角逐渐增大的时候，它满足一个公式tan(fai/2)=-√((sinθ)^2-n^2)/cosθ),从公式可以看出相位会随着入射角的变化而渐变，当θ=π/2时，tan(fai/2)为无穷，所以fai=π，由此可以推断出。 ． 思考题
1. 如何确定入射面?
答：入射光与反射光以及法线共同构成的平面即入射面 2.什么是临界角?临界角是光疏到光密，还是光密到光疏时发生?
答：临界角就是全反射角，他指的是光线由光密介质入射到光疏介质时正好发生全反射时的入射角。
3.利用全反射现象能否产生圆偏振光?
答；利用全反射现象可以产生圆偏振光，一个偏振光在一定角度上经过两次全反射可以产生圆偏振光，菲涅耳棱镜就是利用这个原理所制成的。
4.解释反射系数及透射系数的概念。
答：当电磁波由一个磁导率为μ1、介电常数为ε1的均匀介质，进入另一个具有磁导率为μ2、介电常数为ε2的均匀介质时，一部分电磁波在界面上被反射回来，另一分电磁波则透射过去。反射波与透射波的振幅同入射波振幅之比，分别称之为反射系数与透射系数。
5.根据仿真曲线解释反射及透射光的相位变化规律。
答：图中反应了他们的相位的变化规律，例如图三所示在布儒斯特角处它的相位发生了π的跃变，而根据一个确定的波的表达式来看它是由余弦函数的的变化来确定的，而rp在菲涅耳表达式中是两个确定的余弦函数之比，所以rp由正变为负的时候，其中有一个余弦函数肯定相位发生了变化（奇变偶不变，符号看象限），且在布儒斯特角处，而在全反射角处也会发生变化，而且是逐渐变化的，这是因为当入射角逐渐增大的时候，它满足一个公式tan(fai/2)=-√((sinθ)^2-n^2)/cosθ),从公式可以看出相位会随着入射角的变化而渐变，当θ=π/2时，tan(fai/2)为无穷，所以fai=π。
6.试说明布儒斯特角的概念。
12 / 42
答：布儒斯特角，又称偏振角，是自然光经电介质界面反射后，反射光为线偏振光所应满足的条件。
7.试分析布儒斯特角与临界角哪个大。
答：临界角大于布儒斯特角，我们从它们的公式可以简单的推导出来，布儒斯特角为arctan(n2/n1),全反射角为arcsin(n2/n2), 假设n2/n1=x,因为有光密入射到光疏，所以n2>n1，因此x>1，此时布儒斯特角为arctan(x),全反射角为arcsin(x),我们对它两个同时求导得到：（arctan(x)）’=1/(1+x^2),而（arcsin(x)）’=1/√(1+x^2),由此我们可以得出全反射角公式的倒数大，也就是说，在相同变量的情况下它的数值大，从而我们也就说明了临界角大于布儒斯特角。
8.布儒斯特角都有哪些应用？
答：布儒斯特角可以用于生产墨镜的技术中，可以利用一定的技术设计出让p分量尽量的多的透过墨镜，不至于反射光过强导致视线不明确等。 五．实验心得
这是我们所做的第一个实验，也是遇见困难最多的一个实验，开始上手时我发现自己对MATLAB还并不是那么的熟，前两天发现自己几乎没有做什么，时间过得很快，我觉得周围人都做出来至少一个了，我发现自己似乎很笨，我们的实习第一个星期是在上午，第三天的上午我又没弄出来，下午两点多我灰溜溜的和大家一起回到了宿舍，我觉得自己不应该比人差，于是我休息了一个小时，带着电脑又去了实验室，实验室里只有三班同学和我，也相对于早上安静了许多，于是我找了一个人少的地方坐了下来，静静的一个人在那儿检查自己的程序，一个一个的检查，运行，几乎每次都会有错误，语法错误较少，但是有很多MATLAB潜在的知识我知道的还有些少，晚上大家都吃饭去了，我的程序还没好，我请教了刘老师，刘老师给我指点了很多，我晚上回宿舍时也上网查相关的知识，在这期间我也懂得的去查workspace的重要性，其实我的错误很明显，只是我没有想到而已，在if语句中，我本来应该用数组变量空间，而我却用了一个数空间，这就导致了后边在执行赋值语句是，后边的值会覆盖前边的值，导致出来的图像是一个单值图像，我后来在workspace中发现了这点，最后正确的图像终于出来了，我也觉得收获很大，不仅在光学的知识上同时也在MATLAB的知识上，可以说是双重收获，也让我懂得了遇到不会的只要用心去钻研，一定可以收获到很多。
13 / 42
双光束干涉的仿真
一、实验目的
1.掌握光的相干条件；
2.掌握分波阵面双光束干涉的特点。
二、实验原理
1. 两束光的干涉现象
光的干涉是指两束或多束光在空间相遇时， 在重叠区内形成稳定的强弱强度分布的现象。
例如，图5-1所示的两列单色线偏振光
E1?E01cos(?1t?k1?r??01)E2?E02cos(?2t?k2?r??02)
图5-1 两列光波在空间重叠
在空间P点相遇，E1与E2振动方向间的夹角为θ，则在P点处的总光强为
I?I1?I2?2I1I2cos?cos??I1?I2?2I12
式中，I1、I2是二光束的光强；φ是二光束的相位差，且有
??k2?r?k1?r??01??02???t ????1??2I12?I1I2cos?cos?由此可见，二光束叠加后的总强度并不等于这两列波的强度和，而是多了一项交叉项I12，它反映了这两束光的干涉效应，通常称为干涉项。干涉现象就是指这两束光在重叠区内形成的稳定的光强分布。所谓稳定是指，用肉眼或记录仪器能观察到或记录到条纹分布，即在一定时间内存在着相对稳定的条纹分布。显然，如果干涉项I12远小于两光束光强中较小的一个，就不易观察到干涉现象；如果
14 / 42
两束光的相位差随时间变化，使光强度条纹图样产生移动，且当条纹移动的速度快到肉眼或记录仪器分辨不出条纹图样时，就观察不到干涉现象了。 在能观察到稳定的光强分布的情况下，满足
1, ±2，… ??2m? m=0, ±
的空间位置为光强极大值处，且光强极大值IM为
IM?I1?I2?2I1I2cos?
满足
φ=(2m+1)π m=0, ±1, ±2，
的空间位置为光强极小值处，且光强极小值Im为
Im?I1?I2?2I1I2cos?
当两束光强相等，即I1=I2=I0时，相应的极大值和极小值分别为
IM=2I0(1+cosθ) Im=2I0(1-cos θ)
2. 产生干涉的条件
首先引入一个表征干涉效应程度的参量——干涉条纹可见度，由此深入分析产生干涉的条件。
1) 干涉条纹可见度(对比度) 干涉条纹可见度定义为
V??IM?Im
IM?Im当干涉光强的极小值Im=0时，V=1，二光束完全相干，条纹最清晰；当IM=Im时，V=0，二光束完全不相干，无干涉条纹；当IM≠Im≠0时，0＜V＜1，二光束部分相干，条纹清晰度介于上面两种情况之间。
2) 产生干涉的条件
由上述二光束叠加的光强分布关系可见，影响光强条纹稳定分布的主要因素是：二光束频率；二光束振动方向夹角和二光束的相位差。
(1) 对干涉光束的频率要求
由二干涉光束相位差的关系式可以看出，当二光束频率相等，Δω=0时，干涉光强不随时间变化，可以得到稳定的干涉条纹分布。当二光束的频率不相等，Δω≠0时，干涉条纹将随着时间产生移动，且Δω愈大，条纹移动速度愈快，当Δω大到一定程度时，肉眼或探测仪器就将观察不到稳定的条纹分布。因此，为了产生干涉现象，要求二干涉光束的频率尽量相等。
(2) 对二干涉光束振动方向的要求
15 / 42
西安邮电大学光学实验matlab仿真结果分析与程序
