Gröbner基方法入门第I部分:项序和约化

Gröbner基方法入门第I部分:项序和约化

理想的准素分解与矩阵分解是计算代数的核心问题,它们在计算机代数、计算代数几何、代数编码和密码学、多维系统理论等学科都有非常重要的理论意义与应用价值.Gröbner基理论是研究理想的准素分解与矩阵素分解的重要工具之一.更简单地概括:Gröbner基方法是从任意一个多项式理想的一组给定生成元,计算另一组性质良好的生成元,并称为该理想的Gröbner基.它能用来判定任意多项式是否属于该理想.

Gröbner基的代数算法是Buchberger在1960年代开发的,并进一步被他和许多其他研究人员进行了改进.Gröbner基理论及其前身(如Hironaka的标准基)在代数几何和交换代数(尤其是多项式理想理论)方面提供了丰富的理论框架.Gröbner基的代数算法在理论物理学、应用科学和工程学中具有广泛的应用;Gröbner基成功的主要原因是,数学、科学和工程学中的许多问题都可以用多元多项式方程组表示(例如,理想,模块和矩阵),其中Gröbner基的作用类似于一元多项式Euclidean环的算法中Euclidean除法的作用.

项序和约化

先在集合 S S S上引入“项序”概念:

项序
  • △ \triangle 定义1.1:
    全序: S S S中两个不同元 a , b a,b a,b,要么 a < b a<b a<b,要么 b < a b<a b<a;
    良序: S S S上不存在无穷递减的链,即 S S S的任意非空集有极小元;

然后我们把多项式和理想假定在 K [ x ] \mathcal{K}[\boldsymbol{x}] K[x]中, K \mathcal{K} K是域. K [ x ] \mathcal{K}[\boldsymbol{x}] K[x]中所有的的全体记为 T = { x 1 α 1 , x 2 α 2 , … , x n α n } T=\left\{x_{1}^{\alpha 1},x_{2}^{\alpha 2},\ldots,x_{n}^{\alpha n}\right\} T={x1α1,x2α2,,xnαn},也记为 x α x^{\alpha} xα.项序是一种排序的方式,细分有纯字典序 < plex  <_{\text {plex }} <plex 、全幂序或分次字典序 < tdeg <_{\text {tdeg}} <tdeg、分次反字典序 < dileg <_{\text {dileg}} <dileg.

  • △ \triangle 定义1.2:
    项集合 T \mathcal{T} T上的序分别称为纯字典序 < plex  <_{\text {plex }} <plex 全幂序 < tdeg <_{\text {tdeg}} <tdeg分次反字典序 < dileg <_{\text {dileg}} <dileg,如果对 α = ( a 1 , ⋯   , a n ) , β = ( b 1 , ⋯   , b n ) ∈ T \alpha=\left(a_{1}, \cdots, a_{n}\right), \beta=\left(b_{1}, \cdots, b_{n}\right) \in \mathcal{T} α=(a1,,an),β=(b1,,bn)T有:
    i) α < plex  β ⟺ ∃ i ( 1 ⩽ i ⩽ n ) \alpha<_{\text {plex }} \beta \Longleftrightarrow \exists i(1 \leqslant i \leqslant n) α<plex βi(1in),使得 a i < b i a_{i}<b_{i} ai<bi,而 a j = b j ( 1 ⩽ j < i ) a_{j}=b_{j}(1 \leqslant j<i) aj=bj(1j<i);
    ii) α < tdeg  β ⟺ ∑ i = 1 n a i < ∑ i = 1 n b i \alpha<_{\text {tdeg }} \beta \Longleftrightarrow \sum_{i=1}^{n} a_{i}<\sum_{i=1}^{n} b_{i} α<tdeg βi=1nai<i=1nbi ∑ i = 1 n a i = ∑ i = 1 n b i \sum_{i=1}^{n} a_{i}=\sum_{i=1}^{n} b_{i} i=1nai=i=1nbi,且 ∃ i ( 1 ⩽ i ⩽ n ) \exists i(1 \leqslant i \leqslant n) i(1in),使得 a i < b i a_{i}<b_{i} ai<bi,而 a j = b j ( 1 ⩽ j < i ) a_{j}=b_{j}(1 \leqslant j<i) aj=bj(1j<i);
    iii) α < dilex  β ⟺ ∑ i = 1 n a i < ∑ i = 1 n b i \alpha<_{\text {dilex }} \beta \Longleftrightarrow \sum_{i=1}^{n} a_{i}<\sum_{i=1}^{n} b_{i} α<dilex βi=1nai<i=1nbi ∑ i = 1 n a i = ∑ i = 1 n b i \sum_{i=1}^{n} a_{i}=\sum_{i=1}^{n} b_{i} i=1nai=i=1nbi,且 ∃ i ( 1 ⩽ i ⩽ n ) \exists i(1 \leqslant i \leqslant n) i(1in),使得 a i > b i a_{i}>b_{i} ai>bi,而 a j = b j ( 1 ⩽ j < i ) a_{j}=b_{j}(1 \leqslant j<i) aj=bj(1j<i);

项序的重要性质是良序.在多项式环 K \mathcal{K} K中,若理想 J \mathcal{J} J由一组项 x α x^{\alpha} xα生成,则 J = ⟨ x α ⟩ \mathcal{J}=\langle x^{\alpha}\rangle J=xα为项理想.这里 α ∈ A \alpha \in A αA, A ∈ N n A \in \mathbb{N}^n ANn.于是我们得出这个结论:

  • ♡ \heartsuit 定理1.3:
    i) x β ∈ J ⟺ ∃ α ∈ A , x α ∣ x β x^{\beta} \in J \Longleftrightarrow \exists \alpha \in A, x^{\alpha} \mid x^{\beta} xβJαA,xαxβ;
    ii) ∀ F ∈ K [ x ] , F ∈ J ⟺ F \forall F \in \mathcal{K}[\boldsymbol{x}], F \in \mathcal{J} \Longleftrightarrow F FK[x],FJF 的每个项都属于 J ⟺ F \mathcal{J} \Longleftrightarrow F JF J \mathcal{J} J中项的 K \mathcal{K} K线性组合.

下面介绍Dickson引理:每个项理想都是有限生成的,即都存在有限生成元;

  • ♢ \diamondsuit 引理1.4(Dickson引理):
    对于 K [ x ] \mathcal{K}[\boldsymbol{x}] K[x]中的任意项理想 a = ⟨ x α : α ∈ S ⟩ \mathfrak{a}=\left\langle x^{\alpha}: \alpha \in S\right\rangle a=xα:αS,均存在 α ( 1 ) , . . . , α ( s ) ∈ S \alpha(1),...,\alpha(s) \in S α(1),...,α(s)S使得 a = ⟨ x α ( 1 ) , . . . , x α ( s ) ⟩ \mathfrak{a}=\langle x^{\alpha(1)},...,x^{\alpha(s)} \rangle a=xα(1),...,xα(s).

介绍本节中常用的记号.设 F ∈ K [ x ] F \in \mathcal{K}[\boldsymbol{x}] FK[x]在给定项序下按降序排列为:

F = f 1 x α 1 + ⋯ + f k x α k F=f_{1} \boldsymbol{x}^{\alpha_{1}}+\cdots+f_{k} \boldsymbol{x}^{\alpha_{k}} F=f1xα1++fkxαk

那么 lt ⁡ ( F ) = x α 1 , lm ⁡ ( F ) = f 1 x α 1 , lc ⁡ ( F ) = f 1 \operatorname{lt}(F)=\boldsymbol{x}^{\alpha_{1}}, \operatorname{lm}(F)=f_{1} \boldsymbol{x}^{\alpha_{1}}, \operatorname{lc}(F)=f_{1} lt(F)=xα1,lm(F)=f1xα1,lc(F)=f1分别为首项首单项式首项系数,我们常用 T ( F ) = { x α 1 , ⋯   , x α k } \mathcal{T}(F)=\left\{\boldsymbol{x}^{\alpha_{1}}, \cdots, \boldsymbol{x}^{\alpha_{k}}\right\} T(F)={xα1,,xαk} F F F 中出现的项的集合;

多项式约化

多项式约化就是一种多元多项式的除法,是域上一元多项式Euclid除法在多元多项式的推广.

  • △ \triangle 定义1.5:
    给定 T \mathcal{T} T 上项序 < < <,并设 F , P ∈ K [ x ] F,P \in \mathcal{K}[\boldsymbol{x}] F,PK[x],且 P ≠ 0 P \neq 0 P=0.如果存在 t ∈ T ( F ) , s ∈ T t \in \mathcal{T}(F), s \in \mathcal{T} tT(F),sT 使得 t = s ⋅ lt ⁡ ( P ) t=s \cdot \operatorname{lt}(P) t=slt(P),则称 F F F 是模 P P P 可约化的. 这时,令:

G = F − c lc ⁡ ( P ) ⋅ s ⋅ P G=F-\frac{c}{\operatorname{lc}(P)} \cdot s \cdot P G=Flc(P)csP

其中 c c c F F F 关于 t t t 的系数,并说 F F F 通过模 P P P 消去 t t t 一步约化 G , G, G, 记作 F → t P G F \xrightarrow[t]{P} G FP tG,或简记为 F ⟶ P G F \stackrel{P}{\longrightarrow} G FPG;

P \mathbb{P} P K [ x ] \mathcal{K}[\boldsymbol{x}] K[x] 中的多项式组. 如果存在多项式 P ∈ P P \in \mathbb{P} PP 使得 F ⟶ P G F \stackrel{P}{\longrightarrow} G FPG, 则称
F F F P \mathbb{P} P 可约化的,记为 F ⟶ P G . F \stackrel{\mathbb{P}}{\longrightarrow} G . FPG. 如果没有这样的 P P P 存在,则称 F F F P \mathbb{P} P约化的或为范式.

如果 F F F 可由 P \mathbb{P} P中的多项式经多步约化到范式 R R R,则记为 F → ∗ P R F \xrightarrow[*]{\mathbb{P}} R FP R或者 F ‾ P = R \overline{F}^{P}=R FP=R;

  • ⋆ \star 例1.6:
    考虑 F = x 2 y 2 + 2 x y + y 2 + 1 , P = { P 1 , P 2 } F=x^{2} y^{2}+2 x y+y^{2}+1, \mathbb{P}=\left\{P_{1}, P_{2}\right\} F=x2y2+2xy+y2+1,P={P1,P2};其中 P 1 = x y + 1 , P 2 = y P_{1}=x y+1,P_{2}=y P1=xy+1,P2=y,并取定纯字典项序.那么:

F → x 2 y 2 P 1 x y + y 2 + 1 → x y P 1 y 2 → y 2 P 2 0 F \xrightarrow[x^{2} y^{2}]{P_{1}} x y+y^{2}+1 \xrightarrow[x y]{P_{1}} y^{2} \xrightarrow[y^{2}]{P_{2}} 0 FP1 x2y2xy+y2+1P1 xyy2P2 y20

F F F P \mathbb{P} P的一个约化过程.另一个约化过程是:

F → x 2 y 2 P 1 x y + y 2 + 1 → x y P 2 y 2 + 1 → y 2 P 2 1 F \xrightarrow[x^{2} y^{2}]{P_{1}} x y+y^{2}+1 \xrightarrow[x y]{P_{2}} y^{2}+1 \xrightarrow[y^{2}]{P_{2}} 1 FP1 x2y2xy+y2+1P2 xyy2+1P2 y21

由此看出不同的约化过程,得到的范式可以不同.因为:

  • △ \triangle 定义1.7:
    F , G ∈ K [ x ] . F, G \in \mathcal{K}[\boldsymbol{x}] . F,GK[x]. 在给定项序 < r <_{r} <r 下,将 F , G F, G F,G 按降序排列如下:

F = f 1 x α 1 + ⋯ + f l x α l G = g 1 x β 1 + ⋯ + g m x β m \begin{array}{l} F=f_{1} x^{\alpha_{1}}+\cdots+f_{l} \boldsymbol{x}^{\alpha_{l}} \\ G=g_{1} \boldsymbol{x}^{\beta_{1}}+\cdots+g_{m} \boldsymbol{x}^{\beta_{m}} \end{array} F=f1xα1++flxαlG=g1xβ1++gmxβm

定义 K [ x ] \mathcal{K}[\boldsymbol{x}] K[x] 上的二元关系 ≺ r \prec_r r(称为由项序 < r <_r <r 诱导的 K [ x ] \mathcal{K}[\boldsymbol{x}] K[x] 上的多项式序) 如下:

F ≺ r G ⟺ F \prec_{r} G \Longleftrightarrow FrG i) ∃ k ( 1 ⩽ k ⩽ l ) \exists k(1 \leqslant k \leqslant l) k(1kl), 使得 α k < r β k \alpha_{k}<_{r} \beta_{k} αk<rβk, 而 α i = β i ( 1 ⩽ i < k ) \alpha_{i}=\beta_{i}(1 \leqslant i<k) αi=βi(1i<k)ii) l < m , l<m, l<m, α i = β i ( 1 ⩽ i ⩽ l ) \alpha_{i}=\beta_{i}(1 \leqslant i \leqslant l) αi=βi(1il);

  • ♡ \heartsuit 定理1.8:
    二元关系 ≺ r \prec_r r K [ x ] \mathcal{K}[\boldsymbol{x}] K[x] 上的良序.

  • ♡ \heartsuit 定理1.9:
    给定 T \mathcal{T} T上项序 <, 并设 ≺ \prec < < <诱导的多项式序,对任意 F , P ∈ K [ x ] F, P \in \mathcal{K}[\boldsymbol{x}] F,PK[x], 若 F ⟶ P G , F \stackrel{P}{\longrightarrow} G, FPG, G ≺ F G \prec F GF;

因为 ≺ \prec K [ x ] \mathcal{K}[\boldsymbol{x}] K[x] 上的良序,所以对任意的 F F F P P P, F F F P P P 的约化过程必然终止.在有限不约化过程后所得的多项式 R R R P P P是约化的.我们就得到如下的范式形式:

F = ∑ i = 1 s Q i P i + R F=\sum_{i=1}^{s} Q_{i} P_{i}+R F=i=1sQiPi+R

其中 P i ∈ P , Q i , R ∈ K [ x ] P_{i} \in \mathbb{P}, Q_{i}, R \in \mathcal{K}[\boldsymbol{x}] PiP,Qi,RK[x], 而 R R R P \mathbb{P} P 是约化的. 称多项式 R R R F F F P \mathbb{P} P范式余式, 记作 nform ⁡ ( F , P ) \operatorname{nform}(F, \mathbb{P}) nform(F,P). 又称由 F F F求得 R R R 的过程为 F F F P \mathbb{P} P约化. 对任意 Q ⊂ K [ x ] \mathbb{Q} \subset \mathcal{K}[\boldsymbol{x}] QK[x], 定义

nform ⁡ ( Q , P ) : = { nform ⁡ ( Q , P ) ∣ Q ∈ Q } \operatorname{nform}(\mathbb{Q}, \mathbb{P}):=\{\operatorname{nform}(Q, \mathbb{P}) \mid Q \in \mathbb{Q}\} nform(Q,P):={nform(Q,P)QQ}

  • ⋆ \star 例1.10:
    考虑下列多项式:

P 1 = x 1 x 4 + x 3 − x 1 x 2 P 2 = 2 x 4 2 − 2 x 3 x 4 + 5 x 1 x 2 x 4 − 5 x 1 x 2 x 3 F = x 1 x 4 2 + x 4 2 − x 1 x 2 x 4 − x 2 x 4 + x 1 x 2 + 3 x 2 \begin{aligned} P_{1} &=x_{1} x_{4}+x_{3}-x_{1} x_{2} \\ P_{2} &=2 x_{4}^{2}-2 x_{3} x_{4}+5 x_{1} x_{2} x_{4}-5 x_{1} x_{2} x_{3} \\ F &=x_{1} x_{4}^{2}+x_{4}^{2}-x_{1} x_{2} x_{4}-x_{2} x_{4}+x_{1} x_{2}+3 x_{2} \end{aligned} P1P2F=x1x4+x3x1x2=2x422x3x4+5x1x2x45x1x2x3=x1x42+x42x1x2x4x2x4+x1x2+3x2

它们中的单项式已按 x 4 ≻ ⋯ ≻ x 1 x_{4} \succ \cdots \succ x_{1} x4x1 的纯字典序排列.用符号表示,我们有:

lt ⁡ ( P 1 ) = x 1 x 4 , lt ⁡ ( P 2 ) = x 4 2 , lt ⁡ ( F ) = x 1 x 4 2 lc ⁡ ( P 1 ) = lc ⁡ ( F ) = 1 , lc ⁡ ( P 2 ) = 2 \begin{array}{l} \operatorname{lt}\left(P_{1}\right)=x_{1} x_{4}, \quad \operatorname{lt}\left(P_{2}\right)=x_{4}^{2}, \quad \operatorname{lt}(F)=x_{1} x_{4}^{2} \\ \operatorname{lc}\left(P_{1}\right)=\operatorname{lc}(F)=1, \quad \operatorname{lc}\left(P_{2}\right)=2 \end{array} lt(P1)=x1x4,lt(P2)=x42,lt(F)=x1x42lc(P1)=lc(F)=1,lc(P2)=2

P = { P 1 , P 2 } . F \mathbb{P}=\left\{P_{1}, P_{2}\right\} . F P={P1,P2}.F P \mathbb{P} P 明显可约,譬如:

F → x 1 x 2 x 4 P 1 H , F = b ⋅ λ ⋅ P 1 + H F \xrightarrow[x_{1} x_{2} x_{4}]{P_{1}} H, F=b \cdot \lambda \cdot P_{1}+H FP1 x1x2x4H,F=bλP1+H

其中

b = − 1 , λ = x 2 H = x 1 x 4 2 + x 4 2 − x 2 x 4 + x 2 x 3 − x 1 x 2 2 + x 1 x 2 + 3 x 2 \begin{array}{l} b=-1, \lambda=x_{2} \\ H=x_{1} x_{4}^{2}+x_{4}^{2}-x_{2} x_{4}+x_{2} x_{3}-x_{1} x_{2}^{2}+x_{1} x_{2}+3 x_{2} \end{array} b=1,λ=x2H=x1x42+x42x2x4+x2x3x1x22+x1x2+3x2

对于上面的约化,消去的项关于项序不是极大. 若选极大项,我们需要先约化 F F F 中的首单项式 x 1 x 4 2 . x_{1} x_{4}^{2} . x1x42. P , F \mathbb{P}, F P,F 到其范式的约化如下:

F → x 1 x 4 2 P 1 H 1 → x 4 2 P 2 H 2 → x 1 x 2 x 4 P 1 H 3 F \xrightarrow[x_{1} x_{4}^{2}]{P_{1}} H_{1} \xrightarrow[x_{4}^{2}]{P_{2}} H_{2} \xrightarrow[x_{1} x_{2} x_{4}]{P_{1}} H_{3} FP1 x1x42H1P2 x42H2P1 x1x2x4H3

其中:

H 1 = F − x 4 P 1 = x 4 2 − x 3 x 4 − x 2 x 4 + x 1 x 2 + 3 x 2 H 2 = H 1 − 1 2 P 2 = − 5 2 x 1 x 2 x 4 − x 2 x 4 + 5 2 x 1 x 2 x 3 + x 1 x 2 + 3 x 2 H 3 = H 2 + 5 2 x 2 P 1 = − x 2 x 4 + 5 2 x 1 x 2 x 3 + 5 2 x 2 x 3 − 5 2 x 1 x 2 2 + x 1 x 2 + 3 x 2 \begin{array}{l} H_{1}=F-x_{4} P_{1}=x_{4}^{2}-x_{3} x_{4}-x_{2} x_{4}+x_{1} x_{2}+3 x_{2} \\ H_{2}=H_{1}-\frac{1}{2} P_{2}=-\frac{5}{2} x_{1} x_{2} x_{4}-x_{2} x_{4}+\frac{5}{2} x_{1} x_{2} x_{3}+x_{1} x_{2}+3 x_{2} \\ H_{3}=H_{2}+\frac{5}{2} x_{2} P_{1}=-x_{2} x_{4}+\frac{5}{2} x_{1} x_{2} x_{3}+\frac{5}{2} x_{2} x_{3}-\frac{5}{2} x_{1} x_{2}^{2}+x_{1} x_{2}+3 x_{2} \end{array} H1=Fx4P1=x42x3x4x2x4+x1x2+3x2H2=H121P2=25x1x2x4x2x4+25x1x2x3+x1x2+3x2H3=H2+25x2P1=x2x4+25x1x2x3+25x2x325x1x22+x1x2+3x2

现在 H 3 H_{3} H3 P \mathbb{P} P 已约化,因而不可能有进一步的约化. 所以

R = nform ⁡ ( F , P ) = H 3 = F + ( 5 2 x 2 − x 4 ) P 1 − 1 2 P 2 R=\operatorname{nform}(F, \mathbb{P})=H_{3}=F+\left(\frac{5}{2} x_{2}-x_{4}\right) P_{1}-\frac{1}{2} P_{2} R=nform(F,P)=H3=F+(25x2x4)P121P2

(本节完)


参考资料

A Tutorial on Gröbner Bases With Applications in Signals and Systems
<<多项式代数>>王东明/牟晨琪(2011)ISBN: 9787040316988
Ideals, Varieties, and Algorithms (4th ed.) [Cox, Little & O’Shea 2015-06-14]
介绍一下Grobner基的概念和应用?谢谢[知乎问题回答:来自"匿名用户"]

  • 5
    点赞
  • 23
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值