Курсовая работа: Приближённые методы решения алгебраического уравнения
|an+ 1 |<qn ·|a1 | (2.4)
В самом деле, из равенства (1.4) имеем:
|a2 |=|j ¢(c1 )|·|a1 |
Но точка c1 лежит на отрезке [a, b] (рис.1.4), и потому:
|j ¢(c1 )|<q
Отсюда следует, что:
|a2 |<q·|a1 |
Точно так же получаем, что:
|a3 |=|j ¢(c1 )|·|a2 |<q·|a2 |< q2 ·|a1 |
и вообще:
|an+ 1 |=qn ·|a1 |
Тем самым наше утверждение доказано.
Так само при 0<q<1 последовательность чисел q, q2 , q3 , … , qn , … стремится к нулю, то и погрешность an+ 1 стремится к нулю с возрастанием n. Иными словами, при указанных выше предположениях числа x1 , x2 , … , xn , … приближаются к числу x, причём разность |x-xn | убывает быстрее, чем qn ·|a1 |.
Точно так же можно доказать, что если на отрезке [a, b] выполнено неравенство:
|j ¢(x)|>1,
то процесс итераций расходится.
Особенно быстро сходится процесс последовательных приближений, если в точке x производная функции j(x) обращается в нуль. В этом случае по мере приближения к x, значение j ¢(x) стремится к нулю. Так как:
|an+ 1 |=|j ¢(cn )|·|an |
то сходимость процесса ускоряется по мере приближения к точке x.
Однако то же самое можно наблюдать в методе Ньютона, при замене f(x)=0 на имеем: и её производная: в точке x: f(x)=0 - в методе Ньютона наблюдается ускорение сходимости процесса приближений.
5. Метод касательных (метод Ньютона)
Метод касательных, связанный с именем И. Ньютона, является одним из наиболее эффективных численных методов решения уравнений. Идея метода очень проста. Возьмём производную точку x0 и запишем в ней уравнение касательной к графику функции f(x):
y=f(x0 )+ f ¢(x) (x-x0 ) (1.5)
Графики функции f(x) и её касательной близки около точки касания, поэтому естественно ожидать, что точка x1 пересечения касательной с осью Ox будет расположена недалеко от корня c (рис. 1.5)
Для определения точки имеем уравнение:
f(x0 )+ f ¢(x0 ) (x1 -x0 )=0
таким образом: x1 =x0 – f (x0 )/ f ¢(x0 ) (2.5)
Повторим проделанную процедуру: напишем уравнение касательной к графику функции f(x) при x=x1 и найдём для неё точку пересечения x2 с осью Ox (см. рис.1.5) x2 =x1 – f (x1 )/ f ¢(x1 ). Продолжая этот процесс, получим последовательность {xn }, определён- ную с помощью рекуррентной формулы:
xn + 1 =xn – f (xn )/ f ¢(xn ), n=0, 1, 2, … (3.5)