Ilin_Poznyak_-_Matanaliz
.pdf§ 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- |
О |
SШ |
|
ep-Т г(Р:1)· |
|
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) |
||||
|
vл |
|
|
|
6 |
ЛVЛ |
Л2 VЛ |
|
|
Подставим равенство (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 < Е. |
Теорема доказана.
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], ().