- •ВВЕДЕНИЕ
- •1 ЧИСЛЕННЫЕ МЕТОДЫ РЕШЕНИЯ УРАВНЕНИЙ
- •1.1 Методы отделения корней
- •1.1.1 Постановка задачи
- •1.1.2 Табличный метод отделения корней
- •1.1.3 Графический метод отделения корней
- •1.1.4 Метод интервалов отделения корней
- •1.2.2 Оценка погрешности приближенного корня
- •1.2.3 Метод половинного деления
- •1.2.3.1 Алгоритм метода половинного деления
- •1.2.4 Метод итераций
- •1.2.4.1 Алгоритм метода итераций
- •1.2.5 Метод Ньютона
- •1.2.5.1 Алгоритм метода Ньютона
- •1.2.6 Метод хорд
- •1.2.6.1 Алгоритм метода хорд
- •1.2.7 Комбинированный метод
- •1.2.7.1 Алгоритм комбинированного метода
- •1.2.8 Пример решения уравнения
- •1.2.8.1 Метод половинного деления
- •1.2.8.2 Метод итераций
- •1.2.8.3 Метод Ньютона
- •1.2.8.4 Метод хорд
- •1.2.8.5 Комбинированный метод
- •1.2.9 Уточнение корней уравнений в Excel с помощью циклической ссылки
- •1.2.9.1 Метод половинного деления
- •1.2.9.2 Метод итераций
- •1.2.9.3 Метод Ньютона
- •1.2.9.4 Метод хорд
- •1.2.9.5 Комбинированный метод
- •1.2.10 Решение уравнений средствами MathCAD
- •ПРИЛОЖЕНИЕ
13
Рис. 5. Пример, что при xn − xn-1 < ε имеет место неравенство ξ − xn > ε .
4. Упрощенная оценка (16) может быть использована для случая ϕ '(x) < 0 и |ϕ '(x)| < 1. Если же ϕ '(x) > 0 и ϕ '(x) < 1, то, вообще говоря, оценку (16) использовать нельзя. Но учитывая, что в этом случае последовательные приближения
xn = ϕ(xn–1), n = 1, 2, …, x0 [a; b]
сходится к корню ξ монотонно, можно организовать сужение интервала [a; b], последовательно получая приближения корня с недостатком и избытком. Тогда возможно применение оценки приближения (16).
1.2.4.1 Алгоритм метода итераций
Вычислительная схема решений имеет вид
xi = x0, x0 [a; b]
xi+1 = ϕ(xi), i = 1, 2, …
и сводится к организации циклического процесса с неизвестным числом повторений до выполнения заданного уровня точности, который может быть оценен неравенством (12) или применением упрощенной оценки (16).
Выполним уточнение корня на отрезке [a; b], применив следующий алгоритм.
1. Установить значения a, b, ε – границы отрезка отделения корня и принятую точность приближения.
2. Определить точку (η, M = max| f ' ( x )|).
x [α ; β ]
3.Установить начальное значение x0 [a; b].
4.Вычислить λ = sgn f '(η) M2 .
5.Начать цикл уточнения корня.
Вычислить очередное приближение x1 = x0 – λf(x0). Вычислить d = |x1 – x0|.
Принять x0 = x1.
6.Конец цикла, если d < ε .
7.Вывод результата x1.
8.Конец алгоритма.
Замечание. Выполнение п.4 алгоритма обеспечивает удовлетворение требований ϕ '(x) < 0 и |ϕ '(x)| < 1, при которых может быть применена упрощенная оценка (16).