A и a – bE имеют общую систему собственных векторов. Если X



Скачать 54.14 Kb.
Дата11.10.2012
Размер54.14 Kb.
ТипЛекция
Лекция 8

30 октября 2006 года

Метод обратной итерации.

Рассмотрим квадратную матрицу A. Для простоты считаем ее самосопряженной невырожденной матрицей, среди собственных значений которой нет кратных. Выбираем действительный параметр b. Матрицы A и A – bE имеют общую систему собственных векторов. Если xk — собственный вектор матрицы A, то верно равенство



Введем обозначение Рассмотрим последовательность равенств



Применим степенной алгоритм для нахождения максимального по абсолютной величине собственного числа матрицы B. Мы уже знаем, что он сойдется. Таким образом, мы можем найти собственное число матрицы A, наиболее близкое к b.

Схема метода, таким образом, есть.

  1. Фиксируем параметр b (для определения нижней границы спектра, очевидно, необходимо положить b = 0).

  2. Фиксируем произвольный ненулевой вектор

  3. Решаем систему линейных уравнений

Это и есть обратная итерация.

После того, как найдена последовательность векторов, считаем



Таким образом, мы можем найти собственное число матрицы, ближайшее к b.

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

5.1 Сжимающие отображения. Итерации. Метод простых итераций (МПИ)

Рассмотрим системы нелинейных алгебраических уравнений, записанные в векторном виде.

Система нелинейных алгебраических уравнений

f(u) = 0 (5.1)

может быть также представлена в равносильном виде

u F(u), (5.2)

где (Ln n-мерное евклидово пространство). Как правило, для нелинейной системы переход от формы записи (5.1) к равносильному виду (5.2) осуществляется не единственным образом.

Поставим в соответствие системе (5.2) итерационный процесс, определяющий последовательность итераций (последовательных приближений к решению). Соответствующий итерационный процесс записывается в форме

gif" name="object8" align=absmiddle width=88 height=28> u0 = a, k = 0,1,… (5.3)

Для дальнейшего изложения потребуется понятие отображения. Отображением называется закон, по которому каждому элементу x некоторого множества X однозначно сопоставляется определенный элемент y множества Y (X может совпадать с Y). Это соотношение между элементами и записывается как y = f(x) или f: x → y. Говорят, что отображение f действует из X в Y (f: X → Y). Отображение f: X → X называют преобразованием множества X, это отображение f преобразует множество X в себя. В функциональном анализе и линейной алгебре вместо термина отображение часто употребляется термин «оператор», в случае, если X и Y — числовые множества, употребляется термин «функция».

Определение. Область называется выпуклой, если наряду с любыми двумя точками и она включает все точки отрезка [ab], т.е. точки с координатами u = a + t(b  a), где .

Определение. Отображение называется сжимающим в замкнутой выпуклой области Ω, если существует такое число 0 < q < 1, что



при любых u1, u2, принадлежащих области Ω, здесь — расстояние между элементами. В линейном нормированном пространстве .

Если отображение f: x → y переводит точку x* в себя, то эта точка называется неподвижной точкой отображения (или стационарной точкой).

Приведем одну из основных теорем функционального анализа.

Теорема (принцип сжимающих отображений). Всякое сжимающее отображение имеет в Ω одну и только одну неподвижную точку u*Ω.
Теорема (о сжимающем отображении).

Последовательность k = 0,1,… элементов n-мерного евклидова пространства, порожденная итерационным процессом



сходится к решению U системы нелинейных алгебраических уравнений = F(u), если отображение



является сжимающим; при этом выполнено



Доказательство.

По определению сжимающего отображения



В таком случае получим цепочку неравенств при р > k:



В соответствии с критерием Коши существования предела последовательности, последовательность {uk} стремится к пределу U, поскольку правая часть неравенства стремиться к нулю при .

Напомним критерий Коши сходимости числовой последовательности: последовательность является сходящейся, если для любого положительного числа ε существует номер N такой, что при всех k > N и любых натуральных p расстояние между членами последовательности uk и uk + 1 меньше ε, т.е. .

Напомним критерий Коши для последовательности элементов метрического пространства: последовательность является сходящейся, если для любого существует номер N такой, что при всех k > N и любом натуральном p расстояние

Продолжим доказательство. Переходя в последнем неравенстве к пределу при , получим



Покажем, что U есть корень уравнения (5.2)





Поскольку k выбрано произвольно, а левая часть от k не зависит, то , или U F(U).

В случае скалярного уравнения имеем ( 0  t  1)



откуда следует условие сходимости итерационного метода Отрезок Δ включает в себя всю последовательность {uk},

В случае решения системы нелинейных уравнений достаточным условием сходимости итерационного процесса будет где — матрица Якоби.

Теорема (без доказательства)

Пусть область выпуклая, а компоненты Fi(u) вектор-функции имеют равномерно непрерывные производные первого порядка. Положим, что норма матрицы Якоби



не превосходит некоторого числа , т.е. для всех

В этом случае отображение F(u) является сжимающим в области G, т.е. или

Геометрическая интерпретация метода простой итерации для скалярного случая приведена на рис. 1. Алгоритм метода простых итераций таков.

1. Локализуем корень, приближенно определяем, на каком отрезке он находится. Вопрос локализации корня не решается алгоритмически, это скорее вопрос искусства вычислителя, хотя во многих случаях локализовать корень достаточно легко.
Рис. 1. Геометрическая интерпретация метода простых итераций.

2. Выбираем точку u0 на оси 0u

3. Вычисляем F(u0)

4. Определяем точку u1 по значению F(u0):

4.1 Пересечение горизонтальной прямой AA’ с прямой v = u есть точка С (ОА = v1, AC = u1)

4.2. Очевидно, что горизонтальная координата точки С и есть u1 (так как F(u0) = u1).

4.3. Опустим перпендикуляр из С на 0u. Поскольку ОА = u1, то u1 — значение на первой итерации.

5. Аналогично строим точки u2, u3, … Получившаяся диаграмма носит название лесенка Ламерея.
Метод релаксации. Без ограничения общности рассмотрим скалярный случай. Положим и построим итерационный процесс

Тогда и τ выбирается из условия причем, чем меньше значение тем быстрее будет сходиться итерационный процесс. В частности, если положить , то, а формулы итерационного процесса будут u0 = a.

Похожие:

A и a – bE имеют общую систему собственных векторов. Если X iconСеминару №1. Определить, является ли данная совокупность векторов линейно зависимой
Показать на примере матрицы А, что, если матрица имеет n попарно различных собственных чисел, ее ранг равен числу отличных от нуля...
A и a – bE имеют общую систему собственных векторов. Если X icon"Нахождение собственных значений и собственных векторов матриц"
Численное решение задачи нахождения собственных значений матриц методом Леверрье-Фаддеев
A и a – bE имеют общую систему собственных векторов. Если X iconПаркеты Автор работы: Холодова Оксана Место выполнения работы
В каждом из этих замощений любые два мно­гоугольника имеют либо общую сторону, либо только общую вершину, либо вовсе не имеют общих...
A и a – bE имеют общую систему собственных векторов. Если X iconТеорема Хелли
На плоскости даны четыре выпуклые фигуры, каждые три из которых имеют общую точку. Докажите, что все четыре фигуры имеют хотя бы...
A и a – bE имеют общую систему собственных векторов. Если X iconТест Вопросы I
Если в четырехугольнике диагонали равны, взаимно перпендикулярны и имеют общую середину, то этот четырехугольник
A и a – bE имеют общую систему собственных векторов. Если X iconОпределение. Система векторов называется линейно зависимой, если существуют вещественные числа не равные одновременно нулю такие, что (*) Если же равенство (*)
Линейная зависимость и независимость системы векторов. Коллинеарные и компланарные векторы
A и a – bE имеют общую систему собственных векторов. Если X iconНахождение собственных чисел и собственных векторов матриц Теорема 19. 1
Доказательство. Пусть столбец собственный вектор матрицы с собственным числом. Тогда, по определению. Это равенство можно переписать...
A и a – bE имеют общую систему собственных векторов. Если X iconВариант 1 №1. Вычислить определитель матрицы. №2. Выполнить действия A(B–C)-1D, где. №3. Решить матричное уравнение №4. Исследовав систему на совместность, найти ее общее решение методом Гаусса №5
Найти собственные значения и собственные и присоединенные векторы матрицы линейного оператора и найти вид этой матрицы в базисе из...
A и a – bE имеют общую систему собственных векторов. Если X iconВведение в общую топологию
Курс ориентирован не на саму общую топологию, а на те ее базовые понятия и конструкции, которые имеют приложения в смежных разделах...
A и a – bE имеют общую систему собственных векторов. Если X icon75 процентов от указанных выше затрат (без учета ндс для получателей субсидии, применяющих общую систему налогообложения), но не более 5 млн рублей
Ндс для получателей субсидии, применяющих общую систему налогообложения, но не более 5 млн рублей одной действующей малой инновационной...
Разместите кнопку на своём сайте:
ru.convdocs.org


База данных защищена авторским правом ©ru.convdocs.org 2016
обратиться к администрации
ru.convdocs.org