2.3. Решение слау методом квадратных корней.
Метод квадратных корней используется для решения линейной системы вида АХ = В, у которой матрица А симметрическая, т.е. для элементов выполняется условие aij = aji (i, j = 1, 2, …, n).
В нашем случае:
Условие симметричности не выполняется, т.к.
a13 = 5, a31 = 1, a13 a31 ,
a23 = - 4, a32 = - 2, a23 a32 .
Следовательно, метод квадратных корней в данном случае не может быть применен.
2.4. Решение слау методом обратных матриц.
Систему уравнений можно представить в матричном виде:
А ∙ X = В,
где
А-1
-
обратная матрица к квадратной матрице
А =
Х = А-1 ∙ В - решение системы уравнений.
Формула для вычисления обратной матрицы
А--1=
,
где
-
определитель матрицы А.
Имеем:
, , .
Т.к. определитель матрицы А не равен 0, то матрица А невырожденная и для нее существует обратная матрица A-1.
Аij
называется
алгебраическим дополнением к элементу
и равно
Aij
= (-1)1+jМ
Минор
элемента
aij
– это определитель, полученный из
данного определителя вычеркиванием
i-ой
строки и j-гo
столбца.
|
|
Отсюда
А-1
=
Согласно формуле Х = А-1 ∙ В, получаем:
Получили: x1 = 3 , х2 = 2 , х3 = 0.
2.5. Решение слау методом релаксаций.
Имеем:
, , .
Основная вычислительная формула метода релаксаций имеет вид:
Преобразуем данное выражение:
Исходная матрица имеет размер 3*3, тогда получим:
Произведем расчеты
значений
при параметре
и ограничив число итераций с точностью
0,1 с контролем по параметру невязки:
На
начальной (нулевой) итерации принимаем
первичное приближение
.
1 итерация:
2 итерация:
3 итерация:
4 итерация:
5 итерация:
6 итерация:
7 итерация:
8 итерация:
9 итерация:
10 итерация:
11 итерация:
12 итерация:
13 итерация:
14 итерация:
15 итерация:
16 итерация:
17 итерация:
18 итерация:
19 итерация:
20 итерация:
Занесем результаты в таблицу:
№ итерации |
х1 |
х2 |
х3 |
1 |
2 |
3 |
макс |
1 |
0 |
0 |
0 |
5 |
12 |
1 |
12 |
2 |
0,833 |
1,722 |
0,268 |
4,605 |
6,24 |
0,807 |
6,24 |
3 |
1,601 |
2,506 |
0,536 |
2,529 |
3,424 |
0,803 |
3,424 |
4 |
2,023 |
2,936 |
0,743 |
1,088 |
2,118 |
0,62 |
2,118 |
5 |
2,204 |
3,29 |
0,914 |
0,398 |
1,378 |
0,634 |
1,378 |
6 |
2,27 |
3,498 |
1,078 |
0,204 |
1,278 |
0,492 |
1,278 |
7 |
2,236 |
3,722 |
1,24 |
0,464 |
1,322 |
0,488 |
1,322 |
8 |
2,159 |
3,968 |
1,416 |
0,621 |
1,442 |
0,529 |
1,442 |
9 |
2,056 |
4,243 |
1,613 |
0,747 |
1,611 |
0,591 |
1,611 |
10 |
1,932 |
4,553 |
1,836 |
0,87 |
1,821 |
0,666 |
1,821 |
11 |
1,787 |
4,905 |
2,088 |
0,991 |
2,063 |
0,759 |
2,063 |
12 |
1,622 |
5,304 |
2,375 |
1,133 |
2,344 |
0,861 |
2,344 |
13 |
1,433 |
5,758 |
2,701 |
1,288 |
2,664 |
0,98 |
2,664 |
14 |
1,218 |
6,274 |
3,072 |
1,466 |
3,03 |
1,114 |
3,03 |
15 |
0,974 |
6,86 |
3,494 |
1,672 |
3,448 |
1,264 |
3,448 |
16 |
0,695 |
7,528 |
3,974 |
1,899 |
3,922 |
1,439 |
3,922 |
17 |
0,378 |
8,287 |
4,52 |
2,16 |
4,463 |
1,636 |
4,463 |
18 |
0,018 |
9,151 |
5,141 |
2,457 |
5,075 |
1,861 |
5,075 |
19 |
-0,392 |
10,133 |
5,847 |
2,793 |
5,773 |
2,117 |
5,773 |
20 |
-0,858 |
11,25 |
6,65 |
3,176 |
6,566 |
2,408 |
6,566 |
21 |
-1,387 |
12,521 |
7,563 |
3,612 |
7,463 |
2,74 |
7,463 |
Таким образом,
видно, что метод релаксаций при
для данной системы не применим, происходит
отдаление от
.
