colcompr

Arguments

A

polynomial matrix

Y

square polynomial matrix (right unimodular basis)

rk

normal rank of A

Ac

Ac=A*Y, polynomial matrix

Description

column compression of polynomial matrix A (compression to the left)

Examples

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⁴
     

See also

• rowcompr — row compression of polynomial matrix