Билет № 22.
1). Метод Зейделя решения систем линейных алгебраических уравнений.
Запись итерационного стационарного процесса в канонической форме. Рассмотрим произвольную квадратную матрицу:
Разложим её на сумму трех матриц , где - диагональная часть матрицы , которая содержит элементы , стоящие на главной диагонали; - нижняя треугольная матрица, - верхняя треугольная матрица.
В классическом методе Зейделя, записанном в канонической форме, полагают
В результате формула принимает вид:
,
или
(*)
Перейдем от векторной формы записи рекуррентной формулы (*) к построчной:
Эти уравнения позволяют последовательно рассчитать компоненты вектора - ой итерации подобно тому, как это делалось во время обратного хода в методе Гаусса:
, .
Формула предполагает, что , . Если матрица удовлетворяет условиям теоремы Самарского: , то, все ее диагональные элементы должны быть строго положительными и, тем самым, не могут обращаться в ноль.
Алгоритм в методе Зейделя прост и удобен для вычислений. Он не требует никаких действий с матрицей . Ранее вычисленные на текущей итерации компоненты сразу же участвуют в расчетах наряду с компонентами и, таким образом, не требуют дополнительного резерва памяти, что существенно при решении больших систем.
Сходимость метода Зейделя в случае, когда матрица удовлетворяет условию теоремы Самарского, т.е. является самосопряженной и положительно определенной. Метод Зейделя сходится для любой системы , в которой матрица обладает свойством диагонального преобладания.
В методе Зейделя при вычислении (k+1)-го приближения неизвестных х1, х2,..., хi-1 вычисления осуществляются следующим образом. Выбираем произвольно начальные приближения корней х1(0), х2(0),..., хn (0)и подставляем в первое уравнение системы
(1):
х1(1) = 1 + 11х1(0) + 12х2(0) + … + 1nхn(0)
полученное первое приближение х1(1) подставляем во второе уравнение системы (1):
х2(1) = 2 + 21х1(0) + 22х2(0) + … + 2nхn(0)
Полученные первые приближения х1(1) и х2(1) подставляем в третье уравнение системы (1):
Х3(1) = 3 + 31х1(0) + 32х2(0) + … + 3nхn(0)
и т.д.
Наконец,
Хn(1) = n + n1х1(0) + n2х2(0) + … + nn-1хn-1(1) + nnхn(0) Аналогично строим вторые, третьи т.д. итерации. Таким образом, предполагая, что k-e приближение корней х1(k) известны, по методу Зейделя строим (k+l)-e приближение по следующим формулам:
…
где k= 0, 1, 2, ..., n. Объясним откуда они появляются. Зейдель рассмотрел Метод простых итераций: x(k+1) = x(k) +
Зейдель предложил разбить = В + С на сумму двух матриц. Тогда итерационная формула Метода простых итераций примет вид: x(k+1) = Cx(k) +
Переносим слагаемое x(k) в левую часть: x(k+1) - x(k) = Cx(k) +
Зейдель предложил поменять индекс k на индекс k + 1.
x(k+1) - x(k+1) = Cx(k) +
Тогда получилось соотношение: (E – B)x(k+1) = Cx(k) +
Окончательно имеем:
x(k+1) = (E – B)(-1)Cx(k) + (E – B)(-1)
zeidel = (E – B)(-1)C
zeidel = (E – B)(-1)
Откуда следуют формулы приведенные выше. Преимущество метода Зейделя заключается в том, что области сходимости не совпадают, так как м. Зейделя может сходиться там, где МПИ расходиться. Число итераций м. Зейделя, меньше так как данные лучше подготовлены (). Скорость сходимости в м. Зейделя такая же как и в МПИ (не хуже).
Сходимость метода Зейделя
В=D+TH; =1; A=TH+D+ТB
Матрица будет положительно определена, если диагональные элементы будут положительны.
Требования к выбору
-
Сходимость стационарной двухслойной схемы
-
Наиболее быстрая сходимость
Требования к матрице В
Матрица-В не должна быть матрицей общего вида. Она должна быть легко обратимой, чтобы трудоемкость не превышала n2 (диагональной, треугольной, их комбинацией).
Точность = kn2, где (k<<n) и k - количество итераций
2). Оценка погрешностей выполнения арифметических операций с множеством чисел.
Выполнение операций над m числами
Ведем обозначения (ab) = (ab)(1+ε), где => +, -, *, /
- произведение чисел
y0 = 1; yi = xi·yi-1;
пренебрегаем
Для сложения: n2ε
В заключение:
-
Математический анализ и линейная алгебра – это абстрактные модели для задач вычислительной математики (все операции не ассоциативны)
-
Программировать оценку погрешности нереально из-за сложности и увеличения затрат. Поэтому используют оценку сверху, на практике данную задачу решают изменением входных параметров итд.