Метод средней точки

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

Суть метода

Пусть требуется вычислить определённый интеграл: \[ I = \int_a^b f(x) \, dx. \]

Отрезок интегрирования \([a, b]\) разбивается на \(n\) равных частей (элементарных отрезков) длиной \(h = \frac{b-a}{n}\). Точки разбиения: \(x_i = a + i h\), где \(i = 0, 1, \dots, n\). На каждом отрезке \([x_{i-1}, x_i]\) выбирается средняя точка \(c_i = \frac{x_{i-1} + x_i}{2} = a + \left(i - \frac12\right) h\).

Площадь под графиком функции на отрезке \([x_{i-1}, x_i]\) заменяется площадью прямоугольника с основанием \(h\) и высотой, равной значению функции в средней точке: \[ \int_{x_{i-1}}^{x_i} f(x) \, dx \approx h \cdot f(c_i). \]

Суммируя по всем отрезкам, получаем приближённое значение интеграла: \[ I \approx I_{\text{ср}} = h \sum_{i=1}^n f\left(a + \left(i - \frac12\right) h\right). \]

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

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

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

Для функции \(f(x)\), дважды непрерывно дифференцируемой на отрезке \([a, b]\), погрешность метода средней точки на одном элементарном отрезке \([x_{i-1}, x_i]\) оценивается как: \[ \left| \int_{x_{i-1}}^{x_i} f(x) \, dx - h f(c_i) \right| \le \frac{h^3}{24} \max_{[x_{i-1}, x_i]} |f''(x)|. \]

Глобальная погрешность на всём отрезке \([a, b]\) при \(n\) отрезках разбиения: \[ |I - I_{\text{ср}}| \le \frac{(b-a) h^2}{24} \max_{[a,b]} |f''(x)|. \]

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

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

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

• Метод левых прямоугольников: высота прямоугольника берётся в левой точке отрезка \([x_{i-1}, x_i]\), то есть \(f(x_{i-1})\).
• Метод правых прямоугольников: высота берётся в правой точке \(f(x_i)\).
• Метод средних прямоугольников (метод средней точки): высота берётся в середине отрезка \(f(c_i)\).

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

Пример вычисления

Вычислим интеграл \(\int_0^1 x^2 \, dx\) методом средней точки с \(n = 4\) (шаг \(h = 0.25\)).

Средние точки: \(c_1 = 0.125\), \(c_2 = 0.375\), \(c_3 = 0.625\), \(c_4 = 0.875\).

Значения функции: \(f(0.125) = 0.015625\), \(f(0.375) = 0.140625\), \(f(0.625) = 0.390625\), \(f(0.875) = 0.765625\).

Сумма: \(0.015625 + 0.140625 + 0.390625 + 0.765625 = 1.3125\).

Умножаем на \(h = 0.25\): \(I_{\text{ср}} = 0.328125\).

Точное значение интеграла: \(\int_0^1 x^2 \, dx = \frac13 \approx 0.333333\). Погрешность: \(0.005208\).

Для сравнения, метод левых прямоугольников с тем же разбиением даёт \(0.21875\) (погрешность \(0.114583\)), метод правых — \(0.46875\) (погрешность \(0.135417\)).

Применение

Метод средней точки применяется в задачах, где требуется быстрое приближённое интегрирование с умеренной точностью. Он используется:

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

Метод также лежит в основе более сложных квадратурных формул, таких как составные формулы Гаусса — Лежандра, где узлы выбираются не равномерно, а по определённым правилам для повышения точности.

Реализация в вычислительной технике

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

Пример реализации на псевдокоде: `` function midpoint_rule(f, a, b, n): h = (b - a) / n sum = 0 for i from 1 to n: x_mid = a + (i - 0.5) h sum = sum + f(x_mid) return h sum ``

Ограничения

Метод средней точки имеет ограничения, общие для всех методов прямоугольников:

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

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

Историческая справка

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

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

• Метод средней точки является частным случаем квадратурной формулы Гаусса — Лежандра с одним узлом, где узел расположен в середине отрезка.
• В некоторых учебниках метод называют «формулой средних прямоугольников» или «методом центральных прямоугольников».
• Для линейных функций метод даёт точное значение интеграла, так как площадь прямоугольника со средней высотой в точности равна площади трапеции.

Источники

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