orth - calcul d'une base orthogonale
Q=orth(A) renvoie Q , une base orthogonale de l'image de A . Im( Q ) = Im( A ) et Q'*Q = I .
Le nombre de colonnes de Q est égal au rang de A , comme déterminé par l'algorithme QR.
A=rand(5,3)*rand(3,4); [X,dim]=rowcomp(A);X=X'; svd([orth(A),X(:,1:dim)])