colcompr

Paramètres

A

matrice de polynômes

Y

matrice carrée de polynômes (base unimodulaire à droite)

rk

rang normal de A

Ac

: Ac=A*Y, matrice de polynômes

Description

Compression de colonnes de la matrice de polynômes A (compression à gauche)

Exemples

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⁴
     

Voir aussi

• rowcompr — compression de lignes d'une matrice de polynômes