[論文レビュー] An application of Bell polynomials in numerical solving of nonlinear differential equations

Josef Rebenda|arXiv (Cornell University)|Jan 27, 2019

Iterative Methods for Nonlinear Equations被引用数 4

ひとこと要約

本稿では、非線形微分方程式の数値的解法を簡略化するため、部分的通常ベル多項式の新規応用を提案する。ファウ・ディ・ブルーノの公式をベル多項式を用いて再定式化することで、非線形項の記号的微分を回避し、再帰的算術演算による効率的な計算が可能となり、正確な級数解が得られる。2つの初期値問題において、パターン認識により正確な解が特定されたことが示された。

ABSTRACT

Partial ordinary Bell polynomials are used to formulate and prove a version of the Fa\\`{a} di Bruno's formula which is convenient for handling nonlinear terms in the differential transformation. Applicability of the result is shown in two examples of solving the initial value problem for differential equations which are nonlinear with respect to the dependent variable.

研究の動機と目的

• 非線形微分方程式の変換法における非線形項の取り扱いの計算複雑性、特に未知解の合成関数を含む項に対して、解決を図ること。
• ベル多項式の組み合わせ的性質を活用した、非線形項の変換の体系的かつ計算効率の良いアプローチの開発。
• アドミアン分解などの既存の半解析的手法（記号的微分や複雑な積分計算を要する）に対する実用的代替手段の提供。
• 提案手法が再帰的算術演算により正確な級数解を生成することの実証、および閉形式での解の同定可能性の可能性の提示。

提案手法

• 本稿では、合成関数の高階微分を表現するため、部分的通常ベル多項式を用いたファウ・ディ・ブルーノの公式の修正版を導入する。
• 関数 $ f(u(t)) $ の微分変換を、変換係数 $ U(k) $ のベル多項式を含む和として定式化し、$ f $ の記号的微分の必要性を排除する。
• 計算の効率性を確保するため、$ i < k $ の $ U(i) $ を事前に計算した結果を用いて $ U(k) $ を再帰的に計算する。
• 逆微分変換を用いて、係数がベル多項式展開から得られるテイラー型級数として解を再構成する。
• 数値的積分、微分、反復解法を一切避けており、基本的な算術演算のみに依存する。
• 2つの初期値問題に対して手法を検証し、係数のパターンが識別可能であることを示し、正確な解への収束を確認した。

実験結果

リサーチクエスチョン

• RQ1 部分的ベル多項式を用いることで、非線形微分方程式における微分変換の非線形項の処理を記号的微分なしに簡略化できるか？
• RQ2 ファウ・ディ・ブルーノの公式をベル多項式を用いて再定式化することで、半解析的手法における計算効率がどのように向上するか？
• RQ3 この手法を非線形初期値問題に適用した場合、閉形式での正確な解をどの程度まで得られるか？
• RQ4 アドミアン分解やホモトピー摂動法などの既存手法と比較して、このアプローチの計算的利点（必要な演算数）は何か？

主な発見

• 非線形項 $ f(u(t)) $ の微分変換が、部分的通常ベル多項式を用いて成功裏に表現され、記号的微分を必要とせず、変換係数 $ U(k) $ からの直接計算が可能となった。
• 最初の例では、係数が $ U(k)[0] = 0 $（$ k $ が偶数の場合）、$ U(k)[0] = (-1)^{(k-1)/2} \frac{2^k}{k!} $（$ k $ が奇数の場合）というパターンに従い、正確な解 $ u(t) = \sin 2t $ が得られた。
• 2番目の例では、係数が $ \sin 2t $ のテイラー級数と一致し、パターン認識により正確な解が確認された。
• アルゴリズムは計算的に効率的であり、算術演算と再帰的計算のみに依存し、積分や微分、初期推定値の必要がない。
• 係数のパターンが既知の関数と一致する場合、閉形式での正確な解の同定が可能となり、標準的な数値近似を上回る精度が得られる。
• 本手法は他の問題、特に境界値問題へも一般化可能であり、展開点近傍での解析性を保つ関数的近似を提供する。