КАТЕГОРИИ: Архитектура-(3434)Астрономия-(809)Биология-(7483)Биотехнологии-(1457)Военное дело-(14632)Высокие технологии-(1363)География-(913)Геология-(1438)Государство-(451)Демография-(1065)Дом-(47672)Журналистика и СМИ-(912)Изобретательство-(14524)Иностранные языки-(4268)Информатика-(17799)Искусство-(1338)История-(13644)Компьютеры-(11121)Косметика-(55)Кулинария-(373)Культура-(8427)Лингвистика-(374)Литература-(1642)Маркетинг-(23702)Математика-(16968)Машиностроение-(1700)Медицина-(12668)Менеджмент-(24684)Механика-(15423)Науковедение-(506)Образование-(11852)Охрана труда-(3308)Педагогика-(5571)Полиграфия-(1312)Политика-(7869)Право-(5454)Приборостроение-(1369)Программирование-(2801)Производство-(97182)Промышленность-(8706)Психология-(18388)Религия-(3217)Связь-(10668)Сельское хозяйство-(299)Социология-(6455)Спорт-(42831)Строительство-(4793)Торговля-(5050)Транспорт-(2929)Туризм-(1568)Физика-(3942)Философия-(17015)Финансы-(26596)Химия-(22929)Экология-(12095)Экономика-(9961)Электроника-(8441)Электротехника-(4623)Энергетика-(12629)Юриспруденция-(1492)Ядерная техника-(1748) |
Метод Зейделя (модификация метода итераций)
ЛЕКЦИЯ № 5
При решении системы линейных алгебраических уравнений вида (1) методом итераций значение вычисляется по значениям с предыдущей итерации , ,..., путем подстановки в правую часть системы (4). Можно ожидать, что приближения будут быстрее сходиться к решению системы, если сразу же после вычисления , при вычислении последующих использовать , а не в правой части (4). Процедура вычисления через , ,..., , , называется методом Зейделя и записывается в развернутой форме в виде: (9) Запишем метод Зейделя в векторной форме, для этого представим матрицу a в виде суммы двух треугольных матриц L и U, где , Тогда систему (9) можно записать в виде матричного равенства: (10) Матрица (E-L) - неособенная, т.е. имеет обратную (E-L)-1, следовательно, можно выразить х(k+1) из (10) Из (10) получаем, что
Обозначим тогда (11) Следовательно, метод Зейделя для системы (1) эквивалентен методу простой итерации x = ax + b для системы x = Px + Q, где матрица P и вектор Q определены выше. Теперь для сходимости (11) достаточно, чтобы ||P||1 < 1 или ||P||2 < 1. Используя собственные значения матрицы P можно дать необходимое и достаточное условие сходимости процесса итераций для системы (11): |l(P)| < 1 Здесь в качестве матрицы a выступает матрица P, а в качестве вектора b - вектор Q. Если для одной и той же системы методы итерации и Зейделя сходятся, то метод Зейделя предпочтительнее. Достаточное условие сходимости процесса Зейделя. ТЕОРЕМА: Если для линейной системы х = aх + b (2) выполнено условие , где , то процесс (9) для системы (2) сходится к единственному решению при любом выборе начального приближения. Доказательство: ………………………………………………………………………………
Оценим погрешности приближений по методу Зейделя. Пусть и - две последовательные итерации процесса Зейделя. Применяя к этим итерациям преобразования, получим:
Выполним аналогичные (как для МПИ) преобразования для разности между (k+m)- м и k -м членами последовательных приближений по Зейделю при некотором mÎN:
Рассматривая итоговое равенство при , переходя к пределу получим утверждения теоремы: , где
Тогда условие окончания итерационного процесса Зейделя будет иметь вид: или
Тогда или ,
Дата добавления: 2014-12-07; Просмотров: 492; Нарушение авторских прав?; Мы поможем в написании вашей работы! Нам важно ваше мнение! Был ли полезен опубликованный материал? Да | Нет |