
- •Глава 8. Численное интегрирование
- •8.1. Постановка задачи численного интегрирования. Квадратурные формулы Ньютона – Котеса Постановка задачи приближенного интегрирования. Квадратурные формулы
- •Квадратурные формулы Ньютона – Котеса
- •Необобщенные формулы прямоугольников, трапеций, Симпсона
- •Обобщенные формулы прямоугольников, трапеций, Симпсона
- •Оценка погрешности и определение порядка точности для формулы трапеций
- •Оценка погрешности и определение порядка точности обобщенной формулы трапеций. Преобразуем разность между точным и приближенным значениями интеграла, вычисленным с помощью обобщенной формулы трапеций:
- •Оценки погрешности и порядки точности обобщенных квадратурных формул прямоугольников, трапеций и Симпсона
- •8.2. Квадратурные формулы Гаусса. Метод неопределенных коэффициентов Квадратурные формулы Гаусса
- •Метод неопределенных коэффициентов
- •8.3. Метод Монте-Карло
- •Первая схема метода Монте-Карло
- •Вторая схема метода Монте-Карло
- •8.4. Вычисление первообразных, несобственных и кратных интегралов
- •Вычисление первообразных
- •Вычисление несобственных интегралов
- •Вычисление кратных интегралов. Кубатурные формулы
- •8.5. Обусловленность квадратурных формул интерполяционного типа
- •Контрольные вопросы и задания
Вычисление несобственных интегралов
Для вычисления несобственных интегралов существует несколько методов. Мы ограничимся рассмотрением только двух видов одномерных несобственных интегралов и одного из методов их приближенного вычисления.
Интегралы
с бесконечным пределом интегрирования.
Пусть требуется вычислить приближенное
значение несобственного интеграла
с погрешностью, не превышающей заданного
положительного числа
.
Согласно определению,
.
(8.4.9)
Поэтому
в качестве приближенного значения
несобственного интеграла можно выбрать
значение
при достаточно большом значении величиныb.
Осталось найти это значение.
Представим искомый несобственный интеграл в виде суммы:
.
(8.4.10)
Выберем значение величины b исходя из требования, чтобы
,
(8.4.11)
а
затем используем для вычисления
приближенного значения определенного
интеграла
какую-нибудь квадратурную формулу.
Обозначим черезQ
приближенное значение этого определенного
интеграла, вычисленное с помощью
квадратурной формулы с погрешностью,
не превышающей
.
ТогдаQ
будет также представлять собой и
приближенное значение несобственного
интеграла
с погрешностью,
не превышающей
.
В самом деле,
(8.4.12)
Пример 1
Требуется
найти приближенное значение несобственного
интеграла
с погрешностью, не превышающей
.
Оценим сверху левую часть неравенства (8.4.11):
.
Для
выполнения неравенства (8.4.11) достаточно
потребовать, чтобы выполнялось
неравенство
.
Таким образом, мы можем выбрать
.
Осталось вычислитьQ
– приближенное значение определенного
интеграла
с погрешностью, не превышающей
.
Это можно сделать с помощью любой из
квадратурных формул, как показано выше
(мы не будем этого делать). Полученное
значениеQ
можно использовать в качестве искомого
приближенного значения несобственного
интеграла с погрешностью, не превышающей
.
Интегралы от
неограниченных функций.
Пусть имеется функция
,
не ограниченная в окрестности точки
и такая, что существует несобственный
интеграл
.
Требуется
вычислить приближенное значение
несобственного интеграла
с погрешностью,
не превышающей заданного положительного
числа
.
Согласно определению,
.
(8.4.13)
Поэтому
в качестве приближенного значения
несобственного интеграла можно выбрать
значение
при значении величиныc,
большем a
и достаточно близком к a.
Осталось найти это значение.
Представим искомый несобственный интеграл в виде суммы
.
(8.4.14)
Выберем значение величины c исходя из требования, чтобы
,
(8.4.15)
а
затем используем для вычисления
приближенного значения определенного
интеграла
какую-нибудь квадратурную формулу.
Обозначим черезQ
приближенное значение этого определенного
интеграла, вычисленное с помощью
квадратурной формулы с погрешностью,
не превышающей
.
ТогдаQ
будет также представлять собой и
приближенное значение несобственного
интеграла
с погрешностью,
не превышающей
.В самом деле,
(8.4.16)
Пример 2
Требуется
найти приближенное значение несобственного
интеграла
с погрешностью, не превышающей
.
Подынтегральная
функция не ограничена в окрестности
точки
.
Оценим сверху левую часть неравенства
(8.4.15):
.
Для
выполнения неравенства (8.4.15) достаточно
потребовать, чтобы выполнялось неравенство
.
Таким образом, мы можем выбрать
.
Осталось вычислитьQ
– приближенное значение определенного
интеграла
с погрешностью, не превышающей
.
Это можно сделать с помощью любой из
квадратурных формул, как показано выше
(мы не будем этого делать). Полученное
значениеQ
можно использовать в качестве искомого
приближенного значения несобственного
интеграла с погрешностью, не превышающей
.