[论文解读] Polynomial Time Nondimensionalisation of Ordinary Differential Equations via their Lie Point Symmetries

Évelyne Hubert, Alexandre Sedoglavic|arXiv (Cornell University)|Apr 13, 2006

Nonlinear Waves and Solitons参考文献 11被引用 5

一句话总结

本文提出了一种针对常微分方程（ODE）的多项式时间算法，利用李点对称性（特别是缩放对称性和平移对称性）实现无量纲化，从而减少参数数量。通过计算这些对称性的有理不变量，该方法将 ODE 系统重写为参数更少的无量纲坐标形式，计算复杂度被限制在 O(ℓ³)，其中 ℓ 是参数的数量。

ABSTRACT

Lie group theory states that knowledge of a $m$-parameters solvable group of symmetries of a system of ordinary differential equations allows to reduce by $m$ the number of equation. We apply this principle by finding dilatations and translations that are Lie point symmetries of considered ordinary differential system. By rewriting original problem in an invariant coordinates set for these symmetries, one can reduce the involved number of parameters. This process is classically call nondimensionalisation in dimensional analysis. We present an algorithm based on this standpoint and show that its arithmetic complexity is polynomial in input's size.

研究动机与目标

• 开发一种系统性算法，用于减少物理和生物模型中出现的 ODE 系统的参数数量。
• 将无量纲化过程形式化为基于对称性的变换，利用李点对称性。
• 证明有理不变量的计算与系统重写过程在输入规模上具有多项式时间复杂度。
• 提供一种系统化、算法化的参数约简框架，避免依赖人为的、非系统性的手动操作。
• 将量纲分析扩展为基于李群理论和不变量理论的计算代数框架。

提出的方法

• 识别使 ODE 系统保持不变的一参数群缩放对称性和平移对称性。
• 通过求解偏微分方程组的判定系统，构造这些对称群的无穷小生成元。
• 通过线性代数方法（特别是对称作用系数矩阵进行高斯消元）计算对称群的有理不变量。
• 将原始 ODE 系统用这些有理不变量重写，得到参数更少的无量纲系统。
• 利用对称性生成的李代数的可解性，在算法中优先处理平移对称性，再处理缩放对称性。
• 通过在对称指数上应用线性代数技术，确保算法的算术复杂度保持多项式，其上界为 O(ℓ³)。

实验结果

研究问题

• RQ1 能否将 ODE 系统的无量纲化过程形式化为一种具有保证多项式时间复杂度的基于对称性的算法？
• RQ2 如何系统地计算由李点对称性导出的有理不变量，以减少 ODE 中的参数数量？
• RQ3 计算这些不变量以及将系统重写为不变量坐标下的计算复杂度是多少？
• RQ4 为最大化参数约简，应以何种顺序应用不同类型的对称性（平移与缩放）？
• RQ5 该方法能否扩展以处理超出缩放和平移的更一般对称类型，例如非线性对称？

主要发现

• 该算法使用高斯消元法从李点对称性计算有理不变量，确保计算复杂度为 O(ℓ³) 的多项式时间上界。
• 该方法通过利用对称群成功减少了 ODE 系统的参数数量，已在 Verhulst 逻辑斯蒂模型和捕食者-猎物系统中得到验证。
• 对于具有线性捕食的 Verhulst 模型，系统通过不变量（t̃ = (a−c)t, x̃ = b/(a−c) x）从三个参数减少到一个参数，得到无量纲形式 dx̃/dt̃ = x̃(1−x̃)。
• 在捕食者-猎物模型中，通过选择缩放变换消除了三个参数（r, k₁, k₂），最终仅剩下三个参数：s, e, h。
• 该方法对对称基的选择具有鲁棒性，通过在高斯消元过程中合理安排变量顺序，可保持用户对哪些参数应被消除的控制权。
• 该框架可扩展至非标准对称，例如 FitzHugh-Nagumo 模型中的一个参数群 ∂/∂y + b∂/∂a + ∂/∂d，通过不变量 ỹ = y−d 和 â = a+bd 实现参数约简。