Re: full QR decomposition
- To: mathgroup at smc.vnet.net
- Subject: [mg51687] Re: [mg51645] full QR decomposition
- From: Daniel Lichtblau <danl at wolfram.com>
- Date: Fri, 29 Oct 2004 03:39:13 -0400 (EDT)
- References: <200410280343.XAA09703@smc.vnet.net>
- Sender: owner-wri-mathgroup at wolfram.com
Jan Schmedes wrote:
> Dear group,
>
> i have the following problem. For a given a NxM (e.g. 3000x600) matrix A i
> get a matrix Q (MxN) and R (MxM) after the QRDecomposition. But i would
> like to get the full decomposition with a matrix QF (NxN) and RF (NxM) where RF
> contains only zero lines below the line with index M. I'm especially
> interrested in the last N-M lines in QF which create the zero lines in RF.
> At the moment i use QRDecomposition and append then N-M random lines to
> the matrix Q. With this new matrix i go into the routine GramSchmidt and extract the last N-M
> lines from the result. But this is very slow and i'm not sure if it will
> always work? Are there any suggestions or hints?
>
> best regards and thank you
>
> jan schmedes
I'll discuss this with respect to real valued matrices so as to avoid
requiring Conjugate in various places. The important defining relation
between a given real matrix mat and the decomposition into q and r is
that Transpose[q].r==mat. From this and your requirement for the full
decomposition it is not hard to see that extra rows in q are simply any
orthonormal set for the null space of q. To make the null space
orthogonal we simply use another invocation of QRDecomposition and
extract the "Q" part.
This may be obtained as in the example below.
SeedRandom[1111];
randomMatrix[n_,m_] := Table[Random[], {n}, {m}]
rows = 7;
cols = 3;
mat = randomMatrix[rows, cols];
{q,r} = QRDecomposition[mat];
First we check the QR decomposition.
In[64]:= Transpose[q] . r == mat
Out[64]= True
In[65]:= Chop[q . Transpose[q]] == IdentityMatrix[Length[q]]
Out[65]= True
Now we form the full decomposition.
qextra = NullSpace[q];
{q2,r2} = QRDecomposition[Transpose[qextra]];
qf = Join[q,q2];
rf = Join[r, Table[0, {rows-Length[r]}, {cols}]];
Now check the full decomposition:
In[66]:= Transpose[qf] . rf == mat
Out[66]= True
In[67]:= Chop[qf . Transpose[qf]] == IdentityMatrix[Length[qf]]
Out[67]= True
In[68]:= Dimensions[qf] == {rows,rows}
Out[68]= True
Daniel Lichtblau
Wolfram Research
- References:
- full QR decomposition
- From: Jan Schmedes <schmedes@rz.uni-potsdam.de>
- full QR decomposition