
- •Глава 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
разбиений интервала интегрирования, и
неустранимой погрешности, зависящей
от погрешности вычисления подынтегральной
функции
.
Контрольные вопросы и задания
Как ставится задача численного интегрирования? Что такое квадратурные формулы? Как получаются квадратурные формулы Ньютона-Котеса?
Получите формулы прямоугольников, трапеций, Симпсона (простые и обобщенные). Каков их геометрический смысл?
Получите оценку погрешности формулы трапеций (простой и обобщенной).
Найдите порядок точности формулы трапеций (простой и обобщенной).
Запишите оценки погрешности и порядки точности обобщенных формул прямоугольников, трапеций и Симпсона. Как используется эта информация для вычисления интеграла с заданной точностью?
Получите квадратурные формулы Гаусса. Запишите оценку погрешности.
В чем смысл метода неопределенных коэффициентов? Приведите пример его использования.
Опишите первую схему метода Монте-Карло.
Опишите вторую схему метода Монте-Карло.
Как вычисляются первообразные функций и несобственные интегралы с помощью формул численного интегрирования?
Как вычисляются кратные интегралы с помощью кубатурных формул?
Как влияют погрешности вычисления подынтегральной функции
на погрешность квадратурной формулы? Из каких составляющих складывается погрешность вычисления интеграла?