
- •«Петербургский государственный университет путей сообщения» (фгбоу впо пгупс)
- •Глава 1. Методы решения нелинейных уравнений…………..............15
- •Глава 2. Задача приближения функций………………………….…...23
- •Введение
- •Цели и задачи дисциплины
- •Программа дисциплины
- •Особенности современных инженерных задач
- •О культуре вычислений
- •Правила записи приближенных чисел
- •Глава 1. Методы решения нелинейных уравнений
- •Постановка задачи. Отделение корней.
- •Метод бисекции
- •Метод простой итерации
- •Метод Ньютона
- •Глава 2. Задачи приближения функций
- •Метод интерполирования
- •Интерполяционный многочлен в форме Лагранжа
- •Интерполяционный многочлен в форме Ньютона
- •Нелинейная задача метода наименьших квадратов
- •Глава 3. Численные методы решения задачи Коши для обыкновенных дифференциальных уравнений и систем дифференциальных уравнений
- •Задача Коши для дифференциального уравнения первого порядка
- •Численные методы решения задачи Коши
- •Метод сеток. Основные понятия и определения.
- •Примеры разностных схем и соответствующих им численных методов
- •Явная схема Эйлера
- •Неявная схема Эйлера
- •Симметричная схема
- •Схемы Рунге – Кутты
- •Многошаговые схемы (схемы Адамса)
- •Схемы Адамса – Моултона (неявные многошаговые схемы)
- •Сходимость. Погрешность аппроксимации. Порядок точности
- •Оценка погрешности. Правило Рунге
- •Библиографический список
Метод простой итерации
Чтобы применить
метод простой итерации для решения
нелинейного уравнения
,
необходимо преобразовать это уравнение
к следующему виду:
.
Функцию
будем называть итерационной функцией.
Выберем каким-либо образом приближенное
значение корня
и подставим его в правую часть последнего
уравнения. Получим значение
.
Подставляя теперь
в правую часть уравнения
,
имеем
.
Продолжая этот процесс, получим
последовательность приближений к корню,
вычисляемых по формуле
.
Итерационный процесс следует продолжать
до тех пор, пока после некоторой итерации
с номером
не будет выполняться неравенство
.
После этого с погрешностью
полагают:
.
Неравенство
является критерием окончания, который
на практике заменяет собой априорные
и апостериорные оценки погрешности.
Пример: С
помощью метода простой итерации вычислить
положительный корень уравнения
с точностью
.
Решение:
Результат предыдущего примера дает
отрезок локализации
.
Преобразуем уравнение к виду
,
где
.
Заметим, что
,
.
Так как
на
,
то производная
монотонно убывает и
.
Следовательно, условие сходимости
выполнено. Возьмем
,
и будем вести итерации до выполнения
критерия
.
В таблице приведены соответствующие
приближения с 10 знаками мантиссы.
|
|
|
0 |
0.7000000000 |
|
1 |
0.7046714292 |
|
2 |
0.7029968319 |
|
3 |
0.7035984939 |
|
4 |
0.7033824994 |
|
5 |
0.7034600632 |
|
Критерий окончания
выполняется при
.
После округления значения
до четырех значащих цифр получим
.
Метод Ньютона
Этот метод является
одним из наиболее эффективных методов
решения самых разных нелинейных
уравнений. Расчетная формула метода
Ньютона
.
При выборе начального приближения из
достаточно малой окрестности корня
метод Ньютона сходится квадратично.
Практический критерий окончания итераций
метода Ньютона
.
Пример: Используя
метод Ньютона, найти с точностью
положительный корень уравнения
.
Решение: Корень
уравнения локализован на отрезке
.
Для
имеем
.
Очевидно, что
,
то есть
– простой корень. Возьмем начальное
приближение
(середина отрезка) и будем выполнять
итерации метода Ньютона по формуле
.
Результаты первых итераций с 10 знаками
мантиссы приведены в таблице.
|
|
|
0 |
0.5000000000 |
|
1 |
0.7392185177 |
|
2 |
0.7042444088 |
|
3 |
0.7034399951 |
|
4 |
0.7034395712 |
|
При
вычисления следует прекратить и после
округления получим
.
В качестве недостатка
метода Ньютона следует отметить
необходимость вычисления значения
производной
на каждом шаге итерации. Приведем
некоторые расчетные формулы методов
решения нелинейного уравнения
,
которые являются, по существу, модификациями
метода Ньютона. Эти методы используют
на каждой итерации некоторую процедуру
линеаризации нелинейного уравнения.
Упрощенный метод
Ньютона:
,
;
Метод ложного
положения:
,
,
где
– фиксированная точка из окрестности
корня;
Метод секущих:
,
;
Метод Стеффенсена:
,
.
Модифицированный метод Ньютона для поиска кратных корней:
,
.
Контрольные вопросы
Из каких этапов состоит численное решение уравнения
?
Как графически можно произвести отделение корней?
Приведите пример отделения корней нелинейного уравнения.
Какая теорема лежит в основе локализации корней уравнений.
Назовите известные Вам методы решения нелинейных уравнений.
Опишите одну итерацию метода бисекции.
Во сколько раз сократится отрезок локализации корня после
- ой итерации метода бисекции.
С какой скоростью сходится метод бисекции?
Сформулируйте критерий окончания метода бисекции.
К какому виду необходимо привести уравнение для применения метода простой итерации.
Опишите алгоритм метода простой итерации.
Сформулируйте критерий окончания метода простой итерации.
По какой формуле рассчитывается очередное приближение к корню нелинейного уравнения по методу Ньютона.
Опишите алгоритм метода Ньютона.
С какой скоростью сходится метод Ньютона?