- •Тема 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) целесообразно выбирать так, чтобы выполнялось условие .
Для скорости сходимости метода справедливы оценки
,
,
где , . Эти оценки указывают на квадратичную сходимость метода Ньютона.
Геометрически метод Ньютона означает, что в качестве точек приближения к корню берутся точки пересечения с осью касательной к кривой .
Последовательные приближения сходятся к действительному корню уравнения монотонно со стороны .