Уточнение корней до заданной точности.
То есть сужение отрезка локализации корня [a,b]. Рассмотрим несколько методов.
Метод половинного деления (дихотомии).
Пусть корень отделён и принадлежит отрезку . Находим середину отрезка по формуле (рис.3). Если , то с – искомый корень.
|
Рис. 3.
Рис.3
Вычисляем , выбираем отрезок и т.д. Как только будет выполнено , то в качестве приближенного значения корня, вычисленного с точностью , можно взять . После каждой итерации отрезок, на котором расположен корень уменьшается вдвое, то есть после n итераций он сокращается в 2n раз. Таким образом, число итераций n в данном методе зависит от предварительно заданной точности ε и от длины исходного отрезка и не зависит от вида функции f(x). Это является важным преимуществом метода половинного деления по сравнению с другими методами. Метод, однако, медленно сходится при задании высокой точности расчёта.
Метод хорд.
Пусть на отрезке [a,b] функция f(x) непрерывна и принимает на концах отрезка значения разных знаков, а производные f ′(x) и f ″(x) сохраняют постоянный знак на интервале (a,b). Тогда возможны четыре случая расположения дуги кривой (рис.4).
|
|
Рис.4.
В методе хорд за очередное приближение берём точку пересечения с осью Х прямой (рис.5), соединяющей точки (a,f(a)) и (b,f(b))
Причём одна из этих точек фиксируется − та, для которой знаки f(x) и f ″(x) одинаковы.
Для рис.5 неподвижным концом хорды является х =a.
Уравнение хорды АВ:
Точка пересечения хорды с осью Х (у=0): .
Теперь корень находится на отрезке [a,c1]. Заменяем b на с1.
Рис.5. Иллюстрация метода хорд.
Применяя метод хорд к этому отрезку, получим:
.
Продолжим и т.д., получим: (2)
Условие окончания вычислений:
│сn+1 − cn│< ε или │f(cn)│< ε1.
Для оценки погрешности можно пользоваться общей формулой:
, где
Итак, если f (x)∙f″(x) > 0, то приближённое значение корня находят по формуле (2), если f′(x)∙f″(x) < 0 (т.е. фиксируется х = b), то по формуле:
. (3)
|
|
Метод Ньютона (касательных).
Пусть на отрезке [a,b] функция f(x) непрерывна и принимает на концах отрезка значения разных знаков, а производные f ′(x) и f ″(x) сохраняют постоянный знак на интервале (a,b).
Геометрический смысл метода касательных состоит в том, что дуга кривой
y = f(x) заменяется касательной к этой кривой.
Рис.7.Иллюстрация метода касательных.
Выберем в качестве начального приближения х0 = a и проведём в точке А0(a,f(a)) касательную к графику функции f(x). Абсцисса пересечения касательной с осью Ох (у = 0) является первым приближением к корню (рси.7):
или х0 = .
Через точку А1(х1;f(x1)) снова проведём касательную, абсцисса точки пересечения которой даст второе приближение х2 корня ξ и т.д. Очевидно, что в точке Аn(xn;f(xn)):
y − f(xn) = f ′(xn)(x−xn)
и алгоритм метода Ньютона запишется так:
(4)
Заметим, что в нашем случае, если положить х0 = b и провести касательную к кривой у = f(x) в точке b, то первое приближение не принадлежит отрезку [a,b].
Таким образом, в качестве начального приближения х0 выбирается тот конец интервала [a,b], для которого знаки f(x) и f ″(x) одинаковы.
|
|
Условие окончания вычислений: │сn+1 − cn│< ε или │f(cn)│< ε1.
Для оценки погрешности можно пользоваться общей формулой
, где
Дата добавления: 2018-08-06; просмотров: 339; Мы поможем в написании вашей работы! |
Мы поможем в написании ваших работ!