
книги из ГПНТБ / Дьяченко, А. Н. Интегральное исчисление функций нескольких переменных. Теория вероятностей и элементы математической статистики учебное пособие
.pdfсоответствующее этому значению, является криволинейной трапе цией АВСЕ (см. рис. 13), образованной функцией г = f [х, у). При фиксированном у эта функция является функцией только одного аргумента х, который меняется от х х (у) до х 2 (у). Площадь кри волинейной трапеции будет зависеть от фиксированного у, т. е. будет являться функцией от у. Обозначим эту функцию 1Р (у), тогда
*2 (У) |
|
xP(y)= I fix, у) dx. |
(1.26) |
*1 (У) |
|
Формула (1.26) верна для любого у из интервала [с, d]. Воспользуемся, как и в предыдущем случае, формулой для вы
числения объема тела, заданного площадями сечений:
v —I ^ (У) dy.
С
Подставив в это равенство формулу (1.26), получим, учитывая,
что J — V. |
*2J(У) /(х, |
|
|
|
|
■М |
у) dx |
dy. |
|
||
|
_ *1 (У) |
|
|
I |
|
Эту формулу обычно записывают в виде |
|||||
|
|||||
|
d |
х, (у) |
|
|
|
Ш (* » y)dxdy = \dy J |
f(x, у) dx. |
(1.27) |
|||
D |
с |
х, (у) |
|
|
Если область D является простой (как в направлении оси х, так и в направлении оси у), то справедливы оба равенства: (1.25) и (1.27). Каждое из выражений, стоящих в правых частях этих ра венств, называется двукратным (или повторным) интегралом от функции / (х, у) по области D. Оба двукратных интеграла полу
чаются в |
результате |
последовательного интегрирования сначала |
по одной |
переменной, |
затем — по другой, и отличаются Друг от |
друга порядком интегрирования. Представление двойного интег рала двукратным иногда называют р а с с т а н о в к о й п р е д е л о в в двойном интеграле.
Теорема. Если функция f (х, у) интегрируема в простой области.
D, то двойной интеграл от этой функции по области D равен каж дому из двукратных интегралов от функции f (х, у) по области D,
т. е.
|
Ь |
у, (х) |
d |
х , (у) |
Ш (* . |
y)dxdy = \dx |
| |
f(x, у) dy = l dy |
J f{x, y)dx. |
D |
“ |
!/i (x) |
c |
x , (y) |
Приведенные выше рассуждения доказывают эту теорему для частного случая, когда функция f (х, у) в области D непрерывна
иположительна. Доказательство теоремы для более общего случая, т. е. для любой функции, интегрируемой в D, выходит за рамки на шего курса.
Таким образом, формулы (1.25) и (1.27) могут служить для вы числения двойного интеграла.
20
Замечание. Если область D не является простой, то формулы (1.25) и (1.27) применимы в том случае, когда область D может быть разбита на конечное число простых областей.
Пример 3. Пусть дан двойной интеграл
f | / (*> У) dxdy,
'd
где D — область D, ограниченная эллипсом
Требуется расставить пределы интегрирования двумя способами т. е. двумя способами представить двойной интеграл в виде двукрат ного.
Р е ш е н и е . Для того, чтобы расставить пределы интегрирова ния, необходимо задать границу области D в одной из двух возможных канонических форм. Это было сделано в примере 1.
1. Пусть граница области задана в следующей форме:
х = — а, х — а\
у = -----— Y a2>— х2, у = - — а2— х2.
Тогда, согласно формуле (1.25), имеем:
•Уа1-
у) dxdy = j' dx |
| |
f(x, y)dy. |
|
------Va-— x- |
|
Это есть один из двух возможных способов расстановки пределов в двойном интеграле.
2. Пусть граница области задана в следующей форме:
У = |
— Ь, |
у = |
Ь; |
|
|
|
|
|
- — |
V b 2 — y2\ x = — Y b 2 — y2. |
|||||
|
b |
|
|
|
|
b |
|
Тогда, согласно формуле (1.27), имеем: |
|||||||
f(x, |
у) dxdy = |
Jьdy |
|
|
j" |
f(x, y)dx. |
|
D |
|
|
-Ь |
- |
а |
/-ТЗ---- |
|
|
|
|
|
T |
V b-y- |
||
Это равенство представляет собой второй способ расстановки пре |
|||||||
делов в двойном интеграле. |
|
|
|
|
|
||
Пример 4. Пусть дан двойной |
интеграл |
j*j"/ |
(дг, у) dxdy, |
D
где D — область, ограниченная линиями
У= х2+ 1,
У— х -)- 3.
Требуется расставить пределы интегрирования двумя способами.
Р е ш е н и е . Воспользуемся результатом примера 2, в котором граница рассматриваемой области была представлена в канонических формах.
1. Пусть граница области представлена в первой канонической форме:
С х = — 1, |
х = 2; |
\ </ = х2+ 1 , |
у = х + 3. |
Тогда, согласно формуле (1.25), можно представить двойной ин теграл двукратным:
2 |
я-ЬЗ |
|
JJ f(x, у) dxdy = J dx |
J f{x, y)dy. |
(1.28) |
D—1 x2+l
2.Пусть граница области представлена во второй канонической
форме:
у= 1 - , у — 5
\ —VУ— 1 при 1 < ( / < 2 , |
------ - |
(1.29) |
X — < |
X = у у — 1. |
|
{у — 3 при 2 < у < 5;
Согласно формуле (1.27), имеем:
j*J / (х, y)dxdy = |
5 |
У у- |
1 |
|
§dy |
j |
f {x, y) d x, |
(1.30) |
|
D |
1 |
x, (у) |
|
|
где x, (у) — линия АВС на рис. |
11. |
Учитывая, что линиях, (у) состоит |
из частей различных линий, имеющих разное аналитическое задание, двукратный интеграл, стоящий в правой части формулы (1.30), удоб нее представить в виде суммы двух интегралов:
5 |
|
2 |
У у ^ л |
5 |
тТ рГ |
\dy |
J |
/ (х, у) dx — ^dy |
J |
f ( x , y ) d x + § d y |
j' f (x, y) dx . |
1 |
X, (у) |
1 |
_ |
2 |
if—3 |
(1.30а)
Очевидно, что для вычислений формула (1.28) удобнее, чем (1.30а).
Пример 5. Вычислить двойной интеграл JJx2ydxdy, где D — область, ог-
раниченная линиями у = |
х, |
у = |
D |
|
|
2х, х = 3. |
|
||||
Р е ш е н и е . |
Представим |
границу области D в канонической |
|||
форме. Это можно сделать двумя способами (рис. |
14): |
||||
|
Г |
х = |
0; |
х = 3, |
|
и |
I |
У = |
х; |
у — 2х |
|
|
|
|
|
|
|
( у = |
0; |
у = 6 , |
|
|
|
х __ _У_ . |
|
Г У п р и 0 < у < 3 , |
|||
|
2 |
|
| |
3 при 3 < г/ < |
6. |
Соответственно, исходный двойной интеграл можно двумя спо собами представить в виде двукратного:
з 2х |
|
J | x2ydxdy — Чdx j x2ydy |
(1.31) |
D |
0 x |
и
3 |
у |
6 |
3 |
(1.31a) |
С( x2ydxdy = I dy f x2ydx + |
[ dy |
Г x2ydx. |
||
D |
J L |
3 |
y_ |
|
|
2 |
|
2 |
|
Легко видеть, что для |
вычислений более удобна формула |
(1.31), |
с помощью которой мы и вычислим исходный двойной интеграл. Вы числим «внутренний» интеграл, входящий в двукратный интеграл формулы (1.31):
2х |
2х |
„ |
2х |
i l |
3^ •4 |
j" |
x2ydx = х2 f ydy = |
x2 У— |
|
||
X |
X |
2 |
x |
2 |
2 |
Подставив этот результат в формулу (1.31), получим:
' |
3 |
* |
1 |
3 |
3^ |
3= |
729 |
|
|
|
£5 |
72,9. |
|||||||||
J | x2ydxdy = — |
j x4dx |
|
2 |
5 |
10 |
|||||
d |
2 |
о |
2 |
5~ 0 |
|
|
||||
Советуем читателю самостоятельно убедиться в том, что вычисле |
||||||||||
ния по формуле (1.31а) приводят к тому же результату. |
|
|
||||||||
Пример 6. |
Вычислить двойной интеграл И(*2 + |
У2) dxdy, |
где |
D — об- |
||||||
|
|
|
|
|
Ь |
|
|
|
|
|
ласть, ограниченная |
линиями |
у2 = |
2рх\ |
у2 = |
— 2рх; |
у = |
Л (/г)> 0, |
р> 0 ) .
Ре ш е н и е . Представим границу области D во второй кано
нической форме (рис. 15):
' У = |
0; |
у = |
h, |
х = |
у2 |
х = |
у2 |
------— ; |
— . |
||
|
2р |
|
2р |
В этом случае расстановку в двойном интеграле следует выпол нить по формуле (1.27), получаем следующую формулу:
|
|
У2 |
|
ft |
2р |
[ j (х2 + |
у2) dxdy = \ dy |
j (х2 + у2) dx. |
D |
о |
_ j/y_ |
|
|
2Р |
23
Вычислим «внутренний» интеграл:
у- |
|
JL |
2р |
_ |
2р |
|
(х2+ у2) dx = (— + У2х |
2у6 |
|
24р3 |
|
jr_ |
|
|
|
|
|
2Р |
|
2р |
Подставляя этот результат в предыдущую
2yi
2р 12р3
формулу, получаем:
(х2+ у2) dxdy ■ |
( Ув |
dy = |
5р |
|
\ 12р3 |
||||
|
84р3 |
|||
|
И? |
А5 |
|
|
|
84р3 |
5р |
|
Заметим, что если расставить пределы интегрирования другим способом (изменить порядок интегрирования), то формула для вычис ления двойного интеграла окажется более сложной. Действительно, границу области в этом случае необходимо будет представить в первой канонической форме:
х — |
К2 |
h2 |
|
|
------; |
х = — |
|
|
|
|
2р |
2р |
|
|
|
Y — 2рх при |
----- < х < О |
|
|
У = |
|
|
2р |
У = h, |
|
при О < |
к2 |
||
|
У 2рх) |
|
||
|
х < ------ |
|
||
|
|
|
2Р |
|
и формула для вычисления двойного интеграла примет вид: |
||||
|
|
|
hr |
|
|
О |
Ь |
2р |
h |
j J (х2+у2) dxdy = j dx j (x2 + у2) dy + j dx |
| (x2+y2) dy. |
|||
D |
ъ? |
Y — 2px |
0 V 2px |
|
|
2P |
|
|
|
1.4.Замена переменных в двойном интеграле
Постановка задачи. Пусть требуется вычислить двойной интег рал J от непрерывной функции / (х, у) по области D:
J = JI/ (х, |
П |
|
у) ds= lim 2 / (xk, yk) &sk. |
(1.32) |
|
D |
%n^0 k=l |
|
Рассмотрим кроме системы координат Оху, еще систему коорди нат Огт (рис. 16). Пусть заданы функции, связывающие перемен ные х и у с новыми переменными и и о:
х = х(и, и); у = у(и, v). |
(1.33) |
Будем считать, что эти функции удовлетворяют следующим ус ловиям:
24
1)в плоскостй (uv) существует такая область D ', что между точками D' и D функции (1.33) устанавливают взаимно-однознач ное соответствие;
2)в области D' функции х (и, v) и у (и, v) имеют непрерывные частные производные первого порядка;
3)каждой линии области D' соответствует одна и только одна линия области D\
4)определитель
дх ду
ди |
ди |
/ = |
(1.34) |
дх |
ду |
dv |
dv |
называемый я к о б и а н о м , отличен от нуля для всех точек об ласти D '. При этих условиях требуется преобразовать интеграл J
As'
Ри с. 16
кдвойному интегралу по новым переменным и, следовательно, по области D'.
Геометрический смысл якобиана. Рассмотрим произвольное дробление области D и соответствующее ему дробление области D'.
Пусть криволинейный четырехугольник А 1А 2А3А4 (рис. 16) с пло щадью As есть часть области D, а соответствующий ему криволи
нейный четырехугольник А\А2АЪ4 с площадью As' — часть об
ласти D'. Пусть координаты точек А{ будут (uh vt), а координаты точек Аг— (хг- г/г), где i = 1, 2, 3, 4. Согласно ранее принятым допущениям, х, = х (щ, ог); yt = у (Щ, vt).
Рассмотрим матрицы:
щ — щ i>2—v± \ _
и3---U% Vs--- /
25
Вычислим произведение этих матриц:
АВ = |
U2 |
Ui |
щ — |
|
и 3 |
и 2 |
щ — |
||
|
||||
— Uj) |
дх |
, |
|
|
dv |
(U2 — Ul); |
|||
|
|
|
||
~ ( u 3 — u2) + |
~ |
(v3— vt) ; |
||
du |
dv |
|
|
( |
дх |
ду |
|
du |
|||
|
ди |
||
|
дх |
ду |
|
\ |
dv |
dv |
|
ду |
(и 2— |
“ i H |
|
ди |
|
|
dy_ (u3— tl2)
du
ду_ (Vz—Vi)
dv
(V3— V2)
Заметим, что все элементы полученной матрицы — полные диф ференциалы функций х (и, v) я у (и, и).* Считая, что приращения аргументов и и v малы, заменим дифференциалы функций их при ращениями ** и получим следующее приближенное равенство:
|
|
/ |
дх |
ду |
\ |
|
|
/ |
Ujj—- Ux |
V2 — V! \ |
ди |
ди |
| ^ |
Хх |
у2— |
|
|
V3— V2 ) |
|
|
|||
[ |
U3 U2 |
дх |
ду |
I ~ |
V *з — *2 |
Уз — У2 ! ‘ |
|
|
|
\ |
dv |
dv |
У1 |
|
|
По теореме об определителе произведения матриц *** имеем
-и г |
v2—v1 |
I ж |
х2—хг у2—ух |
и$ ——11% |
v3—v2 |
Х3---Х 2 Уз--- Уз |
Будем считать, что рассматриваемые криволинейные четырех угольники — параллелограммы. (При малом ранге дробления об ласти D погрешность этого допущения мала.) Площадь паралле лограмма может быть вычислена как удвоенная площадь треуголь ника, а площадь последнего — по формуле (2.77) [3], т. е.
«2 — «1 Ц2 — щ
;As'
х3 —х2 Уз Уз |
|
|
и3— «2 |
|
Подставляя эти формулы в равенство (1.35), получим: |
|
|||
As : |
Л As'. |
(1.36) |
||
Равенство становится точным при стремлении к нулю ранга |
||||
дробления области D, т. е. |
|
|
|
|
|/ |= |
П т |
. |
(1.37) |
|
1 |
jl -о |
As' |
v |
' |
Таким образом, абсолютная величина якобиана может быть рассмотрена как коэффициент деформирования площадей частич ных областей при преобразовании по формулам (1.33).
*См. [3], формула (4.12).
**См. [3], формула (4.15).
***См. [1], формула (1.82).
26
Формула замены переменных в двойном интеграле. Для вычис ления двойного интеграла J составим интегральную сумму:
Vn = 2 f(xk, yk) Ask. |
(1.38) |
ft=i |
|
Применяя формулу (1.36) к каждой частичной области, получим:
Ask ж 11 1Ask (k= 1, 2, . . . , п).
Подставив это соотношение в (1.38), получим приближенное равенство
|
П |
|
|
vn tt 2 / ( * ( “ *> Vk), y(uk, vk))\I\Ask. |
(1.39) |
v |
fe=i |
|
|
|
Перейдем в этом равенстве к пределу при стремлении к нулю ранга дробления области D, заметив, что при этом стремится к нулю и ранг дробления области D ' . В левой части равенства (1.39) содержится интегральная сумма по области D от функции f (х, у), а в правой части равенства — интегральная сумма по области D' от функции
f (х(и, v), у (и, v) ) 11 1.
При переходе к пределу в равенстве (1.39) приближенное равен ство станет точным; получим:
Ш (* , y)ds = l\f(x{u, |
v), |
у {и, |
v))\I\ds' |
|
|
D |
D’ |
|
|
|
|
ИЛИ |
|
|
|
|
|
j f f(x, |
у) dxdy=\\ f (х(и, |
v), |
у (и, |
v))\I\dudv. |
(1-40) |
DЪ!
Спомощью этой формулы и осуществляется замена переменных
вдвойном интеграле.
Переход в двойном интеграле к полярным координатам. Переход от декартовых координат к полярным осуществляется, как известно, с помощью формул х == р cos ср; у = р sin ф. Вычислим якобиан этого преобразования:
дх |
ду |
|
р cos ф |
дер |
дф |
— р э т ф |
|
/ = |
ду |
COS ф |
— р sin2 ф— р cos2 ф = — р. |
дх |
ЭШф |
||
др |
др |
|
|
Подставив этот результат в формулу (1.40), получим:
Ш ( Х , |
у) dxdy = J f f(pcoscp, р sin ф) pdcpdp. |
(1-41) |
D |
D’ |
|
Это и есть формула перехода в двойном интеграле от декартовых координат к полярным.
27
Пример 7. Двойной интеграл J = J [ xydxdy, где D — четверть круга
D
х2 + у2 = R2, вычислить двумя способами: 1) в декартовой системе координат, 2) в полярной системе координат.
1. Расставим пределы интегрирования:
R |
Yr2—x- |
|
J = j dx |
[ xydy. |
(1.42) |
оо
Вычислим внутренний интеграл:
VR*—х1 |
Y R*-x°- |
xydy = |
= — (xR2 — x3). |
|
2 |
Воспользуемся формулой (1.41), имеем:
sin 2<ppd<pdpl = |
fp3 sin 2фйфф. |
||
D ' |
|
|
|
Заметим, что область |
D ограничена дугой |
окружности р = R |
|
и отрезками координатных |
осей |
= 0, ф = |
В системе коорди |
нат О'фр области D соответствует область D' — прямоугольник, изо браженный на рис. 17. Расставляя пределы по области D' , имеем:
|
|
|
л |
|
|
|
ТС |
|
|
R |
2 |
|
R |
1 |
~2 |
J |
— |
J |
dp j |
р3 sin 2cpdcp = — |
f p^ |
|
|
-------- cos 2ф |
|
||||||
|
2 |
о |
о |
2 |
о |
2 |
0 |
|
|
|
|
|
|
|
I P3rfP |
„ |
8 ’ |
о |
2 |
28
1.5.Приложения двойных интегралов
Вычисление площадей плоских фигур. Как уже [известно [см. (1.16)], площадь области D может быть выражена двойным интег ралом по этой области от функции f (х, у) = 1, т. е.
|
S = ttdxdy. |
|
|
|
(1.43) |
|
|
D |
|
|
|
|
|
Пример 8. Вычислить площадь круга х2 + |
y2< R 2. |
По формуле (1.43) |
||||
имеем: |
|
|
|
|
|
|
|
S = J J dxdy. |
|
|
|
|
|
|
D |
|
|
|
|
|
Перейдя к полярным координатам, получим: |
|
|||||
S = J J |
pd(pdp = |
R |
2л |
dtp = |
|
|
[ pdp f |
|
|||||
D' |
|
0 |
6 |
|
|
|
R |
r |
|
|
P2 |
R |
|
= .f P |
n |
pdp = |
2я |
2 |
о |
nR2. |
|
2n |
|
= |
Вычисление объемов. Пусть функция f (х, у) непрерывна в об ласти D. Рассмотрим цилиндроид, образованный функцией f (х, у) над областью D. Если в обла
сти D f(x, у)>>0, то цилиндроид имеет вид, изображенный на рис. 2, а объем цилиндроида, как было показано выше (1.5), может быть определен по фор муле
0 = |
Л 7 (* . У) dxdy. |
(1.44) |
|
D |
|
Если в области Df(x, у)<СО, |
||
то цилиндроид расположен ниже |
||
плоскости |
Оху и объем его |
мо |
жет быть определен по формуле
Р и с. 18
o = f Л fix, у) |dxdy. (1.45)
D
Заметим, что формула (1.45) является обобщением формулы (1.44) , так как при f (х, у) > 0 из формулы (1.45) следует формула
(1.44) .
Пусть в области D функция / (х, у) меняет знак и пусть, область D можно разбить на конечное число областей D lt D %, . . . , Dn, в каждой из которых функция / (х, у) имеет один и тот же знак. Для каждой из областей Dk применим формулу (1.45) и просуммируем результаты. Получим:
Л |
Я |
Г |
Т |
0 = 2 ] » * = |
2 |
1Л/(*> |
y)\dxdy\ . |
А=1 |
|
L Dk |
J |
29