1.2.3 Матричная запись методов Якоби и Зейделя

Исходную матрицу системы (1.11) представим в виде суммы трёх матриц

A=A₁+D+A₂,

где D - диагональная матрица;

D =diаg[а₁₁а₂₂…а_mm];

A₁ - нижняя треугольная матрица;

A₂ - верхняя треугольная матрица.

Пример: Дана матрица размерности (33):

.

А₁ А₂ D

Тогда исходную систему (1.11) можно записать в виде

x=-D^-1A₁ x – D^-1A₂ x+D^-1 b.

Тогда метод Якоби можно записать в виде:

или

. (1.17)

В матричной форме метод Зейделя будет выглядеть:

или

. (1.18)

Преобразуем формулы (1.17) и (1.18):

, (1.19)

. (1.20)

Из (1.19) и (1.20) видно, что если итерационный метод сходится, то он сходится к точному решению. Иногда при решении задач большой размерности, в итерационные методы вводятся числовые параметры, которые могут зависеть от номера итерации.

Пример для метода Якоби.

,

где t – числовой параметр.

Возникают вопросы:

При каких значениях t сходимость будет наиболее быстрой?

При каких значениях t метод сходится?

На примере двух методов просматривается вывод о том, что одни и те же методы можно записывать несколькими способами. Поэтому вводят каноническую (стандартную) форму записи:

. (1.21)

Формула (1.21) получена путем объединения (1.19) и (1.20).

Матрица D_n+1 здесь задает тот или иной метод. Если существует обратная матрица к этой матрице, то из последней системы мы можем найти все неизвестные.

1. Метод (1.21) – явный, если матрица D_n совпадает с единичной матрицей и неявный - в противном случае.

2. Метод (1.21) – стационарный, если матрица D_n+1 = D, и параметр t не зависит от номера итерации и нестационарный - в противном случае.