Sylvester矩阵、子结式、辗转相除法的三者关系(第一部分)

最新推荐文章于 2024-06-19 22:50:03 发布

书中玉

最新推荐文章于 2024-06-19 22:50:03 发布

阅读量231

点赞数 5

文章标签：矩阵算法线性代数

本文链接：https://blog.csdn.net/tianwangwxm/article/details/139106325

版权

设 ${R}$ 为域，非 $0$ 多项式 $a_{0}x^{m} + a_{1}x^{m - 1} + \ldots + a_{m}、G = b_{0}x^{l} + b_{1}x^{l - 1} + \ldots + b_{l}{\in R}\lbrack x\rbrack$ ，并且 $\deg(F) = m、\deg(G) = l、m \geq l > 0$ 。

【Sylvester矩阵定义】

对于一般情况：

$\begin{matrix} \begin{pmatrix} a_{0} & a_{1} & \ldots & a_{m} & 0 & 0 \\ 0 & \ddots & \ddots & \ddots & \ddots & 0 \\ 0 & 0 & a_{0} & a_{1} & \ldots & a_{m} \\ b_{0} & \ldots & b_{l} & 0 & 0 & 0 \\ 0 & b_{0} & .. & b_{l} & 0 & 0 \\ 0 & 0 & \ddots & \ddots & \ddots & 0 \\ 0 & 0 & 0 & b_{0} & \ldots & b_{l} \end{pmatrix} & \begin{matrix} \left. \ \begin{matrix} \ \\ \ \\ \ \end{matrix} \right\} & l行 \\ \left. \ \begin{matrix} \ \\ \ \\ \ \\ \ \end{matrix} \right\} & m行 \end{matrix} \end{matrix}$

以 $m = 8 、 l = 7$ 举例如下：

$\begin{pmatrix} a_{8} & a_{7} & a_{6} & a_{5} & a_{4} & a_{3} & a_{2} & a_{1} & a_{0} & 0 & 0 & 0 & 0 & 0 & 0 \\ 0 & a_{8} & a_{7} & a_{6} & a_{5} & a_{4} & a_{3} & a_{2} & a_{1} & a_{0} & 0 & 0 & 0 & 0 & 0 \\ 0 & 0 & a_{8} & a_{7} & a_{6} & a_{5} & a_{4} & a_{3} & a_{2} & a_{1} & a_{0} & 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & a_{8} & a_{7} & a_{6} & a_{5} & a_{4} & a_{3} & a_{2} & a_{1} & a_{0} & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 & a_{8} & a_{7} & a_{6} & a_{5} & a_{4} & a_{3} & a_{2} & a_{1} & a_{0} & 0 & 0 \\ 0 & 0 & 0 & 0 & 0 & a_{8} & a_{7} & a_{6} & a_{5} & a_{4} & a_{3} & a_{2} & a_{1} & a_{0} & 0 \\ 0 & 0 & 0 & 0 & 0 & 0 & a_{8} & a_{7} & a_{6} & a_{5} & a_{4} & a_{3} & a_{2} & a_{1} & a_{0} \\ b_{7} & b_{6} & b_{5} & b_{4} & b_{3} & b_{2} & b_{1} & b_{0} & 0 & 0 & 0 & 0 & 0 & 0 & 0 \\ 0 & b_{7} & b_{6} & b_{5} & b_{4} & b_{3} & b_{2} & b_{1} & b_{0} & 0 & 0 & 0 & 0 & 0 & 0 \\ 0 & 0 & b_{7} & b_{6} & b_{5} & b_{4} & b_{3} & b_{2} & b_{1} & b_{0} & 0 & 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & b_{7} & b_{6} & b_{5} & b_{4} & b_{3} & b_{2} & b_{1} & b_{0} & 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 & b_{7} & b_{6} & b_{5} & b_{4} & b_{3} & b_{2} & b_{1} & b_{0} & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 & 0 & b_{7} & b_{6} & b_{5} & b_{4} & b_{3} & b_{2} & b_{1} & b_{0} & 0 & 0 \\ 0 & 0 & 0 & 0 & 0 & 0 & b_{7} & b_{6} & b_{5} & b_{4} & b_{3} & b_{2} & b_{1} & b_{0} & 0 \\ 0 & 0 & 0 & 0 & 0 & 0 & 0 & b_{7} & b_{6} & b_{5} & b_{4} & b_{3} & b_{2} & b_{1} & b_{0} \end{pmatrix}$

【子结式定义】

对于 $\leq j < l$ ，称多项式

$S_{j}(x) = subres_{j}(F,G) = \sum_{i = 0}^{j}{\det\left( S_{ij} \right)x^{i}}$

其中 $S_{ij}$ 为通过删除矩阵 $S$ 中 $l$ 行 $F$ 系数中的最后 $j$ 行， $m$ 行 $G$ 系数中的最后 $j$ 行，和最后 $2 j + 1$ 列，但第 $m + l - i - j$ 列除外，所得的子矩阵，这里 $\leq i \leq j < l$ 。

如果 $m > l + 1$ ，则额外拓广定义 $S_{l} = b^{m - l - 1}G、R_{l} = b_{0}^{m - l}$ ，并且当 $l < j < m - 1$ 时，定义 $S_{j}(x) = R_{j} = 0$ 。

如果 $\deg\left( S_{j},x \right) = r < j$ ，则称 $S_{j}$ 为 $r$ 次亏损的；否则，称 $S_{j}$ 为正则的。

容易看出 $S_{0} = R_{0}$ 为 $F$ 和 $G$ 关于 $x$ 的结式。

以 $m = 8 、 l = 7 、 j = 3$ 举例如下：

$S_{30} = \begin{matrix} \begin{matrix} \ & \ & \ & \ & \ & \ & \ & \ & 第m + l - 0 - j = 12列 \end{matrix} & \ \\ \begin{pmatrix} a_{8} & a_{7} & a_{6} & a_{5} & a_{4} & a_{3} & a_{2} & a_{1} & 0 \\ 0 & a_{8} & a_{7} & a_{6} & a_{5} & a_{4} & a_{3} & a_{2} & 0 \\ 0 & 0 & a_{8} & a_{7} & a_{6} & a_{5} & a_{4} & a_{3} & 0 \\ 0 & 0 & 0 & a_{8} & a_{7} & a_{6} & a_{5} & a_{4} & a_{0} \\ b_{7} & b_{6} & b_{5} & b_{4} & b_{3} & b_{2} & b_{1} & b_{0} & 0 \\ 0 & b_{7} & b_{6} & b_{5} & b_{4} & b_{3} & b_{2} & b_{1} & 0 \\ 0 & 0 & b_{7} & b_{6} & b_{5} & b_{4} & b_{3} & b_{2} & 0 \\ 0 & 0 & 0 & b_{7} & b_{6} & b_{5} & b_{4} & b_{3} & 0 \\ 0 & 0 & 0 & 0 & b_{7} & b_{6} & b_{5} & b_{4} & b_{0} \end{pmatrix} & \begin{matrix} \left. \ \begin{matrix} \ \\ \ \\ \ \\ \ \end{matrix} \right\} l - j = 4行 \\ \left. \ \begin{matrix} \ \\ \ \\ \ \\ \ \\ \ \end{matrix} \right\} m - j = 5行 \end{matrix} \end{matrix}$

$S_{31} = \begin{matrix} \begin{matrix} \ & \ & \ & \ & \ & \ & \ & \ & 第m + l - 1 - j = 11列 \end{matrix} & \ \\ \begin{pmatrix} a_{8} & a_{7} & a_{6} & a_{5} & a_{4} & a_{3} & a_{2} & a_{1} & 0 \\ 0 & a_{8} & a_{7} & a_{6} & a_{5} & a_{4} & a_{3} & a_{2} & 0 \\ 0 & 0 & a_{8} & a_{7} & a_{6} & a_{5} & a_{4} & a_{3} & a_{0} \\ 0 & 0 & 0 & a_{8} & a_{7} & a_{6} & a_{5} & a_{4} & a_{1} \\ b_{7} & b_{6} & b_{5} & b_{4} & b_{3} & b_{2} & b_{1} & b_{0} & 0 \\ 0 & b_{7} & b_{6} & b_{5} & b_{4} & b_{3} & b_{2} & b_{1} & 0 \\ 0 & 0 & b_{7} & b_{6} & b_{5} & b_{4} & b_{3} & b_{2} & 0 \\ 0 & 0 & 0 & b_{7} & b_{6} & b_{5} & b_{4} & b_{3} & b_{0} \\ 0 & 0 & 0 & 0 & b_{7} & b_{6} & b_{5} & b_{4} & b_{1} \end{pmatrix} & \begin{matrix} \left. \ \begin{matrix} \ \\ \ \\ \ \\ \ \end{matrix} \right\} l - j = 4行 \\ \left. \ \begin{matrix} \ \\ \ \\ \ \\ \ \\ \ \end{matrix} \right\} m - j = 5行 \end{matrix} \end{matrix}$

$S_{32} = \begin{matrix} \begin{matrix} \ & \ & \ & \ & \ & \ & \ & \ & 第m + l - 2 - j = 10列 \end{matrix} & \ \\ \begin{pmatrix} a_{8} & a_{7} & a_{6} & a_{5} & a_{4} & a_{3} & a_{2} & a_{1} & 0 \\ 0 & a_{8} & a_{7} & a_{6} & a_{5} & a_{4} & a_{3} & a_{2} & a_{0} \\ 0 & 0 & a_{8} & a_{7} & a_{6} & a_{5} & a_{4} & a_{3} & a_{1} \\ 0 & 0 & 0 & a_{8} & a_{7} & a_{6} & a_{5} & a_{4} & a_{2} \\ b_{7} & b_{6} & b_{5} & b_{4} & b_{3} & b_{2} & b_{1} & b_{0} & 0 \\ 0 & b_{7} & b_{6} & b_{5} & b_{4} & b_{3} & b_{2} & b_{1} & 0 \\ 0 & 0 & b_{7} & b_{6} & b_{5} & b_{4} & b_{3} & b_{2} & b_{0} \\ 0 & 0 & 0 & b_{7} & b_{6} & b_{5} & b_{4} & b_{3} & b_{1} \\ 0 & 0 & 0 & 0 & b_{7} & b_{6} & b_{5} & b_{4} & b_{2} \end{pmatrix} & \begin{matrix} \left. \ \begin{matrix} \ \\ \ \\ \ \\ \ \end{matrix} \right\} l - j = 4行 \\ \left. \ \begin{matrix} \ \\ \ \\ \ \\ \ \\ \ \end{matrix} \right\} m - j = 5行 \end{matrix} \end{matrix}$

$S_{33} = \begin{matrix} \begin{matrix} \ & \ & \ & \ & \ & \ & \ & \ & 第m + l - 3 - j = 9列 \end{matrix} & \ \\ \begin{pmatrix} a_{8} & a_{7} & a_{6} & a_{5} & a_{4} & a_{3} & a_{2} & a_{1} & a_{0} \\ 0 & a_{8} & a_{7} & a_{6} & a_{5} & a_{4} & a_{3} & a_{2} & a_{1} \\ 0 & 0 & a_{8} & a_{7} & a_{6} & a_{5} & a_{4} & a_{3} & a_{2} \\ 0 & 0 & 0 & a_{8} & a_{7} & a_{6} & a_{5} & a_{4} & a_{3} \\ b_{7} & b_{6} & b_{5} & b_{4} & b_{3} & b_{2} & b_{1} & b_{0} & 0 \\ 0 & b_{7} & b_{6} & b_{5} & b_{4} & b_{3} & b_{2} & b_{1} & b_{0} \\ 0 & 0 & b_{7} & b_{6} & b_{5} & b_{4} & b_{3} & b_{2} & b_{1} \\ 0 & 0 & 0 & b_{7} & b_{6} & b_{5} & b_{4} & b_{3} & b_{2} \\ 0 & 0 & 0 & 0 & b_{7} & b_{6} & b_{5} & b_{4} & b_{3} \end{pmatrix} & \begin{matrix} \left. \ \begin{matrix} \ \\ \ \\ \ \\ \ \end{matrix} \right\} l - j = 4行 \\ \left. \ \begin{matrix} \ \\ \ \\ \ \\ \ \\ \ \end{matrix} \right\} m - j = 5行 \end{matrix} \end{matrix}$

$S_{3} = \left| \begin{matrix} a_{8} & a_{7} & a_{6} & a_{5} & a_{4} & a_{3} & a_{2} & a_{1} & a_{0} \\ 0 & a_{8} & a_{7} & a_{6} & a_{5} & a_{4} & a_{3} & a_{2} & a_{1} \\ 0 & 0 & a_{8} & a_{7} & a_{6} & a_{5} & a_{4} & a_{3} & a_{2} \\ 0 & 0 & 0 & a_{8} & a_{7} & a_{6} & a_{5} & a_{4} & a_{3} \\ b_{7} & b_{6} & b_{5} & b_{4} & b_{3} & b_{2} & b_{1} & b_{0} & 0 \\ 0 & b_{7} & b_{6} & b_{5} & b_{4} & b_{3} & b_{2} & b_{1} & b_{0} \\ 0 & 0 & b_{7} & b_{6} & b_{5} & b_{4} & b_{3} & b_{2} & b_{1} \\ 0 & 0 & 0 & b_{7} & b_{6} & b_{5} & b_{4} & b_{3} & b_{2} \\ 0 & 0 & 0 & 0 & b_{7} & b_{6} & b_{5} & b_{4} & b_{3} \end{matrix} \right|x^{3} + \left| \begin{matrix} a_{8} & a_{7} & a_{6} & a_{5} & a_{4} & a_{3} & a_{2} & a_{1} & 0 \\ 0 & a_{8} & a_{7} & a_{6} & a_{5} & a_{4} & a_{3} & a_{2} & a_{0} \\ 0 & 0 & a_{8} & a_{7} & a_{6} & a_{5} & a_{4} & a_{3} & a_{1} \\ 0 & 0 & 0 & a_{8} & a_{7} & a_{6} & a_{5} & a_{4} & a_{2} \\ b_{7} & b_{6} & b_{5} & b_{4} & b_{3} & b_{2} & b_{1} & b_{0} & 0 \\ 0 & b_{7} & b_{6} & b_{5} & b_{4} & b_{3} & b_{2} & b_{1} & 0 \\ 0 & 0 & b_{7} & b_{6} & b_{5} & b_{4} & b_{3} & b_{2} & b_{0} \\ 0 & 0 & 0 & b_{7} & b_{6} & b_{5} & b_{4} & b_{3} & b_{1} \\ 0 & 0 & 0 & 0 & b_{7} & b_{6} & b_{5} & b_{4} & b_{2} \end{matrix} \right|x^{2} + \left| \begin{matrix} a_{8} & a_{7} & a_{6} & a_{5} & a_{4} & a_{3} & a_{2} & a_{1} & 0 \\ 0 & a_{8} & a_{7} & a_{6} & a_{5} & a_{4} & a_{3} & a_{2} & 0 \\ 0 & 0 & a_{8} & a_{7} & a_{6} & a_{5} & a_{4} & a_{3} & a_{0} \\ 0 & 0 & 0 & a_{8} & a_{7} & a_{6} & a_{5} & a_{4} & a_{1} \\ b_{7} & b_{6} & b_{5} & b_{4} & b_{3} & b_{2} & b_{1} & b_{0} & 0 \\ 0 & b_{7} & b_{6} & b_{5} & b_{4} & b_{3} & b_{2} & b_{1} & 0 \\ 0 & 0 & b_{7} & b_{6} & b_{5} & b_{4} & b_{3} & b_{2} & 0 \\ 0 & 0 & 0 & b_{7} & b_{6} & b_{5} & b_{4} & b_{3} & b_{0} \\ 0 & 0 & 0 & 0 & b_{7} & b_{6} & b_{5} & b_{4} & b_{1} \end{matrix} \right|x^{1} + \left| \begin{matrix} a_{8} & a_{7} & a_{6} & a_{5} & a_{4} & a_{3} & a_{2} & a_{1} & 0 \\ 0 & a_{8} & a_{7} & a_{6} & a_{5} & a_{4} & a_{3} & a_{2} & 0 \\ 0 & 0 & a_{8} & a_{7} & a_{6} & a_{5} & a_{4} & a_{3} & 0 \\ 0 & 0 & 0 & a_{8} & a_{7} & a_{6} & a_{5} & a_{4} & a_{0} \\ b_{7} & b_{6} & b_{5} & b_{4} & b_{3} & b_{2} & b_{1} & b_{0} & 0 \\ 0 & b_{7} & b_{6} & b_{5} & b_{4} & b_{3} & b_{2} & b_{1} & 0 \\ 0 & 0 & b_{7} & b_{6} & b_{5} & b_{4} & b_{3} & b_{2} & 0 \\ 0 & 0 & 0 & b_{7} & b_{6} & b_{5} & b_{4} & b_{3} & 0 \\ 0 & 0 & 0 & 0 & b_{7} & b_{6} & b_{5} & b_{4} & b_{0} \end{matrix} \right|$

【辗转相除法】

以 $m = 8 、 l = 7$ 为例子，并令 $F_{0} = F、F_{1} = G$ ，可得

$F_{0} = Q_{1} \times F_{1} + F_{2}\ \ \ \ \ \ \ \ \ \ \deg\left( F_{0} \right) = 8\ \ \ \ \ \ \deg\left( F_{1} \right) = 7\ \ \ \ \ \ \deg\left( F_{2} \right) = 6$

$F_{1} = Q_{2} \times F_{2} + F_{3}\ \ \ \ \ \ \ \ \ \ \deg\left( F_{1} \right) = 7\ \ \ \ \ \ \deg\left( F_{2} \right) = 6\ \ \ \ \ \ \deg\left( F_{3} \right) = 5$

$F_{2} = Q_{3} \times F_{3} + F_{4}\ \ \ \ \ \ \ \ \ \ \deg\left( F_{2} \right) = 6\ \ \ \ \ \ \deg\left( F_{3} \right) = 5\ \ \ \ \ \ \deg\left( F_{4} \right) = 4$

$F_{3} = Q_{4} \times F_{4} + F_{5}\ \ \ \ \ \ \ \ \ \ \deg\left( F_{3} \right) = 5\ \ \ \ \ \ \deg\left( F_{4} \right) = 4\ \ \ \ \ \ \deg\left( F_{5} \right) = 3$

$F_{4} = Q_{5} \times F_{5} + F_{6}\ \ \ \ \ \ \ \ \ \ \deg\left( F_{4} \right) = 4\ \ \ \ \ \ \deg\left( F_{5} \right) = 3\ \ \ \ \ \ \deg\left( F_{6} \right) = 2$

$F_{5} = Q_{6} \times F_{6} + F_{7}\ \ \ \ \ \ \ \ \ \ \deg\left( F_{5} \right) = 3\ \ \ \ \ \ \deg\left( F_{6} \right) = 2\ \ \ \ \ \ \deg\left( F_{7} \right) = 1$

$F_{6} = Q_{7} \times F_{7} + F_{8}\ \ \ \ \ \ \ \ \ \ \deg\left( F_{6} \right) = 2\ \ \ \ \ \ \deg\left( F_{7} \right) = 1\ \ \ \ \ \ \deg\left( F_{8} \right) = 0$

注意，上面的举例仅仅是最一般的情况，实际算法运行步骤可能小于 $7$ 步。