
- •Численные методы (Вычислительная математика) Санкт-Петербург
- •Оглавление
- •Глава 1. Элементарная теория погрешностей........................10
- •Глава 2.Численные методы решения уравнений……………24
- •Глава 3. Приближенное представление функций.................50
- •Предисловие
- •Введение
- •Глава 1. Элементарная теория погрешностей
- •1.1. Приближенные числа и их погрешности
- •1.2. Правила приближенных вычислений и оценка погрешностей при вычислениях
- •1.3. Неустранимая погрешность вычисления значения функции
- •1.4. Обеспечение заданной точности результата
- •Глава 2.Численные методы решения уравнений
- •2.1. Итерационные методы решения нелинейных уравнений
- •2.2. Метод половинного деления
- •2.3. Метод хорд
- •2.4. Метод касательных (метод Ньютона).
- •2.5. Комбинированный метод хорд и касательных
- •2.6.Метод итераций
- •2.7. Решение систем линейных уравнений методом итераций.
- •2.9. Метод итераций для нелинейных систем уравнений
- •Глава 3. Приближенное представление функций (аппроксимация, интерполяция)
- •3.1. Интерполирование функции алгебраическими многочленами
- •3.2. Единственность интерполяционного многочлена
- •3.3. Интерполяционная формула Лагранжа
- •3.4. Погрешность интерполяционной формулы Лагранжа
- •3.5. Погрешность линейной интерполяции
- •3.6.Конечные разности и их свойства
- •3.7.Интерполяционная формула Ньютона (для равностоящих узлов)
- •3.10. Эмпирические формулы и метод наименьших квадратов
- •Глава 4. Численное дифференцирование и интегрирование
- •4.1. Дифференцирование на основе интерполяционной формулы Ньютона
- •4.2. Дифференцирование на основе формулы Стирлинга
- •4.3. Дифференцирование на основе интерполяционной формулы Лагранжа
- •Глава 5. Численное интегрирование функций
- •5.1. Квадратурные формулы Ньютона-Котеса
- •5.2. Формула трапеций
- •5.3. Формула Симпсона (формула парабол)
- •5.4. Формулы Ньютона – Котеса высших порядков
- •Глава 6. Численное интегрирование обыкновенных дифференциальных уравнений
- •6.1. Метод Эйлера
- •6.2. Метод Рунге – Кутта
- •6.3. Метод Адамса
- •Глава 7. Численные методы оптимизации
- •7.1. Градиентные методы
- •7.2. Метод Франка – Вулфа
- •7.3. Метод наискорейшего спуска
- •7.4. Метод штрафных функций
- •Литература
2.2. Метод половинного деления
Простейшим итерационным методом нахождения корней нелинейных уравнений является метод дихотомии (метод бисекций, половинного деления). Его алгоритм содержит следующие этапы:
Найти отрезок
, в котором распложен корень уравнения.
В качестве начального приближения взять
.
Вычислить значения функции на концах отрезков
и
:
.
Выбрать отрезок, на концах которого значения функции имеют разные знаки, и использовать его в качестве нового отрезка для следующей итерации.
В результате каждой
итерации отрезок, на котором расположен
корень, уменьшается вдвое, после n
итераций он сокращается в
раз. Точность (погрешность) вычислений
задается некоторым малым числом
.
Процесс вычислений продолжается до
тех пор, пока значение функции
после n-й
итерации не станет меньшим по модулю
этого заданного малого числа
:
,
длина полученного отрезка становится
меньше погрешности. Метод дихотомии
медленно сходится к точному решению,
однако он сходящийся.
Пример. Отделить действительные корни аналитически и уточнить их методом дихотомии с точностью до 0,001.
Обозначим
Найдем производную
Составим таблицу знаков
|
-∞ |
-2 |
0 |
1 |
+∞ |
Signf(x) |
- |
- |
- |
- |
+ |
Существует одна перемена знака функции, поэтому уравнение имеет один действительный корень на промежутке [1;+∞). Уменьшаем промежуток, содержащий корень так, чтобы его длина была равна единице.
|
1 |
2 |
Sign f(x) |
- |
+ |
Следовательно, корень необходимо искать на промежутке [1;2].
Все вычисления приведены в таблице.
|
|
|
|
|
||
0 |
1,0000 |
2,0000 |
1,5000 |
|
-3,3125 |
|
1 |
1,5000 |
2,0000 |
1,7500 |
|
7,8242 |
|
2 |
1,5000 |
1,7500 |
1,6250 |
|
1,3953 |
|
3 |
1,5000 |
1,6250 |
1,5625 |
|
-1,1567 |
|
4 |
1,5625 |
1,6250 |
1,5938 |
|
0,0677 |
|
5 |
1,5625 |
1,5938 |
1,5781 |
|
-0,5571 |
|
6 |
1,5781 |
1,5938 |
1,5859 |
|
-0,2479 |
|
7 |
1,5859 |
1,5938 |
1,5898 |
|
-0,0909 |
|
8 |
1,5898 |
1,5938 |
1,5918 |
|
-0,0118 |
|
9 |
1,5918 |
1,5938 |
1,5928 |
|
0,0279 |
|
10 |
1,5918 |
1,5928 |
1,5923 |
|
0,0080 |
|
11 |
1,5918 |
1,5923 |
1,5920 |
|
-0,0019 |
|
12 |
1,5920 |
1,5923 |
1,5922 |
|
0,0031 |
|
13 |
1,5920 |
1,5922 |
1,5921 |
|
0,0006 |
|
14 |
1,5920 |
1,5921 |
1,5921 |
|
-0,0006 |
|
15 |
1,5921 |
1,5921 |
1,5921 |
|
0,0005 |
Найден корень