matlab球面波衍射,单色点源矩孔菲涅耳衍射光场的计算与模拟

点源照射的衍射较平面光照射的衍射更具一般性,对光束的传输与变换、光学系统的衍射成像以及光信息处理等方面有重要作用.单色点源照射矩孔产生的菲涅耳衍射光场的计算并未见报道,为此对单色点源照射矩孔产生的菲涅耳衍射进行分析,推导它的计算式,并进行计算机模拟试验.图1计算点源P0照射小孔在P点产生的衍射场Fig.1ThecalculationofthediffiractionfieldatpontPbehindasmallapertureirradiatedbyapointsourceP01单色点源照射矩孔衍射的积分表达式1882年基尔霍夫根据麦克斯韦的电磁场理论,从亥姆霍兹方程出发由格林定理导出在基尔霍夫边界条件及其衍射装置给定条件(衍射孔的线度以及点源和考察点到衍射孔的距离远大于波长)下的衍射公式,使建立在假设基础上的惠更斯菲涅耳原理有了较完善的表达式,即菲涅耳基尔霍夫衍射公式[1-4].如图1所示,由基尔霍夫衍射公式知单色点源P0照明平面屏幕在其后P点产生的光场为[1-4]U(P)=i2Aeikrreikss[cos(n,r)-cos(n,s)]d(1)其中,A为常数,为单色点源P0的波长,k=2/,r、s分别为波面上Q点处次波面元d到P0(x0,y0,-z0)和P(x,y,z)的距离,(n,r)和(n,s)分别为矢量r和s跟波面上Q点处法线n的夹角.图2计算单色点源照射矩孔产生的衍射光场的示意图Fig.2Thedrawinforcalculatingthediffractionfieldofa如图2所示,以矩孔中心为坐标原点o,其长和宽分别平行于x轴和y轴,z轴则为矩孔的中垂线.观察屏与孔相距为z1,x1、y1轴分别平行于x、y轴.位于(x0,y0,-z0)处的单色点源P0照射矩孔产生的衍射光场的积分表达式为U(P)=i2A2bb222aaeik(rrs+s)(zr0+zs1)drdy(2)其中r=z02+(x-x0)2+(y-y0)2,s=z21+(x1-x)2+(y1-y)2.2单色点源矩孔菲涅耳衍射的计算式由于甚高频光波的波长很小,相位随光程变化极快,而振幅变化相对较为缓慢.当相位作级数展开k{[(x-x0)2+(y-y0)2]2/8z03+[(x-x1)2+(y-y1)2]2/8z31}max可被舍去时,振幅中可取rz0和sz1,且z0/r+z1/s2cos常量[1-2].则单色点源矩孔菲涅耳衍射的积分表达式可写为[3-4]U(P)iAczo0sz122bb22aaeik[z0+(x-x0)22+z0(y-y0)2]eik[z1+(x-x1)22z+1(y-y1)2]dxdy=Acoseik[z0+z1+x02+y022z0+x12+y122z1iz0z122aaei2k{(z10+z11)[x-(zx00+zx11)/(z10+z11)]2-(xz00+xz11)2/(z10+z11)}dx.2bb2ei2k{(z10+z11)[y-(zy00+zy11)/(z10+z11)]2-(zy00+zy11)2/(z10+z11)}dy(3)令[2,4]u=k(z10+z11)[x-(xz00+xz11)/(z10+z11)],v=k(z10+z11)[y-(yz00+yz11)/(z10+z11)],=2k[(zx00+zx11)2/(z10+z11)],=2k[(zy00+zy11)2/(z10+z11)],则(3)式可改写为[3-4]U(P)=i2(zA1co+sz0)eik(z0+z1+x202+z0y20+x122+z1y21)uu1