基于matlab的talbot效应光栅的计算机模拟[1],基于matlab的talbot效应光栅的计算机模拟.doc...

基于matlab的talbot效应光栅的计算机模拟.doc

30D0618F8C4DC3CCF74DBDD562A550D5PDF1基于MATLAB的TALBOT效应光栅的计算机模拟四川省通江县陈河小学张洪摘要计算机模拟技术广泛应用在教学和科研中，在TALBOT效应实验中引入计算机能更深动和深入揭示光学现象的物理内涵，本文提出利用MATLAB模拟TALBOT效应光学实验方法，例如用条纹光栅和网格光栅模拟的TALBOT效应，模拟的结果很好，还可以简单分析TALBOT效应在不同距离所产生的不同现象。该方法的优点是操作简单灵活，能完成一般光学中较难实现的操作，并给出TALBOT效应的实验结果。关键词计算机模拟；MATLAB；光栅；TALBOT效应；自成像绪论在大学理工科课程的教学中引入计算机模拟技术正日益受到重视,关于物理图形比如光学实验图像的模拟,已经有人做过一些工作，但还较少见到有人利用MATLAB软1,2件做模拟光学实验的工作。光学信息处理具有容量大、速度快、装置简单等优点,但在复杂和精密光路设计过程中为了获得最佳的光学信息处理效果,需要进行纷繁的计算和困难的实验验证。而计算机模拟技术可以为光路设计、相关光学器件的设计以及图像处理提供一条简捷的途径。利用MATLAB模拟光学实验简单灵活,操作者可以在计算机上自由设计图形的形状、尺寸以及实验参数变量,不受实际实验室条件的限制,能完成实际中较难完成的操作,从而加深了对物理原理、概念和图像的理解。因此,在3TALBOT效应光学实验中引入计算机模拟技术是必要的。介于此,本文介绍怎样利用MATLAB来模拟TALBOT效应光学实验,并给出了利用此模拟技术解决的TALBOT效应光学实验的模拟,以及TALBOT效应光学原理的相干光学信息处理应用的实验模拟。第一章TALBOT效应1TALBOT效应的原理1．1TALBOT效应的概念通常的情况下要得到一个物体的光学像，一个光学成像系统是必不可少的，最简单的是用一个透镜将物体发出的或散射的光聚集在像面上而得到物体的像。但是，在1836年TALBOT却发现当用相干光照射光栅时，在离光栅某些特定的距离上，能够形成光栅的像，这一现象称为TALBOT效应。它告诉我们在相干光场中，周期性的物体西南石油学院2005届本科毕业设计2能自成像，称为无透镜自成像或TALBOT效应自成像。TALBOT效应有许多有意义的应用。例如可以用来检验和复制衍射光栅，确定光束的准直性，实现图像相减以及构成TALBOT干涉仪检测位相物体；借助于叠栅条纹技术来确定出TALBOT像面位置，从而由TALBOT成像位置公式计算出光波波长等。第二章计算机模拟ＴALBOT２计算机模拟方法2．1常用的傅立叶变换光路光路系统(4F系统)准直透镜；为傅立叶透镜。1L23,L频谱面是平面，空间频率1XF。(21)1YF物像比例关系 201FDXYM(22)2．2傅立叶变换的定义与MATLAB计算4设序列长度为，则的点离散傅立叶变换对定义为XNMXNNM(210,0,12,KNNXKDFTWN3)(210,01,2NKNKXNITX4)其中，，称为变换区间长度。2JNWEDFT30D0618F8C4DC3CCF74DBDD562A550D5PDF3可以把X(K)这个式子写成矩阵乘法运算(2XKXNW5)其中，为序列行向量，是一阶方阵，通常称为旋转因子矩阵。XNNN(20,12,1XX6)(2001011011NNNNNWWNK7)(3)式用MATLAB矩阵运算表示为(2 01WNKN8)因此，可得到用矩阵乘法计算点的程序如下。DFT用矩阵乘法计算N点DFTCLEARALLXN‘请输入序列X’NLENGTHXNN0N1KNNKN‘KWNEXPJ2PI/NWNKWNNKXKXNWNK只要输入序列，运行该程序，即XN可实现的点，这种计算离散傅立XNDFT叶变换的方法概念清楚，编程简单，TALBOT效应的计算机模拟主要过程就是傅立叶正变换和逆变换，所以要清楚它的计算机实现过程。MATLAB基础部分提供了FFT，IFFT，FFT2和IFFT2等快速傅立叶变西南石油学院2005届本科毕业设计4换函数，它使的运算速度量提高了若干数量级，我们在进行傅立叶正变换和逆变DFT换时，只要调用FFT和IFFT及FFT2和IFFT2，高维则用FFTN和IFFTN就行了。在光线通过图13中的扩束镜C的时候，扩束镜C就是光空间的时域特性变成空间的频域特性，就是把光经过一次二维的傅立叶变换(FFT2)，当光再次通过准直镜L时，就把光的空间的频域特性变成空间的时域特性，也就在经过一次二维反傅立叶变换(IFFT2)，TALBOT就是基于这种特性形成无透镜成像的。23计算机模拟过程根据公式(11)，我们要产生一个周期性光栅，但是公式中的很难确定，即使MC确定了，计算投射系数也很难，我们就想办法避开计算也能达到相同的MC0GX0GX计算效果。计算机模拟TALBOT效应程序流程如图22所示。具体操作过程是物体图像的生成可以直接由矩阵运算生成,也可利用WINDOWS下的画图工具,生成一幅二值图像属性为黑白,并将建立的图像存入MATLAB的工作目录\WORK下然后调用命令函数IMREAD输入图像,输入图像是一个二维矩阵。在此，我们的周期性光路图像直接由矩阵运算生成。例如模拟周期性条纹光栅程序FZEROS1000,1000FORI010F20I22990I,1ENDIMSHOWF在这里，有必要对各个语句作用进行说明FZEROS1000,1000是产生一个的全0矩阵，也就是在MATLAB图形窗1口显示全黑的图形。语句结尾以分号结束，这样抑制屏幕输出，可以节省大量的运算时间。FORI010是一个循环语句，结构形式为FOR，，END的格式，进行循环语句计算。F20I22990I,1在的全0矩阵中，令F的值从20行开始，循环111次，间隔为22行，一直到990I结束的值为1，意思就是光栅缝，让光通过，其余不让光通过，这样就形成周期性光栅。END它和FOR语句形成循环嵌套，要有END语句，程序才能结束循环。IMSHOWF一幅图形的方式显示F图像。30D0618F8C4DC3CCF74DBDD562A550D5PDF5全零矩阵和光栅图像如下2．3．1模拟过程因为TALBOT效应需要周期性光栅，最直接简单的方法就是用MATLAB程序编写一个周期性光栅，如前面所述，先建立一个周期性条纹光栅，检验模拟TALBOT效应是否成立。设一维周期光栅的透射系数是一个矩形波函数，其复振幅透过率为(200EXP2,12,MGXCJMD9)在这里，我们用光栅模拟了，对作傅立叶变换得物分布的空间频F0GX0Y率为(20MGCD10)由菲涅耳其程序如下所示(22EXPHJKZJZ11)在这里，我们对传递函数作了特定的处理，也就是在TALBOT距离看到光栅像西南石油学院2005届本科毕业设计6的条件，当时，，这时传递函数21,3NDZ2EXP1JZ(2EXPHJKZ12)这里为空间角频率，为TALBOT距离，为1。KZN空间频率和传递函数作用