Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

матан

.pdf
Скачиваний:
41
Добавлен:
27.05.2015
Размер:
3.81 Mб
Скачать

§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

a

a

452 Глава 18. Несобственные интегралы, зависящие от параметра

Здесь мы воспользовались также непрерывностью левого интеграла как функции y, вытекающей из теоремы 5.1.

УПРАЖНЕНИЕ 1. Соответствующие утверждения о несобственных интегралах с конечными пределами интегрирования (непрерывность по параметру, изменение порядка интегрирования, дифференцирование по параметру) сформулируйте самостоятельно.

ПРИМЕР 1 (интеграл вероятностей). Вычислим несобственный интеграл

e−x2 dx.

−∞

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

∞ ∞

e−x2 dx = 2 e−x2 dx = 2I.

−∞

0

Положим 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) = τ,

= 2u(t2 + 1)du =

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e−τ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= 0

 

 

 

=

1

 

 

 

 

 

0

e−τ =

 

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

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 =

 

 

 

 

 

 

 

 

0

0

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)

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

§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+b1e(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.