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

UptoLike

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

Рубрика: 

Рассмотрим два варианта существования погрешностей:
1)
δх=0, δλ≠0;
2)
δх0, δλ=0.
Последовательно рассмотрим оба варианта.
Вариант 1. δх=0, δλ≠0.
Тогда из (4.7) имеем
δАх
i
=δλ
i
х
i
,
(y
i
, δАх
i
)=δλ
i
(y
i
, x
i
),
A
)x,y(
xy
ii
ii
i
δδλ .
Отношение
1
cos
1
)x,y(
xy
i
iii
ii
ω=
α
=
называется i-ым коэффициентом перекоса матрицы А,
где
α
i
угол между собственными векторами x
i
и y
i
.
Таким образом,
A
ii
δωδλ .
Следовательно, если погрешность
δА мала и мал i-ый коэффициент перекоса, то мала по-
грешность определяемого i-го собственного значения.
Отметим, что для симметричной матрицы все
ω
i
=1, поэтому задача нахождения собст-
венных значений симметричной матрицы является устойчивой.
Вариант 2. δх0, δλ=0.
Тогда из (4.7) имеем
А
δх
i
+δАх
i
=λ
i
δх
i
,
(А
δх
i
, y
j
)+(δАх
i
, y
j
)= λ
i
(δх
i
, y
j
) , (4.8)
где (А
δх
i
, y
j
)=( δх
i
, A
*
y
j
)=(δх
i
, λ
j
y
j
)=λ
j
(δх
i
, y
j
).
Подставляя последнее в (4.8) получим
λ
j
(δх
i
, y
j
)+(δАх
i
, y
j
)= λ
i
(δх
i
, y
j
),
(
δх
i
, y
j
)=(δАх
i
, y
j
)/(λ
i
-λ
j
), при ij. (4.9)
Пусть
δх
i
=
=
β
n
1j
jij
x
, (4.10)
тогда
(
δх
i
, y
j
)=β
ij
(x
j
, y
j
)
после подстановки этого выражения в (4.9) получим
β
ij
=
))(y,x(
)y,Ax(
jijj
ji
λλ
δ
при ij,
A
x
x
ji
j
j
i
ij
δ
λλ
ω
β
. (4.11)
Из (4.10), (4.11) видно, что если мала погрешность определения элементов матрицы А и ма-
лы все коэффициенты перекоса, то мала погрешность определения i-го собственного векто-
ра, соответствующего i-му собственному значению.
Следствие 4.1. Если матрица А=A
*
и все собственные значения простые, то задача нахож-
дения собственных векторов устойчива.