(完整word版)圆孔矩孔的菲涅尔衍射模拟(matlab实现)-工程光学_第1页
(完整word版)圆孔矩孔的菲涅尔衍射模拟(matlab实现)-工程光学_第2页
(完整word版)圆孔矩孔的菲涅尔衍射模拟(matlab实现)-工程光学_第3页
(完整word版)圆孔矩孔的菲涅尔衍射模拟(matlab实现)-工程光学_第4页
(完整word版)圆孔矩孔的菲涅尔衍射模拟(matlab实现)-工程光学_第5页
已阅读5页,还剩3页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

1、工程光学综合练习-圆孔、矩孔的菲涅尔衍射模拟圆孔和矩孔的菲涅尔衍射模拟由惠更斯菲涅尔原理可知接收屏上的P点的复振幅可以表示为£(P)=c“£(Q)二处K(e)do(1)其中£(Q)为衍射屏上的复振幅分布,K(e)为倾斜因子。根据基尔霍夫对此公式的完善,有1C=iA1+cosOK(0)=-«1设衍射屏上点的坐标为(x1,yI),接收屏上点的坐标为(x,y),衍射屏与接收屏间距离为Z1,当满足菲涅尔近似条件时,即3('xl)2+(7yi)2max«万(2)此时可得到菲涅尔衍射的计算公式户G,y)=exp(ikzjiAzidxrdyr把上式指

2、数项中的二次项展开,并改写成傅里叶变换的形式,可以写成E(x,y)=eX1Zi)=xp/(/+y2)jexp+yf)(4)上式为菲涅尔衍射的傅里叶变换表达式,它表明除了积分号前面的一个与xl、yl无关的振幅和相位因子外,菲涅尔衍射的复振幅分布是孔径平面的复振幅分布和一个二次相位因子乘积的傅里叶变换。相对于夫琅和费衍射而言,菲涅尔衍射的观察屏距衍射屏不太远。在菲涅尔衍射中,输入变量和输出变量分别为衍射孔径平面的光场分布和观察平面的光场以及光强分布,考虑到这三个量都是二维分布,而且Matlab主要应用于矩阵数值运算,所以本程序选择用二维矩阵来存储衍射孔径平面和观察平面的场分布,并分别以矩阵的列数和

3、行数来对应平面的直角坐标值(x,y)以及(xi,yi)o孔菲涅尔衍射用MATLAB分别构造表示衍射屏和接收屏的二维矩阵。注意使两矩阵阶次相同,考虑到运算量的要求,采样点数不能过多,所以每个屏的X和y方向各取200到300点进行运算。根据式(4),选取合适的衍射屏和接收屏尺寸和相距的距离,模拟结果如下:取典型的He-Ne激光器波长A=632.8nm,固定衍射屏和接收屏尺寸和相距的距离,分别取不同的圆孔半径,得到以下三组衍射图样,其圆孔半径分别为12mm,20mm,50mm同孔形状3103回2CD253vn衍封国近1DO2C030口图1(r=12mm)今射回在图2(r=20mm)300250200

4、15010050国孔形状1002DD3CO衍射图标图3(r=50mm)三、矩孔的菲涅尔衍射步骤与上述相同,仅需改变与衍射屏形状对应的矩阵。这里选择矩孔的长宽相等,分别为15mm,20mm,30mm,其衍射图样及强度分布如图4、5、6桁射屏刑状30025020015010050100200300衍射后的图样衍射屏形状30D25020015010D50100200300衍射后的图样-200.200图5(a=b=20mm)衍射所形状100200300衍射后的图样-3J0-200图6(a=b=30mm)四、MATLAB程序%圆孔图像画在2行2列的第一个位置%画衍射屏的形状%颜色以黑白区分%所仃长度单位

5、为毫米lamda=632.8e-6;k=2*pi/lamda;z=1000000;%先确定衍射屏N=300;%圆屏采样点数a=lS;b=15;m,n=meshgrid(linspace(-N/2/N/2-l/N);l=rect(m/(2*a).*rect(n/(2*b);q=exp(j*k*(m.A2+n.A2)/2/z);subplot(2,2,l);imagesc(l)colormap(000;111)axisimagetitle('衍射屏形状')若为圆孔,方框内替换为以下程序r=12;a=l;b=l;I二zeos(N,N);mrn=meshgrid(linspace(-N

6、/2/N/2-l/N);D=(m-a).A2+(n-b).A2).A(l/2);i=find(D<=r);l(i)=l;%孔半径范围内透射系数为1L=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);%接4攵屏%网孔频谱%公式中为卷枳,空间域中相卷相当于频域中相乘%求逆变换,得到更振幅分布矩阵%归一化H=fftshift(fft2(h);B=fftshift(fft2(l);G=H.*B;U=fftshift(ifft2(G);Br=(U/max(U);subplot(2,2,2);imshow(abs(U);%画三维图形axisimage;colormap(hot)%f

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论