Page 68 - 363_

69

Приклад:
>> disp(det(A))
30
>> disp(trace(A))
3
>> disp(null(A))
>> disp(orth(A))
0.3395 0.4082 -0.8474
0.2793 0.8165 0.5053
0.8982 -0.4082 0.1632
>> disp(rref(A))
1 0 0
0 1 0
0 0 1
Функція R=chol(A) знаходить розклад Холецького для дійсних

симетричних та комплексних ермітових матриць. Наприклад:

>> A = [1 2 3; 2 15 8; 3 8 400]
A =
1 2 3
2 15 8
3 8 400
>> dlsp(chol(A))
1.0000 2.0000 3.0000
0 3.3166 0.6030
0 0 19.7645
Функція lu(А) здійснює LU-розкладання матриці А у вигляді добутку

нижньої трикутної матриці L (можливо, з переставленнями) і верхньої
трикутної матриці U так, що A = L * U.

Звернення до цієї функції вигляду [L, U, P] = lu(A) дозволяє отримати

три складові цього розкладу – нижню трикутну матрицю L, верхню трикутну U
і матрицю переставлень P такі, що P*A =L*U. Наведемо приклад:

A =
1 2 3
2 15 8
3 8 400
>> displu(A))
3.0000 8.0000 400.0000
-0.6667 9.6667 -258.6667
-0.3333 0.0690 -148.1724
>> [L,U,P]=lu(A);
>> L
L =
1.0000 0 0
0.6667 1.0000 0
0.3333 -0.0690 1.0000
>> U
U=
3.0000 8.0000 400.0000