[论文解读] Spectral accuracy for the Hahn polynomials

René Goertz, Philipp Öffner|arXiv (Cornell University)|Sep 23, 2016

Mathematical functions and polynomials参考文献 17被引用 4

一句话总结

本文在双曲守恒律的谱方法背景下，建立了Hahn多项式谱精度的理论。通过使用离散正交多项式投影，作者证明了Hahn多项式展开对光滑函数可实现谱收敛，同时通过替代节点分布和重构技术解决了龙格现象。

ABSTRACT

We consider in this paper the Hahn polynomials and their application in numerical methods. The Hahn polynomials are classical discrete orthogonal polynomials. We analyse the behaviour of these polynomials in the context of spectral approximation of partial differential equations. We study series expansions $u=\sum_{n=0}^\infty \hat{u}_n ϕ_n$, where the $ϕ_n$ are the Hahn polynomials. We examine the Hahn coefficients and proof spectral accuracy in some sense. We substantiate our results by numericals tests. Furthermore we discuss a problem which arise by using the Hahn polynomials in the approximation of a function $u$, which is linked to the Runge phenomenon. We suggest two approaches to avoid this problem. These will also be the motivation and the outlook of further research in the application of discrete orthogonal polynomials in a spectral method for the numerical solution of hyperbolic conservation laws.

研究动机与目标

• 在双曲守恒律背景下，建立Hahn多项式展开对光滑函数的谱收敛性。
• 解决在等距网格上Hahn多项式插值时出现的龙格现象。
• 提出替代的节点分布与重构策略，以提升逼近稳定性。
• 为在高阶方法（如通量重构FR与通过重构的校正过程CPR）中使用离散正交多项式奠定基础。

提出的方法

• 在谱展开 $ u = \sum_{n=0}^\infty \hat{u}_n \varphi_n $ 中使用Hahn多项式作为基函数，其定义基于超几何函数 $ Q_n(x;\alpha,\beta,N) = {}_3F_2(-n, n+\alpha+\beta+1, -x; \alpha+1, -N; 1) $。
• 通过求和运算应用离散内积，精确计算系数，避免连续正交多项式方法中常见的积分误差。
• 采用基于离散内积的投影方法，计算系数 $ \hat{u}_n(t) = \frac{1}{\|\varphi_n\|^2} \sum_{x \in \Omega} h(x) u(x,t) \varphi_n(x) $，确保系数计算的精确性。
• 研究通过Vandermonde矩阵系统 $ u(x_j,t) = \mathcal{V}(\hat{u}_n(t)) $ 进行插值，重点关注条件数与节点选择问题。
• 提出使用非等距节点（如切比雪夫型）以缓解截断Hahn级数中的龙格现象。
• 探索Gegenbauer重构技术以减少吉布斯型振荡，分析重构阶数 $ m $ 与截断阶数 $ N $ 的比值。

实验结果

研究问题

• RQ1 在双曲守恒律背景下，Hahn多项式能否在逼近光滑函数时实现谱精度？
• RQ2 龙格现象如何影响等距网格上Hahn多项式逼近的效果？是否可被缓解？
• RQ3 在何种函数 $ u $ 或参数 $ \alpha, \beta $ 条件下，可确保截断Hahn级数保持稳定且准确？
• RQ4 在Gegenbauer型重构中，重构阶数 $ m $ 与截断阶数 $ N $ 的比值如何影响总逼近误差？
• RQ5 离散正交多项式能否扩展至三角形网格，以应用于非结构化高阶方法？

主要发现

• Hahn多项式展开对解析函数可实现谱收敛，系数衰减快于任意 $ n $ 的幂次，证实了谱精度。
• 数值实验表明，当 $ m=10 $ 时，截断Hahn级数表现出显著的点态误差与振荡，尤其在 $ \alpha=\beta=0 $ 时，表明存在龙格现象。
• 参数 $ \alpha, \beta $ 的选择显著影响逼近质量，测试案例中 $ (5,0) $ 的收敛性优于 $ (0,0) $ 或 $ (0.5,0.5) $。
• 在等距网格上插值导致Vandermonde矩阵病态且不稳定，证实了龙格现象的影响。
• 采用切比雪夫型节点或其它节点分布可提升稳定性并降低点态误差，优于等距网格。
• 当控制 $ m/N $ 比值时，Gegenbauer重构方法可减少振荡并改善收敛性，但若 $ m $ 与 $ N $ 增长过快，收敛可能失败。