基于Dini展开的高阶Hankel变换及其在光束传输中的应用.PDFVIP

物理学报 Acta Phys. Sin. Vol. 62, No. 14 (2013) 140203 基于Dini 展开的高阶Hankel 变换及其在 光束传输中的应用* † 游开明 林燕玲 王友文 陈列尊 戴志平 陆世专 ( 湖南衡阳师范学院物理与电子信息科学系, 衡阳 421002 ) ( 2013年2月2 日收到; 2013年3月30 日收到修改稿) 基于Dini 级数展开, 导出了p 0 阶准离散Hankel 变换(the p th-order quasi-discrete Hankel transform based on Dini series expansion, p DQDHT) 算法, 并给出了该算法在光束传输中的应用. 通过不同输入函数分别对p DQDHT 算 法进行测试及光束通过透镜传输的应用实例, 结果表明: p DQDHT 算法不仅精度高于现有的Hankel 变换算法, 可以 进行多次正、逆变换, 能广泛应用于光束分步传输问题, 而且执行速度也与一般的快速Hankel 变换算法相当. 关键词: Dini 级数展开, 高阶Hankel 变换, 光束传输 PACS: 02.60.Gf, 42.30.Kq, 42.40.My, 42.25.Bs DOI: 10.7498/aps.62.140203 算, 使HT 变换在精度上有所提高. 1998 年Yu 等16 1 引言 对径向坐标按零阶Bessel 函数零点分布的比例进 行离散, 将变换函数做Fourier-Bessel 级数展开, 实 由于Hankel 变换(Hankel transform, HT) 在物 现了准离散HT (quasi-discrete HT, QDHT), 2004 年 理学的很多分支(如光学、声学、电磁学、分子生 ´ 17 Guizar-Sicairos 和Gutierrez-Vega 又把它推广到 物学等) 有着广泛的应用 1−4 , 引起了很多科学工 p 阶HT (the p th-order quasi-discrete HT, p QDHT), 作者的关注, 并对它的数值计算算法进行了研究. 这类方法利用变换矩阵使计算变得简单, 因而速 1977 年, Siegman5 首先对径向坐标抽样点用指数 度也较快, 由于能量守恒, 逆变换精度很高, 不依 形式量化, 将HT 化为离散形式, 用快速傅里叶变换 赖输入函数, 能适应反复变换, 具有很高的实用价 (FFT) 进行计算, 得出了准快速Hankel 变换算法. 此 值. 遗憾的是基于Fourier-Bessel 级数展开HT 方法 后, 出现了很多用于物理学各个方面的HT 算法, 如 处于条件收敛, 不能计算空间域和空间频率域端点 文献[6—18]. 文献[9—11] 把HT 看作是一个二变 处的变换值, 同时其0 阶HT 也不能计算对称轴上 量函数的简单积分的傅氏变换, 先用一般方法计算 点的变换值, 不能直接用于研究传输问题中对称轴 简单积分, 再利用FFT 得到HT. 2003 年, Markham 上的特性. 文献[18] 给出了一种基于Dini 级数展 和Conchello12 总结了计算HT 的6 种方法, 并把 开的0 阶准离散HT 算法, 由于Dini 级数展开比 它们用于计