BÀI GIẢNG MATLAB - Trang 56

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