- 1、本文档共12页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
异型孔的夫琅禾费衍射图样的计算机模拟
异型孔的夫琅禾费衍射图样的计算机模拟 第 PAGE 12 页 共 NUMPAGES 12 页
异型孔的夫琅禾费衍射图样的计算机模拟
叶小军(安庆师范学院物理与工程学院 安徽 安庆 246011)
指导教师:张杰
摘要: 衍射是光的波动性重要特征之一,是光在传播过程中遇到障碍物时偏离原来的传播方向,在观察屏上出现光强的不均匀分布,当光源与观察屏离衍射屏无限远的这一类衍射称为夫琅和费衍射。本文用分步积分法推导出夫琅禾费矩孔衍射的光强分布公式,对光强分布进行讨论,并且使用计算机模拟其衍射图样。然后从基尔霍夫衍射积分公式出发,利用曲线坐标关系和椭圆的傅里叶变换,导出了椭圆的夫琅禾费衍射精确解。继而由已知的矩形夫琅禾费衍射,求得了一类椭圆与矩形围成的几何图形的夫琅禾费衍射表达式。作出了几种该类图形的夫琅禾费衍射图样,并进行了分析,接着从衍射积分公式出发,利用贝塞尔函数,求得夫琅禾费环孔衍射积分的精确解,模拟了它的衍射图样,并对衍射图样进行了分析。
关键词:夫琅禾费衍射,异形孔衍射,光强,bessel函数,衍射积分,衍射图样
1.引言
许多的基础光学教材在讨论夫琅禾费衍射时,一般都只是对单缝、圆孔等一些简单规则孔的夫琅禾费衍射进行讨论,而对于象矩孔、椭圆或者由这些简单形状的孔所组成的一些复杂异形孔的夫琅禾费衍射则很少提及。其实,夫琅禾费单缝衍射只不过是夫琅禾费矩孔衍射的一种特例。本文将根据惠更斯-菲涅耳原理,利用分步积分法推导出夫琅禾费矩孔衍射的光强分布公式,而夫琅禾费单缝衍射仅是矩孔衍射的一种特例;然后巧妙地利用曲线坐标关系,从菲涅耳-基尔霍夫衍射积分公式出发,导出了椭圆孔径的夫琅禾费衍射精确解,继而求出了一类椭圆与矩形围成的几何图形的夫琅禾费衍射精确解;最后由衍射积分公式出发,利用贝塞尔函数得出了环形孔夫琅禾费衍射精确解,而夫琅禾费圆形孔的衍射可以看作环形孔衍射的一种特例。同时我们还将利用计算机对描绘出上述各种衍射的衍射光强分布曲线以及模拟其衍射图样。
2.衍射光强的计算及衍射图样的分析
2.1夫琅禾费矩孔衍射的光强
设用于夫琅禾费矩孔衍射实验的光源为单色点光源。实验时,把点光源置于一会聚透镜的物方焦点上,从点光源发出的光经过会聚透镜折射后变成平行光。让平行光垂直入射到矩孔上,在矩孔后置一会聚透镜L2,在L2的象方焦平面上放置观察屏E,则在屏上会观察到夫琅禾费矩孔衍射图样。实验装置如图1所示,设矩孔沿x、y轴方向的边长分别为a和b,屏上任意一点P的衍射方向分别用与x方向对应的衍射角θ1和与y方向对应的衍射角θ2表示。在矩孔上任取一平行于y轴的带面元,该带面元的面积为bdx,设该带面元上所有次波在P点的合振动为dE,则把dE对整个矩孔进行积分即可得出矩孔上所有次波在P点的合振动E。
。
图1
为了求出,在带面元上任取一小面元,小面元的面积为,设该面元上的次波在P点的振动为,则对进行积分即可得出带面元上所有次波在P点的合振动。为了便于计算,设光波在矩孔面上的初位相为零,小面元到P点的光程为r。根据惠更斯 菲涅耳原理,该面元上的次波在P点的振动为:
(1)
式中,C为常数,在近轴条件下,可以视倾斜因子K(θ1,θ2)为常数,并把振幅分母中的r也视为常数。令C’=CK(θ1,θ2)/r,则(1)式简化为:
(2)
设带面元的中点(x,0)到P点的光程为r0’,则根据几何关系可知,r=r0’+ysinθ2,
(3)
对上式进行积分便可得出带面元上所有次波在P点的合振动,即:
(4)
设矩孔的中心O点到P点的光程为r0,根据几何关系可知, ,把其代入(4)式得:
(5)
对上式进行积分便可得出矩孔上所有次波在P点的合振动E,即:
(6)
式中,、。由(6)式可知,P点的光强为:
(7)
式中,I0=(Cab)2。由于P为观察屏上的任意一点,因此,(7)式即为夫琅禾费矩孔衍射的光强分布公式。由(7)可使用MATLAP软件,作出不同矩孔衍射图样如下:
图2(a=b=0.05) 图3(a=0.03,b=0.08)
a=0.05;
b=0.05;
f=600;
lambda=632;
文档评论(0)