- •31. Теория вероятностей
- •31.1. Основные понятия теории вероятностей.
- •31.2. Теоремы сложения и умножения вероятностей
- •31.3. Формула полной вероятности. Формула Байеса
- •31.4. Повторение испытаний
- •31.5. Дискретные случайные величины. Функция
- •31.6. Непрерывные случайные величины
- •31.7. Основные законы распределения
- •31.8. Функция одного случайного аргумента
- •31.9. Функция двух случайных аргументов
- •31.10. Закон больших чисел
- •32. Математическая статистика
- •32.1. Основные понятия математической
- •32.2. Точечные и интервальные оценки
- •32.3. Статистическая проверка
- •32.4. Критерии согласия Пирсона и Колмогорова
- •32.5. Элементы теории корреляции. Выборочное
- •33. Вычислительная математика
- •33.1. Элементы теории погрешностей
- •33.2. Аппроксимация функций
- •33.3. Приближенное решение нелинейных уравнений
- •33.4. Приближенное вычисление интегралов
- •33.5. Приближенное решение обыкновенных
- •1. Таблица значений функции
- •2. Таблица значений функции лапласа
- •3. Распределение пуассона
- •4. Таблица значений
- •5. Таблица значений
- •6. Критические точки распределения фишера–снедекора
- •7. Критические точки распределения χ2
- •8. Критические точки распределения стьюдента
- •9. Критические значения распределения колмогорова
- •Содержание
33.4. Приближенное вычисление интегралов
Пусть требуется вычислить определенный интеграл

Интегрирование с помощью рядов
Предположим,
что подынтегральная функция
разлагается в степенной ряд
(33.19)
сходящийся
в интервале
который содержит отрезок интегрирования
Тогда искомый интегралI
можно представить в виде числового ряда
(33.20)
Используем
непрерывную аппроксимацию для функции
и в качестве аппроксимирующей функции
возьмемn-ю
частичную сумму ряда (33.19):

Если
ряд (33.20) знакочередующийся, то для оценки
его остаточного
члена
можно воспользоваться известным
следствием
признака Лейбница. В случае знакопостоянного
ряда (33.20) обычно применяется мажорирование
бесконечно убывающей геометрической
прогрессией.
Численное интегрирование
Задача численного интегрирования заключается в вычислении значения определенного интеграла I на основании ряда значений подынтегральной функции в точках промежутка интегрирования. Формулы численного интегрирования называются квадратурными формулами. Рассмотрим простейшие из них, использующие локальную интерполяцию при точечной аппроксимации подынтегральной функции.
С
помощью точек
разобьем отрезок интегрирования
наn
равных частичных отрезков

сшагом
интегрирования
и обозначим
Величина
где
– численное значение выбранной
квадратурной формулы, называетсяостаточным
членом
квадратурной формулы.
Формула левых прямоугольников:

Формула правых прямоугольников:

Если
непрерывна на
то для формул левых и правых прямоугольников
выполняется


Формула средних прямоугольников:
(33.21)
Если
непрерывна на
то для формулы средних прямоугольников
выполняется

(33.22)
Формула трапеций:
(33.23)
Если
непрерывна на
то для формулы трапеций выполняется

(33.24)
Формула Симпсона (формула парабол):
(33.25)
где

Если
непрерывна на
то для формулы Симпсона выполняется

(33.26)
Пример
1. Вычислить
интеграл
с точностью до
путем разложения подынтегральной
функции в степенной ряд.
Решение. Воспользуемся формулой Маклорена


Заменяя
в ней x
на
получим ряд

Этот
ряд сходится при любом
поэтому его можно почленно интегрировать.
Тогда

Полученный
числовой ряд является знакочередующимся.
Так как

но
то отсюда имеем оценку остаточного
члена:

Поэтому для обеспечения заданной степени точности достаточно взять два члена ряда

Следовательно,

Пример
2. Вычислить
интеграл
разлагая подынтегральную функцию в
степенной ряд и используя три члена
этого разложения. Определить абсолютную
погрешность результата.
Решение. Воспользуемся формулой


Получим


Найдем абсолютную погрешность



т.
е.
В окончательном результате оставим
шесть значащих цифр:
Тогда

т.
е.
Следовательно,
Пример
3. Вычислить
с помощью формулы средних прямоугольников
интеграл
приняв
Определить абсолютную погрешность
результата.
Решение.
По заданным
пределам интегрирования
и
найдем шаг
В данном случае точками разбиения
отрезка
будут
и
Тогда по формуле (33.21), где
получим


Находим
абсолютную погрешность метода
Так как

то, используя формулу (33.22), имеем

т. е.
В окончательном результате оставим
четыре значащие цифры:
Тогда

т. е.
Следовательно,
Для
сравнения приведем точное значение
интеграла:

Пример
4. Вычислить
с помощью формулы трапеций интеграл
с точностью до
Решение.
Для
определения необходимого числа разбиений
n
оценим остаточный член
по формуле (33.24). Так как

то

Решая неравенство

получим
примем
Зная
пределы интегрирования, найдем шаг
Точками разбиения служат


Найдем соответствующие значения
подынтегральной функции



Тогда по формуле (33.23) получаем

Следовательно,

Для
сравнения приведем точное значение
интеграла:

Пример
5. Вычислить
с помощью формулы Симпсона интеграл
с точностью до
Решение.
Для
определения n
оценим остаточный член
по формуле (33.26). Можно показать (проверить
самостоятельно!), что

Поэтому

Решая неравенство

получим
примем
Нужно
определить значения подынтегральной
функции при
для следующих значений аргумента:



Найдем соответствующие значения




Тогда по формуле (33.25) получим


Следовательно,

