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