colcompr
多項式行列の列圧縮
呼び出し手順
[Y, rk, ac] = colcompr(A)
パラメータ
- A
多項式行列
- Y
正方多項式行列 (右ユニモジュラ基底)
- rk
Aの通常ランク- Ac
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⁴
参照
- rowcompr — 多項式行列の行圧縮
| Report an issue | ||
| << coffg | Polynomials | degree >> |