Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Численные методы. Учебно-методический комплекс.pdf
Скачиваний:
0
Добавлен:
15.08.2026
Размер:
1 Мб
Скачать

Методические рекомендации

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

Пусть имеем характеристическое уравнение

det( A λE ) = 0.

Корни этого уравнения λ1 , λ2 ,…,λn являются собственными значениями матрицы A . Пусть им соответствуют линейно независимые собственные векторы x(1 ) , x( 2 ) ,…, x( n ) . Укажем некоторые итерационные методы вычисления наибольшего по модулю собственного значения матрицы A , не требующие раскрытия ее векового определителя.

Пусть среди собственных значений матрицы A есть одно, наибольшее по модулю. Для определенности предположим, что

| λ1 |>| λ2 || λ3 |≥K≥| λn | ,

(1)

так что наибольшим по модулю является первое собственное значение. Очевидно, что для действительной матрицы наибольшее по модулю собственное значение λ1 действительно. Заметим, что такой случай имеет место, если матрица A - действительная и элементы ее положительны.

Возьмем произвольный вектор y

Шаг 1. Последовательно

образуем итерации Ay ,

A2 y ,…, Am y

( m =1,2,3,K) .

 

 

 

Шаг 2. При достаточно больших m полагаем

 

 

 

( m+1 )

 

λ1

 

yi

( i =1,2,K,n ) ,

(2)

 

( m )

 

 

yi

 

где yi( m ) - i -я координата вектора y( m ) = Am y .

То есть, взяв достаточно большой номер итерации m, мы сможем с любой степенью точности определить по формуле (2) наибольший по модулю корень λ1 характеристического уравнения матрицы A . Для нахождения этого корня может быть использована любая координата вектора y( m ) . В частности, можно взять среднее арифметическое соответствующих отношений.

Теперь рассмотрим задачу определения собственного вектора, соответствующего только что найденному собственному значению λ1 .

КазНПУ 703-10-11 Учебно-методический комплекс дисциплины. Издание третье.

70

Заметим, что вектор y( m ) = Am y приближенно представляет собой собственный вектор матрицы A , соответствующий собственному значению λ1 . При достаточно большом m с любой степенью точности будем иметь:

Am y c1λ1m x( 1 ) ,

т.е. вектор Am y = y( m ) лишь числовым множителем отличается от собственного вектора x(1 ) и, следовательно, также является собственным вектором, соответствующим тому же самому собственному значению λ1 .

Пример. Найти наибольшее собственное значение матрицы

4

1

0

A = 1

2

1

0

1

1

и соответствующий ему собственный вектор.

1

Решение. Выбираем начальный вектор y = 1 .

1

4 1

0 1

5

Найдем первую итерацию, т.е. вектор y( 1 ) = Ay = 1

2

1 1

= 4 .

 

 

 

 

 

 

 

1

 

 

 

 

0

1 1

2

Найдем вторую итерацию,

т.е. вектор y( 2 ) = A2 y = Ay(1 )

4 1 0 5

24

= 1

2 1 4

= 15

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и т.д. y( m ) = Am y = Ay( m1 ) ( m =1,2,3,K) .

 

 

 

 

 

 

 

 

 

 

0 1 1 2

6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Составляем таблицу

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

Ay

 

A2 y

 

 

 

 

A3 y

 

A4 y

 

 

A5 y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

5

 

24

 

 

 

 

 

111

 

504

 

 

2268

 

 

 

 

 

 

 

1

 

4

 

15

 

 

 

 

 

60

 

 

252

 

 

1089

 

 

 

 

 

 

 

1

 

2

 

6

 

 

 

 

 

 

21

 

 

81

 

 

333

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A6 y

 

A7 y

 

 

 

A8 y

 

 

 

 

A9 y

 

 

A10 y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4 038939

 

 

 

 

 

 

10161

 

45433

 

202833

 

905238

 

 

 

 

 

 

 

4779

 

21141

 

93906

 

417987

 

1 862460

 

 

 

 

 

 

1422

 

6201

 

 

27342

 

121248

 

539235

 

 

 

 

 

 

Остановившись на итерациях A9 y = y( 9 ) и

A10 y = y( 10 ) ,

получаем зна-

чения

 

 

 

 

 

 

y(10 )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

4 038939

= 4,462 ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y( 9 )

905238

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(10 )

 

 

1862460

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y2

 

=

= 4,456 ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( 9 )

 

 

417987

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y2

 

 

 

 

 

 

 

 

 

 

 

 

 

КазНПУ 703-10-11 Учебно-методический комплекс дисциплины. Издание третье.

71

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]