§7. Дифференцирование несобственного интеграла по параметру |
451 |
Таким образом, второй из интегралов правой части (3) стремится к 0 при A → ∞ и, переходя к пределу при A → ∞, выводим, что второй из несобственных интегралов в (2) существует и равен первому несобственному интегралу в (2). 
§7. Дифференцирование несобственного интеграла по параметру
ТЕОРЕМА 7.1. Пусть f(x, y) — функция двух переменных, определенная на [a, ∞) × [c, d] и непрерывная по x на [a, ∞) для любого y [c, d]. Предположим, что существует производная fy(x, y) непрерывная по совокупности пе-
ременных на полуполосе [a, ∞) × [c, d], и, кроме того, интеграл
∞ |
|
I(y) = a |
f(x, y)dx |
сходится при всех y [c, d], а интеграл |
∞ |
|
Φ(y) = a |
fy(x, y)dx |
сходится равномерно относительно y [c, d]. Тогда для всех y [c, d] справедлива формула:
I (y) = Φ(y).
Доказательство. Для всех y [c, d] по теореме 5.1 имеем:
y |
∞ |
|
∞ |
y |
|
∞ |
c |
dy a |
fy(x, y)dx = a |
dx c |
fy(x, y)dy = a |
[f(x, y) − f(x, c)]dx = |
|
|
∞ |
|
|
∞ |
|
|
|
|
= a |
f(x, y)dx − a |
f(x, c)dx. |
Дифференцируя обе части данного равенства по y находим
∞ |
|
|
∞ |
|
|
fy(x, y)dx = |
d |
|
f(x, y)dx. |
|
dy |
452 Глава 18. Несобственные интегралы, зависящие от параметра
Здесь мы воспользовались также непрерывностью левого интеграла как функции y, вытекающей из теоремы 5.1. 
УПРАЖНЕНИЕ 1. Соответствующие утверждения о несобственных интегралах с конечными пределами интегрирования (непрерывность по параметру, изменение порядка интегрирования, дифференцирование по параметру) сформулируйте самостоятельно.
ПРИМЕР 1 (интеграл вероятностей). Вычислим несобственный интеграл
∞
e−x2 dx.
−∞
Если попытаться вычислить этот интеграл путем поиска первообразной, то данный путь тупиковый, поскольку первообразная в элементарных функциях не выражается. Прежде всего заметим, что
∞ ∞
e−x2 dx = 2 e−x2 dx = 2I.
Положим x = ut, dx = udt, где u > 0 — параметр. Тогда |
∞ |
|
I = 0 |
e−u2t2 udt. |
Умножая обе части равенства на e−u2 du и интегрируя, полу- |
чаем: |
|
|
|
|
|
|
|
∞Ie−u2 du = ∞e−u2 ∞e−u2t2 udt du. |
0 |
|
0 |
|
|
0 |
|
Однако |
∞ |
|
∞ |
|
|
|
|
0 |
Ie−u2 du = I |
0 |
e−u2 du = I2. |
Таким образом, получаем |
|
|
|
|
|
I2 = ∞ue−u2 |
∞e−u2t2 dt du = ∞du ∞ue−u2(t2+1)dt. (1) |
0 |
|
0 |
|
|
|
0 |
0 |
§7. Дифференцирование несобственного интеграла по параметру |
453 |
Если бы удалось нам поменять местами порядок интегрирования в правой части (1), то мы бы пришли к равенству
|
|
|
|
|
∞ ∞ |
|
|
|
|
|
|
I2 = |
|
dt |
ue−u2(t2+1)du. |
(2) |
|
|
|
|
|
0 |
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
Далее заметим, что |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∞ue−u2(t2+1)du = u2(t2 + 1) = τ, |
dτ = 2u(t2 + 1)du = |
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∞ |
|
|
|
|
|
|
|
|
|
|
|
|
∞ |
|
|
|
e−τ dτ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
= 0 |
|
|
|
= |
1 |
|
|
|
|
|
0 |
e−τ dτ = |
|
1 |
. |
2(t2 + 1) |
|
|
2(t2 + 1) |
2(t2 + 1) |
Пользуясь (2) выводим |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∞ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
I2 = 0 |
|
|
|
|
dt |
1 |
|
arctg t|0∞ = |
π |
|
|
|
|
|
= |
|
|
|
|
. |
|
|
|
|
2(t2 + 1) |
2 |
4 |
|
Отсюда |
|
|
|
|
|
|
|
|
√ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
I = |
π |
, |
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
и, тем самым,
∞
e−x2 dx = √π.
−∞
Нам остается обосновать переход от (1) к (2), связанный с изменением порядка интегрирования.
Попробуем воспользоваться теоремой 6.1, что кажется наиболее естественным. Проверим выполнение ее условий. Обо-
значим
f(t, u) = e−u2(t2+1)u.
Во-первых, эта функция непрерывна как функция двух переменных на [0, ∞) × [0, ∞).
Во-вторых, интеграл
∞
0 |
2 |
2 |
|
1 |
1 |
ue−u |
(t |
+1)du = |
|
|
|
|
2 |
1 + t2 |
непрерывен по переменной t.
454 Глава 18. Несобственные интегралы, зависящие от параметра
Однако интеграл
∞ |
|
u2 |
(t2+1) |
|
u2 |
∞ |
|
u2t2 |
dt = |
e−u2 I, |
при u > 0, |
0 |
ue− |
|
|
dt = e− |
|
u 0 |
e− |
|
0, при |
u = 0. |
является разрывной функцией в точке u = 0. Поэтому воспользоваться теоремой 6.1 напрямую нельзя. Применим эту теорему к множеству [u0, ∞) × [0, ∞) и перейдем к пределу при u0 → +0. В результате получаем
∞
e−x2 dx = √π.
−∞
Заметим, что существование одного из повторных интегралов, что необходимо для применения теоремы 6.1, нами было доказано (мы его просто посчитали). 

§8. Интеграл Эйлера первого рода
ОПРЕДЕЛЕНИЕ 8.1. B("Бета")-функцией (интегралом Эйлера первого рода) называется функция
1
B(a, b) = xa−1(1 − x)b−1dx,
0
где a > 0, b > 0.
УПРАЖНЕНИЕ 1. Докажите, что данный интеграл сходится тогда и только тогда, когда a > 0 и b > 0.
ЗАМЕЧАНИЕ. При a < 1, либо при b < 1, данный интеграл является несобственным интегралом, зависящим от параметра.
ТЕОРЕМА 8.1. Для любых a > 0, b > 0 справедливо равенство
B(a, b) = B(b, a).
Доказательство.
1
B(a, b) = xa−1(1 − x)b−1dx = |1 − x = y, −dx = dy| =
0 1
= − (1 − y)a−1yb−1dy = yb−1(1 − y)a−1dy = B(b, a).
1 0
|
ТЕОРЕМА 8.2. Для любых a > 0, |
b > 1 справедливо |
|
равенство |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
B(a, b) = |
|
b − 1 |
|
B(a, b 1). |
|
|
|
|
|
|
|
|
|
|
|
|
|
a + b |
− |
1 |
|
|
|
|
− |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Доказательство. |
|
|
|
|
|
− |
|
|
|
|0 |
|
|
|
0 |
|
|
− |
0 |
1 |
|
− |
|
|
|
a |
a |
|
|
|
|
|
a |
|
|
|
B(a, b) = (1 x)b−1d( |
xa |
) = |
xa |
(1 x)b−1 1 + |
b − 1 |
|
xa(1 x)b−2dx = |
|
|
|
|
|
|
|
|
|
|
|
a |
0 |
|
|
|
− |
|
|
|
|
− |
|
|
− |
|
|
|
|
|
|
|
|
b − 1 |
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
= |
|
xa−1(1 x)b−2(1 (1 x))dx = |
|
|
|
|
|
|
|
a |
|
|
|
|
|
− |
|
|
− |
|
|
a |
|
|
|
|
|
− |
|
|
|
b − 1 |
1 |
|
|
|
|
|
|
|
|
|
|
b − 1 |
|
1 |
|
|
|
|
|
|
|
= |
|
|
xa−1(1 x)b−2dx |
|
|
|
|
|
xa−1(1 x)b−1dx = |
|
|
|
|
|
|
|
|
B(a, b) = b −a 1B(a, b − 1),
456 |
Глава 18. Несобственные интегралы, зависящие от параметра |
= b −a 1B(a, b − 1) − b −a 1B(a, b).
Таким образом
1 + b −a 1
откуда получаем нужное
ТЕОРЕМА 8.3. Для любых a > 0 и n N справедливо равенство
B(a, n) = |
n − 1 |
|
n − 2 |
|
|
1 1 |
= |
|
(n − 1)! |
. |
|
|
|
|
|
|
· · · a + 1 a |
|
|
a + n − 1 a + n − 2 |
|
k4 |
|
1) |
|
|
|
n (a + k |
− |
|
|
|
|
|
|
|
|
|
|
=1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Доказательство. Применяя последовательно теорему 8.2 получаем:
B(a, n) = |
n − 1 |
|
n − 2 |
|
B(a, n |
− |
2) = ... |
|
a + n |
− |
1 a + n |
− |
2 |
|
|
|
|
|
|
|
|
|
|
a + n − 1 a + n − 2 |
· · · |
|
|
|
|
... = |
n − 1 |
|
n − 2 |
|
|
B(a, 1). |
|
|
|
|
|
Однако |
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
B(a, 1) = 0 |
|
|
|
|
|
|
1 |
|
xa−1(1 − x)0dx = |
|
, |
a |
откуда получаем нужное.
ТЕОРЕМА 8.4. Для любых натуральных m, n справедливо равенство
B(m, n) = (m − 1)!(n − 1)!. (n + m − 1)!
Доказательство.
B(m, n) =
Согласно теореме 8.3 имеем:
(n − 1)! (m − 1)!(n − 1)! 4n = (n + m − 1)! .
(m + k − 1)
§9. Интеграл Эйлера второго рода |
457 |
|
ТЕОРЕМА 8.5. Для любых a > 0, b > 0 справедливо |
|
равенство |
∞ |
|
|
|
|
|
|
|
|
|
ya−1 |
|
|
|
|
|
B(a, b) = 0 |
|
|
|
|
|
|
|
|
|
|
|
dy. |
|
|
|
|
|
(1 + y)a+b |
|
|
|
Доказательство. Мы имеем |
|
|
|
|
|
|
|
|
|
1 |
|
x = |
|
|
|
|
|
dy |
|
|
B(a, b) = |
xa−1(1 − x)b−1dx = |
|
y |
1 |
= |
|
|
|
, dx = |
|
|
1 + y |
(1 + y)2 |
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∞ |
|
y |
a |
1 |
|
y |
dy |
∞ |
ya−1 |
= 0 |
− |
|
|
= 0 |
|
|
|
(1 − |
|
)b−1 |
|
|
dy. |
1 + y |
|
|
1 + y |
(1 + y)2 |
(1 + y)a+b |
УПРАЖНЕНИЕ 2. Из предыдущей теоремы следует, что при 0 < a < 1
∞ |
ya |
B(a, 1 − a) = 0 |
|
dy. |
1 + y |
Попробуйте доказать, что |
|
|
B(a, 1 − a) = |
π |
|
. |
sin aπ |
§9. Интеграл Эйлера второго рода
ОПРЕДЕЛЕНИЕ 9.1. Γ("Гамма")-функцией (интегралом Эйлера второго рода) называется функция
∞
Γ(a) = xa−1e−xdx,
0
где a > 0.
УПРАЖНЕНИЕ 1. Докажите, что данный интеграл сходится тогда и только тогда, когда a > 0.
ТЕОРЕМА 9.1. Для всех a > 0 справедливо равенство
Γ(a + 1) = aΓ(a).
458 |
Глава 18. Несобственные интегралы, зависящие от параметра |
Доказательство.
∞ |
|
|
|
∞ |
|
xa |
|
xa |
∞ xa |
Γ(a) = 0 |
xa−1e−xdx = |
0 |
e−xd( |
|
) = |
|
e−x|0∞ + |
0 |
|
e−xdx = |
a |
a |
a |
|
|
|
∞ |
|
|
|
|
|
|
|
|
|
|
|
1 |
0 |
|
1 |
|
|
|
|
|
|
|
= |
|
xae−xdx = |
|
Γ(a + 1). |
|
|
|
|
a |
a |
|
|
|
ТЕОРЕМА 9.2. Для всех натуральных n справедливо равенство
Γ(n + 1) = n!.
Доказательство.
Γ(n + 1) = nΓ(n) = n(n − 1)Γ(n − 1) = ... = n(n − 1)...1Γ(1).
Учитывая, что
∞
Γ(1) = e−xdx = − e−x ∞0 = 1,
0
Γ(a1 + b)
§9. Интеграл Эйлера второго рода |
|
|
459 |
|
|
|
|
|
ТЕОРЕМА 9.3. Для любых a > 0, |
b > 0 справедливо |
|
|
|
|
|
|
|
|
равенство |
Γ(a)Γ(b) |
|
|
|
|
|
|
|
B(a, b) = |
. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Γ(a + b) |
|
|
|
|
|
|
|
|
|
|
Доказательство. Подстановкой x = ty, t > 0 преобразуем
∞ |
|
|
∞ |
|
∞ |
|
Γ(a) = 0 |
xa−1e−xdx = |
0 |
(ty)a−1e−tytdy = ta |
0 |
ya−1e−tydy, |
или |
|
|
|
∞ |
|
|
|
|
|
|
|
|
|
|
|
|
Γ(a) |
= 0 |
ya−1e−tydy. |
|
|
|
|
|
|
|
|
|
ta |
|
|
Данное равенство верно для всех a > 0 и всех t > 0. Полагая здесь a = a1 + b, t = t1 + 1, получаем:
∞
(t1 + 1)a1+b =
или
Γ(a + b)
(t + 1)a+b =
ya1+b−1e−(t1+1)ydy
0
∞
ya+b−1e−(t+1)ydy. (1)
0
Умножая обе части данного равенства на ta−1 и интегрируя, находим
∞ |
ta−1 |
|
|
∞ |
|
|
∞ |
|
0 |
|
|
|
0 |
|
0 |
Γ(a + b) |
|
|
dt = |
ta−1 |
|
ya+b−1e−(t+1)ydy dt. |
(t + 1)a+b |
Но в силу теоремы 8.5 предыдущего пункта |
|
∞ |
|
|
ta−1 |
|
|
|
|
|
|
0 |
|
|
|
dt = B(a, b) |
|
|
|
(t + 1)a+b |
и, следовательно, |
|
|
∞ ta−1 |
∞ya+b−1e−(t+1)ydy dt. |
Γ(a + b)B(a, b) = |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
0 |
|
|
460 Глава 18. Несобственные интегралы, зависящие от параметра
Переставляя интегралы по y и по t, находим
∞ |
∞ |
|
Γ(a + b)B(a, b) = 0 |
ya+b−1dy 0 |
ta−1e−(t+1)ydt = |
∞ ∞
= ya+b−1e−ydy ta−1e−tydt.
|
|
|
|
0 |
|
|
|
0 |
|
|
|
|
Далее |
|
|
|
|
|
|
|
|
|
|
|
|
∞ |
|
|
|
|
|
|
|
|
|
|
|
|
0 |
ta−1e−tydt = |ty = x, |
y > 0, ydt = dx| = |
|
|
∞ |
|
x |
a−1 |
dx |
1 |
∞ |
|
Γ(a) |
|
= 0 |
|
e−x |
|
= |
|
0 |
xa−1e−xdx = |
|
. |
(2) |
y |
y |
ya |
ya |
Подставляя данное выражение в предыдущий интеграл, приходим к равенству
∞ |
|
|
|
Γ(a + b)B(a, b) = 0 |
ya+b−1e−y |
Γ(a) |
dy = |
|
ya |
∞
= Γ(a) yb−1e−ydy = Γ(a)Γ(b),
0
и окончательно находим
B(a, b) = Γ(a)Γ(b). Γ(a + b)
Для завершения доказательства нам осталось обосновать перемену порядка интегрирования в повторном интеграле
∞ ∞
ta−1dt ya+b−1e−(t+1)ydy.
00
Сделаем это сначала в предположении a > 1, b > 1. Для функции
f(t, y) = ta−1ya+b−1e−(t+1)y
выполнены все условия теоремы 6.1.