- •Численные методы
- •Введение
- •1. Решение нелинейных уравнений
- •Постановка задачи
- •1.2. Основные этапы отыскания решения
- •1.3. Метод половинного деления
- •1.4. Метод простой итерации
- •Приведение уравнения к виду для обеспечения выполнения неравенства
- •1.5. Метод Ньютона (метод касательных)
- •1.6. Видоизменённый метод Ньютона
- •1.7. Метод хорд
- •1.8. Комбинированный метод
- •Решение систем линейных алгебраических уравнений
- •2.1. Постановка задачи
- •2.2. Метод простой итерации
- •2.3. Метод Зейделя
- •3. Решение систем нелинейных уравнений
- •3.1. Постановка задачи
- •3.2. Метод Ньютона для системы нелинейных уравнений
- •3.3. Метод итерации для нелинейной системы уравнений
- •3.4. Метод скорейшего спуска решения нелинейных систем
- •3.5. Метод скорейшего спуска для случая линейной системы
- •4. Приближение функций
- •4. 1. Метод наименьших квадратов
- •4.2. Построение интерполяционных многочленов
- •Многочлен Лагранжа
- •Многочлен Ньютона с конечными разностями
- •5. Вычисление собственных значений матрицы Методом Данилевского
- •6. Вычисление определённых интегралов. Метод симпсона (метод парабол)
- •7. Численное решение дифференциальных уравнений
- •7.1. Постановка задачи Коши
- •7.2. Метод Эйлера
- •7.3. Модифицированные методы Эйлера
- •7.4. Метод Рунге – Кутта
- •Решение краевой задачи для линейного дифференциального уравнения второго порядка методом прогонки
- •Приложение
- •Библиографический список
- •Оглавление
- •3.2. Метод Ньютона для системы нелинейных уравнений 26
- •3.3. Метод итерации для нелинейной системы уравнений 30
6. Вычисление определённых интегралов. Метод симпсона (метод парабол)
Заменим график функции на отрезке , , , параболой, проведенной через точки , , где – середина отрезка . Эта парабола есть интерполяционный многочлен второй степени с узлами . Нетрудно убедиться, что уравнение этой параболы имеет вид:
,
где .
Проинтегрировав эту функцию на отрезке , получим
.
Суммируя полученные выражение по , получим квадратурную формулу Симпсона (или формулу парабол):
.
Оценка погрешности. Для оценки погрешности формулы Симпсона воспользуемся следующей теоремой.
Теорема. Пусть функция имеет на отрезке непрерывную производную четвертого порядка . Тогда для формулы Симпсона справедлива следующая оценка погрешности: , где
.
Замечание. Если число элементарных отрезков, на которые делится отрезок , четно, т.е. , то параболы можно проводить через узлы с целыми индексами, и вместо элементарного отрезка длины рассматривать отрезок длины . Тогда формула Симпсона примет вид: , а вместо последней оценки будет справедлива следующая оценка погрешности:
.
Правило Рунге практической оценки погрешности.
Оценка погрешности зависит от длины элементарного отрезка , и при достаточно малом справедливо приближенное равенство: , где приближенное значение интеграла. Если уменьшить шаг в два раза, то получим: .
Вычитая одно из другого, получим: , или .
Это приближенное равенство дает оценку погрешности. Вычисление этой оценки называется правилом Рунге. Правило Рунге – это эмпирический способ оценки погрешности, основанный на сравнении результатов вычислений, проводимых с разными шагами . Для формулы Симпсона , и оценка принимает вид: . Используя правило Рунге, можно построить процедуру приближенного вычисления интеграла с заданной точностью . Нужно, начав вычисления с некоторого значения шага , последовательно уменьшать это значения в два раза, каждый раз вычисляя приближенное значение . Вычисления прекращаются тогда, когда результаты двух последующих вычислений будут различаться меньше, чем на .
Пример. Вычислить .
Решение. Возьмём , тогда .
|
|
|
|
|
0 |
0 |
|
|
10)0,54) |
1 |
0,125 |
0,984625 |
|
|
2 |
0,250 |
|
0,9411761) |
|
3 |
0,375 |
0,876712 |
|
|
4 |
0,5 |
|
0,82) |
|
5 |
0,625 |
0,7191 |
|
|
6 |
0,750 |
|
0,643) |
|
7 |
0,875 |
0,566389 |
|
|
8 |
1 |
|
|
|
|
|
3,45955 |
1,62818 |
1,5 |
.
.
.
Следовательно, значение интеграла можно счесть .