Mercurial > hg > ishara

view general/numerical/matrix/orthogonalise.m @ 19:1eb0ea29ec40

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Doesn't belong here.
author samer
date Thu, 17 Jan 2013 13:32:19 +0000
parents db7f4afd27c5
children
line wrap: on
line source
function B1=orthogonalise(B),
% orthogonalise - Orthogonalise a basis matrix
%
% orthogonalise :: [[N,M]] -> [[N,M]].
%
% Works using SVD.

[U,S,V] = svd(B,0);
B1 = U*V';

% alternative method, seems to be slower
% B1 = B*real((B'*B)^(-0.5));