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

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 23.10.2024

Просмотров: 64

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

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

у-

 

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