[论文解读] On the numerical solution of second order differential equations in the high-frequency regime

James Bremer|arXiv (Cornell University)|Sep 21, 2014

Numerical methods for differential equations参考文献 21被引用 4

一句话总结

该论文提出了一种快速、与频率无关的算法，用于求解高频二阶线性常微分方程，通过利用非振荡相位函数实现。该方法通过采用Kummer方程的对数形式重新表述问题，计算避免了昂贵的振荡解析的相位函数，从而实现运行时间与高频参数λ无关，如在贝塞尔函数、勒让德函数和拉长旋转椭球函数上的实验所示，该方法具有高精度和高效率。

ABSTRACT

We describe an algorithm for the numerical solution of second order linear differential equations in the highly-oscillatory regime. It is founded on the recent observation that the solutions of equations of this type can be accurately represented using nonoscillatory phase functions. Unlike standard solvers for ordinary differential equations, the running time of our algorithm is independent of the frequency of oscillation of the solutions. We illustrate the performance of the method with several numerical experiments.

研究动机与目标

• 为解决标准ODE求解器在高频区域效率低下的问题，其中解的振荡导致与λ成比例的计算成本过高。
• 开发一种适用于广泛二阶ODE类别的通用算法，包括定义特殊函数的方程。
• 通过使用源自Kummer方程对数形式的非振荡相位函数，实现与λ无关的解精度。
• 实现对贝塞尔函数、勒让德函数和拉长旋转椭球函数等特殊函数的快速评估，达到高精度且计算开销极小。

提出的方法

• 该方法构造一个非振荡相位函数α(t)，使得ODE y'' + λ²q(t)y = 0 的解可通过 u(t) = cos(α(t))/|α'(t)|¹ᐟ² 和 v(t) = sin(α(t))/|α'(t)|¹ᐟ² 表示。
• 通过令 r'(t) = λ exp(r(t)/2)，将Kummer的三阶非线性ODE转化为更易处理的对数形式，得到方程 r'' - ¼(r')² + 4λ²(exp(r) - q) = 0。
• 该算法使用在奇点附近采用非均匀网格、子区间上采用高阶切比雪夫谱配点法的数值方法求解该对数形式。
• 解在每个子区间上使用高阶切比雪夫展开表示，从而实现在任意点的高效评估。
• 该方法通过直接处理相位函数而非解析快速振荡，避免了对振荡的显式解析，确保计算成本与λ无关。
• 该方法在贝塞尔函数、勒让德函数和拉长旋转椭球函数方程上进行了验证，与现有求解器如Osipov-Rokhlin进行了比较。

实验结果

研究问题

• RQ1 能否为高频二阶ODE开发一种通用的数值方法，使其与振荡频率λ无关？
• RQ2 对于包括定义特殊函数的方程在内的广泛ODE类，是否存在非振荡相位函数，且能否高效计算？
• RQ3 Kummer方程的对数形式能否以高精度进行数值求解，从而实现解的快速评估？
• RQ4 与现有求解器如Osipov-Rokhlin相比，该方法在拉长旋转椭球函数上的性能和精度如何？
• RQ5 该方法能否实现对高阶勒让德函数等特殊函数的机器精度级精度？

主要发现

• 该算法实现了约 (μλ)⁻¹exp(−μλ) 量级的解精度，其中μ依赖于q但不依赖于λ，确保在大λ值下仍保持高精度。
• 该算法的运行时间与λ无关，而标准求解器的计算成本随λ线性增长。
• 对于阶数 ν = π×10⁵ 的勒让德函数，该方法使用16阶切比雪夫展开在非均匀网格上计算，最大绝对误差低于10⁻¹⁴。
• 在拉长旋转椭球函数的评估中，该方法在速度上优于Osipov-Rokhlin算法，同时保持了相当的精度。
• 该方法准确计算了贝塞尔方程的非振荡相位函数，尽管不存在初等闭式表达式，证明了该方法的通用性。
• 数值实验表明，该方法在包括端点处具有奇异系数的各类ODE中均保持了高精度和高效率。