
- •«Омский государственный технический университет»
- •2. Основные этапы отыскания решения
- •Метод половинного деления
- •Метод простой итерации
- •Приведение уравнения к видудля обеспечения выполнения неравенства
- •Метод Ньютона (метод касательных)
- •6. Видоизменённый метод Ньютона
- •7. Метод хорд
- •8. Комбинированный метод
- •Решение систем линейных алгебраических уравнений
- •1. Постановка задачи
- •2. Метод простой итерации
- •3. Метод Зейделя
- •Решение систем нелинейных уравнений
- •1. Постановка задачи
- •2. Метод Ньютона для системы нелинейных уравнений
- •3. Метод итерации для нелинейной системы уравнений
- •4. Метод скорейшего спуска решения нелинейных систем
- •5. Метод скорейшего спуска для случая линейной системы
- •Приближение функций
- •1. Метод наименьших квадратов
- •2. Построение интерполяционных многочленов
- •Многочлен Лагранжа
- •Многочлен Ньютона с конечными разностями
- •Вычисление собственных значений матрицы методом Данилевского
- •Вычисление определённых интегралов
- •1. Метод Симпсона (метод парабол)
- •Численное решение дифференциальных уравнений
- •1. Постановка задачи Коши
- •2. Метод Эйлера
- •3. Модифицированные методы Эйлера
- •4. Метод Рунге – Кутта
- •Решение краевой задачи для линейного
- •Приложение
- •Библиографический список
Приведение уравнения к видудля обеспечения выполнения неравенства
В общем случае
получить подходящую итерационную форму
возможно, проведя равносильное
преобразование исходного уравнения,
например, умножив его на коэффициент
:
.
Прибавив затем к обеим частям уравнения
и обозначив
можно потребовать выполнения достаточного
условия
.
Отсюда определяется необходимое значение
.
Так как условие
должно выполняться на всем отрезке
,
то для выбора
следует использовать наибольшее
значение
на этом отрезке, т.е.
.
Это соотношение определяет диапазон
значений коэффициента
,
изменяющий величину
в пределах
.
Обычно принимают
.
На рис. 3- 6 показаны
четыре случая взаимного расположения
линий
и
и соответствующие итерационные процессы.
Рис. 3 и 4 соответствуют случаю
,
и итерационный процесс сходится. При
этом, если
(рис. 3), сходимость носит односторонний
характер, а если
(рис. 4), сходимость носит двусторонний,
колебательный характер. Рис. 5 и 6
соответствуют случаю
– итерационный процесс расходится. При
этом может быть односторонняя (рис. 5) и
двусторонняя (рис. 6) расходимость.
|
Рис. 3
Рис. 4
Рис. 5
Рис. 6
Погрешность метода. Оценка погрешности была доказана (5).
Критерий окончания.
Из оценки (5) следует, что вычисления
надо продолжать до выполнения неравенство
.
Если же
,
то оценка упрощается:
.
Пример 1. Используем
метод простой итерации для решения
уравнения
с точностью
.
Преобразуем уравнение к виду:
,
т. е.
.
Нетрудно убедиться,
что корень уравнения находится на
отрезке
.
Вычислив значения
на концах
отрезка, получим:
,
а
,
т. е. функция на концах отрезка имеет
разные знаки,
поэтому внутри отрезка есть корень. Расположение корня наглядно иллюстри-
рует рис. 7.
Рис. 7
Подсчитаем, первую
и вторую производные функции
:
.
Так как
на отрезке
,
то производная
монотонно возрастает на этом отрезке
и принимает максимальное значение на
правом конце отрезка, т. е. в точке
.
Поэтому, справедлива оценка:
.
Таким образом,
условие выполнено,
,
и можно воспользоваться критерием
окончания вычислений. В табл. 2 приведены
приближения, полученные по расчетной
формуле. В качестве начального приближения
выбрано значение
.
Таблица 2
-
0
1
2
3
4
5
1
0,8415
0,8861
0,8712
0,8774
0,8765
Критерий окончания
выполняется при
,
.
Сходимость двусторонняя, качественный
характер такой сходимости представлен
на рис. 4. Приближенное значение корня
с требуемой точностью
.
Пример 2.
Решить методом простой итерации уравнение
на отрезке
с точностью 0,025. Для решения исходное
уравнение приводится к виду
.
Для выбора величины
используем приведенную выше формулу
.
Тогда расчетная формула имеет вид
.
В качестве начального приближения можно
выбрать верхнюю границу заданного
отрезка
.
-
0
1
2
1
0,8
0,78
Так как
,
то
.