- •Итерационные методы решения слау
- •1. Общая характеристика итерационной схемы.
- •2. Каноническая форма одношаговых итерационных методов.
- •3. Метод простой итерации.
- •4. Метод Якоби.
- •5. Метод Зейделя.
- •Метод верхней релаксации.
- •Сходимость стационарных итерационных методов
- •Сходимость конкретных итерационных методов
3. Метод простой итерации.
Если ,то стационарный итерационный метод называется методом простой итерации и в матричном виде имеет вид
. (6)
В координатном виде он запишется в виде
.
Метод простой итерации является явной двухслойной схемой с постоянным параметром . Матрица перехода имеет вид .
4. Метод Якоби.
Запишем систему (1) в виде
. (7)
Если и задано - k-ое приближение к решению системы, то i-я компонента (k+1)-ого приближение вычисляется из i-го уравнения
. (8)
То есть
. . (10)
В выражении (10) слагаемое является i-ой компонентой столбца , где - диагональная матрица. Тогда (10) можно записать в виде
. (11)
Или в матричном виде
. (12)
Хотя формально эта схема неявная ( ), но т.к. D – диагональная матрица, то получаются явные формулы.
5. Метод Зейделя.
Если и задано - k-ое приближение к решению системы, то в итерационном методе Зейделя i-я компонента (k+1)-ого приближение вычисляется из i-го уравнения, но при этом используются уже вычисленные компоненты (k+1)-ого приближения.
Используются два вида формул
, (13)
. (14)
Из формул (13) значения вычисляются последовательно: , а из формул (14) вычисляются в обратном порядке.
Чтобы понять, как вычисляются значения запишем подробнее систему (13)
, (15)
, (16)
……………………………………
. (17)
Первая компонента вектора находится из уравнения (15) явным образом, для ее вычисления нужно знать вектор и значение . При вычислении из уравнения (16) используется только что вычисленное значение и известные значения , с предыдущей итерации. Таким образом, компоненты вектора вычисляются из уравнения (17) последовательно, начиная с i=1.
Запишем метод Зейделя в матричной форме. Для этого представим матрицу А в виде суммы трех матриц
,
где - диагональная матрица с той же главной диагональю, что и матрица А,
- нижняя треугольная (поддиагональная) матрица с нулями на главной диагонали ( ),
- верхняя треугольная (наддиагональная) матрица с нулями на главной диагонали ( ).
Тогда уравнения (13) можно записать в виде
, (18)
или
(19)
Здесь - треугольная матрица и, следовательно, находятся по явным формулам.
Метод верхней релаксации.
Обобщением метода Зейделя является метод верхней релаксации, который в матричном виде имеет вид
(20)
где - заданный числовой параметр.
Для получения расчетных формул перепишем (20) в виде
. (21)
В покомпонентной форме получим
.
Отсюда последовательно, начиная с i=1, находим все
, (22)
, (23)
………………………………………………………….
. (24)