多項式行列の列圧縮
[Y, rk, ac] = colcompr(A)
多項式行列
正方多項式行列 (右ユニモジュラ基底)
A
の通常ランク
Ac=A*Y
, 多項式行列
多項式行列 A
の列圧縮 (左への圧縮)
s = %s; p = [s ; s*(s+1)^2 ; 2*s^2+s^3] [Y, rk, ac] = colcompr(p*p'); Y clean(ac) | ![]() | ![]() |
--> p = [s ; s*(s+1)^2 ; 2*s^2+s^3] p = s s +2s² +s³ 2s² +s³ --> [Y, rk, ac] = colcompr(p*p'); --> Y Y = 1 +2s +s² 2s +s² 1 -1 0 0 0 -1 0 --> clean(ac) ans = 0 0 s² 0 0 s² +2s³ +s⁴ 0 0 2s³ +s⁴