matlab进行波动光学的计算机模拟..doc

Matlab进行波动光学的计算机模拟 简述：在教学和科研中，用Matlab进行计算机模拟正越来越被重视。分析讨论了夫琅禾费衍射以及双光束,多缝,牛顿环和干涉平面光栅的衍射等理论，用Matlab编写出相应程序后，进行了计算机模拟，这有助于理解和研究衍射和干涉理论。 1功能介绍：光的波动性常常表现为干涉、衍射、偏振等现象。虽然关于这些现象的描述及其阐述有好多，但是未能配以精彩的直观形象图形。运用MATLAB获得模拟图形，将定性的语言描述和抽象的数学表示变为生动直观的表现，可以进一步分析和描述有关波动光学的现象和规律等理论。 光的衍射现象研究,不仅具有重要的理论意义,而且在光学仪器研制和成象分析等均有重要价值.针对衍射计算困难的问题,选择合适的计算公式,并基于科学计算软件MATLAB5.3编写计算程序,使计算得变得简洁,大大提高了实验的效率.用MATLAB编程得出的计算结果与用数码相机记录的衍射图样进行比较,两者吻合得很好。 2 功能实现：运用MATLAB获得模拟图形，将定性的语言描述和抽象的数学表示变为生动直观的表现，可以进一步分析和描述有关波动光学的现象和规律等理论。使用M文件形式。M文件有两个形式：一种称为命令文件（scriptfile）就好像dos下的批处理文件一样。这类程序包含了一连串的MATLAB命令，执行时依次执行。另一种称为函数文件（functionfile）。它的第一句可执行语句是以function引导的定义语句，在函数文件中的变量都是局部变量。 2.1 夫琅禾费衍射： 夫琅禾费衍射，是认为光源和观察屏离衍射屏（孔）处于无穷远处衍射现象。实验装置如图? S为单色点光源，放置在透镜L1的物方焦点处，所得平行光垂直入射到障碍物，借助于透镜L2将无穷远处的衍射图样移至L2的像方焦面上观察。? 将衍射屏制作成输入图像，用imread()函数读入，然后利用傅里叶变换函数fft2()对其进行傅里叶变换，得到其傅里叶频谱。由函数fft2()实现的傅里叶变换频谱的直流分量位于图像的左上角，而由透镜实现的光学傅里叶变换的直流分量位于图像中心。因此，为了得到模拟的光学傅里叶变换，需调用函数fftshift()将零频移到频谱中心。? 2.2双光束干涉: 双光束干涉，有分波面类型的杨氏干涉、劳埃镜干涉、菲涅耳双棱镜和双面镜干涉等，还有分振幅类型的麦克耳孙干涉等。 下图所示的是双缝干涉装置，是分波阵面的双光束干涉的典型实例。下面分析它的干涉原理与干涉条纹特点。 如图7.1.1所示，两狭缝间距为d，双缝所在平面与屏幕平行，两者之间的垂直距离为D，O为屏幕上的坐标原点且与两狭缝对称。两个狭缝光源满足振动方向相同、频率相同、相位差恒定的相干条件。故两列光在空间相遇，将产生干涉现象。屏幕上将出现明暗相间的干涉条纹。设OP=y，则由几何关系可知，两个相干光源到达屏幕上任意点P的距离分别为：? r1=?D2+?y?d2^2 r2=?D^2+?y+d2?^2 ?这样两列相干光的光程差为?r=r2?r1，相位差为?φ=2π?r/λ。设两狭缝光源的光波单独到达屏幕P点处的振幅分别为E1和E2，光强分别为I1和I2。则两 列光波叠加后的振幅为E=?E1^2+E2^2 +2E1*E2cos?φ：，而叠加后的光强为： I=I1+I2+2?I1*I2cos?φ。设两列光波在屏幕上相遇点振幅相等，则P点光强为：I=4I0?cos ?φ2^2 ?φ=2kπ,k=0,1,2,3?时为干涉明条纹。??φ=?2k+1?π,k=0,1,2,3?时为干涉暗条纹。?? 其实绝对的单色光只是理论上存在，实际的“单色光”都是有一定谱线宽度的准单色光，这种准单色光入射到干涉装置后，其中的每一种波长成分都将产生自己的干涉条纹，由于波长不同，除零级明条纹外，其他级次的条纹将彼此错开，并发生不同级条纹的重叠。在重叠处总的光强为各种波长的条纹的光强的非相干相加。随着级次的增大，干涉条纹的明暗对比减小，当级次增大到某一值后，干涉条纹就消失了。对于谱线宽度为?λ的准单色光，干涉条纹消失的位置应是波长?λ+?λ2的成分的k级明条纹与波长为λ??λ2的成分的k+1级明纹重合的位置。由于两成分在此位置上有同一光程差，根据光程差与明纹级次的关系可知，条纹消失时，应满足：(λ+?λ2)k=?λ??λ2?k+1?。由于?λ?λ，于是可得：k=λ ?λ由此可知，?λ愈大，即光的单色性愈差，能观察到得干涉条纹级次就愈小。 clearall a=-4*pi:0.01*pi:4*pi;%通过改变PI的倍数而改变条纹数P=1-sin(2*a).^2./sin(a).^2; %当要求P的曲线分布图时P=sin(2*a).^2./sin(a