[論文レビュー] 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

ひとこと要約

本稿では、天体物理学に現れる高次非線形Emden-Fowler方程式を解くために、ベルンシュタイン多項式に基づく正確な作用行列（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）を、強力な非線形性と天体物理学的意義を有するEmeden-Fowler方程式に適用すること。
• OOMsが以前に正確な結果を出せなかった方程式に対して、EOMsの包括的ベンチマークを提供すること。
• 高精度数値解（Mapleのdverk78を用いて）を用いて手法を検証し、複数の多項式次数における誤差および残差ノルムを比較すること。
• さまざまな非線形項におけるノルム1誤差および残差収束の定量的比較を通じて、EOMsの優位性を示すこと。

提案手法

• 非線形項 $ g(y(x)) $ を、系列作用行列を介して $ m \cdot N $ 次のベルンシュタイン多項式で表現することで、係数の正確な近似が可能となる。
• ベルンシュタイン基底関数における積分、微分、積の作用に対して、正確な作用行列（EOMs）が導出され、行列演算における数学的正確性が保証される。
• ガラーキン法をEOMsを用いて適用し、残差関数を基底関数に対して直交させる条件を満たすことで、微分方程式を代数方程式の連立に変換する。
• 定義域の切り詰めと変数変換技術を用いて、Emeden-Fowler方程式の無限区間 $[0, \infty)$ を、ベルンシュタイン多項式が定義される $[0,1]$ に対応する $[0, M]$ に変換する。
• 収束特性を評価するために、$ L^1 $-ノルムを用いて残差および誤差ノルムを計算し、$ m $ の増加に伴いプロットする。
• 検証のための基準解として、高次連続ルンゲ・クッタ法（Mapleのdverk78）を用い、真の誤差を計算する。

実験結果

リサーチクエスチョン

• RQ1 正確な作用行列（EOMs）は、通常作用行列（OOMs）と比較して、高次非線形Emeden-Fowler方程式の解法において顕著に高い精度を達成できるか？
• RQ2 EOMsとOOMsを用いた場合に、多項式次数 $ m $ の増加に伴う誤差および残差ノルムの収束はどのように異なるか？
• RQ3 EOMsは、$ y^{3/2} $、$ y^{5/2} $、$ \sinh(y) $、$ \sin(y) $、$ e^y $ のような強い非線形性を有する項を、Emeden-Fowler方程式でどれほど適切に扱えるか？
• RQ4 EOMsの使用により、定義域全域でほぼゼロに近い残差関数が得られ、解の忠実性が高まるか？
• RQ5 解析解が存在しない方程式に対しても、EOMsアプローチは、数値ベンチマークによる検証のみを用いて信頼性を持って適用可能か？

主な発見

• EOMに基づく手法は、すべてのテスト非線形性において、OOMsよりも顕著に低い $ L^1 $-ノルム誤差および残差を達成し、誤差収束がより速やかに進行した。
• $ f(x) = 1 $、$ g(x) = y^{3/2}(x) $、$ N = 12 $ の場合、EOM法はすべての $ m $ 値において、OOMsと比較して著しく低い誤差および残差ノルムを示した。
• $ f(x) = 1 $、$ g(x) = \sin(y(x)) $、$ N = 12 $ の場合のノルム1誤差および残差プロットは、EOMsが残差をゼロに近づけ、誤差がより速やかに収束することを示した。
• 本手法は、OOMsが行列演算における近似誤差により達成できなかった、ほぼ正確な残差関数を達成した。
• 高次ルンゲ・クッタ法（dverk78）を基準解として用いた結果、EOMの結果は常にOOMの結果よりも基準解に近いことが確認された。
• 系列作用行列により、$ g(y(x)) $ が高次ベルンシュタイン多項式の線形結合として正確に表現可能となり、強力な非線形性を扱う上で極めて重要であった。