[论文解读] Jacobi's bound and normal forms computations. A historical survey

François Ollivier|arXiv (Cornell University)|Nov 13, 2009

Polynomial and algebraic computation参考文献 22被引用 6

一句话总结

本文对微分代数中的雅可比界（Jacobi's bound）进行了全面的历史与技术综述——这是关于常微分方程组阶数的一个基本结果。文章重构了雅可比的原始方法，包括其计算该界值的多项式时间算法、等式成立的截断雅可比行列式条件，以及最短约化方法，同时追溯了这些方法在算法微分代数与计算机代数系统中的重新发现过程及其现代相关性。

ABSTRACT

Jacobi is one of the most famous mathematicians of his century. His name is attached to many results in various fields of mathematics and his complete works in seven volumes have been available since the end of the XIXth century and are very often quoted in many papers. It is then surprising that some of his results may have fallen into oblivion, at least in part. We will try to describe some of Jacobi's results on ordinary differential equations and the available, published or unpublished material he left. We will then expose the selective interests of his followers and their own contributions. There are in fact many interrelated results: a bound on the order of a differential system, a necessary and sufficient condition, given by a determinant, for the bound to be reached, an algorithm to compute the bound in polynomial time, and processes for computing normal forms using as few derivatives as possible. We give for all of them the form under which they could have been proved or rediscovered, sometimes independently of Jacobi's findings. In conclusion, we give the state of the art and suggest some possible applications of Jacobi's bound to improve some algorithms in differential algebra.

研究动机与目标

• 重构并阐明雅可比关于微分系统（特别是常微分方程组阶数的界）的原始成果。
• 记录雅可比未发表手稿的历史轨迹及其对后世数学家的影响。
• 分析雅可比算法及截断雅可比行列式条件在判断界是否紧致时的计算意义。
• 评估雅可比最短约化方法在最小化标准型计算中导数数量方面的相关性。
• 将历史发展与现代微分代数及计算机代数系统中的算法改进联系起来。

提出的方法

• 使用最大横截和 $ J = \max_{\sigma \in S_n} \sum_{i=1}^n a_{i,\sigma(i)} $ 重构雅可比界，其中 $ a_{i,j} $ 表示 $ u_i $ 在 $ x_j $ 中的阶数。
• 应用雅可比算法通过一个典范 $ \lambda \in \mathbb{N}^n $ 在 $ \mathcal{O}(n^3) $ 时间内计算 $ J $，避免对 $ n! $ 个排列进行枚举。
• 引入截断雅可比矩阵 $ \nabla $，其元素为 $ \partial u_i / \partial x_j^{(\alpha_i + \beta_j)} $，其中 $ \alpha_i = \Lambda - \ell_i $，$ \beta_j = \max_i a_{i,j} - \alpha_i $。
• 建立截断雅可比行列式条件：在正则性假设下，当且仅当 $ |\nabla| \neq 0 $ 时，系统阶数等于界 $ J $。
• 分析最短约化方法，断言 $ \ell_i $ 阶导数对计算标准型而言既必要又充分，且该数量不可再减少。
• 综述现代对雅可比界的证明，包括 Kondratieva 等人（1982, 2008）以及 Volevich（1960）的成果，其在独立性或正则性条件下成立。

实验结果

研究问题

• RQ1 雅可比如何推导出常微分方程组阶数的界？最大横截和在该推导中起到什么作用？
• RQ2 为何雅可比计算该界值的算法在一百多年间被忽视，后来又如何被 Kuhn 和 Egerváry 重新发现？
• RQ3 在何种条件下，截断雅可比行列式 $ |\nabla| \neq 0 $ 能够推出系统阶数等于界 $ J $？该判据存在哪些局限性？
• RQ4 最短约化方法在多大程度上最小化了计算标准型所需的导数数量？与其他约化策略相比有何差异？
• RQ5 雅可比成果对改进微分代数与计算机代数系统中算法的现代意义是什么？

主要发现

• 雅可比界 $ J $ 定义为矩阵 $ A = (a_{i,j}) $ 的最大横截和，为 $ n $ 个未知数的 $ n $ 个常微分方程组的阶数提供了上界。
• 雅可比算法通过寻找一个典范 $ \lambda $，在 $ \mathcal{O}(n^3) $ 次运算内计算 $ J $，避免了对 $ n! $ 个排列的指数级搜索。
• 在强界与正则性假设下，截断雅可比行列式条件 $ |\nabla| \neq 0 $ 是系统阶数等于 $ J $ 的充要条件。
• 最短约化方法确保 $ \ell_i $ 阶导数对计算标准型而言既必要又充分，且该数量无法进一步减少。
• 该界及其相关算法被 Kuhn（1955）、Volevich（1960）以及 Kondratieva 等人（1982, 2008）独立重新发现，后者在独立性假设下提供了证明。
• 现代应用包括微分代数中复杂度界限的改进与算法设计，特别是在特征集与形式幂级数解的计算中。