On Monday, March 31, 2014 8:56:45 PM UTC-5, Gustavo Lacerda wrote: > > Do we have something like Matlab's 'orth'? i.e. given a matrix M, find an > orthonormal > basis of the column-space of M? >
If M has full column rank the simplest method is orth(M) = qr(M)[1] or, equivalently, qrfact(M)[:Q] although this returns an object of type QRCompactWYQ. If M can be singular then you are probably better off using svdfact and taking the columns of the [:U] matrix corresponding to non-zero values in the [:S] vector.
