[论文解读] 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 个一阶常微分方程（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.

研究动机与目标

• 识别具有多项式非线性的特殊一阶 ODE 系统类，使其能获得显式解析解。
• 提出一种新方法，同时约束系数与初值，而非仅约束系数。
• 通过复时间变换与频率调制，将可解框架扩展至周期解。
• 证明即使对于任意的 N 和 M，此类系统也可精确求解，并给出动力学的显式公式。
• 突出这些系统在物理与工程背景中的适用性与潜力，尤其在初值可控的场景中。

提出的方法

• 建立一个包含 N 个一阶 ODE 的系统，其右端项为 N 个自变量中关于次数 M 的齐次多项式。
• 引入形如 zn(t) = zn(0)(1 + Kt)^{1/(1−M)} 的试探解，以探索显式的时间演化行为。
• 推导出代数约束条件 (3b)，将参数 K、系数 cnm1…mN 与初值 zn(0) 联系起来，以确保与 ODE 的一致性。
• 应用复时间变换 τ = (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 个一阶 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) 可显式求解出任意两个系数 cnm，以其余八个系数与初值表示。
• 周期性变体通过变换 (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}，其中复系数约束 (8c) 联系了 KR、KI、初值与系数。
• 当求解系数或 K 时，约束条件 (3b) 在未知量上为线性，从而可实现显式代数求解；非线性情形则需数值处理。
• 该方法适用于任意 N 与 M，所导出的解为精确解，对所有 t ≥ 0 有效，且周期性由复时间映射自然产生。