Эйлеровы интегралы

2 раздела
от теории до практики
1 пример
Примеры решения задач
видео
Примеры решения задач
Содержание
  1. Гамма-функция Эйлера.
    Начать изучение
  2. Бета-функция Эйлера.
    Начать изучение

Гамма-функция Эйлера.

Определение.

Гамма-функция Эйлера \(\Gamma(x)\) (эйлеров интеграл второго рода) определяется как несобственный интеграл
$$
\Gamma(x) = \int\limits_{0}^{+\infty} t^{x-1}e^{-t} dt,\ x > 0,\label{ref1}
$$
с двумя особыми точками, \(t = 0\) и \(t = +\infty\).

Для того чтобы можно было применить теоремы об интегралах, зависящих от параметра, представим интеграл \eqref{ref1} в виде суммы двух интегралов
$$
\Gamma(x) = \int\limits_{0}^{1} t^{x-1}e^{-t} dt+\int\limits_{1}^{+\infty} t^{x-1}e^{-t} dt.\label{ref2}
$$
Оба интеграла сходятся равномерно по параметру \(x\) на любом конечном отрезке \([a, b] \subset (0, +\infty)\) по признаку Вейерштрасса. Действительно, пусть \(0 < a < 1\), \(b > 1\). Тогда \(0 \leq t^{x-1}e^{-t} \leq t^{a-1}\) при \(x \geq a\), \(0 \leq t \leq 1\) и \(\int\limits_{0}^{1} t^{a-1} dt = a^{-1}\). Следовательно, интеграл \(\displaystyle\int\limits_{0}^{1} t^{x-1}e^{-t} dt\) сходится равномерно на \([a, b]\).

Аналогично \(0 \leq t^{x-1}e^{-t} \leq t^{b-1}e^{-t}\) при \(x \leq b\), \(t \geq 1\), интеграл \(\displaystyle\int\limits_{1}^{+\infty} t^{x-1}e^{-t} dt\) сходится, а интеграл \(\displaystyle\int\limits_{1}^{+\infty} t^{x-1}e^{-t} dt\) сходится равномерно на \([a, b]\).

Так как подынтегральная функция \(t^{x-1}e^{-t}\) непрерывна при \(t > 0\), \(x > 0\), то в силу теоремы 4 отсюда оба интеграла в формуле \eqref{ref2} будут непрерывными функциями параметра \(x\) на произвольном отрезке \([a, b] \subset (0, +\infty)\), а поэтому \(\Gamma(x)\) есть непрерывная функция при \(x > 0\).

При \(x > 0\) функция \(\Gamma(x)\) непрерывно дифференцируема, причем
$$
\Gamma'(x) = \int\limits_{0}^{1} t^{x-1}e^{-t} \ln t\ dt+\int\limits_{1}^{+\infty} t^{x-1}e^{-t} \ln t\ dt = \int\limits_{1}^{+\infty} t^{x-1} \ln t\ e^{-t} dt.\label{ref3}
$$
Дифференцирование под знаком интеграла законно, так как оба интеграла в формуле \eqref{ref3} сходятся равномерно по параметру \(x\) на любом отрезке \([a, b] \subset (0, +\infty)\).

По индукции можно доказать, что \(\Gamma(x)\) есть бесконечно дифференцируемая функция при \(x > 0\) и
$$
\Gamma^{(n)}(x) = \int\limits_{0}^{+\infty} t^{x-1}e^{-t} (\ln t)^{n} dt.\label{ref4}
$$
В частности,
$$
\Gamma″(x) = \int\limits_{0}^{+\infty} t^{x-1}e^{-t} (\ln t)^{2} dt > 0.\nonumber
$$
Поэтому \(\Gamma(x)\) — выпуклая вниз функция при \(x > 0\) и имеет единственный положительный минимум. Нетрудно было бы показать, что формула \eqref{ref3} имеет место и для комплексных \(x\) при \(Re\ x > 0\), и поэтому \(\Gamma(x)\) есть регулярная функция комплексной переменной \(x\) в правой полуплоскости \(Re\ x > 0\).

Выведем теперь основное функциональное соотношение для гамма-функции. Пусть \(x > 0\). Интегрируя по частям, находим
$$
\begin{array}{cc}
& \Gamma(x+1) = \left.\int\limits_{0}^{+\infty} t^{x}e^{-t} dt = -e^{-t}t^{x}\right|_{0}^{+\infty}+x \int\limits_{0}^{+\infty} t^{x-1}e^{-t} dt = x\Gamma(x),\\
& \Gamma(x+1) = x\Gamma(x),\ x > 0.
\end{array}\label{ref5}
$$

Это и есть основное функциональное соотношение для гамма-функции, найденное Эйлером. На нем в значительной мере основана теория гамма-функции.

Прежде всего заметим, что если \(x \in (0, 1]\), то \(x+1 \in (1, 2]\). Поэтому, зная значения \(\Gamma(x)\) на промежутке \((0,1]\), можно при помощи формулы \eqref{ref5} найти значения \(\Gamma(x)\) на промежутке \((1,2]\), а следовательно, и на любом отрезке \([n, n+1]\), \(n = 1, 2, \ldots\) Это существенно облегчает вычисление значений гамма-функции.

Далее, формула \eqref{ref5} позволяет исследовать поведение \(\Gamma(x)\) при \(x \rightarrow +0\). Имеем
$$
\Gamma(x) = \frac{\Gamma(x+1)}{x} \sim \frac{\Gamma(1)}{x}\ \mbox{при}\ x \rightarrow +0,\nonumber
$$
так как \(\Gamma(x+1)\) — непрерывная функция при \(x = 0\), а \(\Gamma(1) = \displaystyle\int\limits_{0}^{+\infty} e^{-t} dt = 1\). Таким образом, \(\Gamma(x) \rightarrow +\infty\) при \(x \rightarrow +0\).

Из формулы \eqref{ref5} находим
$$
\Gamma(n+1) = n\Gamma(n) = n(n-1) \ldots 1 \cdot \Gamma(1) = n!.\nonumber
$$
Функция \(n!\) определена для натуральных \(n\). Гамма-функция \(\Gamma(x)\) непрерывна для всех \(x > 0\) и \(\Gamma(n+1) = n!\).

Формула \eqref{ref5} позволяет продолжить функцию \(\Gamma(x)\) с сохранением ее свойств на отрицательные значения \(x\), не равные \(-1, -2, \ldots, -n, \ldots\)

Положим по определению
$$
\Gamma(x) = \frac{\Gamma(x+1)}{x},\ -1 < x < 0.\label{ref6}
$$

Так как при \(x \in (-1, 0)\) имеем \(x+1 \in (-1, 0)\), то определение \eqref{ref6} корректно. Исследуем поведение \(\Gamma(x)\) при \(x \rightarrow -1+0\). Полагая \(y = x+1\), получаем, что \(x \rightarrow -1+0\) эквивалентно \(y \rightarrow +0\). Поэтому при \(y \rightarrow +0\), используя \eqref{ref6}, получаем
$$
\Gamma(y-1) = \frac{\Gamma(y)}{y-1} \sim -\Gamma(y) \sim -\frac{1}{y} = -\frac{1}{x+1}.\nonumber
$$
Итак,
$$
\Gamma(x) \sim -\frac{1}{x+1},\ x \rightarrow -1+0.\nonumber
$$

По индукции теперь можно определить \(\Gamma(x)\) на любом интервале \((-(n+1), -n)\), где \(n \in N\), формулой \(\Gamma(x) = \displaystyle\frac{\Gamma(x+1)}{x}\), \(x \in (-(n+1), -n)\), причем \(\Gamma(x) \sim \displaystyle\frac{(-1)^{n}}{x+n}\) при \(x \rightarrow -n\).

График \(\Gamma(x)\) изображен на рис. 73.1.

Рис. 73.1
Рис. 73.1

Бета-функция Эйлера.

Определение.

Интеграл
$$
B(x, y) = \int\limits_{0}^{1} t^{x-1}(1-t)^{y-1} dt,\label{ref7}
$$
зависит от двух параметров и называется бета-функцией Эйлера или эйлеровым интегралом первого рода.

У интеграла две особых точки, \(t = 0\) и \(t = 1\). Записывая интеграл \eqref{ref7} в виде
$$
B(x, y) = \int\limits_{0}^{1/2} t^{x-1}(1-t)^{y-1} dt+\int\limits_{1/2}^{1} t^{x-1}(1-t)^{y-1} dt,\nonumber
$$
получаем, что первый интеграл сходится при \(x > 0\), а второй при \(y > 0\), так что бета-функция определена при \(x > 0\), \(y > 0\).

Опишем свойства бета-функции.

Свойство 1.

\(B(x, y) = B(y, x)\).

Доказательство.

\(\circ\) Делая замену переменной \(\tau = 1-t\), получаем
$$
B(y, x) = \int\limits_{0}^{1} t^{y-1}(1-t)^{x-1} dt = \int\limits_{0}^{1} (1-\tau)^{y-1}\tau^{x-1}d\tau = B(x, y).\ \bullet\nonumber
$$

Свойство 2.

Справедливы формулы
$$
B(x, y) = \int\limits_{0}^{+\infty} \frac{u^{x-1}du}{(1+u)^{x+y}} = \int\limits_{0}^{1} \frac{u^{x-1}+u^{y-1}}{(1+u)^{x+y}}\ du.\label{ref8}
$$

Доказательство.

\(\circ\) Первая из формул \eqref{ref8} получается, если в интеграле \eqref{ref7} сделать замену переменной \(t = \displaystyle\frac{u}{1-u}\). Вторая формула \eqref{ref8} получается из первой, если разбить интеграл на два: по отрезку \([0,1]\) и интервалу \((1, +\infty\)), и во втором интеграле сделать замену переменной \(\displaystyle\frac{1}{u} = v\). \(\bullet\)

Свойство 3.

Справедлива формула
$$
B(x, 1-x) = \int\limits_{0}^{+\infty} \frac{u^{x-1}}{1+u}\ du = \frac{\pi}{\sin x\pi},\ 0 < x < 1.\label{ref9}
$$

Доказательство.

\(\circ\) Полагая в формуле \eqref{ref8} \(y = 1-x\) и пользуясь тождеством
$$
\frac{1}{1+u} = \sum_{\substack{k = 0}}^{\substack{n}}(-1)^{k}u^{k}+(-1)^{n+1} \frac{u^{n+1}}{1+u},\nonumber
$$
получаем
$$
B(x, 1-x) = \int\limits_{0}^{+\infty} \frac{u^{x-1}}{1+u}\ du = \int\limits_{0}^{1} \dfrac{u^{x-1}+u^{-x}}{1+u}\ du =\\= (-1)^{n+1} \int\limits_{0}^{1} \frac{u^{n+1}}{1+u} (u^{x-1}+u^{-x})\ du+\sum_{\substack{k = 0}}^{\substack{n}}(-1)^{k} \left(\frac{1}{k+x}+\frac{1}{k-x+1}\right).\label{ref10}
$$

Так как
$$
0 \leq \int\limits_{0}^{1} \dfrac{u^{n+1}(u^{x-1}+u^{-x})}{1+u}\ du \leq \int\limits_{0}^{1} (u^{n+x}+u^{n+1-x})\ du = \dfrac{1}{n+x+1}+\dfrac{1}{n+2-x},\nonumber
$$
то, переходя в формуле \eqref{ref10} к пределу при \(n \rightarrow +\infty\), получаем
$$
B(x, 1-x) = \sum_{\substack{k = 0}}^{\substack{n}}(-1)^{k} \left(\frac{1}{k+x}+\frac{1}{k-x+1}\right) =\\=\frac{1}{x}+\sum_{\substack{k = 0}}^{\substack{n}}(-1)^{k} \left(\frac{1}{k+x}+\frac{1}{x-k}\right) = \frac{\pi}{\sin \pi x}.\nonumber
$$

Последняя формула получается при \(z = \pi x\) из формулы, задающей разложение \(\displaystyle\frac{1}{\sin z}\) на элементарные дроби. \(\bullet\)

Свойство 4.

\(B(x, y)\) выражается через гамма-функцию, а именно
$$
B(x, y) = \frac{\Gamma(x)\Gamma(y)}{\Gamma(x+y)}.\label{ref11}
$$

Доказательство.

\(\circ\) В интеграле \(\Gamma(x) = \displaystyle\int\limits_{0}^{+\infty} t^{s-1}e^{-t} dt\) сделаем замену переменной \(t = (1+u)v\), \(v > 0\), и положим \(s = x+y\). Тогда
$$
\dfrac{\Gamma(x+y)}{(1+u)^{x+y}} = \int\limits_{0}^{+\infty} v^{x+y-1} e^{-(1+u)v} dv.\nonumber
$$
Умножим это равенство на \(u^{x-1}\) и проинтегрируем по \(u\) от 0 до \(+\infty\). В левой части, пользуясь формулой \eqref{ref8}, получим произведение \(\Gamma(x+y)B(x, y)\), а в правой — изменим порядок интегрирования. Тогда
$$
B(x, y)\Gamma(x+y) =\\= \int\limits_{0}^{+\infty} du \int\limits_{0}^{+\infty} u^{x-1} v^{x+y-1} e^{-(1+u)v} dv = \int\limits_{0}^{+\infty} dv \int\limits_{0}^{+\infty} u^{x-1} v^{x+y-1} e^{-v} e^{-uv} du.\nonumber
$$
Сделаем еще замену переменной \(uv = t\). Тогда
$$
B(x, y)\Gamma(x+y) = \int\limits_{0}^{+\infty} dv \int\limits_{0}^{+\infty} t^{x-1} v^{y-1} e^{-v} e^{-t} dt =\\= \int\limits_{0}^{+\infty} v^{y-1} e^{-v} dv \int\limits_{0}^{+\infty} t^{x-1} e^{-t} dt = \Gamma(x) \Gamma(y).\nonumber
$$

Обоснование изменения порядка интегрирования производится при помощи теоремы 7 отсюда аналогично тому, как это делалось в примере при вычислении интеграла вероятностей (интеграл Эйлера-Пуассона). \(\bullet\)

Следствие.

(Формула дополнения для гамма-функции).

При \(x \neq \pm k\), \(k = 0, 1, \ldots, n, \ldots\), справедлива формула
$$
\Gamma(x)\Gamma(1-x) = \frac{\pi}{\sin x\pi}.\label{ref12}
$$

\(\circ\) Используя \eqref{ref9} и \eqref{ref11}, получаем
$$
\frac{\pi}{\sin x\pi} = B(x, 1-x) = \dfrac{\Gamma(x)\Gamma(1-x)}{\Gamma(1)} = \Gamma(x)\Gamma(1-x).\nonumber
$$
При \(x = \frac{1}{2}\) из формулы дополнения следует, что \(\Gamma(\frac{1}{2}) = \sqrt{\pi}\). \(\bullet\)

Многие интегралы могут быть выражены через эйлеровы интегралы.

Пример 1.

Выразить через гамма-функцию следующий интеграл:
$$
I = \int\limits_{0}^{\pi/2} \sin^{a-1} \varphi \cos^{b-1} \varphi d\varphi,\ a > 0,\ b > 0.\nonumber
$$

Решение.

\(\vartriangle\) Сделаем замену переменной \(x = \sin^{2} \varphi\). Тогда
$$
dx = 2\sin \varphi \cos \varphi d\varphi,\ d\varphi = \frac{1}{2} x^{-1/2}(1-x)^{-1/2} dx,\nonumber
$$
$$
I = \frac{1}{2} \int\limits_{0}^{1} x^{(a/2)-1}(1-x)^{(b/2)-1} dx = \frac{1}{2} B\left(\frac{a}{2}, \frac{b}{2}\right) = \frac{1}{2} \dfrac{\displaystyle\Gamma\left(\frac{a}{2}\right)\Gamma\left(\frac{b}{2}\right)}{\displaystyle\Gamma\left(\frac{a+b}{2}\right)}.\nonumber
$$
Полагая \(a = 1+c\), \(b = 1-c\) и пользуясь формулой дополнения \eqref{ref12}, получаем при \(|c| < 1\) формулу
$$
\int\limits_{0}^{\pi/2} \operatorname{tg}^{c} \varphi\ d\varphi = \frac{1}{2} \Gamma\left(\frac{1}{2}+\frac{c}{2}\right)\Gamma\left(\frac{1}{2}-\frac{c}{2}\right) =\\= \frac{1}{2} \dfrac{\pi}{\displaystyle\sin \pi\left(\frac{c}{2}+\frac{1}{2}\right)} = \frac{1}{2} \dfrac{\pi}{\displaystyle\cos \frac{c\pi}{2}},\quad |c| < 1.\ \blacktriangle\nonumber
$$

Оставить комментарий