[论文解读] Subresultants in Recursive Polynomial Remainder Sequence

Akira Terui|arXiv (Cornell University)|Jun 3, 2008

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

一句话总结

本文引入了递归多项式余数序列（递归PRS）与递归结式——这是经典结式理论在具有非平凡最大公因式多项式上的扩展。通过递归地对最大公因式及其导数应用PRS，该方法能够借助递归Sturm序列稳定计算实根的重数，而递归结式则以行列式形式表示系数与原始多项式系数之间的关系，从而在近似代数计算中提升数值稳定性。

ABSTRACT

We introduce concepts of "recursive polynomial remainder sequence (PRS)" and "recursive subresultant," and investigate their properties. In calculating PRS, if there exists the GCD (greatest common divisor) of initial polynomials, we calculate "recursively" with new PRS for the GCD and its derivative, until a constant is derived. We call such a PRS a recursive PRS. We define recursive subresultants to be determinants representing the coefficients in recursive PRS by coefficients of initial polynomials. Finally, we discuss usage of recursive subresultants in approximate algebraic computation, which motivates the present work.

研究动机与目标

• 将经典结式理论扩展至处理由重复计算最大公因式与导数所产生的递归多项式余数序列（递归PRS）。
• 将递归结式定义为行列式，用以表示递归PRS中元素与原始多项式系数之间的关系。
• 改进近似代数计算中的数值稳定性，特别是检测递归Sturm序列中微小或近零首项系数的能力。
• 建立类似于结式基本定理的递归结式理论基础。
• 通过将递归结式与Sylvester矩阵条件关联，支持浮点数运算中的鲁棒零识别。

提出的方法

• 将递归PRS定义为：递归地对最大公因式及其导数进行PRS计算，直至得到常数项为止。
• 将递归结式引入为由初始多项式系数构成的行列式，用于表示递归PRS中的元素。
• 建立递归结构：每一新的PRS序列均由前一序列中最后一个非常数项及其导数计算得出。
• 利用矩阵变换与列交换，将递归结式行列式与后续PRS层级的结式关联起来。
• 利用Sylvester矩阵条件，并通过首项系数缩放进行归一化，以稳定数值计算。
• 将递归结式框架应用于递归Sturm序列，以计数具有重数的实根。

实验结果

研究问题

• RQ1 经典结式理论如何被扩展以处理由重复计算最大公因式与导数所产生的递归多项式余数序列？
• RQ2 递归结式的代数结构是什么？它们如何与递归PRS中多项式系数相关联？
• RQ3 递归结式能否用于改进近似代数计算的数值稳定性，特别是检测微小或近零首项系数？
• RQ4 Sylvester矩阵需满足何种条件，才能确保在递归Sturm序列中实现可靠的零识别？
• RQ5 如何利用递归结式计算具有重根或根位置相近的多项式的实根数量（含重数）？

主要发现

• 递归结式被定义为行列式，用于表示初始多项式系数在递归PRS中系数的表达形式。
• 本文证明了递归结式的“基本定理”类比，表明递归结式决定了递归PRS的结构。
• 递归结式可表示为低层级结式的缩放乘积，缩放因子涉及中间最大公因式首项系数的幂次。
• 证明了递归结式矩阵的行列式与连续PRS层级中结式的乘积成正比，且给出了显式的缩放常数。
• 该框架通过分析递归结式结构，实现了在递归Sturm序列中对微小首项系数的稳定检测。
• 该方法通过将结式行列式与Sylvester矩阵条件关联，支持浮点数运算中的可靠零识别，从而提升了实根计数的鲁棒性。