- •СИЛЛАБУС
- •1. Информация о дисциплине
- •2. Краткое описание дисциплины.
- •3. Пререквизиты дисциплины
- •4. Постреквизиты дисциплины
- •7. Критерии оценки
- •8. Требования преподавателя.
- •Вопросы
- •Абсолютная и относительная погрешности
- •Значащая цифра. Число верных знаков
- •Тема. Интерполирование функций.
- •1. Постановка задачи.
- •2. Интерполяционная формула Лагранжа.
- •3. Первая и вторая интерполяционные формулы Ньютона.
- •Вопросы
- •Тема. Методы решения СЛАУ.
- •1. Нормы векторов и матриц.
- •4.1. Случай нормальной системы.
- •5. Итерационные методы вариационного типа.
- •5.1. Метод наискорейшего спуска.
- •5.2. Метод минимальных невязок.
- •7. Метод прогонки.
- •Вопросы
- •Постановка задачи.
- •1. Развертывание вековых определителей.
- •2. Метод А. М. Данилевского.
- •3. Метод А. Н. Крылова.
- •Вопросы
- •Тема. Методы решения нелинейных уравнений.
- •1. Метод хорд.
- •2. Метод Ньютона (метод касательных).
- •3. Метод итерации.
- •4. Метод итераций решения систем нелинейных уравнений.
- •Вопросы
- •Тема. Численное интегрирование.
- •1. Постановка задачи.
- •2. Формула средних прямоугольников.
- •3. Формула трапеций.
- •4. Формула Симпсона.
- •5. Квадратурные формулы Ньютона-Котеса.
- •Вопросы
- •1. Постановки задач для ОДУ.
- •2. Задача Коши.
- •3. Метод Эйлера. Усовершенствованный метод ломаных.
- •7. Многошаговые методы.
- •7.1. Метод Адамса.
- •7.2. Метод Милна.
- •Вопросы
- •а) Разностные схемы для уравнения колебания.
- •б) Трехслойные схемы для уравнения теплопроводности.
- •Вопросы
- •1. Исходная задача.
- •3. Разностные схемы.
- •Вопросы
- •2. Разностная схема «крест».
- •3. Погрешность аппроксимации.
- •Вопросы
- •Задания для СРСП
- •4. Тематика и краткое описание лабораторных и студийных работ
- •Методические рекомендации
- •Тема. Интерполяционная формул Ньютона
- •Методические рекомендации
- •Тема. Методы решения СЛАУ. Метод итерации.
- •Методические рекомендации
- •Тема. Метод Зейделя.
- •Методические рекомендации
- •Методические рекомендации
- •Тема. Методы решения нелинейных уравнений.
- •Методические рекомендации
- •Тема. Численное интегрирование.
- •Методические рекомендации
- •Методические рекомендации
- •ГЛОССАРИЙ
Методические рекомендации
Изложим метод, позволяющий вычислять наибольшее по модулю собственное значение матрицы и соответствующий собственный вектор при помощи вычисления последовательности итераций произвольного вектора. Излагаемый метод называется степенным и является простейшим итерационным процессом для решения частичной проблемы собственных значений.
Пусть имеем характеристическое уравнение
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( m−1 ) ( 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
