¿Cómo se calcula una base ortogonal de una matriz usando hmatrix?

Question

orth(somematrix) está integrado en MATLAB pero no parece estar disponible en la Biblioteca Haskell Hmatrix.

Answer 1

import Numeric.LinearAlgebra    

orth :: Field a => Matrix a -> [Vector a]
orth m = toColumns $ fst $ qr m

o sin punta

orth = toColumns . fst . qr

Wikipedia tiene una explicación.

Answer 2

Quizás esto es lo que necesitas:

orth m = toColumns u
        where (u,_,_) = compactSVD m

https://github.com/albertoruiz/hmatrix/issues/10#issuecomment-4077403