Файл: Дьяченко, А. Н. Интегральное исчисление функций нескольких переменных. Теория вероятностей и элементы математической статистики учебное пособие.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 23.10.2024
Просмотров: 64
Скачиваний: 0
Вычислим «внутренний» интеграл:
у- |
|
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