Численное решение нелинейных уравнений

Пусть дано нелинейное уравнение

Численное решение нелинейных уравнений - student2.ru , (20)

корни которого необходимо найти.

Допустим, что известен интервал [a,b], на котором непрерывная функция f(x) меняет знак (рис. 6). В этом случае можно считать, что известно приближенное решение уравнения x=(a+b)/2 с погрешностью, оценкой которой служит полуширина интервала неопределенности D(0)=|b-a|/2.

Численное решение нелинейных уравнений - student2.ru

Рис. 6. Решение нелинейного уравнения

Метод бисекций

Метод бисекций (дихотомии, половинного деления) заключается в вычислении значения функции в середине отрезка [a,b], сравнении его знака со знаком функции, например, в точке а и отбрасывании той половины отрезка [a,b], на концах которой функция имеет одинаковые знаки. Далее это повторяется до тех пор, пока оценка погрешности D(k)=|b(k)-a(k)|/2 не станет меньше заданного числа e.

Алгоритм решения задачи может выглядеть следующим образом:

Численное решение нелинейных уравнений - student2.ru

Метод простых итераций

Приведем (20) к виду

Численное решение нелинейных уравнений - student2.ru , (21)

например, выбрав Численное решение нелинейных уравнений - student2.ru в виде: Численное решение нелинейных уравнений - student2.ru .

Выберем начальное приближение x(0). Тогда можно организовать итерационный процесс по правилу

Численное решение нелинейных уравнений - student2.ru ,

Численное решение нелинейных уравнений - student2.ru ,

………………

Численное решение нелинейных уравнений - student2.ru , (22)

………………

Условие сходимости итерационного процесса записывается в виде неравенства Численное решение нелинейных уравнений - student2.ru (предполагая существование производной Численное решение нелинейных уравнений - student2.ru в области поиска).

Процесс сходимости приближенного решения к точному иллюстрируется на рис. 7,а.

Условие сходимости a<1 является существенным, т.е. при a³1 имеет место расходимость (см. рис. 7,б).

Численное решение нелинейных уравнений - student2.ru Численное решение нелинейных уравнений - student2.ru
а б

Рис. 7. Решение нелинейного уравнения методом простых итераций:

а сходимость при Численное решение нелинейных уравнений - student2.ru , б расходимость при Численное решение нелинейных уравнений - student2.ru

Метод Ньютона (касательных)

Допустим, что функция f(x) непрерывно дифференцируема. Метод Ньютона заключается в построении касательной к графику функции y=f(x) в точке x=x(0), определению точки пересечения касательной с осью абсцисс, которая считается следующим приближением x(1), построению касательной в этой точке и т.д.
(рис. 8).

Используем уравнение касательной, проведенной в точке Численное решение нелинейных уравнений - student2.ru

Численное решение нелинейных уравнений - student2.ru (23)

и условие пересечения графиком касательной оси абсцисс (y=0), найдем

Численное решение нелинейных уравнений - student2.ru . (24)

Численное решение нелинейных уравнений - student2.ru

Рис.8. Метод Ньютона

Метод хорд (секущих)

Недостаток метода Ньютона связан с необходимостью вычисления производной. Однако производную можно приближенно вычислять конечно-разностным способом. На графике это можно изобразить, отметив на кривой y=f(x) две точки при x=x(0), x=x(1) и проведя через них секущую до пересечения с осью абсцисс (рис. 9). Если x(0)и x(1) расположены с разных сторон от корня x*, то ось абсцисс пересекает хорда, соединяющая две точки кривой (рис. 9).

Уравнение секущей, проходящей через две точки (x(0),f(x(0))), (x(1),f(x(1)))

Численное решение нелинейных уравнений - student2.ru .

Точка пересечения секущей с осью абсцисс (y=0) является следующим приближением решения уравнения

Численное решение нелинейных уравнений - student2.ru

Численное решение нелинейных уравнений - student2.ru . (25)

Численное решение нелинейных уравнений - student2.ru

Рис. 9. Метод секущих (хорд)

Наши рекомендации