[论文解读] The application of the exact operational matrices for solving the Emden-Fowler equations, arising in astrophysics

Kourosh Parand, Hossayni, Sayyed A.+1 more|arXiv (Cornell University)|Oct 7, 2013

Fractional Differential Equations Solutions参考文献 72被引用 8

一句话总结

本论文提出基于切比雪夫多项式的精确运算矩阵（EOMs），用于求解天体物理学中出现的高非线性艾米达-福勒方程，相较于传统的普通运算矩阵（OOMs）展现出更优的精度。通过精确的系数逼近实现近乎精确的残差计算，该方法实现了快速的误差收敛，并在多种具有不同非线性的测试案例中，于残差和误差范数方面均优于OOMs。

ABSTRACT

The objective of this paper is to apply the well-known exact operational matrices (EOMs) idea for solving the Emden-Fowler equations, illustrating the superiority of EOMs versus ordinary operational matrices (OOMs). Up to now, a few studies have been conducted on EOMs and the differential equations solved by them do not have high-degree nonlinearity and the reported results are not regarded as appropriate criteria for the excellence of the new method. So, we chose Emden-Fowler type differential equations and solved them by this method. To confirm the accuracy of the new method and to show the preeminence of EOMs versus OOMs, the norm1 of the residual and error function of both methods are evaluated for multiple $m$ values, where $m$ is the degree of the Bernstein polynomials. We reported the results in form of plots to illustrate the error convergence of both methods to zero and also to show the primacy of the new method versus OOMs. The obtained results have demonstrated the increased accuracy of the new method.

研究动机与目标

• 为解决普通运算矩阵（OOMs）在求解高阶非线性微分方程时的局限性，特别是针对具有复杂非线性的方程。
• 将近期提出的精确运算矩阵（EOMs）应用于艾米达-福勒方程，该类方程以强非线性和天体物理重要性著称。
• 通过在OOMs先前未能获得准确结果的方程上测试EOMs，提供EOMs的全面基准评估。
• 使用高精度数值解（通过Maple的dverk78）验证该方法，并在多个多项式阶次下比较误差和残差范数。
• 通过定量比较不同非线性项下的范数1误差和残差收敛性，展示EOMs的优越性。

提出的方法

• 该方法采用次数为 $ m \cdot N $ 的切比雪夫多项式，通过级数运算矩阵表示非线性项 $ g(y(x)) $，实现对系数的精确逼近。
• 推导出切比雪夫基函数上积分、微分和乘积运算的精确运算矩阵（EOMs），确保矩阵运算的数学精确性。
• 应用伽辽金法结合EOMs，通过强制残差与基函数正交，将微分方程转化为代数方程组。
• 采用域截断与变量映射技术，将艾米达-福勒方程的无限定义域 $[0, \infty)$ 转换为 $[0, M]$，以适配定义在 $[0,1]$ 上的切比雪夫多项式。
• 通过计算 $ L^1 $-范数并绘制随 $ m $ 值增加的残差和误差范数，评估收敛行为。
• 使用高阶连续Runge-Kutta方法（Maple的dverk78）作为参考解，以计算真实误差并验证结果。

实验结果

研究问题

• RQ1 与普通运算矩阵（OOMs）相比，精确运算矩阵（EOMs）是否能显著提高求解高非线性艾米达-福勒方程的精度？
• RQ2 当使用EOMs与OOMs时，随着多项式阶次 $ m $ 的增加，误差和残差范数如何收敛？
• RQ3 EOMs在处理如 $ y^{3/2} $、$ y^{5/2} $、$ \sinh(y) $、$ \sin(y) $ 和 $ e^y $ 等强非线性项时，其能力如何？
• RQ4 使用EOMs是否能生成在整个定义域上几乎为零的残差函数，从而表明解的高保真度？
• RQ5 在无已知解析解的方程中，是否可仅通过数值基准可靠地应用EOM方法进行验证？

主要发现

• 在所有测试的非线性情形下，基于EOM的方法均显著低于OOMs的 $ L^1 $-范数误差和残差，且误差收敛至零的速度更快。
• 对于 $ f(x) = 1 $、$ g(x) = y^{3/2}(x) $ 且 $ N = 12 $ 的情形，EOM方法在所有 $ m $ 值下均表现出显著降低的误差和残差范数。
• 对于 $ f(x) = 1 $、$ g(x) = \sin(y(x)) $ 且 $ N = 12 $ 的情形，范数1误差与残差图显示，EOMs产生的残差更接近零，且误差收敛更快。
• 该方法实现了近乎精确的残差函数，而这是OOMs由于矩阵运算中的近似误差而无法实现的。
• 使用高阶Runge-Kutta方法作为参考解，证实EOM结果始终比OOM结果更接近基准解。
• 级数运算矩阵使得 $ g(y(x)) $ 能够被准确表示为高阶切比雪夫多项式的线性组合，这对处理强非线性至关重要。