- •Глава 8. Численное интегрирование
- •8.1. Постановка задачи численного интегрирования. Квадратурные формулы Ньютона – Котеса Постановка задачи приближенного интегрирования. Квадратурные формулы
- •Квадратурные формулы Ньютона – Котеса
- •Необобщенные формулы прямоугольников, трапеций, Симпсона
- •Обобщенные формулы прямоугольников, трапеций, Симпсона
- •Оценка погрешности и определение порядка точности для формулы трапеций
- •Оценка погрешности и определение порядка точности обобщенной формулы трапеций. Преобразуем разность между точным и приближенным значениями интеграла, вычисленным с помощью обобщенной формулы трапеций:
- •Оценки погрешности и порядки точности обобщенных квадратурных формул прямоугольников, трапеций и Симпсона
- •8.2. Квадратурные формулы Гаусса. Метод неопределенных коэффициентов Квадратурные формулы Гаусса
- •Метод неопределенных коэффициентов
- •8.3. Метод Монте-Карло
- •Первая схема метода Монте-Карло
- •Вторая схема метода Монте-Карло
- •8.4. Вычисление первообразных, несобственных и кратных интегралов
- •Вычисление первообразных
- •Вычисление несобственных интегралов
- •Вычисление кратных интегралов. Кубатурные формулы
- •8.5. Обусловленность квадратурных формул интерполяционного типа
- •Контрольные вопросы и задания
8.5. Обусловленность квадратурных формул интерполяционного типа
При вычислении интегралов, как правило, приходится использовать не точные значения подынтегральной функции, а приближенные значения. В предположении, чтодля всех, справедлива оценка
, |
|
указывающая на то, что абсолютное число обусловленности этой задачи равно , то есть длине отрезка интегрирования.
Какова же чувствительность квадратурной формулы (8.1.2) к погрешностям задания функции ? Отметим, что
, |
|
Таким образом, квадратурная формула устойчива к ошибкам задания функции и ее число обусловленности равно.
Все квадратурные формулы интерполяционного типа точны для многочленов нулевой степени, и поэтому . Следовательно, если все весаквадратурной формулы интерполяционного типа положительны, то ее число обусловленностисовпадает с. Чувствительность такой формулы к ошибкам адекватна чувствительности вычисляемого интеграла. Если же среди весовимеются отрицательные, то
.
Известно, что при больших значениях n среди весов квадратурной формулы (8.1.6) появляются отрицательные и значение числа обусловленности становится большим. Например, для формул Ньютона–Котесаприn =10, приn =20, приn =30. Из-за плохой обусловленности эти формулы уже при используются весьма редко.
Таким образом, если подынтегральная функция вычисляется с погрешностью, то погрешность вычисления определенного интеграла складывается из погрешности квадратурной формулы, которую можно уменьшить, увеличив числоn разбиений интервала интегрирования, и неустранимой погрешности, зависящей от погрешности вычисления подынтегральной функции .
Контрольные вопросы и задания
Как ставится задача численного интегрирования? Что такое квадратурные формулы? Как получаются квадратурные формулы Ньютона-Котеса?
Получите формулы прямоугольников, трапеций, Симпсона (простые и обобщенные). Каков их геометрический смысл?
Получите оценку погрешности формулы трапеций (простой и обобщенной).
Найдите порядок точности формулы трапеций (простой и обобщенной).
Запишите оценки погрешности и порядки точности обобщенных формул прямоугольников, трапеций и Симпсона. Как используется эта информация для вычисления интеграла с заданной точностью?
Получите квадратурные формулы Гаусса. Запишите оценку погрешности.
В чем смысл метода неопределенных коэффициентов? Приведите пример его использования.
Опишите первую схему метода Монте-Карло.
Опишите вторую схему метода Монте-Карло.
Как вычисляются первообразные функций и несобственные интегралы с помощью формул численного интегрирования?
Как вычисляются кратные интегралы с помощью кубатурных формул?
Как влияют погрешности вычисления подынтегральной функции на погрешность квадратурной формулы? Из каких составляющих складывается погрешность вычисления интеграла?