版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、工程光学综合练习-圆孔、矩孔的菲涅尔衍射模拟圆孔和矩孔的菲涅尔衍射模拟一、原理由惠更斯-菲涅尔原理可知接收屏上的p点的复振幅可以表示为ep=ceqexpikrrkd (1)其中eq为衍射屏上的复振幅分布,k()为倾斜因子。根据基尔霍夫对此公式的完善,有c=1ik=1+cos21设衍射屏上点的坐标为(x1, y1),接收屏上点的坐标为(x, y),衍射屏与接收屏间距离为z1,当满足菲涅尔近似条件时,即k8z13(x-x1)2+(y-y1)2max2 (2)此时可得到菲涅尔衍射的计算公式ex,y=expikz1iz1-ex1,y1expik2z1x-x12+y-y12dx1dy1 (3)把上式指数
2、项中的二次项展开,并改写成傅里叶变换的形式,可以写成ex,y=exp(ikz1)iz1expik2z1x2+y2fex1,y1expik2z1x12+y12 (4)上式为菲涅尔衍射的傅里叶变换表达式,它表明除了积分号前面的一个与x1、y1无关的振幅和相位因子外,菲涅尔衍射的复振幅分布是孔径平面的复振幅分布和一个二次相位因子乘积的傅里叶变换。相对于夫琅和费衍射而言,菲涅尔衍射的观察屏距衍射屏不太远。在菲涅尔衍射中,输入变量和输出变量分别为衍射孔径平面的光场分布和观察平面的光场以及光强分布,考虑到这三个量都是二维分布,而且matlab主要应用于矩阵数值运算,所以本程序选择用二维矩阵来存储衍射孔径平
3、面和观察平面的场分布,并分别以矩阵的列数和行数来对应平面的直角坐标值(x, y)以及(x1, y1)。二、圆孔菲涅尔衍射用matlab分别构造表示衍射屏和接收屏的二维矩阵。注意使两矩阵阶次相同,考虑到运算量的要求,采样点数不能过多,所以每个屏的x和y方向各取200到300点进行运算。根据式(4),选取合适的衍射屏和接收屏尺寸和相距的距离,模拟结果如下:取典型的he-ne激光器波长=632.8nm,固定衍射屏和接收屏尺寸和相距的距离,分别取不同的圆孔半径,得到以下三组衍射图样,其圆孔半径分别为12mm,20mm,50mm图 1(r=12mm)图 2(r=20mm)图 3(r=50mm)三、矩孔的
4、菲涅尔衍射步骤与上述相同,仅需改变与衍射屏形状对应的矩阵。这里选择矩孔的长宽相等,分别为15mm,20mm,30mm,其衍射图样及强度分布如图4、5、6图 4(a=b=15mm)图 5(a=b=20mm)图 6(a=b=30mm)四、matlab程序%所有长度单位为毫米lamda=632.8e-6;若为圆孔,方框内替换为以下程序r=12;a=1;b=1;i=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);i(i)=1; %孔半径范围内透射系数为1k=2*pi/lamd
5、a;z=1000000;%先确定衍射屏n=300; %圆屏采样点数a=15;b=15;m,n=meshgrid(linspace(-n/2,n/2-1,n);i=rect(m/(2*a).*rect(n/(2*b);q=exp(j*k*(m.2+n.2)/2/z);subplot(2,2,1); %圆孔图像画在2行2列的第一个位置imagesc(i) %画衍射屏的形状colormap(0 0 0; 1 1 1) %颜色以黑白区分axis imagetitle('衍射屏形状')l=300;m=300; %取相同点数用于矩阵运算x,y=meshgrid(linspace(-l/2,l/2,m);h=exp(j*k*z)*exp(j*k*(x.2+y.2)/(2*z)/(j*lamda*z);%接收屏h =fftshift(fft2(h);b=fftshift(fft2(i); %圆孔频谱g=h.*b; %公式中为卷积,空间域中相卷相当于频域中相乘u= fftshift(ifft2(g); %求逆变换,得到复振幅分布矩阵br=(u/max(u); %归一化subplot(2,2,2);imshow(abs(u);axis image;colormap(hot)% figure,imshow(c);
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- -学高一历史必修一教案 人民版
- 全国青岛版信息技术七年级上册专题三第1课三、《搜索引擎》教学设计
- 压力管道系统吹扫清洗工艺工艺规程
- 一年级上册音乐教案 第一单元 快乐的音乐娃苏少版
- 通信工程师:现代电信技术测试题(题库版)
- 通信工程师:同步数字体系(SDH)技术题库知识点三
- 市医联体建设动态监测指标 (试行)
- 七年级生物下册 第四单元 第10章 第2节 人体细胞获得氧气的过程教案 (新版)北师大版
- 2023九年级物理下册 第二十一章 信息的传递第2节 电磁波的海洋教案 (新版)新人教版
- 九年级化学上册 第三章 维持生命之气-氧气 3.1 氧气的性质和用途教案 (新版)粤教版
- 江苏省宿迁市泗阳县2023-2024学年下学期八年级期末学业水平监测数学试卷
- DL-T5024-2020电力工程地基处理技术规程
- 2023年危险废物分类目录
- 2023-2024学年八年级下学期期末考试数学试卷含答案
- 单位食堂厨房设计方案
- DB32T4148-2021燃煤电厂大气污染物排放标准
- 高温合金行业市场前景分析及发展趋势与投资风险研究报告2024-2034版
- 2024年中级轧钢工技能理论考试题库(含答案)
- 银行保洁服务投标方案技术标
- 2024年江苏省网格员招聘理论考试复习题库(含答案)
- (高清版)TDT 1013-2013 土地整治项目验收规程
评论
0/150
提交评论