49
Ví dụ
>> M = [2 3 4;3 5 6;4 6 9];
>> L=chol(M)
L =
1.4142 2.1213 2.8284
0 0.7071 0.0000
0 0 1.0000
>> L'*L-M
ans =
1.0e-015 *
0.4441 0 0
0 0 0
0 0 0
Phân tích QR
Phân tích QR là phép phân tích một ma trận vuông hoặc chữ nhật thành tích của
ma trận trực giao và một ma trận tam giác trên:
=
A
QR
với
T
=
Q Q
E ,
R
là ma trận tam giác trên. Ví dụ
>> A = [2 3 4 5; 3 5 6 8; 4 6 9 10];
>> [Q, R] = qr(A)
Q =
-0.3714 0.2491 -0.8944
-0.5571 -0.8305 -0.0000
-0.7428 0.4983 0.4472
R =
-5.3852 -8.3563 -11.5131 -13.7415
0 -0.4152 0.4983 -0.4152
0 0 0.4472 -0.0000
>> Q*R - A
ans =
1.0e-014 *
0 -0.1332 -0.1776 0
0.0444 -0.0888 -0.0888 0
0 -0.0888 0 0.1776