QR

Berechnen der qr-Faktorisierung einer Matrix.

bool  QR(
  matrix&  Q,     // Matrix mit orthonormalen Spalten
  matrix&  R      // obere Dreiecksmatrix
   );

Parameter

Q

[out]  Eine Matrix mit orthonormalen Spalten. Bei mode = 'complete' ist das Ergebnis eine orthogonale/unitäre Matrix, je nachdem, ob a reel/komplex ist oder nicht. Die Determinante kann in diesem Fall entweder +/- 1 sein. Falls die Anzahl der Dimensionen im Eingabefeld größer als 2 ist, wird ein Stapel von Matrizen mit den oben genannten Eigenschaften zurückgegeben.

R

[out]  Obere Dreiecksmatrix.

Rückgabewert

Gibt bei Erfolg true zurück, andernfalls false.

Beispiel

//---  A*x = b
  matrix A = {{01}, {11}, {11}, {21}};
  Print("A \n"A);
  vector b = {1223};
  Print("b \n"b);
//--- A = Q*R
  matrix qr;
  A.QR(qr);
  Print("q \n"q);
  Print("r \n"r);
  matrix qr=q.MatMul(r);
  Print("qr \n"qr);
  /* 
  A
  [[0,1]
  [1,1]
  [1,1]
  [2,1]]
  b
  [1,2,2,3]
  q
  [[0.4082482904638631,-0.8164965809277259,-1.110223024625157e-16,-0.4082482904638631]
  [0.4625425214347352,-0.03745747856526496,0.7041241452319315,0.5374574785652647]
  [-0.5374574785652648,-0.03745747856526496,0.7041241452319316,-0.4625425214347352]
  [-0.5749149571305296,-0.5749149571305299,-0.09175170953613698,0.5749149571305296]]
  r
  [[-1.224744871391589,-0.2415816237971962]
  [-1.22474487139159,-1.466326495188786]
  [1.224744871391589,1.316496580927726]
  [1.224744871391589,0.2415816237971961]]
  qr
  [[-1.110223024625157e-16,1]
  [1,0.9999999999999999]
  [1,1]
  [2,1]]
  */