Неполная бета-функция

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

Определение

Неполная бета-функция для действительных параметров \( a > 0 \) и \( b > 0 \) и переменной \( x \in [0, 1] \) определяется как:

\[ B(x; a, b) = \int_0^x t^{a-1} (1 - t)^{b-1} \, dt \]

При \( x = 1 \) она совпадает с полной бета-функцией:

\[ B(1; a, b) = B(a, b) = \int_0^1 t^{a-1} (1 - t)^{b-1} \, dt = \frac{\Gamma(a)\Gamma(b)}{\Gamma(a+b)} \]

где \(\Gamma(z)\) — гамма-функция.

Функция определена для всех \( a, b > 0 \), а при \( a = 0 \) или \( b = 0 \) интеграл расходится, если \( x \) не равно 0 или 1 соответственно.

Свойства

Основные свойства

• Непрерывность и монотонность: при фиксированных \( a, b > 0 \) функция \( B(x; a, b) \) непрерывна и строго возрастает по \( x \) на отрезке \([0, 1]\).
• Граничные значения: \( B(0; a, b) = 0 \), \( B(1; a, b) = B(a, b) \).
• Симметрия: \( B(x; a, b) = B(1 - x; b, a) \) (свойство, вытекающее из замены переменной \( t \to 1 - t \)).

Связь с гамма-функцией

Неполная бета-функция выражается через неполную гамма-функцию \(\gamma(s, x)\) при помощи преобразования, но в общем случае не сводится к элементарным функциям. Однако для целых значений \( a \) и \( b \) она может быть выражена через полиномы.

Дифференцирование

Производная по верхнему пределу:

\[ \frac{\partial}{\partial x} B(x; a, b) = x^{a-1} (1 - x)^{b-1} \]

Разложение в ряд

Для \( 0 \le x < 1 \) неполная бета-функция может быть разложена в степенной ряд:

\[ B(x; a, b) = \frac{x^a}{a} \sum_{n=0}^{\infty} \frac{(1-b)_n}{(a+1)_n} x^n \]

где \((q)_n\) — символ Похгаммера (растущий факториал). Этот ряд сходится при \(|x| < 1\).

Регуляризованная неполная бета-функция

На практике чаще используется регуляризованная неполная бета-функция, определяемая как:

\[ I_x(a, b) = \frac{B(x; a, b)}{B(a, b)} \]

Она представляет собой функцию распределения бета-распределения с параметрами \( a \) и \( b \). Свойства:

• \( I_0(a, b) = 0 \), \( I_1(a, b) = 1 \).
• \( I_x(a, b) = 1 - I_{1-x}(b, a) \) (свойство симметрии).
• Для целых \( a \) и \( b \) выражается через биномиальные суммы.

Частные случаи

• При \( a = 1 \): \( B(x; 1, b) = \frac{1 - (1 - x)^b}{b} \).
• При \( b = 1 \): \( B(x; a, 1) = \frac{x^a}{a} \).
• При \( a = b = 1 \): \( B(x; 1, 1) = x \).
• При \( a = 0.5, b = 0.5 \): функция связана с арксинусом: \( I_x(0.5, 0.5) = \frac{2}{\pi} \arcsin(\sqrt{x}) \).

Применение

Статистика и теория вероятностей

Регуляризованная неполная бета-функция является кумулятивной функцией распределения (CDF) бета-распределения, которое используется для моделирования случайных величин, ограниченных отрезком \([0, 1]\), например, вероятностей, долей и концентраций.

Она также применяется для вычисления вероятностей в биномиальном и отрицательном биномиальном распределениях. Например, функция распределения биномиального распределения с параметрами \( n \) и \( p \) выражается через \( I_{1-p}(n - k, k + 1) \).

Математический анализ и специальные функции

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

\[ B(x; a, b) = \frac{x^a}{a} \, {}_2F_1(a, 1-b; a+1; x) \]

Физика и инженерные науки

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

Комбинаторика

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

Вычисление

Для численного вычисления неполной бета-функции и её регуляризованной версии используются различные алгоритмы:

• Ряды и цепные дроби: наиболее эффективный метод — представление \( I_x(a, b) \) в виде непрерывной дроби Ленца (Lentz's continued fraction), которая сходится быстро для большинства значений параметров.
• Интегрирование: для малых \( a \) и \( b \) может использоваться прямое численное интегрирование, например, методом Гаусса-Лежандра.
• Специализированные библиотеки: в большинстве математических пакетов (MATLAB, Mathematica, SciPy, R) есть встроенные функции для вычисления \( I_x(a, b) \). В языке Python функция betainc доступна в модуле scipy.special.

История

Понятие неполной бета-функции было введено в XVIII веке Леонардом Эйлером в контексте изучения бета-интеграла. Систематическое исследование её свойств и приложений было проведено в XIX—XX веках, в частности, в работах Карла Фридриха Гаусса, Адриена Мари Лежандра и Пьера-Симона Лапласа. В современной математической статистике функция стала широко использоваться после работ Джорджа Бокса и других статистиков в середине XX века.

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

• Неполная бета-функция является частным случаем неполной гамма-функции при замене переменной, но не сводится к ней в общем виде.
• Регуляризованная неполная бета-функция \( I_x(a, b) \) является обратной к функции распределения бета-распределения, что позволяет вычислять квантили.
• В некоторых контекстах (например, в физике) используется нестандартное обозначение: \( B_x(a, b) \) вместо \( B(x; a, b) \).

Источники

• Абрамовиц М., Стиган И. Справочник по специальным функциям. — М.: Наука, 1979.
• Бейтмен Г., Эрдейи А. Высшие трансцендентные функции. Том 1. — М.: Наука, 1973.
• Градштейн И. С., Рыжик И. М. Таблицы интегралов, сумм, рядов и произведений. — М.: Физматлит, 1963.
• Press W. H., Teukolsky S. A., Vetterling W. T., Flannery B. P. Numerical Recipes: The Art of Scientific Computing. — 3rd ed. — Cambridge University Press, 2007.