Что такое численное интегрирование
Определённый интеграл — это площадь под графиком функции на отрезке [a; b]. Первообразная есть не у каждой функции (например, у e−x² или sin(x)/x её нет в элементарных функциях), а у таблицы измерений формулы нет вовсе. Тогда отрезок делят на n частей шагом h = (b − a) / n и складывают площади простых фигур.
Формулы методов
- Прямоугольники (левые, правые, средние):
I ≈ h · Σ f(xi), где высота берётся в левой, правой точке или в середине каждого отрезка. Средние прямоугольники точнее левых и правых: их погрешность пропорциональна h², а не h. - Трапеции:
I ≈ h · ( (f0 + fn)/2 + f1 + … + fn−1 ), погрешность порядка h². - Симпсон (параболы):
I ≈ h/3 · ( f0 + 4f1 + 2f2 + 4f3 + … + 4fn−1 + fn ), n — чётное. Погрешность порядка h⁴: обычно самый точный из простых методов, а многочлены до третьей степени он интегрирует точно.
Оценка погрешности по правилу Рунге
Интеграл считается дважды — с шагом h и h/2. Тогда погрешность ε ≈ модуль(Ih/2 − Ih) / (2p − 1), где p — порядок метода (1 для левых и правых прямоугольников, 2 для средних и трапеций, 4 для Симпсона). Если погрешность велика — увеличьте n.
Пример
∫ x²·sin(x) dx от 0 до π. Точное значение π² − 4 ≈ 5,869604. При n = 10 метод трапеций даёт ≈ 5,788 (ошибка 0,08), а Симпсон — ≈ 5,8695 (ошибка около 0,0001): при том же числе точек он точнее почти в тысячу раз.
Интеграл по таблице
Если есть только измерения — расход по часам, мощность по времени, скорость по секундам, — переключитесь на «По таблице точек». Метод трапеций работает с любым шагом, Симпсон — только с равным шагом и нечётным числом точек. Так считают, например, объём воды по графику расхода или энергию по графику мощности.