Path: news.mathworks.com!not-for-mail
From: <HIDDEN>
Newsgroups: comp.soft-sys.matlab
Subject: Re: Uniform random rotation matrix
Date: Sun, 1 Apr 2012 13:24:11 +0000 (UTC)
Organization: The MathWorks, Inc.
Lines: 15
Message-ID: <jl9ktr$f4c$1@newscl01ah.mathworks.com>
References: <idr87f$ho7$1@fred.mathworks.com> <idr8mi$q71$1@fred.mathworks.com> <idrfee$3iu$1@fred.mathworks.com> <idrg69$hrj$1@fred.mathworks.com> <idrgpg$so6$1@fred.mathworks.com> <idri4m$nuc$1@fred.mathworks.com> <idriop$5q2$1@fred.mathworks.com> <idrj41$buh$1@fred.mathworks.com> <idrk4e$t80$1@fred.mathworks.com>
Reply-To: <HIDDEN>
NNTP-Posting-Host: www-01-blr.mathworks.com
Content-Type: text/plain; charset=UTF-8; format=flowed
Content-Transfer-Encoding: 8bit
X-Trace: newscl01ah.mathworks.com 1333286651 15500 172.30.248.46 (1 Apr 2012 13:24:11 GMT)
X-Complaints-To: news@mathworks.com
NNTP-Posting-Date: Sun, 1 Apr 2012 13:24:11 +0000 (UTC)
X-Newsreader: MATLAB Central Newsreader 3414301
Xref: news.mathworks.com comp.soft-sys.matlab:762844

"Bruno Luong" <b.luong@fogale.findmycountry> wrote in message <idrk4e$t80$1@fred.mathworks.com>...
> Just wonder John, does the QR of your Matlab always output diag(R)>=0 ?
> 
> If not, a change should perhaps is nevessary:
> 
> [Q, R] = qr(randn(ndim))
> Q = Q*diag(sign(diag(R)))
> 
> Bruno

Hi Bruno,

How did you come up with this formula for uniform random rotation matrices? Is there a proof or some literature available?

Thanks!