《快速傅里叶变换计算衍射光强的分布本科毕业论文（设计）.doc》由会员分享，可在线阅读，更多相关《快速傅里叶变换计算衍射光强的分布本科毕业论文（设计）.doc（26页珍藏版）》请在三一办公上搜索。
1、快速傅里叶变换计算衍射光强的分布目录快速傅里叶变换计算衍射的光强分布40.引言41. 空域连续函数的离散及延拓52. 离散傅里叶变换与傅里叶变换的关系63.快速傅里叶变换计算衍射光强1131单缝衍射1332 圆孔衍射144. 光强分布曲线154.1单缝衍射的光强分布曲线154.2圆孔衍射的光强分布175.结论20参考文献21附录211.用MATLAB软件模拟单缝衍射和光强分布曲线的程序212.用MATLAB软件模拟圆孔衍射和光强分布曲线的程序22致 谢24河西学院本科生毕业论文（设计）题目审批表25河西学院 学院指导教师指导毕业论文情况登记表26河西学院毕业论文（设计）指导教师评审表27河西学
2、院本科生毕业论文（设计）答辩记录表28河西学院本科生毕业论文诚信声明本人郑重声明：所呈交的本科毕业论文（设计），是本人在指导老师的指导下，独立进行研究工作所取得的成果，成果不存在知识产权争议，除文中已经注明引用的内容外，本论文不含任何其他个人或集体已经发表或撰写过的作品成果。对本文的研究做出重要贡献的个人和集体均已在文中以明确方式标明。本人完全意识到本声明的法律结果由本人承担。 作者签名： 二O 一三 年 五 月 二十六 日（打印） 河西学院本科生毕业论文开题报告论文题目快速傅里叶变换计算衍射的光强分布学生姓名所属学院物理与机电工程学院专业年级指导教师所在单位河西学院职称开题日期选题的根据：衍
3、射现象是波动光学中的重要知识，光的衍射的定义指出光的衍射是一种区别于几何光学规律的光的传播现象.在现代光学发展的今天，学习如何计算衍射的光强分布是必要的.综述有关本选题的研究动态和自己的见解：光的衍射的研究是在十七世纪意大利物理学家和天文学家F.M.格里马尔迪首先发现光遇障碍物时将偏离直线传播，他把此现象起名为“衍射”.胡克和R.玻意耳分别观察到现称之为牛顿环的干涉现象.之后托马斯扬在公元1817年提出了光是横波的假说，法国的土木工程师菲涅尔以此为基础，在1818年提交了一篇应征巴黎科学院悬赏征求阐述光折射现象的论文，在这篇文章中，他提出了一整套高度完善的波动说理论阐明了这一现象。1965年,
4、由库利图基(CooleyTukey)提出的PFT技术彻底改变了这种状况，计算机的普及应用为这种快速计算方法的推广创造了良好的条件因此、利用FFT技术计算衍射的方法逐渐被广泛采用然而，必须指出由于快速傅里叶变换只是离散博里叶变换的一种快速算法，对离散傅里叶变换理论进行研究后很快就能发现，只有当被变换的函数是在频域有限区域存在的“带限函数”时、连续函数的傅里叶变换才能由离散傅里叶变换表述否则，由于频谱的混叠效应，离散傅里叶变换只是连续函数傅里叶变换的一种近似不幸的是，衍射计算问题中所遇到的函数基本上都不是带限函数、因此利用快速博里叶变换对衍射所作的计算也是一种近似只有了解离散傅里叶变换与傅里叶变换
5、的关系，通过合适的离散、较好地将频谱混叠的影响控制在允许的误差范围才能使用这种计算方法得到较好的结果主要内容及其主要的研究方法：本文将要介绍光的衍射现象，根据理论解释，计算衍射光强的分布，让一束光经扩束器或不经扩束器直接投射到一狭缝（或圆孔）上，在狭缝（或圆孔）后面放置一个透镜，那么在透镜的焦平面上放置的屏幕上将产生明暗相间的条纹（明暗相间的圆环）.利用傅里叶变换计算分析光在屏幕上的衍射现象的特点及光强分布情况，并用MATLAB软件仿真模拟衍射现象.论文进度安排和采取的主要措施：完成期限：2013.32013. 6主要措施：1. 依据快速傅里叶变换计算衍射的的光强分布。2. 在研究单缝衍射的基
6、础上进行研究圆孔衍射。3. 在计算的基础上用MATLAB分别模拟仿真了单缝和圆孔衍射现象，及衍射的光强分布。4. 在整个研究过程中，发现问题、解决问题、最终得出结论。主要参考资料和文献：1姚启钧.光学教程M.北京：高等教育出版社，2008:106-109.2陈聪、李定国.基于快速傅里叶变换的衍射现象的数值仿真J.大学物理2004,9:47-623李俊昌、熊秉衡.信息光学理论与计算M.科学出版社M.2009:93-1104董克剑.利用MATLAB模拟光的衍射现象J.物理教师,2008,29(5),30-32.5王竞争等.基于MATLAB的光的干涉和衍射现象的模拟研究J.延边大学学报,2009,3
7、5(4),319-322.6喻力华、赵维义.圆孔衍射光强分布的数值计算J.大学物理2010,1,20(1),17-19指导教师意见：签 名： 年 月 日教研室意见负责人签名：年 月 日学 院 意 见负责人签名：年 月 日快速傅里叶变换计算衍射的光强分布摘要：本文利用快速傅里叶变换计算了光的单缝和圆孔衍射的光强分布，根据计算结果利用MATLAB软件仿真模拟了单缝和圆孔衍射现象.分析表明，衍射图样取决于缝宽或孔径的大小，它反映了障碍物和光波之间限制和扩展的辩证关系，限制范围越小，扩张现象愈显著；在哪个方向上限制，就在该方向上扩展. 且在处理实际问题时应合理选择两种算法SFFT，DFFT关键词：离散
8、傅立叶；快速傅里叶变换；衍射光强的分布；快速卷积算法Abstract: By using fast Fourier transform, this paper calculates the light intensity distribution of the single-slit and circular aperture diffraction. And, according to the calculation results simulates the single slit and circular aperture diffraction phenomenon by using
9、MATLAB software. The simulation analysis showed that the diffraction pattern depends on the size of slit or the width of aperture. It reflects the dialectical relationship of restrictions and extensions between obstacles and light. And the smaller limit rang, the more remarkable expansion phenomenon
10、; it extends in the direction which is the direction of restrictions, and in dealing with practical problems, it should be a reasonable choice to use two kinds of algorithms, S-FFT and D-FFT.Keywords: discrete Fourier transform; fast Fourier transform; the light intensity distribution of diffraction
11、; fast convolution algorithm0.引言1965年,由库利图基(CooleyTukey)提出的PFT技术彻底改变了这种状况，计算机的普及应用为这种快速计算方法的推广创造了良好的条件因此、利用FFT技术计算衍射的方法逐渐被广泛采用然而，必须指出由于快速傅里叶变换只是离散博里叶变换的一种快速算法，对离散傅里叶变换理论进行研究后很快就能发现，只有当被变换的函数是在频域有限区域存在的“带限函数”时、连续函数的傅里叶变换才能由离散傅里叶变换表述否则，由于频谱的混叠效应，离散傅里叶变换只是连续函数傅里叶变换的一种近似不幸的是，衍射计算问题中所遇到的函数基本上都不是带限函数、因此利用
12、快速博里叶变换对衍射所作的计算也是一种近似只有了解离散傅里叶变换与傅里叶变换的关系，通过合适的离散、较好地将频谱混叠的影响控制在允许的误差范围才能使用这种计算方法得到较好的结果为了能够正确使用FPT计算衍射p在具体阐述计算方法之前，有必要对二维空间函数的取样、延拓及离散傅里叶变换与傅里叶变换的关系进行研究1. 空域连续函数的离散及延拓函数作二维离散傅里叶变换时，要求是被变换函数是二维空间的周期离散函数由于实际需要作博里叶变换的函数通常是在空域无限大平面上均有定义的连续函数，于是，必须将函数截断在有限的区域进行取样及延拓通常的取样方法是、先将函数的主要部分通过坐标变换放在第一象限、并沿平行于坐标
13、轴的方向将函数截断在一个的矩形区域内；然后取样周期为，从坐标原点开始将函数离散为从个点的二维离散分布值、图l(a)、(b)描述了上述过程(图中用黑点标注出取样点落在函数定义区域上的位置，用小圆圈表示取样为零的位置)图l(c)是二维周期延拓结果图1.空域连续函数的离散及延拓2. 离散傅里叶变换与傅里叶变换的关系很明显，函数经截断及离散处理后无论在空域还是在频域均会引入误差现以z方向的博里叶变换为例进行研究，以后再将结果推广到二维空间图2示出于某一给定的，函数沿s方向进行离散傅里叶变换的过程图中，左边为一列空域的原函数图像，右边一列图像是它们的频谱的模，符号“”表示它们为傅里叶变换对例如，图2(a
14、1)为空域的原函数，图2(a2)为它的频谱的模 对未经截断函数的取样，等于用图2(b1)的梳状函数乘以图2(a1)的原函数，数学表达式为 （1）由于梳状函数为周期的函数，可以表为博里叶级数： （2）其中,于是, （3）被信号调制的结果(见图2(c1)图2. 函数的离散傅里叶变换过程上式表明，取样信号已经不是原信号，而是无穷多个截波信号 （4）现在，通过博里叶变换来考察信号经取样后的频谱与原信号频谱的关系对上式作博里叶变换得 （5）结果表眼在取样信号频谱中除了包含原信号频谱外，还包含了无穷多个被延拓的频谱、延拓的周期为(见图2(c2)并且、由于原函数的频谱宽度大于延拓的周期，相邻的频谱曲线产生了
15、混叠 根据傅里叶变换中频域的卷积定律，图2(c2)也可以通过原函数的频谱函数(图2(a2)与梳状函数的频谱函数(2(b2)的卷积求出: （6）为强调这个关系，图2(c2)的纵坐标由这个卷积表达式标注 由此可见、连续函数经过周期为的无穷序列取样离散后，其频谱与原函数频谱相比有两点区别：（1）频谱发生了周期为的周期延拓如果原函数的频谱宽度大于时,则产生频谱混叠，引入失真(2)离散信号频谱的幅度是原函数频谱的倍.然而上面对连续函数被无穷序列取样离散的后的频谱研究只是一个理论结果、因为实际上不可能作取样点为无限多的数值计算并是，由于离散傅里叶变换事实上讨论的是在空域及频域均是周期离散函数的傅里叶变换问
16、题还要将离散函数截断及延拓才能满足要求因此，将空域非周期的离散函数(图2(c1)先通过下述矩形宙函数图2(d1)截断: (7)得到具有从个点的的离散分布(图2(e1): (8)然后，再将截断后的部分进行周期为的延拓，形成图2(g1)的周期离散序列: (6)按照傅里叶变换理论，空域中矩形窗函数图今2(d1)与离散序列图2(c1)的乘积的频谱函数、可表为矩形函数的频谱函数(图2(d2)与图2(c1)的频谱函数图2(c2)的卷积： (9)对应的频谱函数曲线示于图2(e2).由图可见由于矩形宙函数的频谱具有较大的起伏变化的伤瓣，卷积运算的结果佼图2(c2)的频谱曲线形状产生了失真(为说明问既图中略有夸
17、大)将图2(c2)与图2(a2)比较不难发现,现在得到的是带有畸变的原函数频谱的周期延拓曲线，延拓周期为离散傅里叶变换是对空域及领域均为周期离散函数的变换，因此, 图2(e2)的曲线还将被周期为的梳状函数(图2(f1)取样其结果是一个周期为的频域的离散函数图2(g2).在频域进行上面频谱函数与梳状函数的乘积取样时，就对应着它们在空域原函数的卷积运算图2(e1)与图2(f1)的函数在空域卷积运算的结果成为一周期为的空域离散函数图2(g1)空域及领域离散函数均以为周期，我们只要分别知道一个周期内的离散值或样本点使可以了解离散函数全貌离散傅里叶变换或其快速算法FFT，便是完成从空域到频域、以及从频域
18、到空域的这个样本点的汁算方法 至此，我们已经知道，离散傅里叶变换是博里叶变换的一种近似计算只要能够将衍射的计算表为卷积的形式，并了解离散傅里叶变换与博里叶变换间的定量关系，采取合适的措施抑制晌曳便能对衍射问题求解将菲涅耳衍射积分的卷积形式: （10）中的二次项展开后得: （11）设为物平面光波复振幅；根据第（10）式，经距离d的衍射到达观测平面的光波复振幅可由下形式的菲涅耳衍射积分表出: （12）式中，为光波波长，.若利用快速傅里叶变换FFT进行计算-式物平面取样宽度为，取样数为,取样间距为，（12）式可写为: (13)式中，是离散傅里叶变换后对应的空域取样间距为确定这个数值，根据前面对离散傅
19、里叶变换的讨论，（13）式的计算结果将是取值范围的点的离教值即:或者， (14)因此 (15) 对（10）式两边作博里叶变换并利用空域卷积定律得: （16）令是频域坐标可以定义菲涅耳衍射传递函数为: （17）3.快速傅里叶变换计算衍射光强光是一种电磁波，按的规律随时间传播，电光源发粗的是一组球面波，设光源位于坐标原点处，以速度在电容率为的介质中传播，当光到达半径为的求面时，光的场强是的函数，可以表示为: （18）其中称为波数，为光矢量.点光源从原点出发的球面波，能量密度为: （19）以表示单位时间内光矢量所在空间的体积，则单位时间内通过整个球面的能量为: （20）而 （21）式中是与光源振动有
20、关的常数，是与介质有关的常数，则 （22）为简便，只考虑某时刻的振动，含时间的项可省去。在光学系统中，光从出射光瞳射出，取光瞳坐标为，观察平面的坐标为，两坐标系相平行，原点在它们的公共垂线上，相距为见图（3）图3 衍射光强分布图光瞳面上任意一点到观察面上的某点p的距离为: （23）由（22）式知，光是从s点以球面波的形式传播到p点的。如果s点振幅为,则在P点光的矢量为: （24）为计算球面上p点的光的场强，需要选取包含s点在内的小面元，则 （25）为便于计算，设光瞳与观察点面相距很远，取,的近似值为: （26）在远场衍射的情况下，即时1且，则 （27）式（27）为远场近似情况下的衍射的公式.那
21、么级光源在p 点的是光瞳函数的二维傅里叶变换式.31单缝衍射图4 单缝衍射原理图衍射装置如图12所示，衍射物体为不透明屏上一条方向平行于轴，长度不限，宽度为的狭缝.用轴上单色点光源S和推直透镜C产生的平行光正入射照明，在透镜L的后焦面上观察.在不考虑透镜L孔径大小的情况下，可认为光波在方向不受限制，所以上述单缝衍射实际是一维的问题，只需要计算沿方向的衍射.设衍射物体的复振幅透射系数，由于照明光波复振幅所以透过衍射物体的复振幅为: （28）应用公式(26)，可以直接计算平面衍射图形的复振幅为此，首先计算的傅里叶变换： （30）将代人公式 （31）并利用公式，给出空间频率与频率平面空间坐标的关系，
22、可得出单缝夫琅禾费衍射的复振幅分布为: （32）最后，利用，可算出观察面上衍射图形的幅照度分布为: （33）32 圆孔衍射根据式: （34）其中透镜焦平面上 点光场振幅经过坐标变换后为: （35）式中是衍射方向与光轴的夹角，称为衍射角.在这里。应用了 的近似关系.计算这个积分，得到: （36）因此点的光场强度为: （37）式中是光轴上 点的光强； 是圆孔边缘与中心点在同一 方向上光线间的相位差.4. 光强分布曲线4.1单缝衍射的光强分布曲线图5（a）是利用MATLAB实现的波长，缝宽为0,2mm的衍射图样，得出了光强分布曲线图，图5（b）是缝宽为1mm的衍射图样，图5（c）是缝宽为2mm的衍射
23、图样，从得到的衍射图样来看，光强的分布取决于缝宽和波长不同的缝宽得到的衍射光强的光强曲线图形有很大区别，利用计算机模拟衍射光强的分布快速便捷，和理论的误差相对也较小.图5(a)缝宽0.1mm单缝衍射图图5(b)缝宽1mm的单缝衍射图图5(c)缝宽2mm的单缝衍射图模拟光斑能量与该分图光斑能量的相对误差.可以看出, 当使用的1 次FFT 计算时, 计算结果与实验比较相距甚远, 利用 的计算结果与第1次计算相比较, 求得相对误为79.1% ,说明2 次计算的结果都还不可信, 有必要继续增加取样数.的计算结果与的计算结果相比,相对误差为67.8%说明取样数仍然未达要求.继续往后增加区域分割数再计算后
24、发现,当与的计算结果比较时, 相对误差已经在6.5% ,取样间隔已经接近合理.再往下增加取样数后, 跟踪误差已经下降到0.1%以下.将的模拟光斑与实验相比不难看出,理论模拟已经非常接近实际.4.2圆孔衍射的光强分布构造透光圆孔如图5中衍射前图样所示，图6（a）（b）分别是圆孔衍射图样和衍射光强分布图及三维光强分布图，该实验中波长，从而很好地得到了圆孔衍射的光强分布图及其三维图光强分布图.图6（c）、（d），模拟了衍射图样,得到了二维光强分布图和三维光强分布图.图6(e)(f)分别是衍射图形及其二维光强分布图，三维光强分布图，该试验中使用，.衍射条纹为明暗相间的圆环,中央为亮斑.由于点数太多显示
25、为黑色,但仍然可以看出中央亮斑的强对于其余亮纹要大得多.衍射前图样10020030050100150200250300衍射后图样图6(a) r=10的圆孔衍射图样图6(b) r=10圆孔衍射光强分布图样?衍射前图样10020030050100150200250300衍射后图样图6(c) r=15圆孔衍射图样图6(d) r=15圆孔衍射光强分布图样衍射后图样图6(e) r=20的圆孔衍射图样10020030050100150200250300衍射前图样图6(f) r=20 的圆孔衍射光强分布图样首先,从理论上大致估算一下满足夫琅禾费衍射所需要的条件. 第一,由点光源发出的球面光波在衍射圆孔表面上
26、形成的光波场可近似认为是平面波,这要求圆孔半径满足远场条件;第二, 由衍射圆孔表面上各次波源发出的球面光波在场点相遇时,可被认为是由无穷远射来的平面波,这一问题也可以反过来考虑,认为点发出的球面光波在圆孔上形成的波场近似为平面波,由于点不在光轴(轴)上,因此要求圆孔半径R满足远场条件和P点满足傍轴条件,在通常情况下, 后者比较容易得到满足,起关键作用的是前者. 因此,主要应该比较,与之间数值的大小关系.5.结论本文利用快速傅里叶变换计算了光的单缝和圆孔衍射的光强分布，根据计算结果利用MATLAB软件仿真模拟了单缝和圆孔衍射现象.鉴于积分可以表示为傅里叶变换及卷积两种形式，存在一次快速博里叶变换
27、及快速卷积算法快速卷积算法需要作一次快速博里叶变换及一次快速傅里叶反变换运算，计算出衍射问题，利用MATLAB模拟出单缝衍射以及圆孔衍射的衍射图样和光强分布图，较为理想的解决了衍射问题.参考文献1姚启钧.光学教程M.北京高等教育出版社,008:106-109.2陈聪、李定国.基于快速傅里叶变换的衍射现象的数值仿真J.大学物理2004,9:47-623李俊昌、熊秉衡.信息光学理论与计算M.科学出版社,09:93-1104董克剑.利用MATLAB模拟光的衍射现象J.物理教师,2008,29(5),30-32.5王竞争等.基于MATLAB的光的干涉和衍射现象的模拟研究J.延边大学学报,2009,35
28、(4),319-322.6喻力华、赵维义.圆孔衍射光强分布的数值计算J.大学物理2010,1,20(1),17-19附录1.用MATLAB软件模拟单缝衍射和光强分布曲线的程序clcclearlambda= input(输入光的波长nm):(500);lambda= lambda* 500e-9;awidth= input(输入缝宽:(0.2,1,2);awidth= awidth* 0.001;l= input(输入缝到屏的距离(1);xmax=3* lambda* l/awidth;ys= xmax;nx= 51;xs=linspace(-xmax,xmax,nx);npoints= 51;
29、ypoint= linspace(-awidth/2,awidth/2,npoints);for j=1: nx r= sqrt(xs(j)- ypoint).2+ l2); phi= 2* pi.*(r- l)./lambda; sumcos= sum(cos(phi); sumsin= sum(sin(phi); b(j)= (sumsin2+ sumcos2)/npoints2;endclf;figure(gcf);nclevels= 225;br= 2* b* nclevels;subplot(2,1,1),image(ys,xs,br);colormap(gray(nclevels)
30、;subplot(2,1,2),plot(xs,b,*,xs,b);gridaxis(-xmax,xmax,0.0,1.0);colormap(gray(nclevels)figure(gcf);2.用MATLAB软件模拟圆孔衍射和光强分布曲线的程序clcclear allN=300;r=10，15，20;a=1;b=1;l=zeros(N,N);m,n=meshgrid(linspace(-N/2,N/2-1,N);D=(m-a).2+(n-b).2).(1/2);i=find(D=r);l(i)=1;subplot(1,2,1);imagesc(l)colormap(gray(N)axis
31、 imagetitle(衍射前图样)L=300;M=300;x,y=meshgrid(linspace(-L/2,L/2,M);lamda=632.8e-6;k=2*pi/lamda;z=1000000;h=exp(j*k*z)*exp(j*k*(x.2+y.2)/(2*z)/(j*lamda*z);H=fftshift(fft2(h);B=fftshift(fft2(l);G=H.*B;U=fftshift(ifft2(G);Br=(U/max(U);subplot(1,2,2);imshow(abs(U);axis image;title(衍射后图样)colormap(gray);figu
32、re(2)colormap(gray);subplot(1,2,1);mesh(x,y,abs(U);subplot(1,2,2)plot(abs(Br)致 谢大学四年学习时光已经接近尾声，在此我想对我的母校，我的父母、亲人们，我的老师和同学们表达我由衷的谢意。感谢我的家人对我大学四年学习的默默支持；感谢我的母校河西学院给了我在大学四年深造的机会，让我能继续学习和提高；感谢河西学院的老师和同学们四年来的关心和鼓励。老师们课堂上的激情洋溢，课堂下的谆谆教诲；同学们在学习中的认真热情，生活上的热心主动，所有这些都让我的四年充满了感动。 这次毕业论文设计我得到了很多老师和同学的帮助，其中我的论文指导老师马老师对我的关心和支持尤为重要。每次遇到难题，我最先做的就是向马老师寻求帮助，而马老师每次不管忙或闲，总会抽空来找我谈，然后一起商量解决的办法。我做毕业设计的每个阶段，无论从选题到查阅资料，论文提纲的确定，中期论文的修改，还是到后期论文格式调整等的各个环节，马老师都给予了我悉心的指导。这几个月以来，马老师不仅在学业上给我以精心指导，同时还在思想给我以无微不至的关怀，在此谨向马老师致以诚挚的谢意和崇高的敬意。感谢在整个毕业设计期间和我密切合作的同学，和曾经在各个方面给予过我帮助的伙伴们，在此，我再一次真诚地向帮助过我的老师和同学表示感谢！
