数值振荡积分_基于奇性和波数的积分区间剖分法.pdfVIP

中国科学 数学 年 第 卷 第 期 自然科学基金项目进展专栏 数值振荡积分 基于奇性和波数的积分区间剖分法 献给林群教授 华诞 马云云 许跃生∗ 中山大学数学与计算科学学院和广东省计算科学重点实验室, 广州 510275 E-mail: mayy007@, yxu06@ 收稿日期: 2014-10-30; 接受日期: 2015-04-09; * 通信作者 国家自然科学基金(批准号:和 和广东省引进计算科学创新科研团队计划 (批准号: 资助项目 摘要 本文针对高频振荡函数的积分提出了复合的 的数值积分公式 当被积函数没有 奇点、振子没有驻点时 复合的数值积分公式是基于对积分区间按振子导数值和波数剖分的原则设计 利用 型数值积分公式计算子区间上的振荡积分 当被积函数有奇点、振子存在驻 点时 根据函数的奇性和波数对积分区间剖分 使其于每个子区间上为非剧烈振荡函数的奇异积分 或 光滑函数没有驻点的振荡积分 然后 采用经典的数值奇异积分公式计算非剧烈振荡函数的奇异积分 采用修改的 型数值积分公式计算振荡积分 与计算振荡积分的已知方法相比 本文 发展的复合 的数值积分方法 既不需要计算振子的反函数 也不必借助特殊函数的积分 关键词 振荡积分 代数奇性 驻点 型方法 区间分划 主题分类 引言 高频振荡函数的积分在雷达、声呐和计算机断层扫描等领域有着广泛的应用 随着振荡频率的增 加 被积函数振荡加剧 除在积分区间端点和振子驻点附近 函数值在积分区间的大部分区域相互抵 消 利用传统的数值积分公式计算振荡函数的积分时 为保证计算精度 计算量随振荡频率的增加呈 线性增长 因此 振荡积分的研究 特别是设计高效的、稳定的数值计算方法是当今计算数学研究的一 个重要的且具有挑战性的课题 本文研究形式为 f 的振荡函数在有界区间上的数值积分公式 其中波数 κ ≫ 振幅 f 具有 奇性 振子g 存在驻点 针对此类振荡函数的积分 文献 提出了复合的 的数值积分公 式 首先根据被积函数振荡的频率和奇性对积分区间设计一套剖分方案 再选取恰当的数值积分公式 计算被积函数在每个子区上的积分 剖分方案的设计准则是保证每个子区间上的数值积分公式具有一 致的收敛阶 借鉴文献 中对积分区间剖分的思想 本文设计了一系列新的剖分方法 并基于相应的 剖分方法提出了两类复合的 的数值积分公式 对于计算复杂振子的振荡积分 新的数值 积分公式不仅拥有文献 中方法无需计算振子反函数或借助特殊函数积分的优点 而且在不降低计 算精度的前提下 减少了计算量 提高了计算效率 英文引用格式 马云云等 数值振荡积分 基于奇性和波数的积分区间剖分法 当被积函数充分光滑、振子在积分区间内存在使得所有小于正整数r 阶的导数都为零的点时 随 着波数 κ → ∞ 振荡积分的衰减速度是 Oκ− 参见文献 当用 来近似积分值时 精度是 Oκ− 所以 当一个数值积分公式的精度低于Oκ− 时 视为无效的方法 近年来 振荡积分 数值计算的研究取得了若干可喜的成果 各种高效数值积分公式的共同特点是 它们的计算精度随着 被积函数振荡频率的增大而提高 目前 文献中计算振荡积分的数值方法主要有四种 它们分别是 渐 进展开法、 型方法、 型方法和数值最速下降法 渐进展开法是指利用波数的负指数形成的 级数逼近积分值 为研究振荡积分的特点和分析数值积分公式的收敛阶提供理论基础 推导振荡积分 渐进展开的数学工具有分部积分公式和 逼近等 参见文献 渐进展开方法采 用振荡积分渐进展开前有限项部分和近似积分值 对于固定的波数 渐进展开法的收敛阶较低