colcomp — compression de colonnes, noyau
[W,rk]=colcomp(A [,flag] [,tol])
matrice réelle ou complexe
chaîne de caractères
nombre réel
matrice carré régulière (matrice de changement de base)
entier (rang de"A)
Compression des colonnes de A : Ac = A*W est à colonnes compressées, c'est à dire
Ac=[0,Af] et Af est de rang maximal
rank(Af) = rank(A) = rk.
flag et tol sont des paramètres optionnels : flag = 'qr'
ou 'svd' ('svd' par défaut).
tol = paramètre de tolérance (de l'ordre de %eps par défaut).
Les ma-rk premières colonnes de W forment une base du noyau de A quand size(A)=[na,ma].