Вопрос 19. Метод Якоби решения систем линейных уравнений. Общая характеристика методов, их сходимость и вычислительная  устойчивость.

 

 

 

Пусть требуется решить систему линейных алгебраических уравнений

,                                               (1)

где -квадратная невырожденная матрица.

 

Перепишем матричное уравнение (1) в виде

,

где, например, матрица C = EA   (в этом случае ).

 

В результате может быть организован итерационный процесс, в ходе которого, начиная с некоторого вектора T , строится последовательность  в соответствии с рекуррентной формулой:

,   0, 1, 2, ...         .                 

 

  Если при этом существует предел , где  решение системы (1), то итерационный процесс считают сходящимся.

 

В отличие от прямых методов решение системы  с помощью итерационных методов можно получить лишь с некоторой погрешностью, которая обусловлена критерием остановки расчетов. При этом обычно используют  норму вектора невязки  или норму разности , т.е., если требуется получить решение с точностью >0, то расчеты выполняются до тех пор, пока не будет выполнено одно из условий:

       или             .

 

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

,

где параметр  может выбираться на каждом шаге либо быть константой.

 

 

 

 

Метод Якоби

 

Дана система линейных уравнений

 ,                                                             (1) 

где, как и раньше,   – квадратная невырожденная матрица. Матрицу  всегда можно представить в виде ,   где

;; .

Тогда матричное уравнение (1) можно переписать в виде:

.                                      (1a)

Уравнение (1a)  преобразуем в реккурентную формулу метода Якоби:

,    k = 0,1,2,….                           (2)

Вычитая из обеих частей (2)  получим представление (2) в виде:

,       k = 0,1,2,….                            (2a)

Обозначив   и    получим следующие расчетные формулы метода Якоби:

,     k = 0,1,2,….                            (3)

Вектор  можно взять любой, например, вектор  или нулевой вектор.

Из (2) можно получить расчетные формулы для  в координатной форме:

;     1, ... n,   0, 1, ... .(4)

В качестве критерия остановки расчетов обычно используют .

 

В матричном виде описание одного шага метода следующее:

1) вычисляем , если ,  то вычисление закончены;

2) находим       ;

3) вычисляем .

Итерационный процесс  сходится, если

                                                                  (5).

 

Для симметричной положительно определенной матрицы для сходимости метода достаточно выполнения условия  ,      1,..., n,

т.е. условия преобладания главной диагонали матрицы .

 

В случае, когда (5) не выполнено, необходимо преобразовать систему (1) таким образом, что бы условие сходимости было выполнено.

 

Задача

 

Решить систему уравнений:

при помощи метода Якоби с точностью .

Решение.

;      .

.

Проверим условие сходимости: 0,756<1. Условие выполнено. Пусть (0; 0; 0)T. Следующие расчеты приведены в таблице.

K

0

3,919

(0,866; 0,775; 0,298)T

(0,866; 0,775; 0,298)T

1

2,449

(–0,242; –0,375; –0,571)T

(0,624; 0,400; –0,273)T

2

1,367

(0,406; 0,143; 0,251)T

(1,030; 0,543; –0,022)T

....

.....

.........

.........

13

0,011

(–0,002; –0,001; –0,002)T

(0,892; 0,418; –0,072)T

14

0,007

 

 

 

Так как , вычисления прекращаем.

Ответ:        x1 = 0,89; x2 = 0,42; x3 = -0,07. 

Примечание. Если использовать модификацию метода Якоби с регулировкой шага: , где 0,8 (получено подбором), то этот же результат будет получен на 9-ой итерации.