非对等抽样的衍射计算研究.docVIP

  非对等抽样的衍射计算研究# 潘力平，袁操今，李重光** （昆明理工大学理学院，昆明 650500） 5 摘要：本文提出两种新的算法来计算光波衍射非对等抽样点数问题。通过对观察屏抽样点的 分块处理运算，可以加快计算机处理速度，并且突破了软件和内存对整个抽样点数的限制。 模拟实验得到的结果证明了理论分析的合理性。 关键词：衍射计算；非对等抽样；菲涅尔衍射；傅里叶变换 中图分类号：O436 10 Diffraction Calculation with Different Samplings PAN Liping, YUAN Caojin, LI Chongguang (Science School, Kunming University of Science and Technology, KunMing 650500) Abstract: In this paper, two different methods to deal with the diffraction calculation with 15 different samplings are presented, when number of sampling points between the object plane and the observe plane are different. By dividing the observe plane into small sections, not only the speed can be improved but also the limitation caused by the software and the memory can be free. Simulation on computer is performed to validate the proposed method theoretically. Keywords: diffraction calculation; different samplings; Fresnel diffraction; Fourier Transform 20 0 引言 随着计算机和高分辨率电荷耦合器件的发展，数字全息、数字干涉测量等技术得到了飞 速的发展，这些技术经常涉及到光衍射的数值计算和模拟，因此，找到一种准确、快速地模 拟光衍射的算法具有重要的意义。目前，常用的光衍射算法有标量衍射算法、传递函数算法、 25 30 S-FFT 算法等[1]。为了加快运算速度，前面提及的算法都用到了快速傅里叶变换。受到快速 傅里叶变换的限制，物平面与接收面的抽样点数必须相同[2,3]。如果接收面是一个 CCD，它 的抽样点数就是它的像素点数，直接采用快速傅里叶变换，就不能很好的模拟 CCD 面上光 强的分布。有人提出补零算法来解决这一问题，但是随着补零点数的增加，会导致计算机的 运算速度降低，进而运算量超出软件和内存的限制[1,2]。 本文提出了两种方法，既解决非对等抽样点数的问题，又可以加快运算速度，同时突破 了软件和内存运算量的限制，基本上能满足衍射计算精度的需求。通过模拟实验，分别对比 两种新方法与补零法的光强分布的差异，分析产生差异的原因。 1 计算原理 假设光的衍射满足菲涅尔衍射条件，由菲涅尔衍射公式知其离散化形式为：  35  U?? m?x, n?y??? exp?? jkz?????? k 2 2 j? z?????? 2z ??????????????????? k  0 0  2  ??? n0?y0??  2  ?  ??? ???  (1)  基金项目：高等学校博士学科点专项科研基金资助课题（20095314120009） 作者简介：潘力平，（1987-），男，硕士研究生，主要研究方向：数字全息技术。 通信联系人：李重光，（1966-），男，教授，主要研究方向：数字全息技术。E-mail: licg66@ -1-   其中 U?? m?x, n?y??  表示观察面复振幅， 0 0 0 0 0 表示物光波面复振幅[4]， m, n, m0 , n0 分别是观察面 x, y 和物光波面 x0 , y0 上的离散点数，??x,??y,??x0 ,??y0 为抽样间 隔。直接采用二维离散傅里叶变换，离散点数必须满足 m?? m0 , n?? n0 。为了解决这一问题， 我们提出了以下两种算法。  40  1.1 交叉分块衍射 假设 m?? am0 , n?? bn0 ，其中 a, b 为整数，以 a 为纵向采样数， b 为横向采样数，