
- •Численные методы
- •Правила округления.
- •Оценка погрешностей арифметических операций и функций.
- •Оценка погрешности функции одной переменной.
- •Лекция № 2. Особенности машинной арифметики. Корректность вычислительной задачи.
- •Принцип равных влияний.
- •Представление данных в эвм.
- •Нелинейные уравнения.
- •Метод простой итерации (мпи).
- •Универсальный алгоритм приведения уравнения к виду,
- •Лекция № 5.
- •Обусловленность задачи решения слау.
- •Лекция № 6. Решение слау прямыми методами. Метод Гаусса и его модификации.
- •Модификации метода Гаусса.
- •Lu-разложение или матричная форма метода Гаусса.
- •Применение прямых методов для решения слау.
- •Метод простой итерации.
- •Лекция № 8.
- •Лекция № 9.
- •Лекция № 10.
- •Итерационный процесс
- •С чебышевским набором параметров.
- •Многочлены Чебышёва.
- •Вспомогательная задача.
- •Чебышёвский процесс для вычисления решения системы.
- •Лекция № 11.
- •Лекция № 13.
- •Нелинейная задача метода наименьших квадратов.
- •Метод Фибоначчи.
- •Градиентные методы.
Применение прямых методов для решения слау.
Метод частичного выбора (однократное вычисление системы).
система с разными правыми частями.
Если нашли
, то решаем фактически
.
метод Холецкого.
A разреженная. Метод прогонки.
Метод простой итерации.
Расчётные формулы:
Пример.
Покоординатная форма записи для МПИ:
Теорема 2.
О достаточном условии сходимости МПИ.
Пусть рассматривается
МПИ
и пусть выполняется условие
.
Тогда МПИ сходится при любом начальном
приближении и справедливы оценки:
–решение системы
.
Следствие.
Критерий окончания МПИ:
,
где
– заданная точность.
Существуют две формы записи МПИ:
Матричная форма.
Покоординатная форма.
Если матрица B
обладает свойством диагонального
преобладания ,
то можно доказать, что
.
Иногда диагональное преобладание можно получить, переставляя строки в исходной системе уравнений.
Лекция № 8.
РЕАЛИЗАЦИЯ
МПИ НА MathCAD.
Алгоритм:
Вычисление
по МПИ с заданной точностью .
Методы Зейделя и релаксации.
Покоординатная форма записи метода Зейделя.
Матричная форма записи метода Зейделя.
Достаточное условие
сходимости МПИ .
Достаточное условие
сходимости МЗ .
Теорема 1.
Пусть рассматривается
метод Зейделя .
Тогда если
,
то метод Зейделя сходится при любом
начальном приближении и справедливы
следующие оценки:
Следствие 1.
Метод Зейделя сходится со скоростью геометрической прогрессии.
Следствие 2.
Критерий окончания
.
Геометрическая интерпретация метода Зейделя.
Замечание.
Метод Зейделя может зациклиться.
SOR метод (верхней релаксации).
Каноническая форма записи итерационных методов.
Определение 1.
Канонической формой итерационного метода общего вида называется:
Если ,
то метод называется неявным,
иначе – явным.
Если
зависят от номера итерации, то метод
называется итерационным,
иначе – стационарным.
B – итерационная матрица.
–итерационный
параметр.
Вектор погрешности
Вектор ошибки
удовлетворяет уравнению
Определение 2.
называется матрицей
перехода от
вектора погрешности от
к
для МПИ
.
Теорема 2.
Если все собственные
числа
матрицы S
,
то метод итерационно сходится при любом
начальном приближении.
Необходимость
Достаточность
Предположим дополнительно, что S имеет собственные векторы, образующие собственный базис.
Лекция № 9.
ИТЕРАЦИОННЫЕ
МЕТОДЫ
ДЛЯ ПОЛОЖИТЕЛЬНО ОПРЕДЕЛЁННЫХ МАТРИЦ.
Итерационные методы в канонической форме.
Критерий Сильвестра:
Запишем наиболее распространённые методы в канонической форме:
Метод Якоби.
Метод релаксации.
–параметр, берётся
всегда положительным.
Метод Зейделя
.
Условия сходимости итерационных методов
для симметричных, положительно определённых матриц.
Теорема 1.
Пусть
и
.
Тогда итерационный метод сходится при
любом начальном приближении.
Докажем, что
.
Рассмотрим
числовую последовательность .
Докажем, что
убывает.
Следствие 1.
Пусть матрица A
обладает свойством диагонального
преобладания, т. е. .
Тогда метод Якоби сходится при любом
начальном приближении.
.
Следствие 2.
Если ,
то метод релаксации
и, в частности метод Зейделя
,
сходится при любом начальном приближении.
Утверждение 1.
Пусть
и
– спектр матрицы. Тогда метод простой
итерации
сходится при
и любом начальном приближении.
Пусть
,
тогда метод сходится при любом начальном
приближении.
Если
– собственные
числа A,
то
.
Если
–собственный
вектор.
минимальная погрешность метода;
метод сходится.
Замечание.
В МПИ необходимым
и достаточным условием выбора параметра
является
.