Некоторые дополнительные вычислительные методы (стр. 3 из 9)
 |
f(b)
a  |
Действительно, уравнение хорды АВ имеет вид
При х = х1 и y = 0, получим
Полагая, что на отрезке [a, b] вторая производная f''(x) сохраняет постоянный знак, метод хорд сводится к двум различным вариантам.
Из рис. 1 видно, что конец а неподвижен и последовательные приближения: x0=b;
образуют ограниченную монотонно убывающую последовательность, причем a<ξ<…<xn+1<xn<…<x1<x0.
Из рис. 2 видно, что неподвижен конец b и последовательные приближения: x0=a;
образуют ограниченную монотонно возрастающую последовательность, причем
x0<x1<x2<…<xn<xn+1<…<ξ<b.
Таким образом, для вычисления корня уравнения имеем две различные вычислительные формулы. За неподвижный конец выбираем тот конец, для которого знак функции f(x) совпадает со знаком второй производной f''(x).
Пример. Найти положительный корень уравнения
с точностью до 0,002.Решение. Прежде всего отделяем корень. Так как
и 
, то искомый корень 
лежит в интервале 
. Полученный интервал велик, поэтому разделим его пополам. Так как 
то 
. Последовательно применяя формулы, будем иметь: 
Так как
и при 
имеем 
, то можно принять: 
Таким образом,
, где 
. Заметим, что точный корень уравнения есть 
.Метод Ньютона (метод касательных)
Пусть корень ξ уравнения f(x)=0, отделен на отрезке [a, b], причем первая и вторая производные f'(x) и f''(x) непрерывны и сохраняют определенные знаки при
. Найдя какое-нибудь n-ое приближение корня 
, мы можем уточнить его по методу Ньютона следующим образом. Пусть ξ=xn+hn, где hn - величина малая. Отсюда по формуле Тейлора получим: f(xn + hn) ≈ f(xn)+hn f¢(xn)=0. Следовательно, 
. Подставив полученное выражение в формулу ξ=xn+hn, найдем следующее значение корня: 
Графическое нахождение корня методом Ньютона (рис. 3).  |
Если в качестве начального приближения выбрать точку х0=В0 , то процесс быстро сходится. Если же выбрать точку х0=А, то х1
[a, b], и процесс нахождения корня расходится. Рекомендуется: в качестве х0 выбирать такую точку, где f(x0)f''(x0)>0.Пример. Вычислить методом Ньютона отрицательный корень уравнения
с пятью верными знаками.Решение. Полагая в левой части уравнение
получим 
. Следовательно, искомый корень находится в интервале 
. Сузим найденный интервал. Так как 
то 
. В этом последнем интервале 
и 
.Так как 
и 
, то можем принять за начальное приближение 
. Последовательные приближения 
вычисляем по следующей схеме:  |  |  |  |  |
| 0 | -11 | 3453 | -5183 | 0,7 |
| 1 | -10,3 | 134,3 | -4234 | 0,03 |
| 2 | -10,27 | 37,8 | -4196 | 0,009 |
| 3 | -10,261 | 0,2 | - | - |
Останавливаясь на
, проверяем знак значения 
. Так как 
, то, 
и любое из этих чисел дает искомое приближение.Метод итерации
Заменим уравнение f(x)=0 эквивалентным x=φ(x). Выберем некоторое начальное приближение x0 и вычислим дальнейшие приближения по формулам x1= φ(x0), x2= φ(x1), …, xn= φ(xn-1). Если последовательность xn имеет предел, то итерационный процесс
xn= φ(xn-1) (n=1, 2, …) называется сходящимся. Пусть функция φ(x) непрерывна. Переходя к пределу в равенстве xn= φ(xn-1), получим
Следовательно, 
является корнем уравнения x=φ(x) и может быть вычислен по формуле xn= φ(xn-1) (n=1, 2, …) с любой точностью. Для данного метода существуют две теоремы:Теорема 1. Пусть корень
уравнения x=φ(x), а также его последовательные приближения x0, xn= φ(xn-1) (n=1, 2, …) содержатся в интервале [a, b] и 
на [a, b]. Тогда справедливы утверждения: