'economy size' qr ?

[Pascal A. Dupuis]

Hello,

I would like to solve a overdetermined set of equation Ax=B, where
A is C m x n
x is R n x 1
B is C m x 1

In order for x to be real, the problem is augmented as:

[ A; conj(A) ] x = [b; conj(b)]

The first step is to compute
[qq, qr] =qr([A , b; conj(A) conj(b)]);

qr is a n+1 x n+1 upper real triangular matrix

then the resolution is performed using various approaches on qr.

Given that
1) qq is not needed (2m x 2m complex matrix !)
2) although A and b are complex, qr is real

are there simpler approaches in order to obtain qr, skipping the qq
computation ?

Thanks in advance

Pascal 
 
-- 
Dr. ir. Pascal Dupuis
K. U. Leuven, ESAT/ELECTA (formerly ELEN):  http://www.esat.kuleuven.ac.be/
Kasteelpark Arenberg, 10; B-3001 Leuven-Heverlee, Belgium
Tel. +32-16-32 10 21 -- Fax +32-16-32 19 85



-------------------------------------------------------------
Octave is freely available under the terms of the GNU GPL.

Octave's home on the web:  http://www.octave.org
How to fund new projects:  http://www.octave.org/funding.html
Subscription information:  http://www.octave.org/archive.html
-------------------------------------------------------------