, (5.2.1)
где . По свойству собственных векторов
… (5.2.2)
Тогда
… … … … … … (5.2.3)
Определим некоторую функцию
(5.2.4)
Составляя линейную комбинацию из y(i) с теми же коэффициентами, получим:
Если предположить, что
(5.2.5)
(остаток от деления характеристического многочлена на один из корней), то очевидно, что и .
Тогда от предыдущего утверждения остается, следовательно
(5.2.6)
Итак, если Ci ≠ 0, то данная комбинация задает собственный вектор x(i) с точностью до числового множителя. Коэффициенты qij могут быть легко определены по схеме Горнера
(5.2.7)
где собственные числа, а коэффициенты характеристического многочлена.
Пример:Найдем собственный вектор, соответствующий = -2 (для предыдущего примера):
Функция значит .
Выражение для собственного вектора
.
Так как собственный вектор определяется с точностью до постоянного множителя, то можно считать
5.3. Итерационный метод одновременного нахождения собственных значений λi и собственных векторов матрицы.
Известно, что если А – симметрична и положительноопределенна (если все расположенные по главной диагонали миноры больше 0), то
1) все ;
2) собственные векторы и ортогональны:
(5.3.1)
Запишем систему, служащую для определения собственного вектора , соответствующего λ1.
(5.3.2)
Или
(5.3.3)
Так как координаты вектора определяются с точностью до константы, то одну из них (например ) можно положить равной 1. И последнее уравнение в системе (5.3.3) примет вид
. (5.3.4)
Данную систему можно представить как итерационный процесс:
(5.3.5)
для
Пример: Найдём собственные числа и векторы матрицы при помощи итерационного процесса.
Покажем положительную определенность матрицы:
.
Составим систему для нахождения λi и .
Положим = 1. Следовательно
Для λ1 за нулевое приближение выбираем вектор , затем выполняем решение по схеме. .
Для λ2 систему можно упростить. Из относительной ортогональности следует, что
Положим , и подставим в итерационную схему.
Третий вектор находится только по ортогональности.
5.4 Итерационный метод отыскания максимального по модулю собственного значения матрицы А.
Для простоты полагаем наличие полной системы собственных векторов :
. (5.4.1)
Зададим некоторый ненулевой вектор . Последовательно будем вычислять .
Разложим по собственным векторам матрицы А: …
Так как выполняется условие (5.4.1), то
. (5.4.2)
Вычислим скалярные произведения
. (5.4.5)
. (5.4.6)
Положим
.
Таким образом получаем итерационную схему для нахождения наибольшего собственного числа
(5.4.7)
Кроме того из (5.4.2) следует, что
. (5.4.8)
Полагаем
, где
Итак, мы получили схему итерационного процесса
. (5.4.9)
в результате которого получаем вектор , соответствующий .
Примервычисления первого собственного числа λ1 матрицы методом скалярных произведений.
Выбираем произвольный вектор и строим последовательность векторов , где.
.
,
5.5 Решение полной проблемы собственных значений
при помощи QR алгоритма.
(Основа наиболее совершенных стандартных программ).
Пусть А – произвольная вещественная матрица.
Лемма: Произвольная квадратная матрица может быть представлена в виде произведения ортогональной и правой треугольной матрицы. (без доказательства).
В среде "MathCad" это разложение выполняет функция М=gr(A).
Аргументом ее есть исходная матрица, а результатом
М =
.
Если размерность , то .
Для выделения ортогональной Q и правой треугольной матрицы R можно воспользоваться встроенной функцией
Итак любую вещественную матрицу А можно представить в виде и матрица подобна матрице А. Составим последовательность матриц
; и , (5.5.1)
тогда все матрицы подобны между собой и подобны исходной матрице А. Представим А в виде
, (5.5.2)
где Λ – правая треугольная форма Жордана. Такая матрица в которой:
, где - собственные числа матрицы.
при j < iи j > i + 1, а элементы или 1.
Всегда можно подобрать матрицу Q так, чтобы элементы λi шли в порядке невозрастания модулей.
.
Теорема: Пусть в разложении матрицы А все диагональные миноры матрицы Q не вырождены, тогда последовательностьматриц Ап при п → ∞ сходится по форме к клеточному правому треугольному виду Ẵ. (Без доказательства).
Имеется в виду, что после некоторой перестановки строк и одновременно такой же столбцов будут выполняться соотношения: если или k = 1,…. s.
6 ЧИСЛЕННЫЕ МЕТОДЫ РЕШЕНИЯ НЕЛИНЕЙНОГО УРАВНЕНИЯ С ОДНИМ НЕИЗВЕСТНЫМ
Пусть дано уравнение f(x)=0, предположим, что f(x) имеет только изолированные корни (т.е. для каждого корня существует некая окрестность, не содержащая других корней этого уравнения).
Приближенное вычисление изолированных корней уравнения состоит из двух этапов:
1) отделение корней – выделение отрезка, принадлежащего области существования функции f(x), на котором расположен только один корень;
2) уточнение корней – вычисление их с заданной точностью.
Процесс отделения корней уравнения основан на свойстве непрерывных функций/
Теорема (Больцано-Коши): Если непрерывная на отрезке функция f(x) принимает на концах [а;b] значения разных знаков f(a)× f(b)< 0, то внутри [а;b] есть такое с, что f(c)=0. Этот корень будет единственным, если существует и имеет постоянный знак на [а;b].
6.1 Графический метод отделения корней
Уважаемый посетитель!
Чтобы распечатать файл, скачайте его (в формате Word).
Ссылка на скачивание - внизу страницы.