1:
2:
3:
4:
5:
6:
7:
8:
9:
10:
11:
12:
13:
14:
15:
16:
17:
18:
19:
20:
21:
22:
23:
24:
25:
26:
27:
28:
29:
30:
31:
32:
33:
34:
35:
36:
37:
38:
39:
40:
41:
42:
43:
44:
45:
46:
47:
48:
49:
50:
51:
52:
53:
54:
55:
56:
57:
58:
59:
60:
61:
62:
63:
64:
65:
66:
67:
68:
69:
70:
71:
72:
73:
74:
75:
76:
77:
78:
79:
80:
81:
82:
83:
84:
85:
86:
87:
88:
89:
90:

>> QRmitH
Es wird eine QR-Zerlegung mit Householder berechnet.
====================================================
---------------------------
| H=I-ß*uu^T              |
| u=x+sign(x1)*||x||_2*e1 |
| ß=2/(u^Tu)              |
---------------------------
 
Matrix A eingeben: [3,-16/5;0,3;4,12/5]
A0 =
    3.0000   -3.2000
         0    3.0000
    4.0000    2.4000
 
Durchgang 1 
================
x =
     3
     0
     4
signx =
     1
e1 =
     1
     0
     0
u =
     8
     0
     4
beta =
    0.0250
 
H =
   -0.6000         0   -0.8000
         0    1.0000         0
   -0.8000         0    0.6000
Qk =
   -0.6000         0   -0.8000
         0    1.0000         0
   -0.8000         0    0.6000
Ak =
   -5.0000    0.0000
         0    3.0000
   -0.0000    4.0000
 
Durchgang 2 
================
x =
     3
     4
signx =
     1
e1 =
     1
     0
u =
     8
     4
beta =
    0.0250
 
H =
   -0.6000   -0.8000
   -0.8000    0.6000
Qk =
    1.0000         0         0
         0   -0.6000   -0.8000
         0   -0.8000    0.6000
Ak =
   -5.0000    0.0000
    0.0000   -5.0000
   -0.0000   -0.0000
 
 
Die QR-Zerlegung mit Householder
--------------------------------
A =
    3.0000   -3.2000
         0    3.0000
    4.0000    2.4000
Q =
   -0.6000    0.6400   -0.4800
         0   -0.6000   -0.8000
   -0.8000   -0.4800    0.3600
R =
   -5.0000    0.0000
    0.0000   -5.0000
   -0.0000   -0.0000