Формула Чебышёва

Формула Чебышёва — это математическое выражение, используемое для приближённого вычисления определённых интегралов. Относится к классу квадратурных формул интерполяционного типа, в которых узлы выбираются таким образом, чтобы минимизировать погрешность при заданном числе точек. Разработана русским математиком Пафнутием Львовичем Чебышёвым в XIX веке.

История

Формула Чебышёва была предложена П. Л. Чебышёвым в 1873 году в работе «О квадратурах» (опубликована в «Записках Императорской Академии наук»). Чебышёв стремился создать метод численного интегрирования, который был бы прост в вычислениях и давал бы точные результаты для многочленов высокой степени. В отличие от формулы Гаусса, где узлы и веса подбираются для точного интегрирования многочленов степени \(2n-1\), Чебышёв предложил использовать равные веса для всех узлов, что упрощает расчёты, но требует специального расположения узлов. Изначально формула была разработана для интегралов на отрезке \([-1, 1]\), но затем обобщена на произвольные отрезки.

Определение

Формула Чебышёва для приближённого вычисления определённого интеграла \(\int_{a}^{b} f(x) \, dx\) имеет вид:

\[ \int_{a}^{b} f(x) \, dx \approx \frac{b-a}{n} \sum_{i=1}^{n} f(x_i), \]

где \(n\) — количество узлов, \(x_i\) — узлы, выбираемые из условия, чтобы формула была точна для многочленов степени \(n\) (или \(n-1\), в зависимости от варианта). Узлы \(x_i\) являются корнями многочлена Чебышёва первого рода \(T_n(x)\), преобразованного к отрезку \([a, b]\). Для стандартного отрезка \([-1, 1]\) узлы определяются как:

\[ x_i = \cos\left(\frac{2i-1}{2n}\pi\right), \quad i = 1, 2, \dots, n. \]

Свойства

Точность

Формула Чебышёва точна для многочленов степени не выше \(n\) (если \(n\) чётно) или \(n-1\) (если \(n\) нечётно). Это означает, что для многочленов указанной степени погрешность равна нулю. Для произвольных функций погрешность зависит от гладкости функции и количества узлов.

Погрешность

Остаточный член формулы Чебышёва для \(n\) узлов на отрезке \([a, b]\) выражается как:

\[ R_n = \frac{(b-a)^{n+2}}{(n+1)!} \cdot \frac{f^{(n+1)}(\xi)}{n} \cdot \int_{a}^{b} \prod_{i=1}^{n} (x - x_i) \, dx, \]

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

Узлы

Узлы формулы Чебышёва не являются равномерно распределёнными на отрезке. Они сгущаются к краям отрезка, что уменьшает эффект Рунге — явление, при котором интерполяция многочленами высокой степени на равномерной сетке приводит к большим колебаниям на концах отрезка. Это делает формулу Чебышёва более устойчивой для функций с особенностями на границах.

Виды и модификации

Классическая формула Чебышёва

Использует узлы, являющиеся корнями многочлена Чебышёва первого рода. Применяется для отрезка \([-1, 1]\). Для произвольного отрезка \([a, b]\) узлы преобразуются линейно:

\[ x_i^{(a,b)} = \frac{a+b}{2} + \frac{b-a}{2} \cos\left(\frac{2i-1}{2n}\pi\right). \]

Формула Чебышёва с равными весами

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

Обобщённая формула Чебышёва

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

Применение

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

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

Пример

Вычислим приближённо интеграл \(\int_{0}^{1} e^x \, dx\) с помощью формулы Чебышёва для \(n=3\). Точное значение: \(e^1 - e^0 \approx 1.71828\).

1. Преобразуем узлы с отрезка \([-1, 1]\) на \([0, 1]\):

• Для \(n=3\) узлы на \([-1, 1]\): \(x_1 = \cos(\pi/6) = \sqrt{3}/2 \approx 0.8660\), \(x_2 = \cos(3\pi/6) = 0\), \(x_3 = \cos(5\pi/6) = -\sqrt{3}/2 \approx -0.8660\).
• Преобразование: \(x_i' = \frac{0+1}{2} + \frac{1-0}{2} x_i = 0.5 + 0.5 x_i\).
• Получаем: \(x_1' \approx 0.5 + 0.5 \cdot 0.8660 = 0.9330\), \(x_2' = 0.5\), \(x_3' \approx 0.5 + 0.5 \cdot (-0.8660) = 0.0670\).

1. Вычисляем сумму значений функции:

• \(f(x_1') = e^{0.9330} \approx 2.542\), \(f(x_2') = e^{0.5} \approx 1.6487\), \(f(x_3') = e^{0.0670} \approx 1.0693\).
• Сумма: \(2.542 + 1.6487 + 1.0693 = 5.260\).

1. Приближённое значение интеграла: \(\frac{1-0}{3} \cdot 5.260 \approx 1.7533\).

Погрешность: \(|1.71828 - 1.7533| \approx 0.035\), что составляет около 2%.

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

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

Критика

Основной недостаток формулы Чебышёва — ограниченное число узлов, для которых узлы являются действительными числами. Для \(n > 7\) некоторые узлы становятся комплексными, что делает формулу неприменимой для вещественных функций. Это ограничение было выявлено самим Чебышёвым и его учениками. Впоследствии были разработаны модификации, например, формулы с усечением или использованием других многочленов. Кроме того, для гладких функций метод Гаусса часто даёт более высокую точность при том же числе узлов, хотя и требует больше вычислительных затрат.

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

• Формула Чебышёва является частным случаем квадратурных формул с равными весами, которые изучались ещё в Древней Греции (например, метод Архимеда для вычисления площади круга).
• Пафнутий Чебышёв разработал эту формулу в рамках своих исследований по теории приближений и механизмам, в частности, для расчёта шарнирных механизмов.
• В честь Чебышёва названы не только квадратурные формулы, но и многочлены, неравенства и другие математические объекты.

Источники

• Чебышёв П. Л. «О квадратурах» // Записки Императорской Академии наук, 1873.
• Березин И. С., Жидков Н. П. «Методы вычислений», том 1, 1966.
• Крылов В. И. «Приближённое вычисление интегралов», 1967.
• Демидович Б. П., Марон И. А. «Основы вычислительной математики», 1970.
• Формальский А. М. «Численные методы», 2006.