[論文レビュー] Construction of Single-valued Solutions for Nonintegrable Systems with the Help of the Painleve Test

S. Yu. Vernov|arXiv (Cornell University)|Jul 27, 2004

Nonlinear Waves and Solitons参考文献 6被引用数 6

ひとこと要約

この論文は、Painlevéテストを用いて非可積分な非線形常微分方程式に対する単値的解（楕円型、三角関数型、有理型）を構成するためのコンピュータ代数アルゴリズムを提示する。ローレンツ級数解を固定パラメータ変数における非線形代数的系に変換することで、1階方程式の次数に依存しない方法を採用し、MapleおよびREDUCEにおけるグレブナー基底を用いた効率的なグローバル解の計算を可能にする。

ABSTRACT

The Painleve test is very useful to construct not only the Laurent-series solutions but also the elliptic and trigonometric ones. Such single-valued functions are solutions of some polynomial first order differential equations. To find the elliptic solutions we transform an initial nonlinear differential equation in a nonlinear algebraic system in parameters of the Laurent-series solutions of the initial equation. The number of unknowns in the obtained nonlinear system does not depend on number of arbitrary coefficients of the used first order equation. In this paper we describe the corresponding algorithm, which has been realized in REDUCE and Maple.

研究の動機と目的

• 非可積分な非線形常微分方程式に対するグローバルな単値的解（楕円型、三角関数型、有理型）を体系的に構成するための手法を開発すること。
• 従来のPainlevéテストの応用における制限を克服し、局所的なローレンツ級数解を用いてグローバルな1階多項式常微分方程式を導出すること。
• 不要な計算を避けるために、関連する項のみを計算するコンピュータ代数アルゴリズムを設計すること。
• 得られる代数的系における未知数の数が、1階常微分方程式の次数に依存しないようにすることにより、スケーラビリティを確保すること。

提案手法

• 可動特異点近傍における局所的ローレンツ級数解の支配的挙動と共鳴を特定するため、Painlevéテストを適用する。
• 一般解が単値的となるように保証するBriot-Bouquet定理を満たす形で、$\sum_{k=0}^{m}\sum_{j=0}^{2m-2k}a_{jk}y^j y_t^k = 0$ の形の1階自己同次多項式常微分方程式を構築し、$a_{0m} = 1$ とする。
• 常微分方程式をローレンツ級数展開 $F(y_t, y) = \sum_{s=-m(p+1)}^{N_{\text{max}}}} K_s t^s$ に変換し、係数 $K_s = 0$ を等式として代入することで、$a_{ij}$ を変数とする過剰決定代数的系を形成する。
• MapleおよびREDUCEにおける記号計算を用いて、ローレンツ展開における必要な項のみを生成し、計算のオーバーヘッドを最小限に抑える。
• この系を、$\lambda_1, \lambda_2, H$ およびローレンツ級数からの2つの任意係数という5つのパラメータでのみの非線形代数的系に簡略化する。この簡略化は $m$ に依存しない。
• 得られた系をグレブナー基底法を用いて解き、有効な単値的解を生成する一貫したパラメータ集合を特定する。

実験結果

リサーチクエスチョン

• RQ1 Painlevéテストを、形式的ローレンツ級数を超えてグローバルな単値的解を体系的に構成するために効果的に応用できるか？
• RQ2 1階常微分方程式の次数 $m$ に依存しないように、代数的系における未知数の数をどのように分離できるか？
• RQ3 ローレンツ展開における必要な項のみを効率的に計算する最適な方法は何か？これにより、無駄な計算を回避できるか？
• RQ4 この手法は、一般化されたHénon–Heiles系のような非可積分系に対し、楕円型、三角関数型、または有理型の解を生成できるか？
• RQ5 MapleやREDUCEのような記号計算ツールをどのように最適化すれば、所定のローレンツ級数解を持つ1階常微分方程式の導出を自動化できるか？

主な発見

• アルゴリズムは、1階常微分方程式の次数 $m$ にかかわらず、5つの未知数のみを含む非線形代数的系に変換することで、非可積分系に対するグローバルな単値的解を効果的に構成する。
• ローレンツ展開における不要な項の計算を回避し、必要な成分のみを生成することで、単純なアプローチに比べて計算効率が著しく向上する。
• C = -16/5 の一般化されたHénon–Heiles系に対して、$x^2$ が楕円関数であることが同定され、4階方程式を導出することなく解を構成可能である。
• アルゴリズムはMapleおよびREDUCEの両方で実装されており、`equlaurlist` や `equalist` といった手続きにより、ローレンツ級数および係数系の自動生成が可能である。
• パラメータにおける得られた非線形系はグレブナー基底法により解けるため、有効な特殊解を生成する一貫したパラメータ集合の回復が可能である。
• 本手法は一般性を有し、ローレンツ級数解が既知であれば、任意の系に適用可能であり、解が退化した楕円関数または有理関数である場合にも適用可能である。