column compression of polynomial matrix
[Y, rk, ac] = colcompr(A);
polynomial matrix
square polynomial matrix (right unimodular basis)
normal rank of A
Ac=A*Y
, polynomial matrix
column compression of polynomial matrix A
(compression to the left)
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⁴