- •Лекция 1
- •Лекция 2
- •Уравнении с одним неизвестным
- •Метод деления отрезка пополам (метод бисекции).
- •Метод Ньютона (метод касательных).
- •Метод простой итерации.
- •Лекция 3 Системы линейных уравнений
- •3.Прямые методы
- •Итерационные методы
- •Метод Гаусса — Зейделя.
- •Системы уравнений.
- •Лекция 4 Аппроксимация функций.
- •Точечная аппроксимация.
- •Многочлены Чебышева.
- •Линейная и квадратичная интерполяция.
- •Сплайны.
- •Многочлен Лагранжа.
- •Лекция 5.
- •Метод Симпсона.
- •Использование сплайнов.
- •Кратные интегралы.
- •Метод Монте-Карло.
- •Лекция 6. Численное дифференцированно.
- •Погрешность численного дифференцирования.
- •Метод неопределенных коэффициентов.
- •Улучшение аппроксимации.
- •Частные производные.
- •Лекция 7.
- •Основные понятия.
- •Задачи оптимизации.
- •Одномерная оптимизация.
- •1. Задачи на экстремум.
- •Метод золотого сечения.
- •Метод градиентного спуска.
- •Лекция 8. Задачи с ограничениями.
- •Линейное программирование.
- •Геометрический метод.
- •Симплекс-метод.
- •Задача о ресурсах.
Лекция 2
НЕЛИНЕЙНЫЕ УРАВНЕНИЯ
Уравнении с одним неизвестным
Вводные замечания. Задача
нахождения корней нелинейных уравнений
вида ![]()
встречается в различных областях научных исследований (здесь— некоторая непрерывная функция). Нелинейные уравнения можно разделить на два класса - алгебраические и трансцендентные. Алгебраическими уравнениями называются уравнения, содержащие только алгебраические функции (целые, рациональные, иррациональные). В частности, многочлен является целой алгебраической функцией. Уравнения, содержащие другие функции (тригонометрические, показательные, логарифмические и др.), называются трансцендентными.
Методы решения нелинейных уравнений делятся на прямые и итерационные. Прямые методы позволяют записать корни в виде некоторого конечного соотношения (формулы). Вам известны такие методы для решения тригонометрических, логарифмических, показательных, а также простейших алгебраических уравнений.
Однако встречающиеся на практике уравнения не удастся решить такими простыми методами. Для их решения используются итерационные методы, т.е. методы последовательных приближений. Алгоритм нахождения корня уравнения с помощью итерационного метода состоит из двух этапов:
а) отыскания приближенного значения корня или содержащего его отрезка;
б) уточнения приближенного значения до некоторой заданной степени точности.
Начальное
приближение может быть найдено различными
способами: из
физических соображений, из решения
аналогичной
задачи
при других исходных
данных, с помощью графических методов.
Если такие априорные
оценки исходного приближения провести
не удастся, то находят две близко
расположенные точки
и
,
в которых непрерывная функция
принимает
значения разных знаков, т.е.
.
В этом
случае между
точками
и
есть
но крайней мере одна точка, в которой
.В качестве начального
приближения
можно принять середину отрезка
, т.е.
.
Итерационный процесс состоит
в последовательном
уточнении начального приближения
.
Каждый
такой шаг
называется
итерацией.
В
результате итераций находится
последовательность приближенных
значений корня
Если эти значения с
ростом
стремятся
к истинному значению корня то говорят,
что итерационный процесс сходится.
Метод деления отрезка пополам (метод бисекции).
Это
один из простейших
методов нахождения корней нелинейных
уравнений. Он состоит
в следующем. Допустим, что нам удалось
найти отрезок
,
на
котором расположено искомое значение
корня
,т.е.
.
В качестве
начального приближения корня
принимаем
середину этого отрезка,
т.е.
.
Далее
исследуем значения функции
на
концах отрезков
и
,
т. е. в точках
.
Тот
из отрезков, на концах которого
принимает
значения разных знаков, содержит искомый
корень; поэтому его принимаем в
качестве нового отрезка . Вторую половину
отрезка
, на которой знак
не
меняется, отбрасываем. В
качестве первой итерации корня принимаем
середину нового отрезка и т. д. Таким
образом, после каждой итерации отрезок,
на котором расположен корень, уменьшается
вдвое, а после
итераций он сокращается
в
раз.
Пусть для определенности
,
.
В
качестве начального приближения корня
примем
.
Поскольку в рассматриваемом случае
,
то
,
и рассматриваем только отрезок
.
Следующее приближение:
.
При этом отрезок
отбрасываем, поскольку
и
т.е.
.
Аналогично находим другие приближения:
и т.д.
Итерационный
процесс продолжаем до тех пор, пока
значение функции
после
-ой
итерации не станет меньшим по модулю
некоторого заданного малого числа
,
т.е.
.
Можно также оценивать длину полученного отрезка: если она становится меньше допустимой погрешности, то счет прекращается.
Здесь
сужение отрезка
производится путем замены границ
или
на
текущее значение корня
.
При
этом значение
вычисляется
лишь один раз, поскольку нам
нужен только знак функции
па
левой границе, а он в процессе итераций
не меняется.
Метод
деления отрезка пополам довольно
медленный, однако он всегда сходится,
т.е. при его использовании решение
получается всегда, причем с заданной
точностью. Требуемое обычно большее
число итераций по сравнению с некоторыми
другими методами не является препятствием
к применению этого метода, если каждое
вычисление значения функции
несложно.
Метод хорд.
Пусть
мы нашли отрезок
,
на
котором функция
меняет
знак. Для определенности
примем
,
.
В данном методе процесс итерации
состоит в том, что в качестве приближений
корню
уравнения
принимаются
значения
точек
пересечения
хорды с осью абсцисс.
Сначала находим уравнение
хорды
:
.
Для точки пересечения её с осью абсцисс
получим уравнение
.
Далее, сравнивая знаки величин
и
для рассматриваемого случая, приходим
к выводу, что корень находится в интервале
,
т.к.
.
Отрезок
отбрасываем. Следующая итерация состоит
в определении нового приближения
как точки пересечения хорды
с осью абсцисс и т.д. Итерационный процесс
продолжается до тех пор, пока значение
не станет по модулю меньше заданного
числа
.
К
ак
видим, алгоритмы метода
деления отрезка пополам и метода хорд
похожи, однако
второй из них в ряде случаев дает более
быструю сходимость итерационного
процесса. При этом успех его применения,
как и в методе деления отрезка пополам,
гарантирован.
Блок
схема метода хорд аналогична предыдущей
с той лишь разницей, что вместо вычисления
приближения корня по формуле
нужно использовать формулу
.
Кроме того, в блок-схему необходимо
ввести операторы вычисления значений
на границах новых отрезков.
