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

книги из ГПНТБ / Дьяченко, А. Н. Интегральное исчисление функций нескольких переменных. Теория вероятностей и элементы математической статистики учебное пособие

.pdf
Скачиваний:
37
Добавлен:
19.10.2023
Размер:
7.54 Mб
Скачать

соответствующее этому значению, является криволинейной трапе­ цией АВСЕ (см. рис. 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.

Соответственно, исходный двойной интеграл можно двумя спо­ собами представить в виде двукратного:

з

 

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):

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р

 

В этом случае расстановку в двойном интеграле следует выпол­ нить по формуле (1.27), получаем следующую формулу:

 

 

У2

 

ft

2р

[ j (х2 +

у2) dxdy = \ dy

j (х2 + у2) dx.

D

о

_ j/y_

 

 

23

Вычислим «внутренний» интеграл:

у-

 

JL

_

2р

 

(х2+ у2) dx = (— + У2х

2у6

 

24р3

jr_

 

 

 

 

Подставляя этот результат в предыдущую

2yi

12р3

формулу, получаем:

(х2+ у2) dxdy ■

( Ув

dy =

\ 12р3

 

84р3

 

И?

А5

 

 

84р3

 

Заметим, что если расставить пределы интегрирования другим способом (изменить порядок интегрирования), то формула для вычис­ ления двойного интеграла окажется более сложной. Действительно, границу области в этом случае необходимо будет представить в первой канонической форме:

х —

К2

h2

 

 

------;

х =

 

 

 

 

 

 

Y — 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\ _

и3---U% Vs--- /

25

х2 —х1 Уз У\

Вычислим произведение этих матриц:

АВ =

U2

Ui

щ —

и 3

и 2

щ —

 

Uj)

дх

,

 

dv

(U2 — Ul);

 

 

 

~ ( u 3 — u2) +

~

(v3— vt) ;

du

dv

 

 

(

дх

ду

du

 

ди

 

дх

ду

\

dv

dv

ду

(и 2—

“ i H

ди

 

 

dy_ (u3tl2)

du

ду_ (Vz—Vi)

dv

(V3— V2)

Заметим, что все элементы полученной матрицы — полные диф­ ференциалы функций х (и, v) я у (и, и).* Считая, что приращения аргументов и и v малы, заменим дифференциалы функций их при­ ращениями ** и получим следующее приближенное равенство:

 

 

/

дх

ду

\

 

 

/

Ujj—- Ux

V2 — V! \

ди

ди

| ^

Хх

у2

 

 

V3V2 )

 

 

[

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

dtp =

 

[ pdp f

 

D'

 

0

6

 

 

 

R

r

 

 

P2

R

 

= .f P

n

pdp =

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

Соседние файлы в папке книги из ГПНТБ