Метод А.Н.Крылова для решения системы линейных алгебраических уравнений

Страницы работы

Содержание работы

32.Метод А. Н. Крылова

В начале тридцатых годов нашего столетия А. Н. Крыловым был предложен достаточно удобный, метод нахождения собственных значений и собственных векторов матриц. Пусть D(λ) ≡det(λ E - A) = λn + p1 λn-1 + p2 λn-2 + …., +  рп  (5.1)

— характеристический полином (с точностью до знака) матрицы А. Согласно тождеству Гамильтона — Кели, матрица А обращает в нуль свой характеристический полином; поэтому An  + p1 An-1 +…+ pn E = 0. (5.2) Возьмем теперь произвольный ненулевойвектор    y(0) =Умножая обе части равенства (5.2) справа на  y(0) получим:

An  y(0)  + p1 An-1 y(0)  +…+ pn y(0)   = 0. (5.3)Положим: Ak  y(0) = y(k) (k=1, 2, ..., п); (5.4) тогда равенство (5.3) приобретает вид    y(n) + p1 y(n-1) +…+ pn y(0)  = 0. (5.5) или   (5.5*)гдеy(k) =, (k=1, 2, ..., п).

Следовательно,   векторное   равенство   (5.5*)  эквивалентно   системе уравнений  p1 yj (n-1) + p2 yj (n-2) + …., +  рп yj(0) =  -yj(n)(j = 1. 2. . . ., n),(5.6) из   которой  можно  определить неизвестные коэф­фициенты p1, p2,....., рп . Так как на основании формулы (5.4)  y(k) =A y(k -1),    где (k= 1, 2, . . ., п),  то   координаты   y1 (к) , y2 (к) , …, yn (к)   вектора    y(k)последовательно вычисляются по формулам

  ,    ,…,   (5.7)     Таким образом, определение коэффициентов   pjхарактеристи­ческого полинома (5.1) методом А. Н. Крылова сводится к решению линейной системы уравнений (5.6), коэффициенты которой вычи­сляются   по   формулам   (5.7),   причем   координаты   начального вектора

y(0) =произвольны. Если система (6) имеет единственное решение, то ее корни р1, р2, ..., рп являются коэффициентами характеристического полинома (5.1). Это решение может быть найдено, например, методом Гаусса. Если система (5.6) не имеет единственного решения, то задача усложняется . В этом случае рекомендуется изменить начальный вектор.

Похожие материалы

Информация о работе