Численное интегрирование: методы и погрешности¶

Численное интегрирование — раздел вычислительной математики, посвящённый приближённому вычислению определённых интегралов в случаях, когда точное значение либо недоступно (отсутствие первообразной в элементарных функциях), либо нецелесообразно из-за сложности аналитического решения. Методы численного интегрирования основаны на замене подынтегральной функции более простой аппроксимирующей функцией (обычно полиномом), интеграл от которой вычисляется точно.

¶Постановка задачи

Требуется вычислить значение определённого интеграла:

\[ I = \int_{a}^{b} f(x)\,dx \]

где \(f(x)\) — непрерывная функция на отрезке \([a, b]\). Численные методы дают приближённое значение \(I_n\), которое зависит от количества узлов разбиения \(n\). Основная цель — получить результат с заданной точностью при минимальных вычислительных затратах.

¶Основные методы

¶Метод прямоугольников

Простейший метод, при котором отрезок \([a, b]\) разбивается на \(n\) равных частей с шагом \(h = (b-a)/n\), а подынтегральная функция заменяется кусочно-постоянной. Различают метод левых, правых и средних прямоугольников в зависимости от того, в какой точке (левом, правом конце или середине подынтервала) вычисляется значение функции.

Формула метода средних прямоугольников:

\[ I_n = h \sum_{i=0}^{n-1} f\left(x_i + \frac{h}{2}\right) \]

Метод средних прямоугольников имеет второй порядок точности, тогда как методы левых и правых прямоугольников — только первый.

¶Метод трапеций

Подынтегральная функция заменяется ломаной линией, соединяющей значения функции в узлах сетки. Формула:

\[ I_n = h \left( \frac{f(a) + f(b)}{2} + \sum_{i=1}^{n-1} f(x_i) \right) \]

Метод трапеций имеет второй порядок точности. Геометрически он означает замену площади под кривой суммой площадей трапеций.

¶Метод Симпсона (парабол)

Более точный метод, при котором подынтегральная функция на каждом двойном подынтервале заменяется параболой, проходящей через три точки. Требует чётного числа подынтервалов \(n\):

\[ I_n = \frac{h}{3} \left( f(a) + f(b) + 4\sum_{i=1,3,5,\dots}^{n-1} f(x_i) + 2\sum_{i=2,4,6,\dots}^{n-2} f(x_i) \right) \]

Метод Симпсона имеет четвёртый порядок точности и является одним из наиболее употребительных на практике.

¶Квадратурные формулы Гаусса

Методы, в которых узлы выбираются не равномерно, а специальным образом (как корни ортогональных полиномов Лежандра), что позволяет достичь максимальной алгебраической степени точности. Формула Гаусса с \(n\) узлами точна для полиномов степени до \(2n-1\) включительно. Наиболее известна формула Гаусса-Лежандра, применяемая на отрезке \([-1, 1]\) с последующим линейным преобразованием к произвольному отрезку.

¶Погрешности методов

¶Остаточный член

Погрешность численного интегрирования определяется остаточным членом \(R_n = I - I_n\). Для основных методов оценки имеют вид:

• метод прямоугольников (средних): \(R_n = \frac{(b-a)h^2}{24} f''(\xi)\)
• метод трапеций: \(R_n = -\frac{(b-a)h^2}{12} f''(\xi)\)
• метод Симпсона: \(R_n = -\frac{(b-a)h^4}{180} f^{(4)}(\xi)\)

где \(\xi\) — некоторая точка отрезка \([a, b]\), а \(f''\) и \(f^{(4)}\) — вторая и четвёртая производные функции.

¶Практическая оценка погрешности

Для оценки погрешности на практике применяют правило Рунге. Если известно приближённое значение \(I_n\) с шагом \(h\) и \(I_{2n}\) с шагом \(h/2\), то для метода порядка \(p\) уточнённая оценка:

\[ I \approx I_{2n} + \frac{I_{2n} - I_n}{2^p - 1} \]

Это правило лежит в основе адаптивных алгоритмов, автоматически выбирающих шаг интегрирования для достижения требуемой точности.

¶Сравнение методов

Методы Гаусса при том же числе вычислений функции дают существенно более высокую точность, однако требуют хранения таблиц узлов и весов и менее удобны при адаптивном уточнении.

¶Области применения

Численное интегрирование применяется в физике (вычисление работы, моментов инерции), механике (центры тяжести), гидродинамике, статистике (функция распределения нормального закона), при решении дифференциальных уравнений, в обработке экспериментальных данных, когда функция задана таблично и аналитическое интегрирование невозможно. Для несобственных интегралов и интегралов с особенностями применяются специальные модификации: замена переменных, выделение особенности, формулы Гаусса-Лагерра и Гаусса-Эрмита.