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

Matlab进行波动光学的计算机模拟 在教学和科研中用Matlab进行计算机模拟正越来越被重视。分析讨论了夫琅禾费衍射以及双光束多缝牛顿环干涉等理论，用Matlab编写出相应程序后，进行了计算机模拟，这有助于理解和研究衍射和干涉理论 S为单色点光源，放置在透镜L1的物方焦点处，所得平行光垂直入射到障碍物，借助于透镜L2将无穷远处的衍射图样移至L2的像方焦面上观察。? 将衍射屏制作成输入图像，用imread()函数读入，然后利用傅里叶变换函数fft2()对其进行傅里叶变换，得到其傅里叶频谱。由函数fft2()实现的傅里叶变换频谱的直流分量位于图像的左上角，而由透镜实现的光学傅里叶变换的直流分量位于图像中心。因此，为了得到模拟的光学傅里叶变换，需调用函数fftshift()将零频移到频谱中心。? 2.2双光束干涉: 双光束干涉，有分波面类型的杨氏干涉、劳埃镜干涉、菲涅耳双棱镜和双面镜干涉等，还有分振幅类型的麦克耳孙干涉等。 下图所示的是双缝干涉装置，是分波阵面的双光束干涉的典型实例。下面分析它的干涉原理与干涉条纹特点。 如图7.1.1所示，两狭缝间距为d，双缝所在平面与屏幕平行，两者之间的垂直距离为D，O为屏幕上的坐标原点且与两狭缝对称。两个狭缝光源满足振动方向相同、频率相同、相位差恒定的相干条件。故两列光在空间相遇，将产生干涉现象。屏幕上将出现明暗相间的干涉条纹。设OP=y，则由几何关系可知，两个相干光源到达屏幕上任意点P的距离分别为：? r1=?D2+?y?d2r2=?D^2+?y+d2?^2 ?这样两列相干光的光程差为?r=r2?r1，相位差为?φ=2π?r/λ。设两狭缝光源的光波单独到达屏幕P点处的振幅分别为E1和E2，光强分别为I1和I2。则两 列光波叠加后的振幅为E=?E12+E2^2 +2E1*E2cos?φ：，而叠加后的光强为： I=I1+I2+2?I1I2cos?φ。设两列光波在屏幕上相遇点振幅相等，则P点光强为：I=4I0?cos ?φ22 ?φ=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).^2plot(a,P) lgray=zeros(256,3);fori=0:255 lgray(i+1,:)=(255-i)/255;endimagesc(P)colormap(lgray) 使用imshow建立灰度图像的使用格式：imshow(I,N)?。其中参数N为正整数，指定灰度的层次，当缺省该参数时，系统默认为256级的灰度级；参数I为数值矩阵，imshow的作用就是将数值矩阵I的元素值用N个灰度级的黑白图像可视化显示出来。实际上是在数值矩阵I和N个灰色调之间建立了一种颜色映射关系：I当中元素值最大者映射为白色（将该元素值作为白色显示），元素值最小者映射为黑色（将该元素值作为黑色显示），元素值介于最大和最小之间的则按照某种约定的规则映射到其它的灰度级（显示为不同灰度的灰色）。? ?????为了方便叙述，假定再现的图像尺寸2mm×2mm,使用上述指令可以很方便的将牛顿环干涉条纹在该区域内再现。 首先，利用(4)式获取干涉光强I在该区 （0.0010.001,0.0010.001?xy假定观察屏是xy平面）的数值分 布?x=linspace(-0.001,0.001,200);?y=linspace(-0.001,0.001,200);? [X,Y]=meshgrid(x,y);?%将xy平