[论文解读] Generating Laguerre expansion coefficients by solving a one-dimensional transport equation

Andrew V. Terekhov|arXiv (Cornell University)|Sep 18, 2018

Fractional Differential Equations Solutions参考文献 34被引用 5

一句话总结

本文提出了一种新颖的方法，通过求解一维输运方程而非直接计算高度振荡的积分来计算拉盖尔展开系数。该方法通过避免使用求积规则和高精度算术，显著提升了数值稳定性、精度和计算效率，尤其适用于非光滑或大规模数据，在地震反演问题和谱方法中已取得成功应用。

ABSTRACT

Spectral methods based on integral transforms may be efficiently used to solve differential equations in some special cases. This paper considers a different approach in which algorithms are proposed to calculate integral Laguerre transform by solving a one-dimensional transport equation. In contrast to the direct calculation of improper integrals of rapidly oscillating functions, these procedures make it possible to calculate the expansion coefficients of a Laguerre series expansion with better stability, higher accuracy, and less computational burden.

研究动机与目标

• 解决通过直接积分快速振荡函数计算拉盖尔变换系数时存在的数值不稳定性和高计算成本问题。
• 克服使用128位算术计算高阶拉盖尔函数时出现的溢出和下溢问题。
• 开发一种鲁棒、稳定且高效的算法，适用于大规模时间序列数据，特别是在地震勘探和反问题中。
• 通过基于输运方程的谱方法替代高阶求积和昂贵的拉普拉斯/傅里叶变换方法，降低计算负担。
• 在求和步骤中使用标准64位精度的同时，通过高精度函数计算和基于输运方程的系数计算保持精度。

提出的方法

• 用求解一维输运方程的初边值问题来替代拉盖尔系数积分的直接计算。
• 使用谱方法求解输运方程，利用其结构避免对振荡函数进行数值积分。
• 基于求解输运方程实施校正程序，以消除傅里叶类谱解中的人为周期性。
• 通过后验能量基分析拉盖尔谱，进一步优化结果并减少伪影。
• 使用BLAS MKL优化的矩阵-向量乘法实现最终系数序列的高效计算。
• 在子区间上采用局部逼近并引入缓冲区域以处理不连续性，随后利用输运方程框架进行全局重构。

实验结果

研究问题

• RQ1 通过求解一维输运方程是否能够替代直接积分来计算拉盖尔变换系数，并实现更高的稳定性和精度？
• RQ2 与高阶求积或基于拉普拉斯/傅里叶变换的反演等传统方法相比，输运方程方法在计算成本和数值稳定性方面表现如何？
• RQ3 当拉盖尔函数以128位精度计算时，最终求和步骤中在多大程度上可以使用标准64位精度而不影响精度？
• RQ4 子区间分解和缓冲区域对近似误差和计算效率有何影响？
• RQ5 在去除谱解中虚假周期性方面，基于输运方程的校正策略与基于能量谱的策略相比，哪种更可靠？

主要发现

• 基于输运方程的方法成功避免了对高度振荡积分的直接计算，显著减少了因溢出和下溢导致的数值误差。
• 即使拉盖尔函数以128位精度计算，该方法在使用64位精度进行求和时仍保持了数值稳定性和精度。
• 测试结果表明，该算法实现了10⁻³至10⁻⁵量级的近似精度，足以满足地震勘探等实际应用需求。
• 尽管该方法在严格意义上并非快速算法，但通过消除对小网格步长或高阶求积规则的需求，降低了计算成本。
• 基于输运方程的校正程序在可靠性、精度和效率方面均优于基于能量谱的方法，尤其在不扩展近似区间的情况下表现更优。
• 使用优化的BLAS MKL例程进行矩阵乘法可实现快速计算，但第二步的向量化受限，影响整体加速效果。