Численные методы линейной алгебры Учебное пособие. Ширапов Д.Ш. - 29 стр.

UptoLike

Составители: 

Рубрика: 

Доказательство
Доказательство следует из того, что если
1S < , то и µ
p
<1. Следовательно, по теореме
3.4 стационарный итерационный процесс сходится.
3.5. Вариационно-итерационные методы
В параграфе 3.4 рассматривались такие итерационные методы решения СЛАУ
Ax=b, (3.40)
в которых для задания итерационных параметров требовалось знание границы собственных
значений матрицы А. Что делать, если такой информации нет? В таком случае можно при-
менять методы, не использующие знания о границах собственных значений матрицы А.
Ниже описываются итерационные методы вида
В
1k
)k()1k(
xx
+
+
τ
+ Ax
(k)
=b, k=0, 1, 2,…
в которых параметры
τ
k+1
выбираются из других условий.
3.5.1. Метод минимальных невязок
Рассмотрим явную схему метода минимальных невязок [8, 9] заданную формулой с сим-
метричной положительно определенной матрицей А.
1k
)k()1k(
xx
+
+
τ
+ Ax
(k)
=b, k=0, 1, 2,… (3.41)
Для невязки R
(k)
= Ax
(k)
-b из (3.41) получим
1k
)k()1k(
RR
+
+
τ
+ AR
(k)
=0, k=0, 1, 2,…
R
(k+1)
=R
(k)
-τ
k+1
AR
(k)
.
Параметр
τ
k+1
выбираем из условия минимума невязки R
(k+1)
по норме
+τ=τ=
++
+
)AR,R(2RARRR
)k()k(
1k
2
)k(
2
)k(
1k
)k(
2
)1k(
2
)k(2
1k
AR
+
τ+
=ϕ(τ
k+1
).
Дальше вычислим производную от функции
ϕ(τ
k+1
) по параметру τ
k+1
и приравняем нулю
ϕ
(τ
k+1
)= -2(R
(k)
,AR
(k)
)+2τ
k+1
2
)k(
AR
=0
или
2
)k(
)k()k(
1k
AR
)R,AR(
=τ
+
, k=0, 1, 2,… (3.42)
При этом значении параметра
τ
k+1
вторая производная ϕ
(τ
k+1
)>0, следовательно, достига-
ется
2
1)(k
τ
Rmin
1k
+
+
.
В этом методе переход от k-ой итерации к (k+1)-ой осуществляется следующим обра-
зом. По найденному значению x
(k)
вычисляется вектор невязки R
(k)
= Ax
(k)
-b и по формуле
(3.42) находится параметр
τ
k+1
. Затем по формуле (3.41) вычисляется вектор x
(k+1)
.
Метод минимальных невязок (3.41), (3.42) сходится с той же скоростью, что и метод
простой итерации с оптимальным параметром
τ .