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

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

Суть метода

Идея метода заключается в разбиении области интегрирования \([a, b]\) на \(n\) равных отрезков (шагов) длиной \(h = \frac{b-a}{n}\). На каждом элементарном отрезке \([x_i, x_{i+1}]\) подынтегральная функция \(f(x)\) заменяется прямой, проходящей через точки \((x_i, f(x_i))\) и \((x_{i+1}, f(x_{i+1}))\). Площадь под графиком функции на этом отрезке приближённо принимается равной площади трапеции с основаниями \(f(x_i)\) и \(f(x_{i+1})\) и высотой \(h\).

Суммируя площади всех элементарных трапеций, получают приближённое значение определённого интеграла:

\[ \int_a^b f(x) \, dx \approx \frac{h}{2} \left[ f(x_0) + 2 \sum_{i=1}^{n-1} f(x_i) + f(x_n) \right] \]

где \(x_0 = a\), \(x_n = b\), \(x_i = a + i \cdot h\) для \(i = 1, 2, \dots, n-1\).

Геометрическая интерпретация

Геометрически метод трапеций представляет собой замену площади под кривой \(y = f(x)\) на отрезке \([a, b]\) суммой площадей \(n\) прямоугольных трапеций, каждая из которых опирается на отрезок \([x_i, x_{i+1}]\) оси абсцисс. В отличие от метода прямоугольников, где используется ступенчатая аппроксимация, метод трапеций даёт более точное приближение, так как линейная интерполяция лучше отражает поведение гладкой функции.

История

Метод трапеций известен с древних времён как один из простейших способов приближённого вычисления площадей. В античной математике, в частности в трудах Архимеда, использовались близкие идеи для вычисления площадей фигур, ограниченных кривыми. Формальное математическое обоснование и систематическое применение метода получили в работах математиков XVII—XVIII веков, таких как Исаак Ньютон и Джеймс Грегори, в рамках развития численного анализа и квадратурных формул.

Погрешность метода

Метод трапеций является точным для линейных функций (полиномов первой степени). Для произвольной дважды непрерывно дифференцируемой функции \(f(x)\) погрешность \(R\) на одном элементарном отрезке \([x_i, x_{i+1}]\) оценивается формулой:

\[ R_i = -\frac{h^3}{12} f''(\xi_i), \quad \xi_i \in (x_i, x_{i+1}) \]

Глобальная погрешность на всём отрезке \([a, b]\) имеет порядок \(O(h^2)\):

\[ |R| \leq \frac{(b-a)h^2}{12} \max_{x \in [a,b]} |f''(x)| \]

Таким образом, уменьшение шага \(h\) (увеличение числа разбиений \(n\)) приводит к квадратичному уменьшению погрешности. Однако на практике чрезмерное увеличение \(n\) может привести к накоплению ошибок округления при вычислениях.

Сравнение с другими методами

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

Метод прямоугольников (левых, правых или средних) имеет погрешность порядка \(O(h)\) для левых и правых прямоугольников и \(O(h^2)\) для средних прямоугольников. Метод трапеций обычно точнее левых и правых прямоугольников, но уступает методу средних прямоугольников при одинаковом числе разбиений, так как константа в оценке погрешности для метода средних прямоугольников в два раза меньше.

Метод Симпсона

Метод Симпсона (парабол) использует квадратичную аппроксимацию подынтегральной функции и имеет погрешность порядка \(O(h^4)\). Для гладких функций он значительно точнее метода трапеций при том же числе разбиений, но требует большего объёма вычислений.

Метод Рунге — Кутты

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

Применение

Метод трапеций широко применяется в инженерных и научных расчётах, когда:

• подынтегральная функция задана таблично (например, результаты эксперимента);
• аналитическое вычисление интеграла невозможно или затруднительно;
• требуется быстрая оценка интеграла с приемлемой точностью.

Метод прост в реализации и не требует вычисления производных, что делает его удобным для использования в программировании. Он реализован во многих математических пакетах (MATLAB, Octave, SciPy и др.) и часто используется как базовая квадратурная формула.

Пример

Вычислить приближённо интеграл \(\int_0^1 x^2 \, dx\) методом трапеций при \(n = 4\).

1. Разбиение: \(a = 0\), \(b = 1\), \(n = 4\), \(h = \frac{1-0}{4} = 0.25\).
2. Узлы: \(x_0 = 0\), \(x_1 = 0.25\), \(x_2 = 0.5\), \(x_3 = 0.75\), \(x_4 = 1\).
3. Значения функции: \(f(0) = 0\), \(f(0.25) = 0.0625\), \(f(0.5) = 0.25\), \(f(0.75) = 0.5625\), \(f(1) = 1\).
4. Вычисление по формуле:

\[ S \approx \frac{0.25}{2} \left[ 0 + 2(0.0625 + 0.25 + 0.5625) + 1 \right] = 0.125 \cdot (0 + 1.75 + 1) = 0.125 \cdot 2.75 = 0.34375 \]

Точное значение интеграла: \(\frac{1}{3} \approx 0.33333\). Погрешность составляет \(0.01042\).

Модификации

Составная формула трапеций

Стандартная реализация для равномерной сетки, описанная выше.

Адаптивный метод трапеций

Автоматическое изменение шага разбиения в зависимости от локальной оценки погрешности. Позволяет сгущать сетку в областях быстрого изменения функции и разрежать — на гладких участках, что повышает эффективность.

Правило Рунге

Используется для апостериорной оценки погрешности и уточнения результата (экстраполяция Ричардсона). Применяется, когда известно, что погрешность имеет порядок \(O(h^2)\).

Интересные факты

• Метод трапеций является частным случаем квадратурных формул Ньютона — Котеса при \(k=1\).
• В некоторых учебниках метод трапеций называют «правилом трапеций» или «формулой трапеций».
• Для монотонных функций метод трапеций даёт систематическую ошибку: если функция выпукла вверх, интеграл занижается; если выпукла вниз — завышается.
• В российской учебной литературе метод трапеций часто изучается в курсе численных методов наряду с методом прямоугольников и методом Симпсона.

Источники

• Бахвалов Н. С., Жидков Н. П., Кобельков Г. М. Численные методы. — М.: Бином. Лаборатория знаний, 2008.
• Самарский А. А., Гулин А. В. Численные методы. — М.: Наука, 1989.
• Калиткин Н. Н. Численные методы. — М.: Наука, 1978.
• Формальное определение и оценка погрешности — материалы курса «Численные методы» МГУ им. М. В. Ломоносова.