周期函数的傅里叶变换 周.docx

周期函数的傅里叶变换与周期角谱与衍射角谱的相似性光电1304班 吴妮珊 U201314111一 要求1836 年泰伯(Talbot)发现了周期性物体的自成像效应,即泰伯效应；当用平面波照明一个具有周期性的透过率函数的物体时,会发现在该透射片后的某些距离上出现该周期函数的像。周期函数不满足狄里赫利条件，如果将周期函数进行级数展开，则级数的系数满足广义傅里叶变换，所以对周期函数傅里叶变换，先将周期函数级数展开，再做广义傅里叶变换。角谱分析的本征函数就是余弦函数；当周期函数展开为傅里叶级数，则角谱为系列脉冲；乘传递函数后即为菲涅尔衍射角谱。某些位置z 处的传递函数为1，则衍射角谱与周期函数角谱相似。通过下面几个问题的讨论，加深对泰伯效应的理解。(1) 建立二元振幅光栅透过率函数为：的菲涅尔衍射光场数学模型。（2）对自成像的可能性及条件进行讨论,给出在用单色平面波照明时其自成像的位置。（3）采用matlab 绘出自成像位置, , 光强分布。二 数学模型根据题设条件，透过率函数为。对该函数进行计算可得对u(x)进行傅立叶级数展开，根据公式可以求出对u(x)求傅里叶变换，可以得到单位振幅单色平面波垂直照射物体，物体后表面的场的角谱为物体后的衍射是菲涅尔衍射，是线性空间不变系统，其传递函数为离物体z的平面上场的角谱为物体角谱与传递函数之积，即根据脉冲函数与普通函数的乘积性，有当,即时，，变成对应的场即为上式的傅里叶变换相应地光强为根据自成像效应出现的条件，可以得到在用单色平面波照明时其自成像的位置三 仿真1 z=时2 z=/2时3 z=/4时四 分析结果根据上述三张图，尤其是在z=、/2、/4时，几张图中的光强分布完全相同，显然并不合理，与作业要求中的泰伯效应描述出的现象并不符合。上述matlab仿真结果说明，在数学建模和推导的过程中，或者是在matlab的仿真程序的编写中，一定在某处存在错误，导致最后结果的不正确。鉴于三张图的光强分布图完全相同，我认为错误应该是存在于的计算结果上，即实际上正确的值并非在上面数学建模过程中求出的值。基于上述的的计算结果是根据书上相关章节的内容求出的，所以目前我并不知道正确的求取过程应该是什么样的。希望老师上课的时候能够讲解一下。五 总结与感悟通过这次电子作业matlab仿真得到的错误结果，说明我在物理光学的傅里叶光学的部分的学习中，在知识的把握和理解上存在一定的疏漏和偏差，在平时这方面的学习上也存在着一些不足和欠缺。这次的电子作业给我在物理光学的学习上敲响了警钟。在今后的学习过程中，我一定会更加认真，更加重视学习中的各种细节，不放过存在的任何问题。附 matlab源代码lamda=550*10^-9;d=0.1;zt=2*d^2/lamda;k=2*pi/lamda;j=sqrt(-1);z=zt;x=0:0.01:1;u1=rectpuls((x-d),d/2)*exp(-j*2*pi);u2=rectpuls((x-2*d),d/2)*exp(-j*2*pi);u3=rectpuls((x-3*d),d/2)*exp(-j*2*pi);u4=rectpuls((x-4*d),d/2)*exp(-j*2*pi);u5=rectpuls((x-5*d),d/2)*exp(-j*2*pi);u6=rectpuls((x-6*d),d/2)*exp(-j*2*pi);u7=rectpuls((x-7*d),d/2)*exp(-j*2*pi);u8=rectpuls((x-8*d),d/2)*exp(-j*2*pi);u9=rectpuls((x-9*d),d/2)*exp(-j*2*pi);u10=rectpuls((x-10*d),d/2)*exp(-j*2*pi);u11=rectpuls((x+d),d/2)*exp(-j*2*pi);u12=rectpuls((x+2*d),d/2)*exp(-j*2*pi);u13=rectpuls((x+3*d),d/2)*exp(-j*2*pi);u14=rectpuls((x+4*d),d/2)*exp(-j*2*pi);u15=rectpuls((x+5*d),d/2)*exp(-j*2*pi);u16=rectpuls((x+6*d),d/2)*exp(-j*2*pi);u17=rectpuls((x+7*d),d/2)*exp(-j*2*pi);u18=rectpuls((x+8*d),d/2)*exp(-j*2*pi);u19=rectpuls((x+9*d),d/2)*exp(-j*2*pi);u20=rectpuls((x+10*d),d/2)*exp(-j*2*pi);U=u1+u2+u3+u4+u5+u6+u7+u8+u9+u10+u11+u12