联用高斯求积与连分式求和计算Hankel变换及在地球物理上的应用.pdfVIP

第 卷 第 期 37 　 1 Vol.37　No.1 年 月 COMPUTING TECHNIQUES FOR GEOPHYSICAL AND GEOCHEMICAL EXPLORATION 2015 1 Jan.2015 　 　 : ( ) 文章编号 1001-1749 2015 01-0001-09 联用高斯求积与连分式求和计算 Hankel 变换 及在地球物理上的应用 , , , , 胡瑞华 林 君 孙彩堂 刘长胜 周逢道 　 ( , ) 吉林大学 仪器科学与电气工程学院 长春 　 　130061 : 、 。 摘 要 光学 电磁学和地震学都需要计算 变换 除少数核函数的 变换有解析 　 Hankel Hankel , 。 , 式 绝大多数的通常采用数值积分的方法对其近似计算 虽然数字滤波法是常规的计算方法 但 , , 。 由于滤波系数的不统一 导致计算结果并不完全一致 甚至在某些情况结果是错误的 为了解决 , 。 : 这个问题 采用直接数值积分的方法计算 Hankel 变换 具体过程是 首先将 Hankel 变换分解 为积分限从零到Bessel 函数的第一个零点值的积分段和后续各相邻零点值作为积分限的积分 , 。 , 段的和 然后用高斯求积计算每段积分 对高斯求积结果组成的序列采用连分式求和 无论是计 。 算结果的精确性还是收敛速度都要优于直接求和 这里首先详细阐述了联用高斯求积与连分式 , 求和计算 Hankel 变换的基本原理和相关算法 然后用 Matlab 编程对比了直接求和与连分式求 , 。 和的计算效果 最后将此方法应用于电偶源在均匀半空间激发的地下场计算的两个实例中 应 : , 用表明 将高斯求积与连分式求和进行联用是计算 Hankel 变换的一种有效方法 可以广泛应用 于地球物理应用中。 : ; ; ; ; 关键词 Hankel 变换 高斯求积 连分式求和 均匀半空间 场计算 中图分类号: 文献标志码: : / O 173.2　　 A　DOI 10.3969 .issn.10