[論文レビュー] Explicitly Solvable Systems of First-order Ordinary Differential Equations with Polynomial Right-hand Sides, and Their Periodic Variants

F. Calogero, Farrin Payandeh|arXiv (Cornell University)|Jun 11, 2021

Advanced Differential Equations and Dynamical Systems被引用数 1

ひとこと要約

本稿では、係数と初期条件の間の制約を課すことにより、任意の次数 M の多項式右辺をもつ N 個の 1 階常微分方程式（ODE）系を明示的に解ける形にしている。主な結果は、パラメータと初期データが特定の代数的条件を満たす場合に、zn(t) = zn(0)(1 + Kt)^{1/(1−M)} という閉形式解が系を満たすものであり、複素時間変換と周波数変調を用いて周期的解が導出されている。

ABSTRACT

In this Letter we identify special systems of (an arbitrary number) N of first-order Ordinary Differential Equations with homogeneous polynomials of arbitrary degree M on their right-hand sides, which feature very simple explicit solutions; as well as variants of these systems--with right-hand sides no more homogeneous--which feature periodic solutions. A novelty of these findings is to consider special systems characterized by constraints involving both their parameters and their initial data.

研究の動機と目的

• 多項式非線形性をもつ 1 階 ODE 系の特殊クラスで、明示的な解析的解が得られるものを同定すること。
• 係数のみでなく初期条件を同時に制約するという新しいアプローチを考案すること。
• 複素時間変換と周波数変調を用いて、可解枠組みを周期的解へと拡張すること。
• このような系が任意の N および M に対して正確に解けること、および動的挙動を明示的な式で記述できることを示すこと。
• 初期データの制御が可能である物理的・工学的文脈におけるこれらの系の応用可能性と潜在的価値を強調すること。

提案手法

• N 個の従属変数における次数 M の斉次多項式として右辺を持つ N 個の 1 階 ODE 系を定式化する。
• 時間発展を明示的に探索するため、zn(t) = zn(0)(1 + Kt)^{1/(1−M)} の形の試行解を導入する。
• ODE と整合性を持つために、パrameter K、係数 cnm1…mN、初期値 zn(0) を結ぶ代数的制約 (3b) を導出する。
• 複素時間変換 τ = (exp(iωt) − 1)/(iω) を適用し、新しい変数 xn(t) + iyn(t) = exp(iωt/(M−1)) zn(τ) を定義して周期的解を生成する。
• 元の系を、実部と虚部を含む 2N 個の実数 ODE からなる新しい自己同次系 (6a) に変換し、修正された非線形項 Zn(t) を含む。
• 導出された制約のもとで、変換後の系が周期 T = 2π/|ω| あるいはその整数倍の周期性を示す解を持つことを確認する。

実験結果

リサーチクエスチョン

• RQ1 次数 M の斉次多項式非線形性をもつ N 個の 1 階 ODE 系が、どのような条件下で明示的に解けるか？
• RQ2 係数と初期条件の両方に及ぶ制約が、このような系において閉形式解をもたらす仕組みは何か？
• RQ3 変数変換を用いることで、明示的解の枠組みを周期的解の生成に拡張できるか？
• RQ4 複素時間と周波数変調が、非周期的可解系を周期的系へと変換する役割は何か？
• RQ5 初期データと係数の相互作用が、系の動的挙動をどのように制御可能にするか？

主な発見

• 方程式 (3b) に示される N 個の制約を満たす限り、系 (1) は明示的解 zn(t) = zn(0)(1 + Kt)^{1/(1−M)} を持つ。この制約は K、cnm1…mN、および zn(0) を結ぶ。
• N=2 かつ M=4 の場合、解は zn(t) = zn(0)[1 + Kt]^{-1/3} となり、制約 (4c) は、10 個の係数 cnm のうち 2 個を残りの 8 個と初期データで明示的に表すことができる。
• 周期的バージョンは変換 (5) を用いて構築され、実時間 t を複素時間 τ に写像し、系 (6) を得る。この系は xn(t) および yn(t) に対して周期 T = 2π/|ω| の解を持つ。
• N=2 かつ M=4 の周期的ケースでは、解は ζn(t) = [xn(0)+iyn(0)] exp(iωt/3) · [1 + (KR+iKI)(exp(iωt)−1)/(iω)]^{-1/3} で与えられ、KR、KI、初期データ、係数を結ぶ複素制約 (8c) が成立する。
• 制約 (3b) は、未知係数または K を求める際、線形であるため明示的な代数的解が得られるが、非線形ケースでは数値的処理が必要となる。
• この手法は任意の N および M に適用可能であり、導出された解はすべての t ≥ 0 に対して正確に成り立ち、周期性は複素時間写像から自然に生じる。