用傅里叶变换计算衍射的光强分布

龙岩学院学年论文（设计）论文题目用傅里叶变换计算衍射的光强分布学院物理与机电工程学院专业物理学（光电子技术方向）年级2011级姓名徐武童学号2011042526指导教师兑自强二0一三年四月十二日用傅里叶变换计算衍射的光强分布物理与机电工程学院11物本2011042526徐武童指导老师：兑自强【摘要】：利用傅里叶变换式计算光的单缝和圆孔衍射的光强分布，根据计算结果利用MATLAB软件仿真模拟单缝和圆孔衍射及光强分布，分析计算和模拟结果得知衍射图样取决于缝宽或孔径的大小【关键词】：傅里叶变换；单缝；圆孔；衍射；光强分布目录前言11.傅里叶变换式11.1一维变换式21.2二维变换式31.3三维傅里叶变换式32.用傅里叶变换计算衍射的光强分布42.1计算圆孔衍射的光强分布62.2计算单缝衍射的光强分布73.光强分布曲线83.1单缝衍射的光强分布曲线83.2圆孔衍射的光强分布曲线94.讨论104.1单缝衍射104.2圆孔衍射10总结11致谢110前言衍射现象是波动光学中的重要知识，光的衍射的定义从广义上说是光在传播过程中，遇到障碍物时产生的偏离几何光学规律从而引起光强重新分布的现象，也称为绕射。该定义指出光的衍射是一种区别于几何光学规律的光的传播现象。当所选光学元件的尺度与波长相当时，光的传播现象明显不同于几何光学所描述的。它也明确给出了产生衍射现象的条件“光波遇到障碍物”，对于任何一束光都会因在空间传播过程中遇到障碍物而使自由波面受损，从而改变波前后振幅，使光表现出衍射行为。而傅里叶变换是一种特殊的积分变换，它能将满足一定条件的某个函数表示成正弦基函数的线性组合或者积分。在不同的研究领域，傅里叶变换具有多种不同的变体形式，如连续傅里叶变换和离散傅里叶变换。在数学领域，尽管最初傅立叶分析是作为热过程的解析分析的工具，但是其思想方法仍然具有典型的还原论和分析主义的特征。从纯粹的数学意义上看，傅立叶变换是将一个函数转换为一系列周期函数来处理的。从物理效果看，傅立叶变换是将图像从空间域转换到频率域，其逆变换是将图像从频率域转换到空间域。在现代光学发展的今天，如何运用傅里叶方法解决干涉、衍射和成像等问题成了至关重要的部分。1傅里叶变换式1.1一维变换式某个空间变量的一维函数xf，可以表示为无穷多个谐波分量的线性组合：00sincos1)(kxdkkBkxdkkAxf（1）其中决定各空间频率k的贡献的权重因子kA和kB分别是xf的傅里叶余弦和正弦变换式，由下式给出：xdxkxfkAcos,xdxkxfkBsin(2)将（2）式代入（1）式中dkxdxkxfkxdkxdxkxfkxxfsinsin1coscos10000（3）由于xkkxxkkxxxksinsincoscoscos，上式可改写为0cos1dkxdxxkxfxf（4）方括号中的量是k的偶函数，因而改变外面的一个积分限得到dkxxkxfxfcos21)(（5）由于0sin2dkxdxxkxfi，把它与上式相加，并应用欧拉公式，得dkexdexfxfikxxik21（6）由于可以写为dkekFxfikx21（7）只要dxexfkFikx（8）式（8）中已令xx，函数kF叫做xf的傅里叶变换式，用下面的记号来表示：}{xfFkF(9)常常把xf和kF称为傅里叶变换式。如果f是时间的函数而不是空间变量的函数，为了在时域中得到相应的变换式偶，我们只是把x换成t，再把空间角频率k换成时间角频率w，即dwekFtfiwt21（10）以及dtetfwFiwt（11）1.2二维变换式光学中一般涉及的信号，例如孔上的光场或者象平面上的通量密度分布，将傅里叶变换式推广到二维情况有yxykxkiyxdkdkekkFyxfyx,21,2（12）或者dxdyeyxfyxFykxkiyx,,（13）其中xk和yk分别是沿坐标轴方向的空间角频率。1.3三维傅里叶变换式将傅里叶变换推广到三维情况有zyxzkykxkizyxdkdkdkekkkFzyxfzyx)3,,21,,（14）以及dxdydzezyxfzyxFykykxkizyx,,,,（15）其中xk，yk和zk分别是沿坐标方向的空间角频率。2用傅里叶变换计算衍射的光强分布光是一种电磁波，按jwte的规律随时间传播，电光源发出的是一组球面波，设光源位于坐标原点处，以速度v在电容率为的介质中传播，当光到达半径为r的求面时，光的场强E是tr,的函数，可以表示为krwtjrErtjwrEtrEexpexp,（16）其中2wk称为波数，trE,为光矢量点光源从原点出发的球面波，能量密度为221Ew（17）以v表示单位时间内光矢量所在空间的体积，则单位时间内通过整个球面的能量为VEW221（18）而0EkrE（19）式中0E是与光源振动有关的常数，k是与介质有关的常数，则krwtjErktrEexp,0（20）为简便，只考虑某时刻的振动，含时间的项jwte可省去。在光学系统中，光从出射光瞳射出，取光瞳坐标为00,yx，观察平面的坐标为yx,，两坐标系相平行，原点在它们的公共垂线上，相距为z。见图（1）图（1）光瞳面上任意一点s00,yx到观察面上的某点pyx,的距离为2122020zyyxxr（21）由（20）式知，光是从s点以球面波jkrerk的形式传播到p点的。如果s点振幅为Eyx,,则在P点光的矢量为jkreyxErkyxE00,,（22）为计算球面上p点的光的场强，需要选取包含s点在内的小面元00dydx，则00001,),(dydxeryxEkyxEjkr（23）为便于计算，设光瞳与观察点面相距很远，取Rr,r的近似值为RyyxxRyxRyyxxyxRzyyxxrx002020210020222122020222（24）在远场衍射的情况下，即R时Ryxk222《1且zR，则dszyyxxjkyxEzkeyxEjkz]exp[,,0000（25）式（25）为远场近似情况下的衍射（也称之为夫琅禾费衍射）的公式。式（25）与式（13）形式完全相同，级光源在p点的Eyx,是光瞳函数E00,yx的二维傅里叶变换式。2.1计算单缝衍射的光强分布当光源为线光源时，式（25）可以演化为000dxexEzkexEzjkxxjkz（26）显然，式（26）和式（11）在形式上完全相同，切有以下对应关系：tfxE0zkxw式（24）可以记为}{0xEFzkexEjkz（27）可见，若要求观察平面上的光强分布，只要把表示出瞳光源的光强分布进行傅里叶变换，其中把傅里叶变换的w用zkx置换即可。也就是说只要计算出出瞳的傅里叶变换，就能求出观察面上的光强分布。为简便，把式（26）前的常数省略，则000dxexExEzjkxx（28）如果狭缝上有均匀照度，其值为A,则长度为a2的狭缝上的振幅为axaxAxE000,0,则wajwAeejwAdxAexEjwajwaaajwxsin00（29）所以，长为a2的狭缝光源在观察面上所形成的光强分布为222222sin4sin4wawaaAwawAxExExI（30）2.2计算圆孔衍射的光强分布圆孔上的夫琅禾费衍射，在光学仪器的研究中具有重要意义。由（25）式知，在远场情况下，一个任意孔在p点所产生的光的场强分布为000000exp,,dydxzyyxxjkyxEzkeyxEjkz（36）对于一个圆孔，由于对称性，所以在孔径平面和观察平面上都采用极坐标。因此，令cos0x，sin0y，cosqx，sinqy（37）因而微分面元现在是dddydx00（38）把（37）式和（38）式代入（36）式中，得ddezkeEazqkjjkz020cos,（39）由于问题是完全对称的，其解一定与无关，这样可以令0求解（39）式。二重积分中与变量有关的部分20cosdezqkj，与零阶的（第一类）贝塞耳函数20cos21dveuJviu，在形式上完全相似，故（39）式可写为dzqkJzkeEajkz002,（40）由贝塞尔函数的递推关系uJuuJududmmm10可知，当1m时uuuJuduJu010(41)其中的u只作为虚拟变量。所以，(39)式中的积分部分可记为azkaqwwdwwJkqzdzqkJ00020则zkaqJkaqzazketEjkz122(42)所以，圆孔衍射的光强分布为212222zkaqzkaqJzAkI（43）3.光强分布曲线3.1单缝衍射的光强分布曲线根据式（30），用MATLAB模拟单缝衍射的光强分布曲线，见图（2）。（图a是单缝衍射图样，图b是与之相应的光强分布曲线，波长都为500nm,缝宽分别是0.1mm﹑0.4mm﹑0.8mm）图（2）3.2圆孔衍射的光强分布曲线根据式（43），用MATLAB模拟圆孔衍射的光强分布曲线，见图（3）。（图c是圆孔衍射图样，图d是与之相对应的光强分布曲线，波长都是500nm,孔径分别为0.1mm﹑-2-1.5-1-0.50图a（1）0.511.52-3-2-1.5-1-0.500.511.52-1.5-1-0.50图b（1）0.511.52x10-300.10.20.30.40.50.60.70.80.91-2（单位：x/m）x10x10-2-1.5-1-0.50图a（2）0.511.52-3-2-1.5-1-0.500.511.52x10-1.5-1-0.50图b（2）0.511.52x1000.10.20.30.40.50.60.70.80.91x10-2（单位：x/m）-3-1.5-1-0.50图a（3）0.511.52x10-3-2-1.5-1-0.500.511.52-2-1.5-1-0.50图b（3）0.511.52x10-300.10.20.30.40.50.60.70.80.91x10-2（单位：x/m）-3-3-30.4mm﹑0.8mm）图（3）-0.8-0.6-0.4-0.20图c（1）0.20.40.60.81x10-1-0.8-0.6-0.4-0.200.20.40.60.81-0.8-0.6-0.4-0.20图d（1）0.20.40.60.81x10-300.10.20.30.40.50.60.70.80.91x10-3-1-3-1（单位：x/m）-1-0.8-0.6-0.4-0.20图c（2）0.20.40.60.81x10-3-1-0.8-0.6-0.4-0.200.20