双重贝塞尔函数积分数值计算.pdf

煤田地质与勘探 V o l. 29 N o. 3 ·58 · COAL GEOLO GY EX PLO RA T ION J un. 200 1 ( ) 　　文章编号: 100 11986 200 1 双 重 贝 塞 尔 函 数 积 分 的 数 值 计 算 华　军, 蒋延生, 汪文秉　 ( 西安交通大学, 陕西 西安　7 10049) 摘要: 双重贝塞尔函数积分由于贝塞尔函数互乘项的强振荡与慢衰减特性而难以应用通常的数值 ) ) 积分算法。本文将被积区间[0, ∞ 划分为[0, 0 ]、[ 0 , ∞ 两部分, 应用贝塞尔函数的汉克尔函数表 ) ( ) 述式及后者的大宗量渐近特性, 区间[0 , ∞ 的双重贝塞尔函数积分可被转化为 Fou r ier 正 余 弦 变换, 并可利用各种快速算法对其进行数值计算; 区间[ 0, 0 ]上的双重贝塞尔函数积分的计算可直 接应用一般的数值积分算法并能获得较高的计算精度; 当需大量计算有共同参量的双重贝塞尔函 数积分时, 其计算效率仍显不足。此时, 可应用贝塞尔函数的导数关系式对[ 0, 0 ] 内的双重贝塞尔 函数积分进行恒等变换, 再用差商近似导数, 将其转化为对贝塞尔函数本身的积分, 而该积分又仅 需计算一次。故本算法对双重贝塞尔函数积分的计算效率有明显提高。 关　键　词: 瞬变电磁测深; H ank el 变换; 双重贝塞尔函数积分; 数值积分 中图分类号: P 63 1　　文献标识码: A 本文根据函数 ( ) 的缓变特性、汉克尔函数 F i 1　引言 的渐近特性和贝塞尔函数的导数关系式引入一算 在大回线圈激励下, 若用谱域法进行一维瞬变 法, 计算效率明显提高。 ( ) ( ) 电磁测深 T EM 正演计算时, 会遇到如下 1 式形 2　算法 式的一系列可称为双重贝塞尔函数积分或称为双阶 ( ( ) ( ) 双 变 量 汉 克 尔 变 换 对于 1 式, 若仅考虑某一 , 可忽略其下 th e doub le o rder doub le F i