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

Ilin_Poznyak_-_Matanaliz

.pdf
Скачиваний:
740
Добавлен:
02.06.2015
Размер:
14.78 Mб
Скачать

§ 5

ИНТЕГРАЛЫ ЭЙЛЕРА

301

1. Вычислим интеграл

00

1 = Jx 1/ 4 (1 + х)-2 dx.

О

Обращаясь к формулам (9.44) и (9.48), очевидно, получим

1 = В (~

~) = Г(5j4)Г(Зj4) =

~г (~) (~).

 

4' 4

Г(2)

4 4 4

2. Вычислим интеграл

 

 

 

 

7Г/2

 

 

1 =

J sinp- 1 epcosq- 1 epdep.

 

 

 

О

 

Полагая х =

sin2 ер,

получим

r(P)r(q)

1

х,-l(1

 

 

I ~ %[

-

х)!-l dx ~ %в(И) ~ %г(~) .

3. Обратимся к интегралу

7Г/2

I p - 1 = J sinp- 1 ер dep.

о

Используя результат, полученный в примере 2 (надо положить q = 1), найдем

7Г/2

 

 

 

 

2 2

(Р)

 

J

 

d

 

г 2" .

(9.49)

о

sinp-l

ер

ер

= ~г(~)

 

 

 

г(Р:1)

 

Имеем далее

г(~) = J00 е-Х~ = 2 J00 е-(у'х)2 dvГx·

о

о

 

 

 

 

00

2

100

2

полу-

Полагая vгx = t и замечая, что Je- t

dt равен -

J e- t

dt,

о

 

2 _ 00

 

 

чим, согласно примеру, рассмотренному в п. 2 § 4 гл. 3 (интеграл Пуассона)

г(~) = VК.

Поэтому формула (9.49)

примет вид

 

 

 

 

7Г/2

 

 

 

 

(Р)

 

1

-

J.p-l

ер

d

- V7Г г

2"

(9.50)

p-l-

О

 

ep-Т г(Р:1)·

 

302

ИНТЕГРАЛЫ, ЗАВИСЯЩИЕ ОТ ПАРАМЕТРОВ

ГЛ.9

§ 5. Формула Стирлинга

Формулой Стирлинга называют следующую асимптотическую

формулу:

(9.51)

где аn ---+ О при n ---+ 00.

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

00

Г(л + 1) = Jt>'e- t dt.

(9.52)

О

 

Для этого мы воспользуемся так называемым

м е т о Д о м

Л а п л а с а, опирающимся на следующее утверждение.

Ле,м,,м,а 1. Пусть фу'Н/х;'Ци-я j(t), интегрируема-я при нек;о­

тором а > О на сегменте [-а, а], представима в виде

2n-l

 

f(t) = L Cktk + O(t2n ).

(9.53)

k=O

Тогда имеет место следующа-я асимnтоти'Ч,еск;а-я формула:

a e->..t2f(t) dt =

~ С2тг(m + ~) + 0(1)

(9.54)

J

~

..\.т+~

..\.n+~·

 

т=О

д о к а з а т е л ь с т во. Подставим соотношение (9.53) в ин­ теграл, стоящий в левой части формулы (9.54), и учтем, что

интегралы, отвечающие нечетным степеням t, обратятся в нуль. Для оценки оставшихся интегралов достаточно убедиться в спра­

ведливости при m ? О следующего равенства:

a

2

!

г(т+ !)

+ О(е->..а

2

(9.55)

t 2m e->..t

dt =

2

).

J

 

2

..\.т+~

 

 

 

О

Представим интеграл в левой части (9.55) в следующем виде:

а

2

00

2

00

2

 

Jt 2m e->..t

dt =

Jt 2m e->..t

dt -

Jt 2m e->..t

dt.

(9.56)

О

 

О

 

а

 

 

впервом интеграле в правой части (9.56) произведем замену

х= лt2 и получим

(9.57)

§ 5

ФОРМУЛА СТИРЛИНГА

303

Далее заметим, что при л > 1 и t ? а справедливо следующее

неравенство:

_>..t2 ./ -(>"-1)а2 _t2

е"'" е е.

Применяя это неравенство, оценим второй интеграл в правой

части (9.56)

(9.58)

аа

Из равенств (9.56), (9.57) и оценки (9.58) вытекает требуемая формула (9.55). Лемма доказана.

Для того чтобы применить эту лемму, произведем в интег­

рале (9.52) замену t =

л(l + х).

В результате этот

интеграл

примет вид

 

 

 

 

 

Г(л + 1) =

 

(х)

 

 

 

л>"+lе->" J e->..[x-ln(l+x)] dx.

(9.59)

 

 

-1

 

 

 

Обозначим через g(x) следующую функцию, определенную

на полупрямой х > -1:

 

 

 

 

 

g(x) =

sgnx· Jx -ln (1 + х).

(9.60)

Тогда равенство (9.59)

можно переписать в виде

 

Г(л + 1) =

 

(х)

 

л>"+lе->"

J e->..g2(x) dx.

(9.61)

 

 

 

-1

 

Нашей целью является изучение асимптотического поведе­

ния при л ----7 +00 следующего интеграла:

 

 

 

(х)

2

 

 

I(л)

= J e->..g

(х) dx.

(9.62)

 

 

-1

 

 

 

Для этого рассмотрим подробнее функцию g(x), определен­

ную равенством (9.60). Поскольку

 

 

 

.!!:...g2(x) =

.!!:...(x-ln(l+x)) = _х_,

(9.63)

dx

dx

 

 

1 + х

 

то функция g2 (х) строго убывает при -1 < х < О и строго воз­

растает при х > О. Отсюда вытекает, что функция g(x) строго

возрастает на полупрямой х > -1, причем областью ее значений

является вся числовая прямая. Далее, так как функция g2(x)

имеет в окрестности точки х = О разложение

304 ИНТЕГРАЛЫ, ЗАВИСЯЩИЕ ОТ ПАРАМЕТРОВ ГЛ.9

то существует строго положительная при х > -1 функция h )

такая, что справедливо равенство

g2(x) = x2h(x).

Функция h(x) является бесконечно дифференцируемой при

х > -1, следовательно, бесконечно дифференцируемой будет и

функция g(x) = xJh(x).

Учитывая вышенаписанное, можно утверждать, что для функ­

ции у = g(x), определенной равенством (9.60), существует обрат-

ная функция х = g-l(y), строго возрастающая и бесконечно

дифференцируемая на всей числовой прямой и удовлетворяю­

щая условию g-l (О) = о.

Обозначим эту обратную функцию символом х = с.р(у). Ис­

пользуя указанные выше ее свойства, найдем асимптотику ин­

теграла (9.62). Справедлива следующая теорема.

Теоре,м,а 9.13. Пустъ фу'Н/х;'Цuя х = с.р(у) является обрат­ н.ой. 1>: фун.1>:'ЦUU У = g(x), оnределен.н.оiJ. равен.ством (9.60). Тогда для uн.теграла (9.62) nри любом фU1>:сuрован.н.ом н.омере n сnра­

ведлuва следующая асuмnтотu'Ч,еС1>:ая формула:

1 л =

~ cp(2m+l) (О) г(m +~) + 0(1)

 

()

~

(2т)!

лm+~

лn+~ .

(9.64)

 

т=О

 

 

 

 

 

Д О К а з а т е л ь с т в о.

Фиксируем произвольное положитель­

ное число а и положим Ь =

с.р(-а),

с = с.р(а).

Это означает, что

а = g(c) = -g(b), и,

следовательно,

-1 < Ь < О и с > о.

 

Оценим следующие два интеграла:

 

 

Ь

 

 

 

00

 

 

11(л) = J е-Лg2 (х) dx,

12(л)

= J е-Лg2 (х) dx.

(9.65)

-1

 

 

 

с

 

 

Для оценки первого интеграла заметим, что при -1 < х < Ь

выполняется неравенство g(x) < -а, т. е. g2(x) > а2 и, следова­

тельно,

е-Лg2 ( х) < е-Ла2 .

в таком случае

11 (л) ~ е-Ла

2

Ь

IbI)e-Ла2 .

(9.66)

 

J dx = (1 -

 

 

-1

 

 

Аналогично оценивается интеграл 12(Л). При х > с выпол­

няется неравенство g(x) > а, т. е. g2(x) > а2 . Следовательно,

при л > 1 их> с имеет место оценка

е-Лg2(х) = e-(Л-1)g2(Х)е-g2 (х) < е-(Л-1)а2 e- g2 (x).

§ 5

ФОРМУЛА СТИРЛИНГА

 

305

Отсюда мы получаем

 

 

 

00

 

 

12(>")

~ е-(л-l)а2 J e-g2 (x) dx =

Сlе-Ла2 .

(9.67)

с

Из оценок (9.66) и (9.67), которым удовлетворяют интегралы (9.65), мы получаем для интеграла (9.62) следующее соотношение:

 

с

+ О(е-Ла2 ).

 

1(>") =

Jе-лg2 (х) dx

(9.68)

 

ь

 

 

Произведем в интеграле (9.68) замену переменной t

= g(x),

т. е. х = <p(t). в результате получим

 

 

 

а

 

 

1(>") =

J e-Лt2<р'(t)dt+О(е-Ла2).

(9.69)

Поскольку функция <p'(t) является бесконечно дифференци­

руемой, воспользовавшись формулой Маклорена, представим ее

в виде

2n-l

<р'(t) = ""'rp (k+l) (О) t k + О(t 2n ).

~ k! k=O

Для получения формулы (9.64) остается применить лемму к

функции J(t) = <p'(t). Теорема 9.13 доказана.

В заключение этого параграфа укажем следующий простой

способ вычисления производных <p(k)(o). Из равенства (9.63) мы

получаем

2

.

' -

_Х_ _

rp(t)

g

g

-

х + 1 -

rp(t) + 1

Из этого равенства вытекает соотношение

<p'(t) = ~ = 2g 1 + rp(t) = 2t 1 + rp(t).

g'

rp(t)

rp(t)

Таким образом, мы получаем следующее равенство:

<p(t)<p'(t) =

2t + 2t<p(t).

(9.70)

Последовательно дифференцируя это равенство и полагая t = о,

мы определим все производные <p(k) (о). Найдем для примера

значения первых трех производных функций <p(t) В нуле. Продифференцировав (9.70), получим

[<p'(t)]2 + <p(t)<p"(t) = 2 + 2(t<p'(t) + <p(t)).

(9.71)

Положим t = О и учтем, что <р(0) = о. Тогда <р,2(0) =

2, т. е.

<р'(о) = V2.

 

306

ИНТЕГРАЛЫ, ЗАВИСЯЩИЕ ОТ ПАРАМЕТРОВ

ГЛ.9

 

После дифференцирования равенства (9.71) получим

 

 

3<р' . <р" + . <р'"

=

2(t<p" + 2<р').

 

Приравняв t нулю, получим 3V2<p"(O) = 4V2, т. е. <р"(О)

= 4/3.

Аналогично из равенства

 

 

 

 

 

 

 

 

з<р,,2 + 4<р' . <р'" + <р. <pIV = 2(t<p'" + 3<р").

 

получим <р"'(О) = V2/3.

 

 

 

 

 

 

 

 

Следовательно формулу (9.64) можно записать в виде

 

I(л) = J2 Г(1/2) +

v'2 Г(3/2) +

0(1) .

(9.72)

 

 

 

 

6

ЛVЛ

Л2

 

 

Подставим равенство (9.72)

 

в (9.61) и учтем, что Г(I/2) =

= yJГ, Г(3/2) = (1/2)Г(I/2) = YJГ/2. В результате получим

 

Г(л+1)=V27Глл

Л

е-

Л

1

 

(9.73)

 

 

 

(I+- +0(1)).

 

 

 

 

 

 

12Л

Л2

 

Выпишем первые пять членов асимптотического разложения

гамма-функции Эйлера:

Г(Л+1)= V27ГЛ лле-л (1+_1_+_1__

139

571

+0(1))

12Л 288Л2

51840Л3

2488320Л4

Л5

Отметим без доказательства, что остаток асимптотического

ряда не превосходит последнего удерживаемого слагаемого.

§6. Кратные интегралы, зависящие от параметров

1.Собственные кратные интегралы, зависящие от па­

раметров. Пустьх = (хl, Х2, ... , хт) - произвольная точка об­

ласти D т-мерного евклидова пространства Ет , а у = (Уl, У2, ...

. .. , Yl) - точка области s1 пространства E l . Обозначим симво­

лом D х s1 подмножество (l + т)-мерного евклидова простран­ ства, состоящее из всех точек Z = (Zl, Z2, ... , Zm+l) таких, что

точка (Zl, Z2, ... , Zm) принадлежит D, а точка (Zm+l, Zm+2, ...

... , Zm+l) принадлежит S1. При этом мы будем часто исполь­ зовать обозначение Z = (х, у) Е D х S1. Замыкание области D

мы будем обозначать символом D. Легко видеть, что замыкание

D х s1 совпадает с D х S1.

Пусть f(x, у) -функция, определенная в D х S1, причем для любого уа Е s1 функция f (х, У) интегрируема по х в области D.

Тогда функцию

 

I(y) = JЛХ, у) dx,

(9.74)

D

определенную в области S1, будем называть интегралом, зави­ сящим от параметра у. Заметим, что параметр у является

§ 6

КРАТНЫЕ ИНТЕГРАЛЫ, ЗАВИСЯЩИЕ ОТ ПАРАМЕТРОВ

307

[-мерным вектором и, следовательно, интеграл (9.74) зависит от

[ числовых параметров Уl, У2, ... , Yl·

В полной аналогии с теоремами 9.9-9.12 доказываются сле­

дующие теоремы.

Теорема 9.14 (о непрерывности интегра.ла (9.74) по параметру). Если фун'Х:'Ция Лх, у) непрерывна по сово'Х:уnности

аргументов в зам'Х:нутоi1 области D х Q, то интеграл (9.74)

представляет собоi1неnрерывную фун'Х:'Цию параметра у в обла­

сти Q.

Теорема 9.15 (об интегрировании интегра.ла (9.74)

по параметру). Если фун'Х:'Ция Лх, у) непрерывна по сово'Х:уnно­

сти аргументов в зам'Х:нутоi1 области D х Q, то фун'Х:'Цию

(9.74) MO;JfCHO интегрировать по параметру под зна'Х:ом интег­ рала, т. е. справедливо равенство

JI(y) dy = J dx JЛх, у) dy.

!1 D!1

Теорема 9.16 (о дифференцируемости интегра.ла (9.74)

по параметру). Если фун'Х:'Ция f(x, у) И ее 'Частная nроизвод-

ная.!!.L непрерывны в D х Q, то интеграл (9.74) имеет в облас­

aYk

дI

ти Q непрерывную 'Частную nроизводную - , nри'Чем сnравед­ aYk

ливо равенство J

~ = af(x, У) dx. aYk aYk

D

2. Несобственные кратные интегралы, зависящие от параметров. Понятие несобственного кратного интеграла, за­ висящего от параметров, можно было бы ввести так же, как

ив предыдущем пункте, для случая, когда подынтегральная

функция f(x, у) определена в D х Q , где D с Ет и Q с E l .

Однако наибольший интерес представляет случай D = Q, кото­ рый и будет нами изучен. Кроме того, мы будем предполагать,

что f(x, у) = F(x, у) g(x), где F(x, у) непрерывна при х # у в D х D, а функция g(x) ограничена в D. Таким образом, мы

рассматриваем интегралы вида

V(y) = JF(x, y)g(x) dx,

(9.75)

D

где подынтегральная функция может иметь особенности лишь при х = у. Нас будет интересовать вопрос о непрерывности ин­

тегралов вида (9.75) по параметру у. В связи с этим введем сле­ дующее определение равномернои сходимости интеграла (9.75)

вто'Ч'Х:е. Символом К(Уо, 6) мы будем обозначать шар радиуса 6

сцентром в точке уо.

З08 ИНТЕГРАЛЫ, ЗАВИСЯЩИЕ ОТ ПАРАМЕТРОВ ГЛ.9

Оnреде.ле1-tuе. Интеграл (9.75) назовем с х о д я Щ и м с я

равномерно по параметру У в то'Чnе Уо Е D, если

для любого Е> О можно уnазатъ 6 > О таnое, 'Что К(Уо, 6) с D

и для любоil. nубируемоil. области w с К(Уо, 6) и всех то'Чеn У Е Е К(уо, 6) выполняется неравенство

I!F(x, y)g(x) dxl < Е.

Теоре,м,а 9.17. Если интеграл (9.75) сходится равномерно

по У в то'Чnе уо Е D, то он непрерывен в то'Чnе уо.

Д о к а з а т е л ь с т в о. Требуется доказать, что для любого

Е > О найдется 6 > О такое, что при Iy - Yol < 6 выполняется

неравенство IV (у) - V (Уо) I < Е. ИЗ определения равномерной

сходимости в точке следует существование такого 61 > О, что

К(Уо, 61) с D и при У Е К(Уо, 61)

11 F(x, y)g(x) dxl < Е/3.

(9.76)

Положим

JF(x, y)g(x) dx,

 

V1 (У) =

(9.77)

 

к

V2 (y) =

J F(x, y)g(x) dx.

 

 

D\K

 

ИЗ неравенства (9.76) следует, что при Iy - Yol < 61

 

IV1 (y)1 < Е/3.

(9.78)

Далее заметим, что при х Е D \ К(Уо, 61) и У Е К(Уо, 61/2)

функция F(x, У) будет равномерно непрерывной по совокупно­

сти аргументов. Следовательно, найдется положительное чис­

ло 6 < 61/2 такое, что при Iy - Yol < 6 будет выполняться нера­

венство

IF(x, Уо) - F(x, y)1 < E/(3MIDI),

где М - константа, ограничивающая функцию g, и IDI- объем области D. В таком случае, при Iy - Yol < 6

IV2 (y) -

V2 (yo)1

~ м

J

IF(x, Уо) - F(x, y)1 dx ~ Е/3.

 

 

 

D\K(yo,Jl)

 

 

 

 

 

(9.79)

Из соотношений (9.77)-(9.79) следует, что при Iy - Yol < 6

IV(y) -

V(yo)1

~ IV1 (y)1 + IV1 (yo)1 + IV2 (y) - V2 (yo)1 < Е.

Теорема доказана.

§ 6

КРАТНЫЕ ИНТЕГРАЛЫ, ЗАВИСЯЩИЕ ОТ ПАРАМЕТРОВ

309

 

Укажем одно достаточное условие равномерной сходимости

интеграла в точке, наиболее часто встречающееся в приложениях.

Теорема 9.18. Пусть фу'Н'Кция F(x, у) 'Неnрерыв'На в D х D при х i- у, а фу'Н'Кция g(x) рав'Номер'Но огра'Ни'Ч,е'На в D. Предпо­

ложим, 'Ч,то существуют nостоя'Н'Ные Л, О < л < т, и с > О та'Кие, 'Ч,то для всех х Е D, у Е D выnол'Няется 'Нераве'Нство

 

(9.80)

Тогда и'Нтеграл (9.75)

сходится рав'Номер'Но по у в 'Каждоi1

то'Ч,'Ке уо Е D.

 

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

Пусть уо - произвольная точка обла­

сти D. Требуется доказать, что для любого Е > О найдется <5 > О

такое, что для любой кубируемой области w с К(Уо, <5) и всех у Е К(уо, <5) выполняется неравенство

I!F(x, y)g(x) dxl < Е.

 

(9.81)

Применяя (9.80) и используя условие ограниченности g(x),

получим

 

 

 

I!F(x, y)g(x) dxl

:::;; cl!lx - уl-Л dx.

 

Фиксируем точку у Е К(Уо, <5) и заметим,

что из условия

w С К(Уо, <5) вытекает включение w С К(у, 2<5).

Следовательно,

I JF(x, y)g(x) dxl :::;; Cl

J Ix - уl-Л dx.

(9.82)

w

J((y,2J)

 

 

Переходя в интеграле в правой части (9.82) координатам с центром в точке у (см. гл. 2, § 5,

к сферическим п. ЗО ), получим

 

 

J

F(X, y)g(x)dx

:::;; С2Jrт-1-Л dr = С22m-л <5т-Л = сз<5т-Л.

 

т-Л

 

 

w

 

о

Отсюда следует, что, выбрав <5 достаточно малым, мы получим

неравенство (9.81). Теорема доказана.

3.Приложение к теории ньютонова потенциала. Пусть

внекоторой точке Ро(х, у, z) помещена масса то. По закону все­

мирного тяготения на массу т, помещенную в точку M(~, ТJ, (),

действует сила

F

 

тто_

- ry

r

 

--/ю'

где R = р(Ро, М) = J(x -

~)2 + - ТJ)2 + (z - ()2, 1'- гравита­

ционная постоянная, r = R/R -

единичный вектор, направление

310

ИНТЕГРАЛЫ, ЗАВИСЯЩИЕ ОТ ПАРАМЕТРОВ

ГЛ.9

которого совпадает с направлением вектора РоМ. Считая I = 1

и массу m = 1, получим силу тяготения

F = _то r.

R2

Отметим, что компоненты этой силы имеют вид

Х =

-

~~ (~ -

х),

у =

-

~~ (7] -

у),

Z =

-

то (( - z).

 

 

R3

 

Очевидно, что потенциал силы тяготения, определяемый как скалярная функция U такая, что F = grad и, равен

то И=-.

R

Если масса сосредоточена не в точке Ро(х, у, z), а распределена

по области D с плотностью р(х, у, z), то для потенциала силы

и для компонент силы мы получим следующие выражения:

и(~, 7], () =

111р(х,:, z) dx dy dz,

(9.83)

Х=- 111

D

 

р(Х~~'Z)(~-х)dхdуdz,

 

D

 

 

у = - 111p(x~~, z) (7] - у)dx dy dz,

(9.84)

D

Z=- 111р(Х~~'Z)((-z)dХdУdz.

D

Нетрудно показать, что интегралы (9.84) представляют со­ бой частные производные потенциала (9.83). Поскольку подын­ тегральные функции в интегралах (9.83) и (9.84) мажорируются

функцией ~, где л = 1 для интеграла (9.83) и л = 2 для интег-

Ю,

ралов (9.84), то в силу теоремы 9.18 указанные интегралы схо­ дятся равномерно в каждой точке M(~, 7], (). Следовательно, по

теореме 9.17 эти интегралы представляют собой непрерывные

функции точки M(~, 7], ().

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]