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

Итерационные методы решения СЛАУ

 

Канонический вид итерационных методов решения СЛАУ.

Итерационные методы Якоби и Зейделя

Рассмотрим СЛАУ

, (6.1)

где матрица , , имеет обратную, , .

Рассмотрим сначала два примера итерационных методов. Для их построения предварительно преобразуем систему (6.1) к виду

, (6.2)

(при этом предполагается, что все отличны от нуля).

Условимся, как обычно, считать значение суммы равным нулю, если верхний предел суммирования меньше нижнего. Так, уравнение (6.2) при имеет вид

.

В дальнейшем верхний индекс будет указывать номер итерации.

В методе Якоби исходят из записи системы в виде (6.2), причем итерации определяются следующим образом:

. (6.3)

Начальные значения , , задаются произвольно. Окончание итераций определяется либо заданием максимального числа итераций , либо условием

,

где заданное число.

Итерационный метод Зейделя имеет вид

. (6.4)

Чтобы понять, как находятся отсюда значения , , запишем подробнее первые два уравнения системы (6.4).

, (6.5)

. (6.6)

Первая компонента вектора находится из уравнения (6.5) явным образом, для ее вычисления нужно знать вектор и значение . При нахождении из уравнения (6.6) используются только что найденное значение и известные значения , , с предыдущей итерации. Таким образом, компоненты вектора находятся из уравнения (6.4) последовательно, начиная с .

 

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

Для исследования сходимости итерационных методов удобнее записывать их не в координатной, а в матричной форме. Представим матрицу системы (6.1) в виде суммы матриц

, (6.7)

где диагональная матрица с той же главной диагональю, что и матрица . Матрица нижняя треугольная матрица и верхняя треугольная матрица с нулевыми главными диагоналями.

Например, при матрицы имеют вид

, , .

Представление системы (6.1) в форме (6.2) эквивалентно ее записи в виде матричного уравнения (доказать)

.

Отсюда видно, что метод Якоби (6.3) в векторной записи выглядит следующим образом:

,

или, что то же самое,

. (6.8)

Метод Зейделя (6.4) записывается в виде

,

или

. (6.9)

Учитывая (6.7), методы (6.8) и (6.9) можно переписать соответственно в виде

, (6.10)

. (6.11)

Из этой записи видно, что если итерационный метод сходится, то он сходится к решению исходной системы уравнений.

Часто для ускорения сходимости в итерационные методы вводятся числовые параметры, которые зависят, вообще говоря, от номера итерации. Например, в методы (6.10), (6.11) можно ввести итерационные параметры следующим образом:

,

.

Способ выбора итерационных параметров выясняется при исследовании сходимости. В теории итерационных методов существует два круга вопросов: а) при каких значениях параметров метод сходится, б) при каких значениях параметров сходимость будет более быстрой (соответствующие параметры называются оптимальными).

Методы Якоби и Зейделя относятся к одношаговым итерационным методам, когда для нахождения требуется помнить только одну предыдущую итерацию . Иногда используются и многошаговые итерационные методы, в которых определяется через значения на двух и более предыдущих итерациях, т.е. .