
- •Тема 2. Численное решение нелинейных уравнений Вопросы для самоподготовки:
- •Отделение корней
- •1.1 Общие понятия
- •Задание 1.
- •1.2 Индивидуальные задания
- •Уточнение корней уравнения методами половинного деления и простой итерации
- •2.1 Уточнение корня. Постановка задачи
- •2.2 Метод половинного деления
- •Метод хорд
- •Метод Ньютона
- •Комбинированный метод
- •2.3 Метод простой итерации (метод последовательных приближений)
- •2.4 Решение уравнений средствами MathCad
- •2.4.1. Функции произвольного вида
- •2.4.2 Нахождение корней полиномов
- •2.4.3 Нахождение корней уравнений путем символических преобразований
- •2.4.4 Поиск корней уравнений в Mathcad
- •2.5 Блок-схемы
- •2.6 Индивидуальные задания
- •Задание 3.
Уточнение корней уравнения методами половинного деления и простой итерации
2.1 Уточнение корня. Постановка задачи
Если корень
уравнения f(х)=0
отделен, т.е. найден отрезок [а; b],
на котором имеется один и только один
корень уравнения, то любую точку этого
отрезка можно принять за приближенное
значение корня. Погрешность такого
приближения не превосходит длины [а; b].
Следовательно, задача отыскания
приближенного значения корня с заданной
точностью
сводится к нахождению отрезка [а; b],
такого что │b
- а│<
.
Эту задачу обычно называют задачей
уточнения
корня.
Итерационный процесс состоит в последовательном уточнении начального приближения х0. Каждый такой шаг называется итерацией. В результате итераций находится последовательность приближенных значений корня х1, х2, ..., хn. Если эти значения с увеличением числа итераций n приближаются к истинному значению корня, то говорят, что итерационный процесс сходится.
2.2 Метод половинного деления
Разделим отрезок [а; b], содержащий единственный корень, пополам и определим знак функции f(х) в точке х=(а+b)/2. Если f((а+b)/2)=0, тогда корень уравнения найден. Если же f((а+b)/2)≠0, то на концах одного из отрезков [а; (а+b)/2] или [(а+b)/2; b] функция будет принимать значения разных знаков.
Обозначим этот
отрезок через [a1; b1].
Если │b1-a1│<,
то любая точка интервала (a1; b1)
может быть принята за приближенное
значение корня. Если же │b1-
a1│≥,
то, положив а=а1
, b=b1
и продолжая процесс деления отрезка
пополам, на каком-то конечном шаге
получим точное значение корня, либо
через конечное число шагов длина
[a; b]станет
меньше .
В последнем случае за приближенное
значение корня можно принять любую
точку отрезка [a; b]
(желательно, его середину).
Метод половинного деления практически удобно применять для грубого нахождения корня данного уравнения, метод прост и надежен, всегда сходится.
Количество шагов, необходимых для определения корня с точностью , можно рассчитать как
Метод хорд
В методе хорд точка
пересечения кривой
с осью
на каждом шаге приближается точкой
пересечения хорды, стягивающей концы
дуги
,
,
с осью
.
Последовательные приближения при этом
вычисляются по формуле
. (2)
В
формуле (2) в качестве
берется тот конец отрезка
,
для которого выполняется условие
,
а начальное приближение
выбирается так, что
.
Геометрическая интерпретация метода хорд:
Описанный метод
называется также методом секущих или
методом линейной интерполяции.
Последовательные приближения в методе
хорд образуют монотонную ограниченную
сверху или снизу корнем
последовательность. При этом справедлива
оценка
,
где
,
.
Метод секущих
можно рассматривать как метод итерации
для эквивалентного уравнения ,
где
и начальное приближение берется так,
что
.
Метод Ньютона
Уточнение корня методом Ньютона опирается на соотношение
.
(3)
Начальное
приближение в (3) целесообразно выбирать
так, чтобы выполнялось условие .
Для скорости сходимости метода справедливы оценки
,
,
где
,
.
Эти оценки указывают на квадратичную
сходимость метода Ньютона.
Геометрически
метод Ньютона означает, что в качестве
точек приближения к корню берутся точки
пересечения с осью
касательной к кривой
.
Последовательные
приближения сходятся к действительному
корню уравнения монотонно со стороны
.