мера жордана. кратный интеграл

advertisement
Е.М. РУДОЙ
МАТЕМАТИЧЕСКИЙ АНАЛИЗ:
МЕРА ЖОРДАНА. КРАТНЫЙ ИНТЕГРАЛ
(теория, примеры, упражнения)
2009
2
Предисловие
Учебное пособие состоит из четырех глав, отражающих основные
теоретические и практические аспекты университетского курса математического анализа, читаемого во втором семестре второго курса для
студентов математического факультета Новосибирского государственного педагогического университета, обучающихся по специальности «информатикаматематика».
Большинство задач, включенных в пособие, содержится в качестве
упражнений и примеров в различных изданиях (см. список литературы). Часть примеров составлена специально для настоящего издания.
Нумерация упражнений сквозная и не зависит от номера главы и
параграфа. Сложные задачи помечены звездочкой.
При составлении упражнений автор попытался собрать задачи, соответствующие теоретическому курсу.
3
1. Мера Жордана.
Ранее мы изучили понятие определенного интеграл для функций
одной переменной. При интегрировании мы имели дело только с отрезками, длины которых определяется элементарно (как разность между
концом и началом отрезка).
Одним из основных приложений определенного интеграла является
измерение площадей плоских фигур. Причем, давая точное определение понятия «площади», мы использовали многоугольники как элементарные (простейшие) фигуры. Затем основные свойства их площадей
перенесли на квадрируемые фигуры (неотрицательность площади, конечная аддитивность, монотонность).
Во второй главе мы введем понятие кратного интеграла – интеграла
для функций n-переменных. Для этого нам надо дать точное определение меры множества в Rn . При этом в случае n = 2 мера должна
совпадать с площадью фигуры.
1.1 Сеть в Rn . Кубы и элементарные множества.
Определение. Множество вида
{x = (x1 , . . . , xn ) ∈ Rn | ai ≤ xi ≤ bi , i = 1, . . . , n},
где ai , bi – заданные числа, ai < bi , i = 1, . . . , n, называется n-мерным
брусом в пространстве Rn .
Если bi − ai = bj − aj для всех i, j = 1, . . . , n, то брус называется
n-мерным кубом.
Определение. Пусть k – натуральное число; разобьем каждую координатную ось на отрезки точками xi = m · 10−k , m = 0, ±1, ±2, . . . .
Таким образом, все пространство Rn «разобьется» на кубы вида
}
{
m+1
n m
Q = x = (x1 , . . . , xn ) ∈ R ≤ xi ≤
, i = 1, . . . , n .
10k
10k
В этом случае будем говорить, что в пространстве Rn построена сеть
порядка k, а кубы Q, входящие в эту сеть, будем называть кубами порядка k.
4
Отметим, что в одномерном случае (при n = 1) кубы Q являются
отрезками, в двумерном (при n = 2) — квадратами, в трехмерном (при
n = 3) — обычными трехмерными кубами.
x2
x1
Сеть в пространстве R2
Перечислим основные свойства сети, непосредственно следующие из
определения.
5
1. Различные кубы одного порядка не имеют общих внутренних точек. В пересечение двух кубов одного порядка могут входить только их
граничные точки.
2. Куб порядка k состоит из 10n кубов порядка k + 1.
3. Если G – ограниченное множество в Rn , то для любого натурального k существует конечное число кубов порядка k, объединение
которых содержит G.
Определение. Любой конечный набор попарно различных кубов
одного порядка называется элементарным множеством.
Замечание. В силу свойства 2 сети понятие элементарного множества не зависит от порядка k. Кроме того, объединение и пересечение
конечного числа элементарных множеств есть элементарное множество.
1.2 Мера элементарных множеств.
По аналогии с формулой для длины (n = 1), площади (n = 2) и
объема (n = 3) определим меру n-мерного куба.
Определение. Пусть Q – куб порядка k. Число
µQ = 10−nk
(1)
называется его мерой.
Определение. Пусть S – элементарное множество, т.е.
∪
S=
Qi ,
i
где Qi попарно различные кубы порядка k. Меру µS элементарного
множества S определим по формуле
∑
µS =
µQi .
i
Меру пустого множества будем считать равной нулю, т.е.
µ∅ = 0.
6
Непосредственно из определения меры элементарных множеств получаем следующие свойства:
1. Неотрицательность меры: µS ≥ 0 для любого элементарного
множества S. Причем µS = 0 тогда и только тогда, когда S = ∅.
2. Полуаддитивность меры: µ(P ∪ S) ≤ µP + µS для любых элементарных множеств P и S. Причем, если P и S не имеют общих кубов, то
µ(P ∪ S) = µP + µS
3. Монотонность меры: Если P ⊂ S, то µP ≤ µQ.
1.3 Мера Жордана произвольного ограниченного множества.
Пусть F – ограниченное множество в пространстве Rn . Выберем
натуральное k. Обозначим через Fk множество всех тех кубов порядка
k, которые целиком содержатся в F , а через F k – множество всех тех
кубов порядка k, имеющих непустое пересечение с F .
7
x2
Fk
Fk
x1
Справедливость следующих включений очевидна:
Fk ⊂ F ⊂ F k .
(2)
Отметим, что множества Fk и F k могут быть пустыми. Например, если
множество F есть точка в пространстве Rn , то не существует ни одного
куба, содержащегося в F . Случай же пустого F k возможен, когда F = ∅.
8
Кроме того, в силу монотонности меры элементарных множеств из
(2) следует, что
µFk ≤ µF k .
(3)
Ясно, что при увеличении k множества Fk возрастают, а F k убывают,
т.е. справедливы включения
F1 ⊂ F2 ⊂ . . . ⊂ Fk . . . ⊂ F ⊂ . . . ⊂ F k ⊂ . . . ⊂ F 1 .
(4)
Тем самым, мы имеем две неотрицательные монотонные числовые последовательности мер множеств Fk и F k , одна из которых {µFk } —
возрастающая, другая {µF k } — убывающая, такие, что
0 ≤ µF1 ≤ µF2 ≤ . . . ≤ µFk . . . ≤
≤ . . . ≤ µF k ≤ . . . ≤ µF 2 ≤ µF 1 .
(5)
Из (5) следует, что последовательность {µFk } ограничена сверху, а
{µF k } – снизу. Поэтому для любого ограниченного множества F существуют конечные пределы
µ∗ F = lim µFk ,
k→∞
µ∗ F = lim µF k .
k→∞
Определение. Числа µ∗ F и µ∗ F называются внутренней и внешней
мерами множества F соответственно.
Определение. Если µ∗ F = µ∗ F , то говорят, что множество F измеримо по Жордану. Общее значение его внешней и внутренней мер
называется n-мерной мерой Жордана и обозначается µF , т.е.
µF = µ∗ F = µ∗ F.
В дальнейшем для краткости будем, как правило, выражение «по
Жордану» опускать. Например, вместо «измеримое по Жордану множество» будем говорить «измеримое множество» и т.п.
Перечислим свойства меры, непосредственно следующие из определения:
9
1. Мера элементарного множества совпадает с введенной ранее мерой.
2. Если µ∗ F = 0, то множество F измеримо и µF = 0.
3. Мера µΠ произвольного бруса
Π = {x ∈ Rn | ai ≤ xi ≤ bi , i = 1, . . . , n}
равна произведению длин его сторон, т.е.
µΠ =
n
∏
(bi − ai ).
i=1
4. Мера µA любого измеримого множества A ⊂ Rn неотрицательна,
т.е.
µA ≥ 0.
Оказывается, что не всякое ограниченное множество является измеримым. Действительно, в пространстве R1 рассмотрим множество Q0
рациональных точек на отрезке [0, 1]. Для любого натурального k ни
один куб порядка k не содержится в Q0 . Следовательно, внутренняя
мера µ∗ Q0 равна нулю. С другой стороны, т.к. Q0 плотно в [0, 1], то
внешняя мера µ∗ Q0 равна 1.
Пример. Доказать, что куб Q ранга k0 в пространстве Rn измерим
по Жордану и его мера Жордана совпадает с введенной в (1) мерой,
т.е. равна 10−k0 n .
Решение. Рассмотрим сеть ранга k. Если k < k0 , то никакой куб
ранга k не содержится в Q. Поэтому, множество Fk тех кубов ранга
k, которые лежат в Q пустое, т.е. Fk = ∅ и, следовательно, µFk = 0.
Если k = k0 , то Fk0 = Q и µFk0 = 10−k0 n . Далее, если k > k0 , то опять
получаем, что Fk = Q и µFk = 10−k0 n . Поэтому внутренняя мера куба
Q равна
µ∗ Q = lim µFk = 10−k0 n .
k→∞
Вычислим теперь внешнюю меру куба Q. Данный куб Q по определению есть множеству
{
}
m
mi + 1
i
≤
x
≤
,
i
=
1,
.
.
.
,
n
,
Q = (x1 , . . . , xn ) i
10k0
10k0
10
где mi – некоторые целые числа, i = 1, . . . , n. Пусть k > k0 . Рассмотрим
объединение F k всех кубов ранга k, имеющих хотя бы одну общую точку
с Q. Очевидно, что
}
{
m
1
mi + 1
1
i
k
F = (x1 , . . . , xn ) − k ≤ xi ≤
+ k , i = 1, . . . , n
10k0
10
10k0
10
– куб с ребром длины
1
1
+2 k,
k
10
100
который содержит (10k−k0 + 2)n кубов ранга k. Поэтому
(
µF k = (10k−k0 + 2)n · 10−kn = 10−k0 n 1 +
и, следовательно,
2
)n
10k−k0
µ∗ Q = lim µF k = 10−k0 n .
k→∞
Здесь мы воспользовались тем, что отбрасывание конечного числа членов числовой последовательности не влияет на ее предел. Поэтому мы
не вычисляли меру F k при k ≤ k0 .
Таким образом, мы показали, что внешняя и внутренняя меры множества Q совпадают. Следовательно, куб Q измерим и его мера равна
10−k0 n , что совпадает с определением (1) меры куба.
Пример. Доказать, пользуясь определением меры Жордана, измеримость отрезка F = [a, b] ⊂ R1 , a < b, и найти его меру.
Решение. Для определенности будем считать, что a > 0. Представим
числа a и b в виде бесконечной десятичной дроби
a = a0 , a1 a2 a3 . . . ak . . . ,
b = b0 , b1 b2 b3 . . . bk . . . .
При этом справедливы неравенства
a0 , a1 a2 a3 . . . ak ≤ a ≤ a0 , a1 a2 a3 . . . ak +
1
,
10k
11
1
.
10k
Выберем натуральное k. Пусть Fk множество всех тех отрезков (кубы в одномерном случае) порядка k, которые целиком содержатся в
отрезке [a, b], а F k – множество всех тех отрезков порядка k, имеющих
непустое пересечение с [a, b]. Очевидно, что справедливы следующие
неравенства
b0 , b1 b2 b3 . . . bk ≤ b ≤ b0 , b1 b2 b3 . . . bk +
b0 , b1 . . . bk − a0 , a1 . . . ak −
2
≤ µFk ≤
10k
≤ µF k ≤ b0 , b1 . . . bk − a0 , a1 . . . ak +
2
.
10k
Устремляя k к бесконечности, получаем, что
b − a ≤ µ∗ F ≤ µ∗ F ≤ b − a.
Откуда следует измеримость отрезка [a, b]. При этом его мера µ[a, b]
равна b − a.
Задачи.
1. Доказать, что открытый куб Q порядка k в пространстве Rn является измеримым по Жордану множеством и его мера равна мере замкнутого куба ранга k, то есть µQ = 10−kn .
2. Найти меру следующих множеств в пространстве R1 :
а) A = {1},
б) X = {1, 2, . . . , n},
в) Y = { 21n | n ∈ N},
г)∗ Z = { 2mn | 1 ≤ m ≤ 2n − 1, k = 1, 2, 3, . . .}.
3. Доказать, что конечное число точек в пространстве Rn имеет меру
нуль.
4. Доказать, пользуясь определением меры Жордана, измеримость
следующих множеств и найти их меру:
а) интервал (a, b) в R;
б) замкнутый прямоугольник Π = {(x1 , x2 ) ∈ R2 | ai ≤ xi ≤ bi , i =
1, 2} в R2 ;
12
в) открытый прямоугольник Π = {(x1 , x2 ) ∈ R2 | ai < xi < bi , i =
1, 2} в R2 .
5∗ . Указать неограниченное множество с конечной внутренней мерой.
6∗ . Доказать, что
а) множество с положительной внутренней мерой имеет внутренние
точки;
б) мера измеримого по Жордану множества, не имеющего внутренних точек, равна нулю.
7. Пусть X1 ⊂ X2 ⊂ Rn . Доказать, что
µ∗ (X1 ) ≤ µ∗ (X2 ),
µ∗ (X1 ) ≤ µ∗ (X2 ).
8. Доказать, что если µ∗ (X) = 0, то множество X измеримо и его
мера равна нулю.
9. Пусть множество M измеримо и B ⊂ M . Следует ли, что B –
измеримое множество?
10. Пусть множество M – измеримо и его мера равна нулю. Может
ли множество M содержать внутренние точки?
11. Пусть множество M – измеримо и его мера равна нулю. Измеримо ли замыкание M множества M ?
12. Пусть замыкание M множества M измеримо. Измеримо ли множество M ?
13. Привести пример неизмеримого по Жордану множества в R2 .
14. Пусть X = {(x, y) | x ∈ [0, 1]∩Q, y = 0}. Измеримо ли множество
X в R2 ? Если – да, то найти его меру.
15. Пусть множество A измеримо в Rn . Измеримы ли проекции A
на оси координат в R1 ?
16∗ . Пусть последовательность {xk } сходится в Rn . Доказать, что
множество {xk ∈ Rn | k = 1, 2, 3 . . . , } измеримо. Найти его меру.
1.4 Критерий измеримости. Цилиндрические множества.
Сформулируем критерий измеримости множества по Жордану в
пространстве Rn .
13
Теорема. Множество измеримо по Жордану тогда и только тогда,
когда мера его границы равна нулю.
Используя эту теорему можно легко доказать, что множество Q0 рациональных точек на отрезке [0, 1] в пространстве R1 неизмеримо. Действительно, границей множества Q0 является весь отрезок [0, 1], мера
которого равна единице.
Определение. Пусть в Rn−1 задана область Ω. Пусть в области Ω
определены две непрерывные функции φ(x1 , . . . , xn−1 ), ψ(x1 , . . . , xn−1 )
такие, что φ(x1 , . . . , xn−1 ) ≤ ψ(x1 , . . . , xn−1 ) для всех (x1 , . . . , xn−1 ) ∈ Ω.
Цилиндрическим множеством с основанием Ω называется множество
C = {x = (x1 , x2 , . . . , xn ) ∈ Rn | (x1 , . . . , xn−1 ) ∈ Ω,
φ(x1 , . . . , xn−1 ) ≤ xn ≤ ψ(x1 , . . . , xn−1 )} .
Теорема. Пусть функция f (x1 , . . . , xn−1 ) непрерывна на измеримом, ограниченном и замкнутом множестве Ω ⊂ Rn−1 . Тогда множество
epi f = {x = (x1 , . . . , xn ) ∈ Rn |
(x1 , . . . , xn−1 ) ∈ Ω, xn = f (x1 , . . . , xn−1 )}
измеримо в Rn и его мера равна нулю.
Определение. Множество epi f ⊂ Rn называется графиком функции f (x1 , . . . , xn−1 ), (x1 , . . . , xn−1 ) ∈ Ω ⊂ Rn−1 .
Теорема. Если множество Ω – ограничено, замкнуто и измеримо в
Rn−1 , функции φ(x1 , . . . , xn−1 ) и ψ(x1 , . . . , xn−1 ) – непрерывны в Ω, то
цилиндрическое множество C с основанием Ω измеримо в Rn .
Задачи.
17. Доказать, что для любых множеств A, B ⊂ Rn справедливы
следующие включения:
а) ∂(A ∪ B) ⊂ ∂A ∪ ∂B,
б) ∂(A ∩ B) ⊂ ∂A ∪ ∂B,
в) ∂(A\B) ⊂ ∂A ∪ ∂B,
где ∂X — граница множества X.
14
18. Пусть мера множеств A и B равна нулю. Доказать, что множества A ∪ B, A ∩ B, A\B — измеримы и µ(A ∪ B) = 0, µ(A ∩ B) = 0,
µ(A\B) = 0.
19. Пусть A ∪ B и A ∩ B — измеримые множества. Верно ли, что A
и B — измеримы?
→
20. Пусть множество X ⊂ Rn – измеримо, −
a = (a1 , . . . , an ) – посто−
→
янный вектор. Определим множество X + a = {y = (y1 , . . . , yn ) | yi =
xi + ai , x = (x1 , . . . , xn ) ∈ X, i = 1, . . . , n}. Доказать, что множество
→
X +−
a – измеримо и его мера равна мере множества X.
21∗ . Доказать, что на числовой прямой существует ограниченное
открытое неизмеримое по Жордану множество.
22. Привести пример несчетного неизмеримого по Жордану множества, замыкание которого измеримо по Жордану.
23∗ . Доказать, что множество всех внутренних точек измеримого по
Жордану множества, измеримо.
24. Будут ли цилиндрическими следующие множества:
а) Ω = {(x, y) ∈ R2 | x2 + y 2 ≤ 1} ⊂ R2 ,
б) Ω = {(x, y, z) ∈ R3 | x2 + y 2 + z 2 ≤ 1} ⊂ R3 ,
в) Ω = {(x, y) ∈ R2 | x2 + y 2 = 1} ⊂ R2 ,
г) Ω = {(x, y) ∈ R2 | x ∈ [0, 1] ∩ Q, y ∈ [0, 1]} ⊂ R2 ,
д) Ω = {(x, y, z) ∈ R3 | x2 + y 2 ≤ 1, 0 ≤ z ≤ 1} ⊂ R3 ,
е) Ω = {(x, y, z) ∈ R3 | x2 + y 2 = 1, 0 ≤ z ≤ 1} ⊂ R3 ?
25. Измеримы ли множества из предыдущей задачи в соответствующих пространствах?
1.5 Основные свойства меры Жордана.
Свойства полуаддитивности, аддитивности и монотонности меры
Жордана для элементарных множеств непосредственно следовали из
определения их меры. Эти же свойства справедливы и для любых измеримых по Жордану множеств.
15
Теорема. Пусть множества A и B измеримы. Тогда множества A ∪
B, A ∩ B, A\B измеримы.
Теорема. Пусть множества A и B измеримы. Тогда имеет место
неравенство
µ(A ∪ B) ≤ µA + µB
(полуаддитивность меры).
Следствие. Пусть множества Ai – измеримы, i = 1, . . . k. Тогда
имеет место неравенство
k
k
) ∑
(∪
µ
Ai ≤
µAi .
i=1
i=1
Теорема. Пусть множества A и B измеримы и не имеют общих
внутренних точек. Тогда имеет место равенство
µ(A ∪ B) = µA + µB
(аддитивность меры).
Теорема. Пусть множества A и B измеримы и A ⊂ B. Тогда имеет
место неравенство
µA ≤ µB
(монотонность меры).
Задачи.
26. Привести пример двух множеств таких, что A ⊂ B, A ̸= B, но в
тоже время µA = µB.
27∗ . Пусть множества A и B измеримы. Доказать, что
а) µ(A ∪ B) = µ(A\B) + µ(B\A) + µ(A ∩ B);
б) µ(A ∪ B) = µ(A) + µ(B) − µ(A ∩ B).
28. Пусть Ak – измеримые по Жордану множества меру нуль, k =
1, . . . , p, и пусть
p
∪
A=
Ak .
k=1
16
Доказать, что µA = 0.
29. Пусть Ak – измеримые по Жордану множества меру нуль, k ∈ N,
и пусть
∞
∪
A=
Ak
k=1
– измеримое по Жордану множество. Доказать, что µA = 0.
30∗ . Привести пример таких измеримых множеств Ak , k ∈ N, мера
µAk которых равна нулю, но в тоже время
A=
∞
∪
Ak
k=1
– неизмеримое множество.
2. Кратные интегралы.
Ранее мы рассмотрели физические и геометрические задачи, приводящие к понятию определенного интеграла для функции одной переменной. Типичной задачей такого рода является задача о вычислении
площади криволинейной трапеции.
Легко указать аналогичные «многомерные» задачи, приводящие к
понятию двойного, тройного или n-кратного интегралов. Например, вычисление массы неоднородного тела T по известной объемной плотности
ρ(x). Для этого разобьем тело на малые участки Ti , i = 1, . . . , k. Приближенно будем считать объемную плотность ρ(x) каждого участка Ti
постоянной и равной ρ(xi ), где xi – некоторая точка, принадлежащая
Ti . Тогда приближенная масса тела будет задаваться формулой
k
∑
ρ(xi )|Ti |,
i=1
где |Ti | – объем участка Ti .
Точное значение массы тела естественно определить как предел указанной суммы при неограниченном уменьшении каждого участка Ti .
17
Этот предел и будет взят за определение тройного интеграла от функции ρ(x) по трехмерной области T . Конечно, для этого следует уточнить
термины «разбиение тела», «неограниченное уменьшение» и пр.
2.1 Разбиение измеримых множеств.
Пусть Ω – измеримое по Жордану множество в Rn .
τ
Определение. Конечная система τ = {Ωi }ki=1
измеримых множеств
Ωi , i = 1, . . . , kτ , называется разбиением множества Ω, если:
1. µ(Ωi ∩ Ωj ) = 0 для любых i ̸= j ;
∪ τ
2. ki=1
Ωi = Ω.
Замечание. Отметим, что множества Ωi разбиения τ могут попарно
пересекаться. Главное, чтобы мера такого пересечения равнялась нулю.
Определение. Число
diam Ω = sup ρ(x, y)
x,y∈Ω
называется диаметром множества Ω.
Определение. Число
|τ | = max {diam Ωi },
i=1,...,kτ
где diam Ωi – диаметр множества Ωi , назовем мелкостью разбиения τ .
k τ1
Определение. Пусть дано два разбиения τ1 = {Ω1j }j=1
и τ2 =
lτ
2
{Ω2i }i=1
множества Ω, и для каждого Ω2i ∈ τ2 существует такое Ω1j ∈ τ1 ,
что Ω2i ⊂ Ω1j . В этом случае будем говорить, что разбиение τ2 вписано
в разбиение τ , и писать τ2 ≻ τ1 .
Перечислим основные свойства разбиений множеств:
1. Если τ1 ≺ τ2 и τ2 ≺ τ3 , то τ1 ≺ τ3 (транзитивность).
2. Для любых двух разбиений τ1 и τ2 множества Ω существует такое
его разбиение τ , что
τ ≻ τ1 ,
τ ≻ τ2 (финальность).
18
τ
Лемма. Пусть τ = {Ωi }ki=1
– разбиение множества Ω. Тогда справедливо равенство
kτ
∑
µΩ =
µΩi .
i=1
Лемма. Для любого измеримого множества Ω и для любого ε > 0
существует разбиение τ множества Ω такое, что мелкость |τ | < ε.
2.2 Интегральные суммы. Определение кратного интеграла.
Пусть на измеримом множестве Ω ⊂ Rn задана функция y = f (x),
τ
x ∈ Ω. Пусть τ = {Ωi }ki=1
– разбиение множества Ω. Пусть ξ i – произвольная точка множества Ωi , i = 1, . . . , kτ .
Определение. Число
στ = στ (f ) = στ (f ; ξ 1 , . . . , ξ kτ ) =
kτ
∑
f (ξ i )µΩi
i=1
называется интегральной суммой Римана функции f (x), соответствующей разбиению τ .
Определение. Функция f (x) называется интегрируемой по Риману
на множестве Ω, если один и тот же предел имеет любая последовательность интегральных сумм
σ τm =
k τm
∑
f (ξ i,m )µΩm
i ,
i=1
k
τm
соответствующих разбиениям τm = {Ωm
i }i=1 множества Ω, у которых
их мелкость |τm | стремится к нулю при m → ∞, а точки ξ i,m выбраны
произвольным образом из множеств Ωm
i , i = 1, . . . , kτm , m = 1, 2, . . .
Если этот предел существует, то он называется интегралом Римана
от функции f (x) по множеству Ω и обозначается
∫
∫
∫
f (x) dx или
. . . f (x1 , . . . , xn ) dx1 . . . dxn .
Ω
Ω
19
Таким образом,
∫
f (x) dx = lim στm (f ; ξ 1,m , . . . , ξ kτm ,m ).
m→∞
Ω
∫
Последнее равносильно тому, что существует число, обозначаемое
f (x) dx, которое удовлетворяет следующему условию: для любого ε >
Ω
τ
0 найдется такое δ > 0, что для любого разбиения τ = {Ωi }ki=1
множеi
ство Ω мелкости |τ | и при любом выборе точек ξ ∈ Ωi , i = 1, . . . , kτ ,
имеет место неравенство
∫
f (x)dx − στ (f ; ξ 1 , . . . , ξ kτ ) < ε.
Ω
Кратко это можно записать в виде:
∫
f (x) dx = lim στ (f ; ξ 1 , . . . , ξ kτ ).
|τ |→0
Ω
Определение. Множество Ω, по которому производится интегрирование, называется областью интегрирования.
∫
Если n > 1, то интеграл f (x) dx называется кратным интегралом.
Ω
В случае n = 2 он называется двойным, в случае n = 3 – тройным, а в
случае произвольного натурального n — n-кратным.
Пример. Пусть Ω – измеримое множество в Rn ; пусть f (x) = 1
для всех x ∈ Ω. Вычислим интеграл от функции f (x) по области Ω.
τ
Для этого рассмотрим произвольное разбиение τ = {Ωi }ki=1
области Ω
и составим для него интегральную сумму
στ =
kτ
∑
i=1
f (ξ i )µΩi ,
20
где ξ i произвольная точка, принадлежащая Ωi , i = 1, . . . , kτ . Учитывая,
что f (ξ i ) = 1 для всех i = 1, . . . , kτ , и свойство разбиений, получаем
στ =
kτ
∑
µΩi = µΩ.
i=1
Таким образом, имеем
∫
1 dx = µΩ.
Ω
Замечание. Пусть Ω – множество меры нуль в пространстве Rn ,
n ≥ 2. Тогда любая функция f (x), определенная на Ω, интегрируема на
нем и
∫
f (x) dx = 0.
Ω
В частности, функция f (x) может быть и не ограничена на Ω. Отметим, что в случае функции одной переменной интегрируемая по Риману
функция должна быть ограничена.
Замечание. Можно показать, что в случае n = 1 и Ω = [a, b] – отрезок на числовой прямой R определения определенного интеграла, изученного ранее (когда рассматриваются интегральные суммы, соответствующие только разбиению отрезка на частичные отрезки) и в смысле
определения разбиения множества, данного в этом параграфе, равносильны, т.е. приводят к одному и тому же понятию интеграла.
Пример. Вычислить интеграл
∫ ∫
xy dxy,
0≤x≤1
0≤y≤1
рассматривая его как предел интегральной суммы, разбивая область
интегрирования на квадраты прямыми
x=
i
j
, y = , i, j = 0, . . . , n.
n
n
21
Задачи.
31. Вычислить диаметр следующих множеств:
а) Ω = {(x, y) ∈ R2 | x2 + y 2 ≤ 1},
б) Ω = {(x, y) ∈ R2 | x2 + y 2 = 1},
в) Ω = {(x, y) ∈ R2 | x2 + y 2 < 1},
г) Ω − треугольник с вершинами O(0,0), A(1,0), B(0,1),
д) Ω = {(x, y) ∈ R2 | x ∈ [0, 1], y = 0},
е) Ω = {(0, 0), (1, 1), (2, 0)} ⊂ R2 .
32. Вычислить диаметр n-мерного куба в пространстве Rn .
33. Пусть f (x) = C, где C − const, x ∈ Ω, Ω – измеримое в Rn
множество. Вычислить
∫
f (x) dx.
Ω
34. Вычислить интеграл
∫ ∫
(x2 + y 2 ) dxy,
1≤x≤2
1≤y≤3
рассматривая его как предел интегральной суммы, разбивая область
интегрирования на квадраты прямыми
x=1+
i
2j
, y = 1 + , i, j = 0, . . . , n.
n
n
35. Привести пример непрерывной, не равной тождественно нулю
функции f (x), x ∈ Ω ⊂ Rn , такой, что
∫
f (x) dx = 0.
Ω
22
36∗ . Пусть функция f (x), x ∈ Ω ⊂ Rn , – непрерывна и неотрицательна на Ω, µΩ > 0, и пусть
∫
f (x) dx = 0.
Ω
Доказать, что f (x) = 0 на Ω.
37∗ . Привести пример неограниченной на измеримом множестве Ω ⊂
Rn , n ≥ 2, функции f (x) такой, что интеграл
∫
f (x) dx
Ω
существует.
2.3 Верхние и нижние суммы Дарбу. Критерий интегрируемости.
По аналогии со случаем интеграла от функции одной переменной
введем понятия верхних и нижних сумм Дарбу, соответствующих данному разбиению τ множества.
Определение. Пусть функция f (x) ограничена на измеримом мноτ
жестве Ω ⊂ Rn ; τ = {Ωi }ki=1
– разбиение множества Ω,
mi = inf f (x), Mi = sup f (x), i = 1, 2, . . . , kτ .
x∈Ωi
x∈Ωi
Тогда суммы
sτ =
kτ
∑
i=1
mi µΩi , Sτ =
kτ
∑
Mi µΩi
i=1
называются нижней и верхней суммами Дарбу для данного разбиения
τ области Ω.
Для сумм Дарбу sτ , Sτ и интегральной суммы στ справедливы очевидные неравенства
sτ ≤ στ ≤ Sτ .
23
Кроме того, аналогично одномерному случаю можно доказать, что для
любых двух разбиений τ1 и τ2 множества Ω выполняется неравенство
sτ1 ≤ Sτ2 .
При заданной на множестве Ω функции f (x) суммы sτ и Sτ являются функциями разбиений τ множества Ω. Поэтому для них определим
понятия пределов
lim sτ , lim Sτ .
|τ |→0
|τ |→0
Определение. Пусть задана функция F (τ ), определенная на множестве всех разбиений τ измеримого множества Ω. Будем говорить, что
число a является пределом функции F (τ ) при |τ | → 0, если для любого ε > 0 существует такое δ > 0, что для всех разбиений τ , мелкость
которых |τ | < δ, имеет место
|F (τ ) − a| < ε.
Теорема. Пусть функция f (x) – ограничена на измеримом множестве Ω. Функция f (x) интегрируема по Риману на множестве Ω тогда
и только тогда, когда
lim (Sτ − sτ ) = 0.
|τ |→0
При выполнении этого условия имеют место равенства
∫
lim sτ = lim Sτ = f (x)dx.
|τ |→0
|τ |→0
Ω
Rn .
Определение. Пусть функция f (x) определена на множестве Ω ⊂
Число
ω(f ; Ω) = sup (f (x′ ) − f (x′′ ))
x′ ,x′′
называется колебанием функции f (x) на множестве Ω.
Используя это определение, последнюю теорему можно переформулировать в эквивалентном виде.
24
Теорема. Пусть функция f (x) – ограничена на измеримом множестве Ω. Функция f (x) интегрируема по Риману на множестве Ω тогда
и только тогда, когда
lim
|τ |→0
kτ
∑
ω(f ; Ωi ) µΩi = 0.
i=1
Сформулируем еще один критерий интегрируемости – критерий Дарбу.
Теорема. Пусть функция f (x) ограничена на измеримом множестве
Ω, тогда f (x) интегрируема на Ω тогда и только тогда, когда для любого
ε > 0 существует разбиение τ множества Ω такое, что
Sτ − sτ < ε.
Кратко последнее утверждение можно записать так:
∀ε > 0 ∃τ : Sτ − sτ < ε.
Отметим, что все эти утверждения доказываются аналогично одномерному случаю.
Одним из важнейших классов интегрируемых функций являются
непрерывные функции.
Теорема. Пусть функция f (x) определена и непрерывна на замкнутом, ограниченном и измеримом множестве Ω. Тогда f (x) интегрируема
на Ω.
Данная теорема является лишь достаточным признаком интегрируемости функции.Следующий пример показывает, что множество интегрируемых функций не ограничивается лишь непрерывными.
Пример. Показать, что функция
{
0, если y ∈ [0, 12 ],
f (x, y) =
1, если y ∈ ( 12 , 1]
интегрируема на квадрате Ω = [0, 1] × [0, 1]. Вычислить интеграл
∫ ∫
f (x, y) dxdy.
Ω
25
Решение. Покажем сначала, что функция f (x, y) интегрируема на
Ω. Для этого воспользуемся критерием Дарбу. Выберем произвольное
ε > 0 и найдем такое разбиение τ множества Ω, чтобы Sτ − sτ < ε.
Разобьем множество Ω на две части:
]}
{
[
1
,
Ω1 = (x, y) | x ∈ [0, 1], y ∈ 0,
2
{
(
]}
1
Ω2 = (x, y) | x ∈ [0, 1], y ∈
,1 ,
2
тогда
τ = {Ω1 , Ω2 }
искомое разбиение. Действительно, для Ω1 имеем M1 = m1 = 0, µΩ =
1
1
1
2 ; для Ω2 имеем M2 = m2 = 1, µΩ2 = 2 . Поэтому, Sτ = sτ = 2 и,
следовательно, 0 = Sτ − sτ < ε, ∀ε > 0.
Далее, так как функция f (x, y) интегрируема на Ω и для разбиения
τ ее верхняя и нижняя суммы Дарбу равны, то
∫ ∫
1
f (x, y) dxdy = .
2
Ω
Приведем пример неинтегрируемой на измеримом множестве функции. Рассмотрим в R2 квадрат Ω = [0, 1] × [0, 1] и положим
{
f (x, y) =
0, если x – рациональное,
1, если x – иррациональное.
Затем для любого натурального n рассмотрим разбиение τ квадрата
Ω на одинаковые квадраты Ωij со стороной длины 1/n, i, j = 1, . . . , n.
Так как множество рациональных чисел плотно в R, то для всех i, j =
1, . . . , n
mij =
inf
(x,y)∈Ωij
f (x, y) = 0, Mij =
sup
(x,y)∈Ωij
f (x, y) = 1.
26
Таким образом, получаем
Sτ − s τ =
n
∑
(Mij − mij )µΩij = 1
i, j=1
для любого натурального n. При n → ∞ мелкость |τ | разбиения τ стремиться к нулю, но разность Sτ −sτ равна единице. Следовательно, функция f (x, y) не интегрируема.
Задачи.
38. Составить верхнюю Sτ и нижнюю sτ суммы Дарбу для функции
f (x, y) = xy в области Ω = {(x, y) | 0 ≤ x, y ≤ 1}, где τ – разбиение
области Ω на квадраты прямыми
x=
i
j
, y = , i, j = 0, . . . , n.
n
n
39. Составить верхнюю Sτ и нижнюю sτ суммы Дарбу для функции
f (x, y) = x2 + y 2 в области Ω = {(x, y) | 1 ≤ x ≤ 2, 1 ≤ y ≤ 3}, где τ –
разбиение области Ω на прямоугольники прямыми
x=1+
i
2j
, y = 1 + , i, j = 0, . . . , n.
n
n
2.4 Свойства кратных интегралов.
На кратные интегралы от ограниченных функций переносятся все
основные свойства определенного интеграла функции одной переменной: линейность, аддитивность и пр.
∫
1. Если Ω – измеримое множество, то 1 · dx = µΩ.
Ω
2. Если функции fi (x) – интегрируемы на множестве Ω, то для любых действительных чисел αi , i = 1, . . . , k, функция
k
∑
i=1
αi fi (x)
27
также интегрируема на множестве Ω и
∫ ∑
k
αi fi (x) dx =
Ω i=1
k
∑
i=1
∫
fi (x) dx (линейность интеграла).
αi
Ω
3. Если Ω1 и Ω2 – измеримые множества, Ω1 ⊂ Ω2 , функция f (x)
ограничена и интегрируема на Ω2 , то она интегрируема и на множестве
Ω1 .
τ
4. Если Ω – измеримое множество, τ = {Ωi }ki=1
– его разбиение,
функция f (x) ограничена на множестве Ω, а ее сужения на Ωi интегрируемы на Ωi , i = 1, . . . , kτ , то функция f (x) интегрируема на Ω и
∫
f (x) dx =
Ω
kτ ∫
∑
f (x) dx (аддитивность интеграла).
i=1 Ω
i
5. Если функции f (x) и g(x) интегрируемы на множестве Ω и для
всех x ∈ Ω выполняется неравенство f (x) ≤ g(x), то
∫
∫
f (x) dx ≤ g(x) dx.
Ω
Ω
6. Если функции f (x) и g(x) интегрируемы и ограничены на некотором множестве Ω, то их произведение f (x) g(x) – интегрируемая на
Ω функция. Если к тому же
inf |g(x)| > 0,
x∈Ω
то и частное f (x)/g(x) интегрируемо на множестве Ω.
7. Если функция f (x) ограничена и интегрируема на Ω, то |f (x)|
интегрируема на Ω и справедливо неравенство
∫
∫
f (x)dx ≤ |f (x)|dx.
Ω
Ω
28
8. Пусть функция f (x) – ограничена, интегрируема и неотрицательна на измеримом множестве Ω, µΩ > 0. Пусть, кроме того, функция
f (x) непрерывна на Ω и существует такая точка x0 , что f (x0 ) > 0. Тогда справедливо неравенство
∫
f (x)dx > 0.
Ω
Из последнего свойства следует, что если f (x) – интегрируемая и
непрерывная на измеримом множестве Ω с положительной мерой функция такая, что
∫
|f (x)|dx = 0,
Ω
то f (x) ≡ 0 в Ω.
9. Если функции f (x) и g(x) интегрируемы на множестве Ω ⊂ Rn ,
m ≤ f (x) ≤ M , x ∈ Ω, и функция g(x) не меняет знак на Ω, то существует такое число λ ∈ [m, M ], что
∫
∫
f (x)g(x)dx = λ g(x)dx.
Ω
Ω
Следствие. Если множество Ω – открытое множество, любые две
точки которого можно соединить непрерывной кривой, функции f (x)
и g(x) – интегрируемы, функция f (x) – непрерывна и ограничена, а
функция g(x) не меняет свой знак на Ω, то в Ω существует такая точка
ξ, что
∫
∫
g(x)dx.
f (x)g(x)dx = f (ξ)
Ω
Ω
Задачи.
40∗ . Вычислить приближенно интеграл
∫ ∫
dx dy
√
,
x2 + y 2 + 24
x2 +y 2 ≤25
29
аппроксимируя область интегрирования системой вписанных квадратов, вершины которых находятся в целочисленных точках, и выбирая
значения подынтегральной функции в вершинах этих квадратов, наиболее удаленных от начала координат.
41∗ . Приближенно вычислить интеграл
∫ ∫
√
x + y dxdy,
Ω
где Ω – треугольник, ограниченной прямыми x = 0, y = 0 и x + y = 1,
разбив область Ω прямыми x = const, y = const, x+y = const на четыре
равных треугольника и выбрав значение подынтегральной функции в
центрах масс этих треугольников.
42. Пусть функции X(x) и Y (y) – непрерывны для всех x ∈ [a, b] и
y ∈ [c, d] соответственно. Доказать равенство
∫ ∫
∫b
X(x)dx ·
X(x)Y (y)dxdy =
a
R
∫d
Y (y)dy,
c
где R = [a, b] × [c, d].
43∗ . Какой знак имеют следующие двойные интегралы:
∫ ∫
а)
ln(x2 + y 2 )dxdy,
∫
|x|+|y|≤1
∫
б)
(
)1
1 − (x2 − y 2 ) 3 dxdy.
x2 +y 2 ≤4
44∗ . Пользуясь теоремой о среднем, оценить интеграл
∫∫∫
dxdydz
√
,
2
(x − a) + (y − b)2 + (z − c)2
x2 +y 2 +z 2 ≤R2
где a, b, c, R – константы такие, что a2 + b2 + c2 > R2 .
30
2.5 Сведение двойного интеграла к повторному.
Одним из способов вычисления кратных интегралов является способ
сведения их к повторным. Сначала рассмотрим случай двойного интеграла. Пусть заданы две функции φ(x) и ψ(x) непрерывные на отрезке
[a, b] такие, что φ(x) ≤ ψ(x) для всех x ∈ [a, b]. Определим множество
Ω = {(x, y) ∈ R2 | a ≤ x ≤ b, φ(x) ≤ y ≤ ψ(x)}.
Множество Ω является цилиндрических множеством в пространстве
R2 с основанием [a, b] и, поэтому, – измеримо.
Определение. Пусть функция f (x, y) задана на множестве Ω и при
каждом x ∈ [a, b] интегрируема на отрезке [φ(x), ψ(x)] по y. Тогда функция
ψ(x)
∫
(6)
F (x) =
f (x, y) dy
φ(x)
называется интегралом, зависящим от параметра x, а интеграл
∫b ( ψ(x)
∫
)
f (x, y) dy dx
a
φ(x)
называется повторным интегралом. Обычно повторный интеграл записывают в виде
ψ(x)
∫b
∫
dx
f (x, y) dy.
a
φ(x)
Теорема. (Случай прямоугольной области) Пусть функция f (x, y)
интегрируема на прямоугольнике P = [a, b] × [c, d]. Пусть для всех x ∈
[a, b] существует интеграл,
∫d
F (x) =
f (x, y)dy,
c
31
зависящий от параметра x. Тогда существует повторный интеграл
∫b
∫b
F (x) dx =
∫d
dx
a
a
f (x, y) dy
c
и справедливо равенство
∫ ∫
∫b
∫d
f (x, y) dxdy =
dx
f (x, y) dy.
a
P
(7)
c
Замечание. В этой теореме можно пометь x и y местами, т.е., если предположить существование интеграла, зависящего от параметра
y для всех y ∈ [c, d], то справедливо равенство
∫ ∫
∫d
f (x, y) dxdy =
∫b
dy
f (x, y) dx.
c
P
a
Используя последнюю теорему, можно доказать теорему о сведении
двойного интеграла к повторному для произвольных цилиндрических
областей.
Теорема. (Случай произвольной области) Пусть область Ω – цилиндрическая с непрерывными функциями φ(x) и ψ(x). Пусть функция
f (x, y) интегрируема на Ω и для любых x ∈ [a, b] существует интеграл,
зависящий от параметра
ψ(x)
∫
F (x) =
f (x, y) dy.
φ(x)
Тогда существует повторный интеграл
∫b
∫b
F (x) dx =
a
ψ(x)
∫
dx
a
φ(x)
f (x, y) dy,
32
и справедливо равенство
∫ ∫
∫b
f (x, y) dxdy =
dx
a
Ω
ψ(x)
∫
f (x, y) dy.
φ(x)
Замечание. В случаях, когда область Ω не является цилиндрической, часто удается разбить эту область на объединение конечного числа областей цилиндрического типа, не имеющих общих внутренних точек. В силу аддитивности интеграла Римана интеграл по всей области
равен сумме интегралов по этим областям.
Пример. Вычислить повторный интеграл
∫2
I=
∫2x
dx
ln y dy.
x
1
Решение. Применим метод интегрирования по частям к внутреннему
интегралу по переменной y:
∫2x
y=2x
ln y dy = (yln y − y) = 2x ln (2x) − 2x − x ln x + x =
y=x
x
= x ln x + 2x ln 2 − x
Следовательно, имеем
∫2
(x ln x + 2x ln 2 − x) dx =
I=
1
(
=
x2 ln x x2 (2 ln 2 − 1)x2
−
+
2
4
2
) x=2
= 5 ln 2 − 9 .
4
x=1
Пример. Вычислить интеграл
∫ ∫
I=
xy dxdy,
Ω
33
где Ω — квадрат [0, 1] × [0, 1].
Решение. Так как область интегрирования — квадрат, то для вычисления интеграла можно применить формулу (7). В результате будем
иметь
∫1
∫1
∫1
∫1
1
I=
dx xy dy = x dx · y dy = .
4
0
0
0
0
Пример. Переменить порядок интегрирования в повторном интеграле
2∫
sin x
∫π
I=
dx
f (x, y) dy.
0
0
Решение. Область Ω интегрирования задается неравенствами
0 ≤ x ≤ π, 0 ≤ y ≤ 2 sin x.
Проекция Ω на ось Oy есть отрезок [0, 2]. Каждая прямая y = const ∈
[0, 2] пересекает область Ω по отрезку с концами φ(y) и ψ(y), которые
находятся как решение уравнения y = 2 sin x из отрезка [0, π]: φ(y) =
arcsin (y/2), ψ(y) = π − arcsin (y/2). Следовательно, область Ω задается
неравенствами
0 ≤ y ≤ 2, arcsin (y/2) ≤ x ≤ π − arcsin (y/2).
Таким образом, имеем
π−arcsin
∫ (y/2)
∫2
I=
dy
0
f (x, y) dx.
arcsin (y/2)
Пример. Вычислить интеграл
∫ ∫
(x − y) dx,
Ω
34
где Ω — область, ограниченная линиями y = 2 − x2 , y = 2x − 1.
Решение. Найдем точки пересечения линий. Для этого необходимо
решить систему уравнений
{
y = 2 − x2 ,
y = 2x − 1,
из которой получаем, что A(−3, −7) и B(1, 1) — искомые точки. Таким
образом, область Ω — цилиндрическая с основанием [−3, 1], лежащем
на оси Ox. Следовательно, находим
∫ ∫
∫1
(x − y) dx =
dx
−1
Ω
2−x
∫ 2
(x − y) dy =
2x−1
) y=2−x2
)
∫1 (
∫1 (
1 2 3
64
1 4
3
2
=
xy − y dx = .
− x − x + 2x + x −
dx =
2
2
2
15
−3
y=2x−1
−3
Пример. Вычислить двойной интеграл
∫ ∫ √
4
I=
1 − y 2 dxdy,
Ω
где Ω — треугольник, ограниченный прямыми x = 0, y = 1, y = x.
Решение. Область Ω является цилиндрической с основанием [0, 1],
лежащем на оси Ox, нижней стороной y = x, верхней стороной y = 1.
Следовательно,
∫1
∫1 √
4
I=
dx
1 − y 2 dy.
0
x
Известно, что внутренний интеграл не является элементарной функцией.
35
В то же время, область Ω — цилиндрическая с основанием [0, 1],
лежащем на оси Oy, верхней и нижней сторонами x = y и x = 0 соответственно. Следовательно,
∫1
I=
dy
0
∫y √
4
1−
y 2 dx
=
0
∫1 √
4
2
1 − y 2 y dy = .
5
0
Данный пример показывает, что выбор порядка интегрирования при
сведении кратного интеграла к повторному может существенно облегчить вычисление интеграла.
Задачи.
45. Вычислить
∫2π
∫1
y cos2 xdy.
dx
46. Вычислить
0
0
∫1
∫3
y(x3 − y)dy.
dx
0
47. Вычислить
2
∫x
∫3
(x − y)dy.
dx
1
48. Вычислить
x3
∫e2y
∫1
dy
0
ln xdy.
ey
49. Вычислить следующие интегралы:
∫1 ∫1
а)
∫1 ∫1
x dxdy,
0
0
б)
xy dxdy,
0
0
36
∫π ∫π
в)
∫1 ∫1
0
0
0
∫π ∫π
д)
0
∫1 ∫1
x sin (xy) dxdy,
1
(x2 + y) dxdy,
г)
sin x cos y dxdy,
1
xp y q dxdy, p, q > 0.
е)
0
0
∫ ∫
50. Вычислить
(x2 + 2y) dxdy,
D
если область D ограничена прямыми x = 0, x = 1, y = 0, y = 2.
51. Вычислить
∫ ∫
dxdy
,
(x + y)2
D
если область D ограничена прямыми x = 3, x = 4, y = 1, y = 2.
52. Вычислить
∫ ∫
xy dxdy,
D
если область D ограничена прямыми x = 0, y = 0, x + y = 1.
53. Вычислить
∫ ∫
(cos 2x + sin y) dxdy,
D
если область D ограничена прямыми x = 0, y = 0, 4x + 4y − π = 0.
54. Вычислить
∫ ∫
sin(x + y) dxdy,
D
если область D ограничена прямыми x = 0, y = π/2, y = x.
55. Вычислить
∫ ∫
x dxdy,
D
если область D – треугольник с вершинами O(0, 0), A(1, 0), B(1, 1).
37
∫ ∫
56. Вычислить
y dxdy,
D
если область D – треугольник с вершинами O(0, 0), A(2, 1), B(−2, 1).
57. Вычислить
∫ ∫
(x + y) dxdy,
D
если область D – трапеция с вершинами O(0, 0), A(1, 0), B(1, 2), C(0, 1).
58. Вычислить
∫ ∫
2y dxdy,
D
√
если область D ограничена параболой y = x и прямыми y = 0, x + y =
2.
59. Вычислить
∫ ∫ 2
x
dxdy,
y2
D
если область D ограничена прямыми x = 2, y = x и гиперболой y = 1/x.
60. Вычислить
∫ ∫
ex dxdy,
D
если область D ограничена прямыми x = 0, y = 1, y = 2 и кривой
x = ln y.
61. Изменить порядок интегрирования в двойных интегралах:
∫1
а)
∫x2
dx
0
∫1
в)
0
∫2
f (x, y)dy,
б)
y
dx
3−2y
∫
∫1
г)
f (x, y)dy,
x
0
x3
2−y
∫
f (x, y)dx,
dy
∫2x
dy
0
f (x, y)dx,
√
y
38
∫1
д)
∫ex
f (x, y)dy,
dx
е)
e−x
0
√
∫2ax
∫2a
dx
0
∫1
1−y
∫
ж)
dy
0
√
−
f (x, y)dy,
√
2ax−x2
∫2
f (x, y)dx,
∫y
з)
dy
1
1−y 2
f (x, y)dx.
ln y
62∗ . Вычислить повторные интегралы:
√
∫1
∫1−x2
∫π ∫π
sin x
dx
(1 − y 2 )3/2 dy, б)
dy
dx,
x
а)
0
0
в)
∫1
1/5
y
∫
∫1
(1 − x )
3 1/2
dy
0
y 1/2
y
0
dx, г)
∫1
y2
dx
0
√
3
√
y 4 − x2 dy.
x
63. Вычислить следующие интегралы:
∫ ∫
а)
x3 y 5 dxdy, где D = {(x, y) | |x| + |y| ≤ 1},
D
∫ ∫
x2 dxdy, где D = {(x, y) | |x| + |y| ≤ 1},
б)
в)
г)∗
∫ ∫
∫ ∫ √
D
|x − y 2 | dxdy, где D = {(x, y) | |y| ≤ 1, 0 ≤ x ≤ 2},
D
max{sin x, sin y} dxdy, где D = {(x, y) | 0 ≤ x ≤ π, 0 ≤ y ≤ π},
D
д)
∗
∫ ∫
sign(x2 − y 2 − 2) dxdy, где D = {(x, y) | x2 + y 2 ≤ 4}.
D
39
2.6 Сведение тройного и n-кратного интеграла к повторным.
Сформулируем теоремы о сведении n-кратного интеграла к повторному в случае n ≥ 3.
Теорема. Пусть функция f (x1 , . . . , xn−1 , xn ) определена и интегрируема на цилиндрическом множестве
Ωn = {(x1 , . . . , xn ) ∈ R2 | φ(x1 , . . . , xn−1 ) ≤ xn ≤ ψ(x1 , . . . , xn−1 )},
где функции φ(x1 , . . . , xn−1 ) и ψ(x1 , . . . , xn−1 ) заданы на ограниченном
и замкнутом множестве Ωn−1 ⊂ Rn−1 и непрерывны на нем. Пусть для
любых точек (x1 , . . . , xn−1 ) ∈ Ωn−1 существует интеграл, зависящий от
параметров (x1 , . . . , xn−1 )
ψ(x1 ,...,x
∫ n−1 )
F (x1 , . . . , xn−1 ) =
f (x1 , . . . , xn−1 , xn )dxn .
φ(x1 ,...,xn−1 )
Тогда существует (n − 1)-кратный интеграл
∫ ∫
∫
...
ψ(x1 ,...,x
∫ n−1 )
dx1 dx2 . . . dxn−1
Ωn−1
f (x1 , . . . , xn−1 , xn )dxn
φ(x1 ,...,xn−1 )
и справедлива формула
∫ ∫
∫
. . . f (x1 , . . . , xn−1 , xn )dx1 . . . dxn =
Ωn
∫ ∫
=
∫
...
Ωn−1
ψ(x1 ,...,x
∫ n−1 )
f (x1 , . . . , xn−1 , xn )dxn .
dx1 dx2 . . . dxn−1
φ(x1 ,...,xn−1 )
Рассмотрим случай n = 3. Пусть функция трех переменных f (x, y, z)
определена на множестве Ω ⊂ R3 и непрерывна на нем. Пусть Ω – цилиндрическое множество вида
Ω = {(x, y, z) ∈ R2 | (x, y) ∈ Ωxy , φ(x, y) ≤ z ≤ ψ(x, y)},
40
где Ωxy – ограниченное и замкнутое множество в R2 , функции φ(x, y) и
ψ(x, y) непрерывны на Ωxy . Тогда к интегралу
∫ ∫ ∫
f (x, y, z)dxdydz
Ω
применима теорема о сведении кратного интеграла к повторному, т.е.
справедливо равенство
∫ ∫ ∫
∫ ∫
f (x, y, z)dxdydz =
Ω
ψ(x,y)
∫
dxdy
Ωxy
f (x, y, z)dz.
(8)
φ(x,y)
Если в свою очередь множество Ωxy тоже является цилиндрическим
множеством вида
Ωxy = {(x, y) ∈ R2 | a ≤ x ≤ b, φ1 (x) ≤ y ≤ ψ1 (x)},
где функции φ1 (x) и ψ1 (x) непрерывны на отрезке [a, b], то, применив
в правой части формулы (8) формулу сведения двукратного интеграла
по множеству Ωxy к повторному, получим
∫ ∫ ∫
ψ∫1 (x)
∫b
dx
f (x, y, z)dxdydz =
a
Ω
ψ(x,y)
∫
dy
φ1 (x)
f (x, y, z)dz.
(9)
φ(x,y)
Если обозначить через Ω(x0 ) сечение множества Ω плоскостью {x =
x0 }, т.е.
Ω(x0 ) = Ω ∩ {(x, y, z) ∈ R3 | x = x0 },
то при условии x ∈ [a, b] включение (x, y, z) ∈ Ω(x) равносильно включениям φ1 (x) ≤ y ≤ ψ1 (x) и φ(x, y) ≤ z ≤ ψ(x, y). Поэтому, объединив
в формуле (9) два внутренних интегрирования по переменным y и z,
получим формулу
∫ ∫ ∫
∫b
f (x, y, z)dxdydz =
Ω
∫ ∫
dx
a
f (x, y, z)dydz.
Ω(x)
41
Например, если f (x, y, z) = 1, то
∫ ∫ ∫
∫ ∫
dxdydz = µ3 Ω,
dydz = µ2 Ω(x),
Ω
Ω(x)
где µ3 – трехмерная мера (объем), µ2 – двумерная мера (площадь).
Тогда получаем
∫b
µ3 Ω = µ2 Ω(x)dx,
a
т.е. объем тела Ω равен одномерному интегралу от площадей сечений
Ω(x).
Пример. Вычислить интеграл
∫ ∫ ∫
z dxdydz,
I=
Ω
где область
√ Ω определяется неравенствами 0 ≤ x ≤ 1/2, x ≤ y ≤ 2x,
0 ≤ z ≤ 1 − x2 − y 2 .
Решение. Область Ω — цилиндрическая с основанием Ωxy , где
Ωxy = {(x, y) ∈ R2 | 0 ≤ x ≤ 1/2, x ≤ y ≤ 2x}.
В свою очередь Ωxy также является цилиндрической. Поэтому можно
применить формулу (9). Таким образом, находим
1
∫2
I=
dx
0
1
=
1
2
∫2
0
√
1
z dz =
2
dy
x
1
2
2
1−x
∫ −y
∫2x
0
∫2
∫2x
dx
0
1
z=√1−x2 −y2
z dy =
2
z=0
x
)
∫2 (
∫2x
1 2 y=2x
1
2
2
2
dx (1 − x − y ) dy =
y − yx − y dx =
2
3
y=x
x
0
42
1
)
∫2 (
1
10 3
7
=
x− x
dx =
.
2
3
192
0
Пример. Вычислить интеграл
∫ ∫ ∫
I=
x2 yz dxdydz,
Ω
где область Ω ограничена плоскостями x = 0, y = 0, z = 0, x+y +z −2 =
0.
Решение. Проекция области Ω на плоскость Oxy есть треугольник,
ограниченный прямыми x = 0, y = 0, y = 2 − x. Сверху область Ω
ограничена плоскостью z = 2 − x − y, снизу — z = 0. Таким образом,
имеем
∫2
2−x−y
2−x
2−x
∫
∫
∫2
∫
1
2
z dz =
x dx
y dx
x dx
y(2 − x − y)2 dy =
2
2
I=
0
0
0
1
=
24
0
∫2
x2 (2 − x)4 dx =
0
16
.
315
0
Пример. Вычислить тройной интеграл
∫ ∫ ∫
I=
y dxdydz,
Ω
где
Ω = {(x, y, z) | 0 ≤ z ≤ 1, |x| ≤ z, z ≤ y, x2 + y 2 + z 2 ≤ 4}.
Решение. Неравенства |x| ≤ z, 0 ≤ z ≤ 1 задают треугольник Ω′ на
плоскости Oxz. Решим исходную систему неравенств относительно y.
43
Из первого и третьего неравенства следует, что y ≥ 0, поэтому третье
и четвертое неравенства равносильны системе
√
z ≤ y ≤ 4 − x2 − z 2 .
√
Эта система имеет решение при условии, что 4 − x2 − z 2 ≥ z. Этому условию удовлетворяют все точки треугольника Ω′ . Следовательно,
область Ω — цилиндрическая
с основанием Ω′ , нижней и верхней сторо√
2
2
нами y = z и y = 4 − x − y соответственно. В свою очередь область
Ω′ на плоскости Oxz тоже является цилиндрической
Ω′ {(x, z) | 0 ≤ z ≤ 1, −z ≤ x ≤ z}.
Следовательно,
dz
I=
0
√
∫z
∫1
−z
4−x
∫ 2 −z 2
1
y dy =
2
dx
z
∫z
∫1
dz
0
(4 − x2 − 2z 2 ) dx =
−z
)
∫1 (
7 3
17
=
4z − z
dz = .
3
12
0
Пример. Пусть Πn — пирамида в Rn ,
Πn = {(x1 , . . . , xn ) | 0 ≤ xn ≤ xn−1 ≤ . . . ≤ x2 ≤ x1 ≤ 1}.
Вычислить n-кратный интеграл
∫
In = x1 x2 . . . xn dx.
Πn
Решение. Из уравнений, задающих пирамиду, следует, что Πn — является цилиндрической областью, которую можно представить в виде
Πn = {x1 , . . . , xn ) | 0 ≤ x1 ≤ 1, 0 ≤ x2 ≤ x1 ,
44
0 ≤ x3 ≤ x2 , . . . , 0 ≤ xn ≤ xn−1 }.
Поэтому, сводя кратный интеграл к повторному и применяя метод математической индукции, получим
∫x1
∫1
In =
x∫n−1
dx1
0
dx2 . . .
x1 . . . xn−1 xn dxn =
0
0
Задачи.
64. Вычислить следующие повторные интегралы:
∫1
∫2
а)
dx
0
0
0
∫2
∫2
∫xy
1
dx
y
0
∫0
г)
dy
д)
∫3
y 2 cos x dx,
dz
0
√
z∫2 −x2
∫z
dz
0
∫zy
y
∫4
dz
,
x(1 + x2 y 2 z 2 )
0
∫0
−1
1
0
0
dy
xyz dz,
∫xy
dy (x + y + z) dz,
∫x
dx
в)
е)
dy
0
∫1
б)
∫3
xy 2 z 2 dz,
dx
−z
0
3−y−z
3−z
∫
∫
dz
dy
1−z
0
dx
.
(x + y + z)2
1
.
(2n)!!
45
65. Вычислить интегралы:
∫1 ∫1 ∫1
а)
xyz dxdydz,
0 −1 0
∫1 ∫1 ∫1
б)
(xy + yz + xz) dxdydz,
0
0 0
π
∫ ∫π ∫π
в)
sin(x + y + z) dxdydz.
0
0
0
66. Вычислить тройной интеграл
∫ ∫ ∫
x2 y 2 z dxdydz,
D
если область D ограничена плоскостями x = 1, x = 3, y = 0, y = 2,
z = 2, z = 5.
67. Вычислить тройной интеграл
∫ ∫ ∫
dxdydz
,
(x + y + z + 1)3
D
если область D ограничена плоскостями x = 0, y = 0, z = 0, x+y+z = 1.
68. Вычислить тройной интеграл
∫ ∫ ∫
xyz dxdydz,
D
если область D ограничена плоскостями x = 0, y = 0, z = 0 и сферой
x2 + y 2 + z 2 = 1.
69. Вычислить тройной интеграл
∫ ∫ ∫
x dxdydz,
D
46
если область D ограничена плоскостями z = 0, z = 3 и цилиндром
x2 + y 2 = 1.
70∗ . Вычислить следующие тройные интегралы:
∫ ∫ ∫
xy 2 z 3 dxdydz,
а)
Ω
если область Ω ограничена поверхностями z = xy, y = x, x = 1, z = 0;
∫ ∫ ∫ (
б)
x2 y 2 z 2
+ 2 + 2
a2
b
c
)
dxdydz,
Ω
если область Ω ограничена поверхностью
x2 y 2 z 2
+ 2 + 2 = 1, a, b, c > 0,
a2
b
c
∫ ∫ ∫ √
x2 + y 2 dxdydz,
в)
Ω
если область Ω ограничена поверхностями
x2 + y 2 = z 2 ,
z = 1.
71. Различными способами расставить пределы интегрирования в
следующих повторных интегралах:
∫1
а)
x+y
1−x
∫
∫
f (x, y, z) dz,
dx
dy
0
0
√
∫1
б)
∫1−x2
dx
−1
0
√
− 1−x2
∫1
dy
√
f (x, y, z) dz,
x2 +y 2
47
∫1
в)
∫1
dx
0
2 +y 2
x∫
f (x, y, z) dz,
dy
0
∫1
0
1−x−y
1−x
∫
∫
dx
dy
f (x, y, z) dz,
г)
0
0
∫1
д)
0
dx
0
0
∫2
dx
0
f (x, y, z) dz,
dy
0
∫1
е)
∫y
∫x
1−|y−1|
∫
dy
f (x, y, z) dz.
0
0
72∗ .
Упростить следующие выражения, заменяя тройные интегралы
однократными:
∫x
а)
∫ξ
dξ
∫η
dη
0
0
∫1
б)
f (ζ) dζ,
0
∫1
dx
0
x+y
∫
dy
0
f (z) dz.
0
73∗ . Пусть f (x1 , . . . , xn ) – непрерывная функция в области
Ω = {(x1 , . . . , xn ) | 0 ≤ xi ≤ x, i = 1, . . . , n}.
Доказать равенство
∫x1
∫x
dx2 . . .
dx1
0
x∫n−1
0
∫x
0
∫x
0
∫x
dxn−1 . . .
dxn
=
f (x1 , . . . , xn ) dxn =
xn
f (x1 , . . . , xn ) dx1 .
x2
48
74∗ . Вычислить следующие многократные интегралы:
∫1 ∫1
а)
∫1
(x21 + x22 + . . . + x2n ) dx1 dx2 . . . dxn ,
...
0
0
0
∫1 ∫1
б)
∫1
(x1 + x2 + . . . + xn )2 dx1 dx2 . . . dxn .
...
0
0
0
75∗ .
Пусть Πn — пирамида в пространстве Rn , определенная ранее.
Вычислить следующие интегралы
∫
а)
1 · dx,
Πn
∫
б)
(x1 + x2 + . . . + xn ) dx.
Πn
2.7 Замена переменных в кратных интегралах.
В этом параграфе рассмотрим теорему о замене переменных в кратных интегралах, которая является одним из важнейших средств вычисления кратных интегралов.
Пусть функция f (y1 , . . . , yn ) интегрируема в Ω, т.е. существует интеграл
∫
∫ ∫
f (y) dy =
Ω
...
f (y1 , . . . , yn ) dy1 . . . dyn ,
(10)
Ω
где Ω – ограниченная замкнутая измеримая область в пространстве Rn
переменных y1 , . . . , yn . Предположим далее, что переменные y1 , . . . , yn
в свою очередь зависят от переменных x1 , . . . , xn , т.е. задано преобразование координат

y1 = ψ1 (x1 , . . . , xn ),



y2 = ψ2 (x1 , . . . , xn ),
(11)
..................,



yn = ψn (x1 , . . . , xn ).
49
Для краткости преобразование (11) будем обозначать символом
y = ψ(x),
понимая, что x = (x1 , . . . , xn ), y = (y1 , . . . , yn ) – точки n-мерного пространства, а ψ – совокупность n функций ψ1 , . . . , ψn .
Обозначим через Ω′ такую область в пространстве переменных x,
которая при действии преобразования (11) переходит в Ω. Кратко это
будем обозначать так: Ω = ψ(Ω′ ). Кроме того, будем считать, что преобразование (11) взаимно однозначно, т.е. существует x = ψ −1 (y) такое,
что Ω′ = ψ −1 (Ω).
Определение. Пусть в каждой точке x = (x1 , . . . , xn ) области Ω′ существуют все частные производные первого порядка от функций ψ1 (x), . . . , ψn (x).
Тогда матрица
 ∂ψ

∂ψ2
n
1
. . . ∂ψ
∂x1
∂x1
∂x1
 ∂ψ1 ∂ψ2
∂ψ 
 ∂x2 ∂x2 . . . ∂xn2 


 ...

∂ψn
∂ψ1
∂ψ2
.
.
.
∂xn
∂xn
∂xn
называется матрицей Якоби, а ее определитель – Якобианом, который
обозначается
D(y)
D(y1 , . . . , yn )
или
.
D(x)
D(x1 , . . . , xn )
Сформулируем теорему о замене переменных в кратных интегралах.
Теорема. Пусть преобразование (11) отображает область Ω′ в Ω взаимно однозначно. Пусть функции ψ1 , . . . , ψn имеют непрерывные частные производные первого порядка. Пусть Якобиан преобразования (11)
отличен от нуля в любой точке области Ω′ . Пусть функция f (x) интегрируема на Ω. Тогда справедлива формула замены переменных
∫
∫
D(y) f (y)dy = f (ψ(x)) (12)
dx.
D(x)
Ω′
Ω
В подробной записи формула (12) имеет вид
∫ ∫
∫
. . . f (y1 , . . . , yn ) dy1 dy2 . . . dyn =
Ω
50
∫ ∫
=
∫
...
f (ψ1 (x1 . . . , xn ), . . . , ψn (x1 , . . . , xn ))×
Ω′
D(y , . . . , y ) 1
n ×
dx1 dx2 . . . dxn .
D(x1 , . . . , xn )
Отметим, что доказательство этой теоремы разбивается на несколько этапов: сначала справедливость формулы (12) доказывается для линейного преобразования координат, а затем сводится к общему случаю.
Замечание. В условиях теоремы о замене переменных в кратном
интеграле можно допустить обращение в нуль якобиана преобразования
на некотором множестве S, принадлежащему Ω′ , при условии, что nмерная мера множества S равна нулю.
Замечание. Величину dy1 dy2 . . . dyn называют элементом объема в
декартовой системе координат Oy1 y2 . . . yn . С помощью преобразования
(11) мы переходим от декартовых координат
Oy1 y2 . . . yn к новым криволинейным координатам Ox1 x2 . . . xn . Интеграл
∫ ∫
∫
I=
...
1 dy1 dy2 . . . dyn
Ω′
равен n-мерному объему области Ω′ . Используя формулу замены переменных в кратных интегралах (12), получаем
∫ ∫
∫ D(y) I=
... dx1 dx2 . . . dxn .
D(x)
Ω
Поэтому величину
D(y) dx1 dx2 . . . dxn
D(x)
называют элементом объема в криволинейной системе координат Ox1 x2 . . . xn .
51
2.8 Криволинейные координаты.
Взаимно однозначное отображение

 y1 = y1 (x1 , . . . , xn ),
..................
y = F (x) =

yn = yn (x1 , . . . , xn ),
x ∈ Ω,
(13)
которое отображает множество Ω ⊂ Rnx в пространство Rny ставит в
соответствие каждой точке x = (x1 , . . . , xn ) набор чисел y1 , y2 , . . . , yn .
Эти числа можно рассматривать как другие координаты той же точки x. Такие координаты называют криволинейными, так как кривые,
получающиеся из уравнений
yi = yi (x1 , . . . , yn ),
i = 1, . . . , n,
при фиксировании всех аргументов, кроме одного, называемые координатными кривыми, не являются, вообще говоря, прямыми, в отличие
случае декартовых координат. Отметим, что в случае декартовых координат уравнения (13) имеют вид
y1 = x1 , y2 = x2 , . . . , yn = xn .
Таким образом, формулы (13) можно рассматривать не только как
формулы, задающие отображение из одного пространства в другое, но
и как переход от координат x1 , . . . , xn к координатам y1 , . . . , yn в том же
пространстве. В этом случае множество Ω′ = F (Ω) представляет собой
множество наборов новых координат точек множества Ω.
Полярная система координат. Из криволинейных координат на плоскости отметим полярные координаты r, φ, связанные с декартовыми x,
y соотношениями
x = r cos φ, y = r sin φ.
(14)
В декартовой системе координат координаты x, y точки (x, y) могут
принимать значения от −∞ до +∞. В полярной системе координат
r ∈ [0, ∞), φ ∈ [0, 2π). В этом случае отображение (14) является взаимно
52
однозначным. Координатными линиями полярных координат являются
концентрические окружности r = const с центром в начале координат
и лучи φ = const, выходящие из начала координат.
Якобиан преобразования (14) равен r, т.е.
D(x, y)
= r.
D(r, φ)
Заметим, что якобиан обращается в ноль в точке (0, 0), но 2-мерная
мера точки равна нулю. Поэтому в силу замечания к теореме о замене
переменных в кратном интеграле формула (12) для замены переменных
(14) верна.
Из криволинейных координат в пространстве R3 отметим цилиндрическую и сферическую.
Цилиндрическая система координат. Цилиндрические координаты
r, φ, h связаны с декартовыми x, y, z соотношениями
x = r cos φ, y = r sin φ, z = h,
(15)
r ≥ 0, 0 ≤ φ < 2π, −∞ < h < +∞.
Фактически цилиндрические координаты – это полярные координаты
на плоскости Oxy и обычная система координат в ортогональном дополнении плоскости Oxy – оси Oz. Легко вычислить якобиан преобразования (15)
cos φ −r sin φ 0 D(x, y, z) = sin φ r cos φ 0 = r.
D(r, φ, h) 0
0
1 Сферическая система координат. Сферические координаты r, φ, ψ
связаны с декартовыми x, y, z соотношениями
x = r cos φ cos ψ, y = r sin φ cos ψ, z = r sin ψ,
r ≥ 0, 0 ≤ φ < 2π, −
π
π
≤ψ≤ .
2
2
(16)
53
Так же, как и полярные координаты (r, φ) на плоскости, сферические
координаты (r, φ, ψ) точки M в пространстве R3 имеют простой геометрический смысл: r – длина радиуса-вектора из начала координат в
точку M , φ – угол проекции радиуса-вектора на плоскость Oxy с положительным направлением оси Ox (долгота), ψ – угол радиуса-вектора
с плоскостью Oxy (широта).
Якобиан преобразования (16) равен
cos φ cos ψ −r sin φ sin ψ −r cos φ sin ψ ∂(x, y, z)
= sin φ sin ψ r cos φ cos ψ −r sin φ sin ψ =
∂(r, φ, ψ) sin ψ
0
r cos φ
= r2 cos ψ.
Пример. Переходя к полярной системе координат, вычислить двойной интеграл
∫ ∫ √
x2 + y 2 dxdy,
I=
Ω
где Ω — четверть круга x2 + y 2 ≤ a2 , a > 0.
Решение. Сделав замену переменных
x = r cos φ, y = r sin φ,
получим, что область Ω взаимно однозначно отобразиться на область
Ω′ = {(r, φ) | r ∈ [0, a], φ ∈ [0, π/2]}. В результате будем иметь
∫ ∫ √
r2 dr
r2 cos2 φ + r2 sin2 φr drdφ =
I=
π
∫2
∫a
Ω′
0
Пример. Вычислить интеграл
∫ ∫
I=
Ω
x2 sin
y
( xy )
2
dxdy,
dφ =
0
πa3
.
6
54
если область Ω ограничена четырьмя параболами
x2 =
πy 2 2πy 2
, x =
, y = 2x, y 2 = 4x.
3
3
Решение. Для вычисления интеграла сделаем замену переменных
√
√
3
3
x = uv 2 , y = u2 v.
Такая замена переменных устанавливает взаимно однозначное соответствие между областью Ω и прямоугольником Ω′ , где
{
}
π
2π
Ω′ = (u, v) ∈ R2 | 2 ≤ u ≤ 4,
.
≤v≤
3
3
Найдем якобиан преобразования
2 2
D(x, y) 13 u− 3 v 3
=
1 1
D(u, u) 23 u− 3 v 3
2 13 − 13
3u v
1 23 − 23
3u v
1
=− .
3
Следовательно, имеем
∫ ∫
I=
2
2
3
u v
Ω′
2π
1
=
3
∫3
∫4
v dv
π
3
sin
( uv )
2
( uv )
4
u 3 v 3 sin
2
1
3
|J| dudv =
√
2π
2
du =
3
∫3
(cos v − cos (2v)) dv =
3
.
3
π
3
2
Пример. Переходя к полярным координатам, свести двойной интеграл к однократному
∫ ∫ ( )
y
f
I=
dxdy,
x
Ω
где Ω — множество, ограниченное петлей декартова листа x3 +y 3 = 3xy.
55
Решение. Перейдем в полярную систему координат
x = r cos φ, y = r sin φ,
в которой уравнение петли будет иметь следующий вид:
r=
3 cos φ sin φ
.
cos3 φ + sin3 φ
При этом угол φ меняется от 0 до π2 , т.к. в этом случае r ≥ 0, а значениям φ = 0 и φ = π2 соответствует точка (0, 0) в декартовой системе
координат. Таким образом, множество Ω является образом множества
Ω′ в полярной системе координат, где
{
}
π
3 cos φ sin φ
′
Ω = (r, φ) | 0 ≤ φ ≤ , 0 ≤ r ≤
.
2
cos3 φ + sin3 φ
Окончательно, получаем
(
∫ ∫
I=
f
Ω′
r sin φ
r cos φ
3 cos φ sin φ
cos3 φ+sin3 φ
π
)
∫2
r drdφ =
∫
f (tg φ) dφ
0
r dr =
0
π
=
9
2
∫2
f (tg φ)
cos2 φ sin2 φ
dφ.
(cos3 φ + sin3 φ)2
0
Пример. Вычислить
∫ ∫ ∫ √
I=
z x2 + y 2 dxdydz,
Ω
где область Ω ограничена цилиндром x2 + y 2 = 2x и плоскостями y = 0,
z = 0, z = 1.
Решение. Перейдем к цилиндрическим координатам, в которых цилиндр описывается следующим образом:
r = cos φ, 0 ≤ z ≤ 1.
56
Следовательно, прообразом области Ω в цилиндрических координатах
есть множество
{
}
π
Ω′ = (r, φ, z) | 0 ≤ r ≤ 2 cos φ, 0 ≤ φ ≤ , 0 ≤ z ≤ 1 .
2
Поэтому
π
∫ ∫ ∫
2
I=
2∫
cos φ
∫2
zr drdφdz =
Ω′
∫1
2
r dr
dφ
0
0
z dz =
0
π
∫2
4
=
3
8
cos2 φ dφ = .
9
0
Пример. Вычислить
∫ ∫ ∫
(x2 + y 2 ) dxdydz,
I=
Ω
где область Ω — верхняя половина шара x2 + y 2 + z 2 ≤ a2 .
Решение. Перейдем к сферическим координатам. Тогда прообразом
множества Ω является множество
Ω′ = {(r, φ, ψ) | 0 ≤ r ≤ a, 0 ≤ φ ≤ 2π, 0 ≤ ψ ≤
π
}.
2
Таким образом, получаем
∫ ∫ ∫
∫a
r4 cos2 ψ drdφdψ =
I=
Ω′
r4 dr
0
π
∫2
∫2π
cos3 ψdψ =
dφ
0
0
π
=
2π 5
a
5
∫2
(1 − sin2 φ) d(sin φ) =
0
57
(
) π
2 4π 5
2π 5
1
3
=
a sin φ − sin φ =
a .
5
3
15
0
Задачи.
76. Переходя к полярной системе ∫координат,
расставить пределы
∫
интегрирования в двойном интеграле
f (x, y)dxdy, если
Ω
а) Ω – круг x2 + y 2 ≤ a2 ,
б) Ω – круг x2 + y 2 ≤ ax, a > 0,
в) Ω – кольцо a2 ≤ x2 + y 2 ≤ b2 ,
г) Ω – треугольник 0 ≤ x ≤ 1, 0 ≤ y ≤ 1 − x,
д)∗ Ω – квадрат с вершинами O(0, 0), A(0, 1), B(1, 0), C(1, 1),
е)∗ Ω = {(x, y) | (x − 1)2 + y 2 ≤ 1, x2 + (y − 1)2 ≤ 1} ,
ж)∗ Ω = {(x, y) | (x − 1)2 + y 2 ≤ 1, 1 ≤ x ≤ y}.
77. Переходя к полярным координатам, вычислить следующие интегралы:
∫ ∫
cos(x2 + y 2 ) dxdy, где Ω = {(x, y) | x2 + y 2 ≤ 1},
а)
Ω
∫ ∫
ln(1 + x2 + y 2 ) dxdy, где Ω = {(x, y) | x2 + y 2 ≤ 1},
б)
Ω
∫ ∫
в)
x2
dxdy, где Ω = {(x, y) | x2 + y 2 ≤ x},
x2 + y 2
Ω
∫ ∫ √
г)
x2 + y 2 dxdy, где Ω = {(x, y) | x2 + y 2 ≤ y},
Ω
∫ ∫
y2
д)
Ω
√
1 − x2 dxdy, где Ω = {(x, y) | x2 + y 2 ≤ 1}.
58
78. Переходя к полярным координатам, свести двойные интегралы
к однократным:
∫ ∫
√
а)
f ( x2 + y 2 ) dxdy, где Ω = {x2 + y 2 ≤ x, x2 + y 2 ≤ y},
Ω
(
∫ ∫
б)
f
)
xy
x2 + y 2
Ω
dxdy, где Ω = {
∫ ∫
в)
f
Ω
Ω = {x + y ≤
2
2
(y)
x
√
|x| ≤ y ≤ 1},
dxdy, где
√
6x, (x2 + y 2 ) ≤ 9(x2 − y 2 )}.
79. Вместо x и y ввести новые переменные u и v и определить пределы интегрирования в следующих двойных интегралах:
∫b
а)
a
∫βx
dx f (x, y) dxdy, (0 < a < b, 0 < α < β), если
αx
u = x,
∫2
б)
0
v=
y
,
x
2−x
∫
dx
f (x, y) dxdy, если u = x + y,
v = y − x.
1−x
80. Ввести новые переменные u и v и вычислить следующие интегралы:
∫ ∫
а)
(x2 y 2 + y 2 ) dxdy,
D
где D = {1/x ≤ y ≤ 2/x, x ≤ y ≤ 3x},
∫ ∫
(x + y)2
б)
dxdy,
x
D
59
где D = {1 − x ≤ y ≤ 3 − x, x/2 ≤ y ≤ 2x},
∫ ∫
в)
xy(x + y) dxdy,
D
где D = {−1 ≤ x − y ≤ 1, 1/x ≤ y ≤ 2x},
∫ ∫
г)
x2 dxdy, где D = {x3 ≤ 2x3 , x ≤ 2y ≤ 6x},
D
∫ ∫
д)
xy(x + y) dxdy,
D
где D = {x − 1 ≤ y ≤ x + 1, −x − 1 ≤ y ≤ −x + 1}.
∫∫∫
81. В интеграле
f (x, y, z) dxdydz перейти к сферическим коорD
динатам и затем записать его в виде повторного, если:
а) D = {a2 ≤ x2 + y 2 + z 2 ≤ 4a2 , y ≥ 0},
б) D = {x2 + y 2 + z 2 ≤ R2 , y ≥ 0, x ≥ 0, x + y ≥ 0},
в) D = {x2 + y 2 + z 2 ≤ R2 , x2 + y 2 ≥ 3z 2 },
г) D = {x2 + y 2 + z 2 ≤ 2az, x2 + y 2 ≥ z 2 }.
∫∫∫
f (x, y, z) dxdydz перейти к цилиндрическим
82. В интеграле
D
координатам и затем записать его в виде повторного, если:
{x y z
}
а) D =
+ + ≤ 1, x ≥ 0, y ≥ 0, z ≥ 0 ,
2 3 4
б) D = {x2 + y 2 ≤ z, 0 ≤ z ≤ H},
в) D = {x2 + y 2 ≤ 2x, x2 + y 2 ≤ 1 − z, z ≥ 0},
√
г) D = {|z| ≤ 5 − 3x2 + 3y 2 , z 2 ≤ x2 + y 2 + 1}.
60
83. Вычислить
∫ ∫ ∫
а)
(x2 + y 2 ) dxdydz, D = {x2 + y 2 ≤ 2z, 0 ≤ z ≤ 2},
D
∫ ∫ ∫
(x2 − 4xy + y 2 ) dxdydz,
б)
D
D = {0 ≤ x ≤ 1, 0 ≤ y ≤ 1, 0 ≤ z ≤ 1, },
∫ ∫ ∫
в)
xyz dxdydz,
D
√
D = {x + y + z ≤ 3x, x2 + y 2 + z 2 ≤ y, z ≥ 0},
∫ ∫ ∫
г)
z dxdydz, D = {x2 + y 2 + z 2 ≤ 1, z ≥ 0},
2
∫ ∫ ∫
2
2
D
z 2 dxdydz, D = {x2 + y 2 + z 2 ≥ 1, x2 + y 2 + z 2 ≤ 2z}.
д)
D
84. Показать, что при переходе к обобщенным сферическим координатам
x = ar cos φ cos ψ, y = br sin φ cos ψ, z = cr sin ψ
якобиан отображения равен J = abcr2 cos ψ.
85. Пусть
{ 2
}
x
y2 z2
Ω=
+ 2 + 2 ≤1 .
a2
b
c
Вычислить интегралы:
)
∫ ∫ ∫ ( 2
x
y2 z2
а)
+ 2 + 2 dxdydz,
a2
b
c
Ω
∫ ∫ ∫
(x2 + y 2 ) dxdydz.
б)
Ω
61
86∗ . Вычислить интеграл
∫ ∫ ∫
xyz dxdydz,
D
где область Ω ограничена поверхностями
yz = a1 x, yz = a2 x, a2 > a1 > 0,
zx = b1 y, zx = b2 y, b2 > b1 > 0,
xy = c1 z, xy = c2 z, c2 > c1 > 0,
сделав замену переменных
u=
yz
zx
xy
, v=
, w=
.
x
y
z
∗
∫ ∫ 87
∫ . Подобрав подходящую замену координат, вычислить интеграл
x dxdydz, где
D
D = {a ≤ xyz ≤ b, cx ≤ yz ≤ dx, my ≤ x ≤ ny}.
2.9 Несобственные кратные интегралы.
Выше мы ввели понятие кратного интеграла для ограниченных областей и, кроме того, основные его свойства сформулировали для ограниченных функций. В приложениях же часто встречаются задачи, когда необходимо рассматривать неограниченные области или функции.
В этом параграфе мы введем понятие несобственного кратного интеграла. В отличие от одномерного случая мы не будем рассматривать
два вида несобственных интегралов — по неограниченным областям и
от неограниченных функций, а определим их сразу.
Определение. Последовательность измеримых по Жордану мноn
жеств {Ωm }∞
m=1 , Ωm ⊂ R , называют исчерпывающей множество Ω ⊂
n
R , если
1) Ω1 ⊂ Ω2 ⊂ . . . ⊂ Ωm ⊂ Ω для всех m ≥ 1,
62
2)
∞
∪
Ωm = Ω.
m=1
Далее будем рассматривать такие функции f (x), x ∈ Ω ⊂ Rn , которые интегрируемы на любом открытом измеримом по Жордану множестве ω таком, что ω ⊂ Ω.
Определение. Пусть для любой последовательности множеств {Ωm }∞
m=1 ,
исчерпывающих Ω, существует предел
∫
lim
f (x) dx,
m→∞
Ωm
не зависящий от выбора последовательности {Ωm }∞
m=1 . Тогда этот предел называют несобственным кратным интегралом от функции f (x) на
множестве Ω и обозначают
∫
f (x) dx.
Ω
При этом функцию f называют интегрируемой в несобственном смысле
на Ω. Таким образом,
∫
∫
f (x) dx = lim
f (x) dx.
(17)
Ω
m→∞
Ωm
Определение. Если предел (17) существует, то говорят, что несобственный интеграл сходится. В противном случае, говорят, что несобственный интеграл расходится.
Замечание. Если функция f (x), x ∈ Ω, интегрируема по Риману
на множестве Ω ⊂ Rn , то предел (17) существует для любой последовательности {Ωm }, исчерпывающей множество Ω и равен интегралу от
функции f (x) на множестве Ω, то есть
∫
∫
lim
f (x) dx = f (x) dx.
m→∞
Ωm
Ω
63
Таким образом, понятие несобственного интеграла является обобщением понятия интеграла Римана.
Определение. Несобственный интеграл
∫
f (x) dx
Ω
называется абсолютно сходящимся, если сходится интеграл
∫
|f (x)| dx.
Ω
Теорема. Несобственный интеграл
∫
f (x) dx сходится тогда и толь-
Ω
ко тогда, когда он сходится абсолютно.
Смысл этой теоремы в том, что понятия сходимости и абсолютной
сходимости совпадают, т.е. отсутствует понятие условной сходимости.
Напомним, что в одномерном случае, вообще говоря, из сходимости
несобственного интеграла не следует его абсолютная сходимость.
Перечислим основные свойства несобственных интегралов.
1. Пусть функции f (x) и g(x) интегрируемы в несобственном смысле
на множестве Ω. Тогда αf (x) ± βg(x) интегрируема в несобственном
смысле на Ω и
∫
∫
∫
(αf (x) ± βg(x)) dx = α f (x) dx ± β g(x) dx (линейность).
Ω
Ω
Ω
2. Пусть f (x) ≥ 0, x ∈ Ω. Тогда, если предел (17) существует хотя
бы для одной последовательности {Ωm }, то предел (17) существует для
любой последовательности, исчерпывающей область Ω и они все равны,
т.е. несобственный интеграл
∫
f (x) dx
Ω
сходится.
64
Замечание. Из этого свойства следует, что для исследования сходимости несобственного интеграла от неотрицательной функции достаточно выбрать одну исчерпывающую последовательность множеств и
исследовать существование предела (17).
3. Пусть функции f (x) и g(x) интегрируемы на одних и тех же подмножествах множества Ω. Пусть
∫ 0 ≤ f (x) ≤ g(x) для всех x ∈ Ω. Тогда из сходимости интеграла g(x) dx следует сходимость интеграла
Ω
∫
∫
f (x) dx; из расходимости интеграла f (x) dx следует расходимость
Ω
Ω
∫
интеграла g(x) dx (признак сравнения).
Ω
4. Пусть множества Ω1 ⊂ Rn , Ω2 ⊂ Rn и отображение
φ : Ω 1 → Ω2
удовлетворяют следующим условиям: а) Ω1 и Ω2 – открытые; б) существуют множества S1 и S2 меры нуль такие, что множества Ω1 \S1 и
Ω2 \S2 – открытые и φ : Ω1 \S1 → Ω2 \S2 – диффеоморфизм, т.е. взаимно
однозначное дифференцируемое отображение. Тогда для любой
∫ неотрицательной функции f (x), x ∈ Ω2 из сходимости интеграла f (x) dx
Ω2
∫
следует сходимость интеграла f (φ(t)) |φ′ (t)| dt и оба интеграла равны
Ω1
(теорема о замене переменных в несобственном интеграле).
Пример. Исследовать сходимость интеграла
∫
2
2
I = e−(x +y ) dxdy
R2
и найти его значение.
Решение. Подынтегральная функция положительна. Рассмотрим последовательность множеств Ωk = {x2 + y 2 ≤ k 2 }, k = 1, 2, . . ., исчерпывающую все пространство R2 :
∫
2
2
Ik = e−(x +y ) dxdy.
Ωk
65
Переходя к полярным координатам, получим
∫2π
Ik =
∫k
dφ
0
e
−r2
0
∫k
rdr = 2π
e−r rdr =
2
0
= π(1 − e−k ) → π при k → ∞.
2
Так как подынтегральная функция неотрицательна, то для любой исчерпывающей пространство R2 последовательности предел числовой
последовательности Ik будет один и тот же. Поэтому интеграл I сходится и равен π.
Пример. Исследовать на сходимость интеграл
∫ ∫ ∫
dxdydz
√
I=
.
2
( x + y 2 + z 2 )α
x2 +y 2 +z 2 ≥1
Решение. Так как подынтегральная функция положительна, то можно выбрать одну исчерпывающую множество интегрирования последовательность. В этом примере удобно рассмотреть последовательность
колец Ωk , где
Ωk = {1 ≤ x2 + y 2 + z 2 ≤ k 2 }.
Переходя к сферическим координатам в интеграле
∫ ∫ ∫
dxdydz
√
Ik =
,
2
( x + y 2 + z 2 )α
Ωk
после преобразований получим
∫k
Ik = 4π
dr
rα−2
1
{
=
ln k (при α = 3,
)
.
1
1
−
1
при α ̸= 3
α−3
3−α k
Отсюда видно, что предел последовательности {Ik } при k → ∞ существует тогда и только тогда, когда α > 3. Таким образом, интеграл I
сходится при α > 3, и расходится при α ≤ 3.
66
Пример. Исследовать сходимость интеграла
∫ ∫
x−y
I=
dxdy,
x2 + y 2
Ω
где Ω = {x + y ≥ 1, x ≥ 0, y ≥ 0}.
Решение. Как было указано выше, несобственный кратный интеграл
сходится тогда и только тогда, когда он сходится абсолютно. Кроме
того, множество Ω симметрично относительно прямой y = x, поэтому
сходимость исходного интеграла эквивалентна сходимости интеграла
∫ ∫
x−y
I1 = 2
dxdy,
x2 + y 2
Ω1
где Ω1 = {x + y ≥ 1, x ≥ 0, 0 ≤ y ≤ x}.
Так как подынтегральная функция непрерывна, то в свою очередь
сходимость интеграла I1 эквивалентна сходимости интеграла
∫ ∫
x−y
I2 =
dxdy,
x2 + y 2
Ω2
где Ω2 = {x2 + y 2 ≥ 1, x ≥ 0, 0 ≤ y ≤ x}. Переходя в последнем
интеграле к полярным координатам, получаем
∫∞
∫π/4
(cos φ − sin φ) dφ dr,
I2 =
0
1
который расходится. Следовательно, интеграл I – тоже расходится.
Задачи.
88. Переходя к полярным координатам, вычислить интегралы:
∫ ∫
2
2
а)
e−(x +y ) cos (x2 + y 2 ) dxdy,
R2
67
∫ ∫
б)
e−(x
2 +y 2 )
sin (x2 + y 2 ) dxdy.
R2
89. Исследовать сходимость и вычислить интегралы (параметры положительны):
∫ ∫
∫ ∫
dxdy
dxdy
а)
, б)
,
p
q
x y
(x + y)p
xy ≥ 1
x+y ≥1
x≥1
0≤x≤1
∫ ∫
∫∫
dxdy
в)
,
г)
e−(x+y) dxdy.
(x2 + y 2 )p
x2 +y 2 ≥1
0≤x≤y
90. Показать, что
∫ ∫
sin (x2 + y 2 ) dxdy = π,
lim
n→∞
|x| ≤ n
|y| ≤ n
∫ ∫
тогда как
lim
n→∞
x2 +y 2 ≤2πn
91. Показать, что интеграл
∫∫
sin (x2 + y 2 ) dxdy = 0.
x2 − y 2
dxdy
(x2 + y 2 )2
x≥1, y≥1
расходится, но повторные интегралы
∫∞
∫∞
dx
1
сходятся.
1
x2 − y 2
dy,
(x2 + y 2 )2
∫∞
∫∞
dy
1
1
x2 − y 2
dx
(x2 + y 2 )2
68
92∗ . Исследовать на сходимость следующие интегралы в зависимости от параметра:
∫ ∫
∫ ∫
dxdy
dxdy
а)
,
б)
.
(1 + x2 + xy + y 2 )α
(1 + x4 + y 4 )α
R2
R2
93. Вычислить интегралы:
∫1 ∫1 ∫1
а)
0
0
dxdydz
,
xp y q z r
∫∞ ∫∞ ∫∞
б)
0
∫ ∫ ∫
в)
x2 +y 2 +z 2 ≤1
e−(x
2 +y 2 +z 2 )
dxdydz,
−∞ −∞ −∞
dxdydz
.
(1 − x2 − y 2 − z 2 )p
94. Исследовать сходимость следующих интегралов:
∫ ∫
y
а)
dxdy, где Ω = {x ≥ 1, −1 ≤ y ≤ 1},
x
Ω
∫ ∫
б)
y
dxdy, где Ω = {x ≥ 1, −1 ≤ xy ≤ 1},
x
Ω
∫ ∫
в)
e−x
√ dxdy, где Ω = {x ≥ 0, −1 ≤ y ≤ 1}.
y
Ω
95. Исследовать сходимость следующих интегралов:
∫ ∫ ∫
1−x−y−z
а)
dxdydz,
xyz
Ω
где Ω = {x + y + z ≤ 1, x ≥ 0, y ≥ 0, z ≥ 0},
∫ ∫ ∫
1−x−y−z
б)
dxdydz,
√
xyz
Ω
69
где Ω = {x + y + z ≤ 1, x ≥ 0, y ≥ 0, z ≥ 0},
∫ ∫ ∫
x−y
√ dxdydz, где Ω = {x ≥ 1, y ≥ 1, z ≥ 1}.
в)
z3
Ω
70
3. Геометрические и физические приложения кратных
интегралов.
3.1 Геометрические приложения кратных интегралов.
1. Пусть Ω ⊂ Rn – измеримое по Жордану множество. Тогда мера
µΩ множества Ω вычисляется по формуле
∫
µΩ = dx.
Ω
В пространстве R2 мера – это площадь S, в R3 – объем V :
∫ ∫
S=
dxdy,
Ω
∫ ∫ ∫
V =
dxdydz.
Ω
Определение. Пусть функцию двух переменных z = f (x, y) – непрерывна и неотрицательна, (x, y) ∈ Ω ⊂ R2 . Тело V ⊂ R3 , ограниченное
сверху графиком функции z = f (x, y), снизу – областью Ω, называется
цилиндрическим.
Цилиндрическое тело является измеримым по Жордану множеством
в пространстве R3 и его объем вычисляется по формуле
∫ ∫
f (x, y) dxdy.
V =
Ω
Таким образом, вычисление площадей плоских фигур и объемов цилиндрических тел сводится к вычислению двойных интегралов.
Пример. Вычислить площадь фигуры, ограниченной линиями
x = 4y − y 2 , x + y = 6.
71
Решение. Найдем координаты точек пересечения заданных линий.
Для этого нужно решить систему уравнений x = 4y − y 2 , x + y = 6. В
результате получим две точки пересечения A(4, 2) и B(3, 3). Следовательно, площадь фигуры равна
∫ ∫
S=
∫3
dxdy =
Ω
2
4y−y
∫
dy
2
∫3
dx =
6−y
1
(−y 2 + 5y − 6) dy = .
6
2
Пример. Найти площадь S области, ограниченной кривыми
x2 + y 2 = 1, x2 + y 2 = 4, y = 2x, y = 3x.
72
x2
D
x1
Решение. Рассмотрим область D в полярной системе координат
x = r cos φ, y = r sin φ.
Тогда части границы области D, соответствующие окружностям в полярной системе координат будут описываться уравнениями r = 1 и
r = 2, а прямым – φ = arctg 2 и φ = arctg 3. Таким образом, в по-
73
лярной системе координат D = {1 ≤ r ≤ 2, arctg 2 ≤ φ ≤ arctg 3}.
Следовательно, площадь S равна
arctg
∫ 3
S=
∫2
dφ
arctg 2
1
rdr = (arctg 3 − arctg 2).
2
1
Отметим, что задачу можно решить и другими способами. Например,
область D разбить на цилиндрические подобласти. Но такой способ трудоемкий, т.к. пришлось бы искать точки пересечения линий, и, в этом
случае, необходимо решать систему нелинейных уравнений.
Пример. Вычислить площадь плоской области Ω, ограниченной
кривой
(x2 + y 2 )2 = 2x3 .
Решение. Введем полярные координаты на плоскости. Подставим
x = r cos φ, y = sin φ в уравнение кривой. В результате получим, что в
полярных координатах кривая имеет вид
r4 = 2r3 cos3 φ или r = 2 cos3 φ.
Заметим, что область Ω симметрична относительно оси Ox. Поэтому
достаточно вычислить площадь половины области Ω, лежащей в первой
четверти и результат удвоить.
Для верхней половины области Ω угол φ меняется от 0 до π/2, а
радиус r – от 0 до 2 cos3 φ. Таким образом, получаем
π
∫ ∫
dxdy = 2
S=
π
0
0
π
∫2
cos6 φdφ =
2 cos3 φ
r2 dφ =
0
π
∫2
0
π
∫2
r dr =
dφ
0
Ω
=4
3
2 cos
∫ φ
∫2
∫2
(1 + cos 2φ)3 dφ =
0
+3 cos2 2φ + cos3 2φ)dφ =
(1 + 3 cos 2φ+
0
5π
.
16
74
Пример. Вычислить объем тела, вырезанного цилиндром x2 + y 2 =
Rx из сферы x2 + y 2 + z 2 = R2 .
Решение. Тело симметрично относительно плоскости Oxy. Поэтому
можно найти объем верхней√
половины и удвоить его. Сверху тело ограничено графиком функции R2 − x2 − y 2 , поэтому
V =2
∫ ∫ √
R2 − x2 − y 2 dxdy.
Ω
Здесь область интегрирования Ω — это проекция цилиндра на плоскость
Oxy, т.е. круг x2 + y 2 ≤ Rx. В полярных координатах граница этого
круга имеет вид r = R cos φ. Таким образом, получаем
π
∫2
V =2
R∫
cos φ
√
dφ
− π2
R2 − r2 r dr =
0
π
=−
2R3
3
∫2
(sin3 φ − 1) dφ =
− π2
2(3π − 4) 3
R .
9
Пример. Найти объем тела D, ограниченного поверхностями
x2 + y 2 = x, x2 + y 2 = 2x, z = x2 + y 2 , z = 0.
Решение. Воспользуемся цилиндрическими координатами, в которых поверхности, ограничивающие область, имеют вид:
r = cos φ, r = 2 cos φ, h = r2 , h = 0.
Поэтому D в полярных координатах имеет вид
{
[ π π]
}
D= φ∈ − ,
, r ∈ [cos φ, 2 cos φ], h ∈ [0, r2 ] .
2 2
75
Следовательно, объем V области D равен
π
∫ ∫ ∫
V =
∫2
dxdydz =
dφ
− π2
D
π
∫2
=
rdr
cos φ
15
r3 dr =
4
dφ
∫r2
dh =
0
π
2∫
cos φ
− π2
2∫
cos φ
cos φ
∫2
cos4 φ =
− π2
45
π.
32
Пример. Найти объем тела D, ограниченного следующими поверхностями:
z = x2 + y 2 , z = 2 (x2 + y 2 ), y = x, y = x2 .
Решение. Проекция области D на плоскость Oxy представляет собой
множество {(x, y) ∈ R2 | 0 ≤ x ≤ 1, x2 ≤ y ≤ x}, поэтому
D = {(x, y, z) ∈ R3 | 0 ≤ x ≤ 1, x2 ≤ y ≤ x, x2 + y 2 ≤ z ≤ 2(x2 + y 2 )}.
Объем тела равен
∫ ∫ ∫
∫1
dxdydz =
V =
∫1
=
0
∫1 (
∫x
2
dx
2
(x + y )dy =
x2
dx
0
D
∫x
2(x∫2 +y 2 )
dy
x2
dz =
x2 +y 2
)
1 6
3
4 3
4
x − x − x dx = .
3
3
35
0
Пример. Найти объем тела D, ограниченного поверхностью
(x2 + y 2 + z 2 )2 = (x2 + y 2 − z 2 ).
Решение. Если заменить x, y, z на −x, −y, −z, то уравнение поверхности не изменится. Следовательно, область симметрична относительно
76
всех координатных плоскостей, поэтому можно найти объем тела, находящегося в первом октанте, и умножить его на 8. Далее, перейдем
к сферическим координатам, тогда уравнение поверхности будет иметь
следующий вид:
r2 = cos 2ψ.
Таким образом, имеем
π
∫ ∫ ∫
V =
dxdydz = 8
∫2
cos ψ dψ
0
D
∫4
r2 dr =
dφ
0
π
4π
=
3
√
cos
∫ 2ψ
π
∫4
0
π
4π
cos ψ (cos 2ψ) dψ =
3
3
2
0
∫4
3
(1 − 2 sin2 ψ) 2 d(sin ψ).
0
Вычисляя последний интеграл с помощью метода замены переменной,
получим
π2
V = √ .
4 2
Пример. Пользуясь подходящей заменой переменных, вычислить
объем тела, ограниченного поверхностями
z = x2 + y 2 , z = 2 (x2 + y 2 ), xy = 1, xy = 2, x = 2y,
2x = y, (x > 0, y > 0).
Решение. Для нахождения объема необходимо вычислить интеграл
∫ ∫ ∫
I=
dxdydz
Ω
где Ω — область, которую занимает тело. Произведем в интеграле замену переменных, положив
u = xy, v =
y
, z = z.
x
77
При указанной замене
x = u 2 v − 2 , y = u 2 v 2 , z = z.
1
1
1
Вычислим якобиан
D(x, y, z) =
D(u, v, z) 1
2
1
2
u− 2 v − 2
1 1
u− 2 v 2
1
1
1
1
3 1
− 21 u 2 v − 2 .
=
1 12 − 12
2v
2u v
Кроме того, прообраз множества Ω есть множество
{
1
′
Ω = (u, v, z) | 1 ≤ u ≤ 2, ≤ v ≤ 2,
2
}
−1
−1
u(v + v ) ≤ z ≤ 2u(v + v ) .
Таким образом, получаем
∫ ∫ ∫
I=
Ω′
1
D(x, y, z)
dudvdz =
D(u, v, z)
2
du
1
1
=
2
∫2
1
∫2
∫2
dv
v
1
2
−1 )
2u(v+v
∫
dz =
u(v+v −1 )
)
∫2 (
1
3
u du
1 + 2 dv = .
v
2
1
2
Задачи.
96. Вычислить площади фигур, ограниченные заданными линиями:
а) x = y 2 − 2y, x + y = 0,
в) y 2 = 4x − x2 , y 2 = 2x,
б) y = 2 − x, y 2 = 4x + 4,
г) 3y 2 = 25x, 5x2 = 9y,
д) y 2 + 2y − 3x + 1 = 0, 3x − 3y − 7 = 0,
е) y = 4x − x2 , y = 2x2 − 5x.
78
97. Переходя в полярную систему координат, найти площади плоской фигуры, ограниченной кривой:
а) (x2 + y 2 )2 = 2(x2 − y 2 ),
б) (x2 + y 2 )3 = x4 + y 4 ,
в) (x2 + y 2 )3 = 4x2 y 2 .
98∗ . Найти площадь плоской фигуры, ограниченной кривыми:
а) x2 + y 2 = 2ax, x2 + y 2 = 2bx, y = x, y = 0, b > a > 0,
√
б) (x2 + y 2 )2 = 2a2 (x2 − y 2 ), x2 + y 2 = a2 , ( x2 + y 2 ≥ a > 0).
99. Вычислить объемы тел, ограниченные следующими поверхностями:
x y z
а) + + = 1, x = 0, y = 0, z = 0,
a
b
c
б) z = x2 + y 2 , x + y = 1, x = 0, y = 0, z = 0,
в) z = x2 + y 2 , y = x2 , z = 0, y = 1.
100. Используя цилиндрические координаты, найти объем тела, ограниченного заданными поверхностями:
а) x2 + y 2 + z 2 = 1, x2 + y 2 ≥ |x|,
б)∗ x2 + y 2 − az = 0, (x2 + y 2 )2 = a2 (x2 − y 2 ), z = 0,
в)∗ z = x2 + y 2 , x4 + y 4 = x2 + y 2 , z = 0,
г) x2 + y 2 = x, z 2 = x2 + y 2 , z ≥ 0.
101. Использую сферические координаты, найти объем тела, ограниченного заданными поверхностями
√
а) x2 + y 2 + z 2 = 4, z 2 = x2 + y 2 , 0 ≤ z ≤ x2 + y 2 ,
б) (x2 + y 2 + z 2 )2 = xyz,
в) (x2 + y 2 + z 2 )2 = z(x2 + y 2 ),
79
в) (x2 + y 2 + z 2 )2 = x.
102. Пользуюсь подходящей заменой переменных, вычислить объемы тел, ограниченных следующими поверхностями:
а) xy = 1, xy = 3, x = 2y, x = 3y, z = x2 + y 2 , z = 0,
б) y 2 = x, y 2 = x, x = y, x = 2y, z =
1
, z = 0,
xy
в) y 2 = 2x, y 2 = 3x, x2 = y, x2 = 2y, z = xy, z = 0,
г) xy = 1, xy = 4, y 2 = x, y 2 = 3x, z 2 = xy, z = 0,
д) x + y + z = 1, x + y + z = 2, x + y = z, x + y = 2z,
y = x, y = 3x.
3.2 Приложения кратных интегралов к геометрии масс.
Одной из физических характеристик материальной области (плоской или пространственной) является плотность – неотрицательная функция ρ, заданная на области Ω ⊂ Rn (n = 2 или n = 3) и интегрируемая
на Ω.
Определение. Массой материальной фигуры Ω с плотностью ρ(x, y)
(при n = 2) или ρ(x, y, z) (при n = 3) называют величину
∫ ∫
ρ(x, y) dxdy при n = 2,
(18)
M=
Ω
∫ ∫ ∫
M=
ρ(x, y, z) dxdydz
при n = 3.
(19)
Ω
Определение. Центром масс тела Ω с плотностью ρ(x, y, z) называют точку C ∈ R3 с координатами
∫ ∫ ∫
1
xC =
xρ(x, y, z) dxdydz,
M
Ω
80
1
yC =
M
∫ ∫ ∫
yρ(x, y, z) dxdydz,
Ω
zC =
1
M
∫ ∫ ∫
zρ(x, y, z) dxdydz,
Ω
где M – масса тела, которую можно вычислить по формуле (19).
Аналогично определяется центр масс плоской фигуры.
Определение. Величины Myz = M xC , Mzx = M yC , Mxy = M zC
называют статистическими моментами тела Ω относительно координатных плоскостей Oyz, Ozx, Oxy.
Определение. Моментом инерции тела Ω относительно оси l называют величину
∫ ∫ ∫
d2l ρ dxdydz,
Il =
(20)
Ω
где dl = dl (x, y, z) – расстояние от точки (x, y, z) тела до оси l, ρ =
ρ(x, y, z) – плотность тела. В частности, если в качестве l взять координатную ось Ox, то формула (20) имеет вид
∫ ∫ ∫
Ixx =
(y 2 + z 2 )ρ(x, y, z) dxdydz.
Ω
Аналогично определяются моменты инерции Iyy и Izz относительной
координатных осей Oy и Oz соответственно.
′ , I ′ , I ′ относительно коорОпределение. Моменты инерции Ixy
zx
yz
динатных плоскостей Oxy, Oyz, Ozx определяют по формулам
∫ ∫ ∫
′
Ixy =
z 2 ρ(x, y, z) dxdydz,
Ω
′
Iyz
∫ ∫ ∫
x2 ρ(x, y, z) dxdydz,
=
Ω
81
′
Izx
∫ ∫ ∫
y 2 ρ(x, y, z) dxdydz.
=
Ω
Определение. Момент инерции относительно точки M0 с координатами (x0 , y0 , z0 ) определяют по формуле
∫ ∫ ∫
d2M0 ρ dxdydz,
IM0 =
Ω
√
где dM0 = (x − x0 )2 + (y − y0 )2 + (z − z0 )2 – расстояние от точки (x, y, z)
тела до точки M0 . В частности, момент инерции IO относительно начала координат
∫ ∫ ∫
IO =
(x2 + y 2 + z 2 )ρ(x, y, z) dxdydz.
Ω
Пример. Найти массу тела, занимающего единичный куб Ω = {0 ≤
x ≤ 1, 0 ≤ y ≤ 1, 0 ≤ z ≤ 1} в R3 , если плотность тела в точке M (x, y, z)
дается формулой
ρ(x, y, z) = x + y + z.
Решение. Масса m тела Ω с плотностью ρ(x, y, z) вычисляется по
формуле
∫ ∫ ∫
m=
∫1
ρ(x, y, z) dxdydz =
∫1
dx
0
Ω
∫1
=
0
∫1
dy
0
(x + y + z) dz =
0
)
∫1 (
∫1
3
1
dy = (x + 1) dx = .
dx
x+y+
2
2
0
0
Пример. Найти координаты центра масс однородной пластинки,
ограниченной параболой ay = x2 и прямой y = 2, a > 0.
82
Решение. В силу симметрии пластинки относительно оси Oy (это
следует из того, что при замене x на −x уравнения границы пластинки
не меняются) имеем, что xC = 0. Найдем массу пластинки
∫ ∫
∫2
m=
dxdy = 2
√
∫2
∫ ay
√
8 2a
√
y dy =
dy
dx = 2 a
.
2
0
Ω
√
0
0
Отсюда получаем, что
1
yC =
m
∫ ∫
3
y dxdy = √
4 2a
∫2
√
∫ ay
∫2
3
3
6
y dy
dx = √
y 2 dy = .
5
4 2
0
Ω
0
0
Пример. Найти момент инерции квадрата со стороной a, поверхностная плотность которого равна y, относительно одной из вершин.
Решение. Поместим начало координат в одну из вершин квадрата,
а координатные оси направим по двум взаимно перпендикулярным его
сторонам. Найдем момент инерции относительно начала координат:
∫ ∫
∫a
2
I0 =
2
(x + y )y dxdy =
∫a
0
Ω
∫a
=
(
(x2 y + y 3 ) dy =
dx
0
)
4
a2 x 2 a
+
2
4
dx =
5a5
.
12
0
Пример. Найти статический момент однородного полукруга радиуса R относительно диаметра.
Решение. Поместим начало координат в середину диаметра полукруга , а ось Oy направим по диаметру. Тогда статический момент
полукруга относительно диаметра будет равен статическому моменту
полукруга относительно оси Oy. Вводя полярные координаты, получим
π
∫ ∫
My =
∫2
x dxdy =
Ω
∫R
r2 cos φ dr =
dφ
− π2
0
83
π
=
R3
3
∫2
cos φ dφ =
− π2
2R3
.
3
Пример. Найти координаты центра тяжести однородного (с единичной плотностью) тела, ограниченного поверхностями
x2 + y 2 = z, x + y + z = 0.
Решение. Из второго уравнения, описывающего границу тела Ω следует, что z = −x − y на этой границе. Поэтому проекция Ω на плоскость
Oxy есть множество D = {x2 + y 2 ≤ −x − y}, т.е. круг
(
)
(
)
1 2
1 2 1
x+
+ y+
≤ .
2
2
2
Кроме того, в силу симметрии тела относительно плоскости {x = y}
имеем xC = yC . Рассмотрим следующую замену переменных:
1
1
x = r cos φ − , y = r sin φ − .
2
2
Якобиан такой замены будет такой же, как и полярной, т.е. равен r.
Итак, масса тела
∫ ∫ ∫
dxdydz =
m=
√1
∫2π
dφ
0
D
1
√
=
dφ
0
rdr
0
∫ 2(
∫2π
1−r(cos
∫ φ+sin φ)
∫2
r2 −r(cos φ+sin
)
1
π
3
r − r dr = .
2
8
0
Координаты центра масс
1
x C = yC =
m
∫ ∫ ∫
x dxdydz =
D
dz =
φ)+ 21
84
=
1
m
∫2π
√1
)(
)
∫2 (
1
1
1
− r2 dr = − ,
dφ r r cos φ −
2
2
2
0
0
1
zC =
m
∫ ∫ ∫
z dxdydz =
D
1
=
m
√1
=
4
π
dr
dφ
rz dz =
r2 −r(cos φ+sin
0
0
∫2π
1−r(cos
∫ φ+sin φ)
∫2
∫2π
φ)+ 21
√1
)
∫2 (
1
2
dφ r 1 − 2r(cos φ + sin φ) + r +
·
2
0
0
(
)
1
5
2
· 1−r −
dr = .
2
6
)
(
Таким образом, получили, что координаты центра масс − 12 , − 21 , 65 .
Пример. Найти массу и момент инерции однородного (с единичной
плотностью) тела D, ограниченного поверхностями
x2 + y 2 = 2z, x2 + y 2 = z 2
относительно прямой l:
x = 0, z = 4.
Решение. Масса тела равна
∫ ∫ ∫
m=
∫2π
dxdydz =
dφ
0
D
∫2
∫r
dr
0
rdz =
r2
2
)
∫2 (
r3
4π
2
= 2π
r −
dr =
.
2
3
0
85
Момент инерции тела относительно прямой l равен
∫ ∫ ∫
Il =
d2l dxdydz,
D
где dl – расстояние от точки (x, y, z) тела D до прямой l. Квадрат этого
расстояния находим по формуле d2l = x2 + (z − 4)2 , следовательно
∫ ∫ ∫
( 2
)
Il =
x + (z − 4)2 dxdydz =
D
∫2π
=
∫2
dφ
0
∫2π
=
∫2
dφ
0
∫r
dr
0
(
)2
r cos2 φ + (z − 4)2 dz =
r2
2

 4
2
r cos φ −
r5
2
(
cos2 φ + r
(r −
3
4)3
−r
r2
2
−4
3
)3 

 dr =
0
∫2π (
=
1 5
r6
(r − 4)5
r cos2 φ −
cos2 φ +
+
5
12
15
0
+
4(r −
12
(
4)4
−
r2
2
)4 )
2
40π
.
dφ =
24
3
−4
0
Задачи.
103. Найти координаты центра масс однородной (с единичной плотностью) плоской фигуры:
{ 2
}
y
а) Ω =
≤ x ≤ 2a − y, a > 0 ,
a2
}
{
a2 y 2
≤ x ≤ 2a, a > 0 ,
б) Ω = y ≤ ,
x 8a
86
{√
}
√
√
x + y ≤ a, x ≥ 0, y ≥ 0 ,
{
}
г) Ω = x2 + y 2 ≤ 1, x ≥ 0, y ≥ 0 .
в)∗ Ω =
104. Найти момент инерции относительно координатных осей и относительно начала координат однородной (с единичной плотностью) плоской фигуры Ω, если
}
{
( π)
, a, α − const ,
а) Ω = x2 + y 2 ≤ a2 , |y| ≤ xtg α, α ∈ 0,
2
{
}
б) Ω = (x − a)2 + (y − a)2 ≥ a2 , a ≥ x ≥ 0, a ≥ y ≥ 0 ,
{x y
}
x y
в) Ω =
+ ≤ 1,
+ ≥ 1, y ≥ 1, a > b > 0, c > 0 ,
a
c
b
c
{
}
∗
2
г) Ω = xy = a , xy = 2a, x = 2y, y = 2x, x > 0, y > 0 .
105. Найти моменты инерции относительно координатных осей однородных (с единичной плотностью) тел:
а) D = {0 ≤ x ≤ a, 0 ≤ y ≤ b, 0 ≤ z ≤ c},
б) D = {x2 + y 2 ≤ R2 , 0 ≤ z ≤ h},
{
}
√
в) D = 0 ≤ z ≤ h(R − x2 + y 2 ) ,
г) D = {x2 + y 2 + z 2 ≤ R2 , z ≥ 0}.
106. Найти координаты центра масс тела Ω с плотностью ρ(x, y, z),
если
а) Ω = {0 ≤ x ≤ a, 0 ≤ y ≤ a, 0 ≤ z ≤ a}, ρ(x, y, z) = (x + y + z)2 ,
б) Ω = {x2 + y 2 + z 2 ≤ R2 , x ≥ 0}, ρ(x, y, z) = √
1
x2 + y 2
,
в) Ω = {R2 ≤ x2 + y 2 + z 2 ≤ 4R2 , y ≥ 0}, ρ = x2 + y 2 + z 2 ,
{√
}
г) Ω =
x2 + y 2 ≤ z ≤ h , ρ(x, y, z) = z 2 ,
87
д)∗ Ω = {x2 + y 2 − z 2 = 0, 0 ≤ z ≤ h}, ρ(x, y, z) = z.
3.3 Приложения кратных интегралов к физике.
Пусть в пространстве R2 задана область Ω и функция µ(ξ, η), (ξ, η) ∈
Ω.
Определение. Логарифмическим потенциалом в точке M (x, y) называется интеграл
∫ ∫
1
u(x, y) =
µ(ξ, η) ln r dξdη,
2π
Ω
√
где r = (ξ − x)2 + (η − y)2 – расстояние от точки (x, y) до точки (ξ, η).
Функцию µ называют плотностью (например, плотность масс, плотность электрического заряда).
Определение. Пусть в пространстве R3 задана область G и функция ρ(ξ, η, ζ) ∈ G. Ньютоновским потенциалом в точке M с координатами (x, y, z) называют интеграл
∫ ∫ ∫
ρ(ξ, η, ζ)
u(x, y, z) = k
dξdηdζ,
r
G
√
где r = (ξ − x)2 + (η − y)2 + (ζ − z)2 – расстояние от точки (x, y, z) до
точки (ξ, η, ζ), k = const. Функцию ρ называют плотностью (например,
плотность масс, плотность электрического заряда). Если ρ – плотность
масс, то ньютонов потенциал – это потенциал гравитационного поля
материального тела G.
Определение. Если ньютонов потенциал u(x, y, z) определен в области G, то говорят, что в G задано поле с ньютоновым потенциалом
u(x, y, z). Напряженностью этого поля в точке M (x, y, z) называют вектор
∫ ∫ ∫ −−→
−
→
MN
E (x, y, z) = k
−−→ 3 ρ(ξ, η, ζ)dξdηdζ,
|M N |
G
88
−−→
где N = N (ξ, η, ζ), M N = (ξ − x, η − y, ζ − z),
√
−−→
|M N | = (ξ − x)2 + (η − y)2 + (ζ − z)2 .
Определение. Пусть потенциал определен плотностью ρ(ξ, η, ζ), заданной в области G, и пусть в поле с этим потенциалом находится область G1 с заданной в ней плотностью ρ1 (ξ1 , η1 , ζ1 ). Силой, действующей
на G1 , называют вектор
∫ ∫ ∫
−
→
→
−
F =
E (ξ1 , η1 , ζ1 )ρ1 (ξ1 , η1 , ζ1 )dξ1 dη1 dζ1 .
G1
Определение. Пусть M (x0 , y0 , z0 ) – материальная точка с массой
−
→
m находится в гравитационном поле напряженностью E (x, y, z). Тогда
силой, действующей на материальную точку M , называется вектор
−
→
−
→
F = m E (x0 , y0 , z0 ).
Аналогично определятся сила, действующая на точечный заряд, находящийся в электростатическом поле.
Пример. Гравитационное поле создано полым шаром G с внутренним и внешним радиусами R1 и R2 , имеющим плотность
ρ(x, y, z) =
ρ0
,
r
ρ0 = const,
где r – расстояние от центра шара до точки (x, y, z). Найти силу, действующую на материальную точку M массы m, удаленную на расстояние R от центра шара, R > R2 .
Решение. Поместим начало системы координат в центр шара, ось Oz
направим через точку M . Тогда точка M имеет координаты (0, 0, R).
Область, занимаемая шаром, описывается уравнениями
√
R1 ≤ ξ 2 + η 2 + ζ 2 ≤ R2 .
−−→
Пусть N (ξ, η, ζ) – произвольная точка шара, тогда M N = (ξ, η, ζ − R),
√
−−→
|M N | = ξ 2 + η 2 + ζ 2 + R2 − 2ζR = r2 + R2 − 2ζR, где r = ξ 2 + η 2 + ζ 2 .
89
−
→
Следовательно, составляющие вектора сил F = (Fx , Fy , Fz ) вычисляются по формулам
∫ ∫ ∫
ξ dξdηdζ
Fx = mρ0
,
2
r(r + R2 − 2ζR)3/2
G
∫ ∫ ∫
Fy = mρ0
r(r2
G
∫ ∫ ∫
Fz = mρ0
G
η dξdηdζ
,
+ R2 − 2ζR)3/2
(ζ − R) dξdηdζ
.
r(r2 + R2 − 2ζR)3/2
Очевидно, что первый и второй интегралы равны нулю, т.е. Fx = Fy = 0.
В третьем интеграле перейдем к сферическим координатам
ξ = r cos φ cos ψ, η = r sin φ cos ψ, ζ = r sin ψ.
В результате получим
π
∫2π ∫2 ∫R2
Fz = mρ0
0 − π2 R1
(r sin ψ − R)r cos ψ dφdψdr
=
(r2 + R2 − 2rR sin ψ)3/2
π
∫R2
∫2
= 2πmρ0
dr
− π2
R1
r(r sin ψ − R) cos ψ dψ
.
(r2 + R2 − 2rR sin ψ)3/2
Во внутреннем интеграле сделаем замену переменных
t = r2 − R2 − 2rR sin ψ,
тогда
πmρ0
Fz =
2R2
∫R2
2
(R+r)
∫
dr
R1
(R−r)2
(
)
(r2 − R2 )t−3/2 − t−1/2 dt =
90
)
2πmρ0 ( 2
R2 − R12 .
2
R
Пример. Найти ньютонов потенциал в точке P (0, 0, z), z > 0 неоднородного шара x2 + y 2 + z 2 = 1, если его плотность ρ(ξ, η, ζ) = kξ 2 ,
k − const, т.е. плотность пропорциональна квадрату расстояния точки
шара до плоскости Oxy.
Решение. Ньютонов потенциал вычисляется по формуле
∫ ∫ ∫
ζ dξdηdζ
√
u(0, 0, z) = k
.
ξ 2 + η 2 + (ζ − z)2
=−
G
Перейдем к сферическим координатам, тогда
π
∫2π
u(0, 0, z) = k
∫2
∫R
− π2
0
0
r4 sin2 ψ cos ψ dr
r2 cos2 ψ + (r sin ψ − z)2
=
π
∫2
∫R
r4 dr
= 2πk
√
dψ
dφ
0
− π2
√
sin2 ψ cos ψ dψ
r2 + z 2 − 2rz sin ψ
.
Сделаем замену переменных t = r2 + z 2 − 2rz sin ψ, тогда получим
u(0, 0, z) =
πk
4
2
(z+r)
∫
∫R
dr
0
r(r2 + z 2 − t2 )2
4kπ
dt =
15z 3
z 3 t1/2
(
)
2 7
R + z 2 R5 .
7
(z−r)2
Задачи.
107. Найти логарифмический потенциал, если Ω = {x2 + y 2 ≤ 1} и
√
а) µ(x, y) = 1, б) µ(x, y) = x2 + y 2 .
108. Найти ньютонов потенциал в точке M0 (x0 , y0 , z0 ), создаваемый
шаром с плотностью ρ0 = const и радиусом R.
91
109. Найти силу, с которой конус с плотностью ρ0 = const, высотой
H и радиусом основания R притягивает точку массой m, расположенной в вершине конуса.
110. Найти силу, с которой цилиндр с плотностью ρ0 = const, высотой H и радиусом R притягивает материальную точку массой m, расположенной
а) в центре основания цилиндра,
б) на границе основания цилиндра.
111. Найти силу притяжения материальной точки массой m с координатами (0, 0, h) материальным кругом
Ω = {x2 + y 2 ≤ R2 , z = 0},
по которому равномерно распределена масса с поверхностной плотностью ρ0 = const.
92
4. Примеры контрольных работ
4.1 Контрольная работа №1.
1. Используя определение меры Жордана, вычислить меру множества F на числовой прямой R, если
{
}
1
1.1. F =
|n∈N ;
10n
}
{
1
1.2. F = 1 − n | n ∈ N ;
10
{
}
1
1.3. F = 1 + n | n ∈ N ;
10
{
}
1
1.4. F = 1 − n+1 | n ∈ N ;
10
{
}
(−1)n
1.5. F =
|n∈N ;
10n
√
1.6.∗ F = [0, 2);
√
1.7.∗ F = [0, 3].
2. Дано множество M ⊂ R2 . Найти двумерную меру M , если
2.1. а) M − точка (0, 1),
б) M = {(x, y) | x ∈ (0, 1], y = 0};
2.2. а) M − точка (1, 0),
б) M = {(x, y) | x ∈ (0, 1), y = 0};
2.3. а) M − точка (1, 1),
б) M = {(x, y) | x ∈ (0, 1], y = 1};
2.4. а) M − точка (0, 0),
б) M = {(x, y) | x = 0, y ∈ [0, 1]};
2.5. а) M − точка (0, 0),
б) M = {(x, y) | x√= 0,
√ y ∈ (0, 1)};
∗
2.6 . а) M − точка ( 2, 2),
93
б) M = {(x, y) | √
x =√0, y ∈ [0, 1] ∩ Q};
2.7∗ . а) M − точка ( 3, 3),
б) M = {(x, y) | x ∈ [0, 1] ∩ Q, y = 0};
3. Пусть Ω — квадрат [0, 1] × [0, 1] на плоскости, τ — разбиение квадрата прямыми x = i/n, y = j/n, i, j = 1, . . . , n. Составить интегральную
сумму и суммы Дарбу функции f (x, y) для разбиения τ , если:
√
3.1. f (x, y) = x2 + y 2 ;
3.2. f (x, y) = (x2 + y 2 )2 ;
2 +y 2
3.3. f (x, y) = ex
;
3.4. f (x, y) = sin (x2 + y 2 );
3.5. f (x, y) = cos (x2 + y 2 );
3.6∗ . f (x, y) = ex
2 −y 2
;
3.7∗ . f (x, y) = sin (x2 − y 2 ).
2
2
∫ ∫ 4. Пусть Ω — круг {(x, y) | x + y ≤ 1}. Оценить значение интеграла
f (x, y) dxdy, если:
Ω
4.1. f (x, y) = x2 + xy + y 2 − 12x − 3y;
4.2. f (x, y) = 3 + 2x − y − x2 + xy − y 2 ;
4.3. f (x, y) = 3x + 6y − x2 − xy + y 2 ;
4.4. f (x, y) = 4x2 − 4xy + y 2 + 4x − 2y + 1;
4.5. f (x, y) = x2 + 2y 2 + 3xy − 4x − 5y;
4.6∗ . f (x, y) = x2 cos π(x2 + y 2 );
4.7∗ . f (x, y) = y 2 sin π(x2 + y 2 ).
4.2 Контрольная работа №2.
94
1. Изменить порядок интегрирования в повторных интегралах:
∫1
1.1.
∫x
dx
0
0
y∫2 +y
∫1
1.2.
f (x, y) dx;
dy
0
0
7x+10
6
∫1
∫
1.3.
dx
f (x, y) dy;
−1
2x
∫π
1.4.
cos
∫ x
dx
f (x, y) dy;
−π
−1
∫2
1.5.
∫3x
dx
1
√
2+ ∫ 4−y 2
∫3
dy
√
− 3
∗
f (x, y) dx;
√
12−y 2
∫π
1.7 .
f (x, y) dy;
ln x
√
1.6∗ .
f (x, y) dy;
sin
∫ x
dx
π
4
f (x, y) dy.
cos x
2. Вычислить двойной интеграл:
∫ ∫
2.1.
(x + y) dxdy, Ω = {(x, y) | x2 + y 2 ≤ 1, y ≥ x};
Ω
95
∫ ∫
xy dxdy, Ω = {(x, y) | x2 + y 2 ≤ 25, 3x + y ≥ 5};
2.2.
Ω
∫ ∫
2.3.
x dxdy,
Ω
Ω = {(x, y) | x + y ≤ 2, x2 − y 2 ≤ 1, x ≥ 0, y ≥ 0};
2
2
∫ ∫
(2y − x) dxdy, Ω = {(x, y) | y(y − x) ≤ 2, x(x + y) ≤ 3};
2.4.
Ω
∫ ∫
xy dxdy, Ω = {(x, y) | x2 + y 2 ≤ 1, (x − 1)2 + y 2 ≤ 1};
2.5.
Ω
∗
∫ ∫
x dxdy,
2.6 .
Ω
Ω = {(x, y) | y ≥ 0, cos x ≥ 0, sin x ≥ 0, 0 ≤ x ≤ π};
∗
∫ ∫
xy dxdy,
2.7 .
Ω
Ω = {(x, y) | x2 + y 2 ≤ 1, (x − 1)2 + y 2 ≤ 1, x2 + (y − 1)2 ≤ 1}.
3. Перейдя к полярным координатам, вычислить двойные интегралы:
∫ ∫
3.1.
cos (π
Ω
√
x2 + y 2 ) dxdy, Ω = {(x, y) | x2 + y 2 ≤ 1};
96
∫ ∫
3.2.
x2
Ω
∫ ∫
dxdy
, Ω = {(x, y) | 9 ≤ x2 + y 2 ≤ 25};
+ y2 − 1
2 +y 2
ex
3.3.
dxdy, Ω = {(x, y) | x2 + y 2 ≤ 1};
Ω
∫ ∫
ln (x2 + y 2 ) dxdy, Ω = {(x, y) | 1 ≤ x2 + y 2 ≤ 4, };
3.4.
Ω
∫ ∫
3.5.
ln (x2 + y 2 )
dxdy, Ω = {(x, y) | 1 ≤ x2 + y 2 ≤ 9};
x2 + y 2
Ω
∫ ∫
3.6∗ .
y 2 dxdy, Ω = {(x, y) | 2x ≤ x2 + y 2 ≤ 6x, y ≤ x};
Ω
∗
∫ ∫
x dxdy, Ω = {(x, y) | x ≤ x2 + y 2 ≤ 2x, y ≥ 0}.
3.6 .
Ω
4. Вычислить тройные интегралы:
∫ ∫ ∫
4.1.
(x2 + y + z) dxdydz,
Ω
Ω = {(x, y, z) | 0 ≤ x ≤ 1, 0 ≤ y ≤ 1, 0 ≤ z ≤ 1};
∫ ∫ ∫
(x2 y + z) dxdydz,
4.2.
Ω
Ω = {(x, y, z) | 0 ≤ x ≤ 1, −1 ≤ y ≤ 1, 0 ≤ z ≤ 1};
∫ ∫ ∫
4.3.
xyz dxdydz,
Ω
Ω = {(x, y, z) | 0 ≤ x ≤ 1, 0 ≤ y ≤ 1, 0 ≤ z ≤ 1};
97
∫ ∫ ∫
(x2 + y 2 ) dxdydz,
4.4.
Ω
Ω = {(x, y, z) | − 1 ≤ x ≤ 1, −1 ≤ y ≤ 1, −1 ≤ z ≤ 1};
∫ ∫ ∫
4.5.
(xy + yz + zx) dxdydz,
Ω
Ω = {(x, y, z) | 0 ≤ x ≤ 1, 0 ≤ y ≤ 1, 0 ≤ z ≤ 1};
4.6∗ .
∫ ∫ ∫
xy dxdydz,
Ω
Ω = {(x, y, z) | 0 ≤ x ≤ y ≤ z ≤ 1};
∗
∫ ∫ ∫
zex dxdydz,
4.7 .
Ω
Ω = {(x, y, z) | 0 ≤ x ≤ 1, 0 ≤ y ≤ 2, 0 ≤ z ≤ 3}.
4.3 Контрольная работа №3.
1. Исследовать на сходимость несобственный двойной интеграл:
∫ ∫
x
1.1.
dxdy, Ω = {(x, y) | x ≥ 1, y ≥ 1};
y
∫ ∫
1.2.
Ω
y
dxdy, Ω = {(x, y) | 0 < x ≤ 1, 0 < y ≤ 1};
x
Ω
∫ ∫
√
1.3.
Ω
1
x2
+ y2
dxdy, Ω = {(x, y) | x ≥ 1, y ≥ 1};
98
∫ ∫
√
1.4.
Ω
1
x2
+ y2
∫ ∫
1.5.
x2
dxdy, Ω = {(x, y) | 0 < x ≤ 1, 0 < y ≤ 1};
1
dxdy, Ω = {(x, y) | x2 + y 2 ≥ 1};
+ y2
Ω
1.6∗ .
∫ ∫
√
Ω
∗
∫ ∫
√
1.7 .
Ω
1
x2
+
y4
1
x2 + y 4
dxdy, Ω = {(x, y) | x2 + y 2 ≥ 1};
dxdy, Ω = {(x, y) | x2 + y 2 ≤ 1};
2. Найти площадь плоской фигуры, ограниченной кривыми:
2.1. 4y = x2 − 4x, x = y + 3;
2.2. y 2 = 10x + 25, y 2 = 9 − 6x;
2.3. y 2 = 2x + 1, y 2 = 1 − 2x;
2.4. x2 + y 2 = 4, y 2 = 4 − 4x, x < 1;
2.5. y = cos x, y = cos 2x, 0 ≤ x ≤
2π
;
3
2.6∗ . 2x2 + 2y 2 = 2x + 1, x2 + y 2 = 1, x2 + y 2 ≥ 1;
2.7∗ . (x + y)2 + x2 = 1.
3. Найти площадь плоской фигуры, ограниченной линиями:
3.1.
x2 y 2
+ 2 = 1;
a2
b
3.2. (x2 + y 2 )2 = y 2 ;
3.3. x2 + y 2 = x + y;
99
3.4.
√
√
x + y = 1, x + y = 1;
3.5. (x2 + y 2 )2 = 2xy;
3.6∗ . x4 + y 4 = 2xy;
3.7∗ . x3 + y 3 = xy.
4. Найти объем тела, ограниченного поверхностями:
4.1. x2 + y 2 = x, z = x2 + y 2 , z = 0;
4.2. y = x2 , y = 1, z = 0, z = x2 + y 2 ;
4.3. x + y + z = 1, 4x + y = 1, 4x + 3y = 3, y = 0, z = 0;
4.4. x2 + y 2 = π, z = sin (x2 + y 2 ), z = 0;
4.5. z 2 + y 2 = 1, x + y + z = 1, x + y + z = −1;
4.6∗ . x2 + y 2 = 1, x2 + z 2 = 1;
4.7. z 2 − x2 = 1, z 2 − y 2 = 1, z =
√
2.
100
Список литературы
1. Ильин, В.А. Основы математического анализа / В.А. Ильин, Э.Г. Позняк.
— М.: Наука, 1973. Ч. 2. 448 с.
2. Кудрявцев, Л.Д. Краткий курс математического анализа. Дифференциальное и интегральное исчисление функций многих переменных.
Гармонический анализ: Учебник. / Л.Д. Кудрявцев. — М.: Физматлит,
2002. Т. 2. 400 с.
3. Демидович, Б.П. Сборник задач и упражнений по математическому анализу: Учеб. пособие для вузов / Б.П. Демидович. — М.: ООО
«Издательство Астрель»: ООО «Издательство АСТ», 2005. 558 с.
4. Кудрявцев, Л.Д. Сборник задач по математическому анализу. Функции нескольких переменных: Учеб. пособие для вузов / Л.Д. Кудрявцев,
А.Д. Кутасов, В.И. Чехлов, М.И. Шабунин. — М.: Наука. Физматлит,
1995. 496 с.
5. Виноградова, И.А. Задачи и упражнения по математическому анализу. В 2 кн. Кн. 2. Ряды, несобственные интегралы, кратные и поверхностные интегралы: Учеб. пособие для университетов, ред. вузов / И.А. Виноградова,
С.Н. Олехник, В.А. Са -довничий; Под ред. В.А. Садовничего — 2-е изд.,
перераб. — М.: Высшая школа, 2002. 712 с.
6. Данко П.Е. Высшая математика в упражнениях и задачах. В 2-х
ч. Ч. 2: Учеб. пособие для втузов. — 5-е изд., испр. / П.Е. Данко, А.Г.
Попов, Т.Я. Кожевникова. — М.: Высшая школа, 1999. 416 с.
7. Ярахмедов, Г.Я. Математический анализ. Многомерный математический анализ / Г.Я. Ярахмедов. Новосибирск, 2001.
101
Приложение.
Таблица неопределенных интегралов элементарных функций
∫
0 · dx = C,
∫
xα dx =
∫
xα+1
+ C, ∀α ∈ R\{−1},
α+1
dx
= ln |x + a| + C,
x+a
∫
ax dx =
∫
ax
+ C, a > 0, a ̸= 1,
ln a
ex dx = ex + C,
∫
sin x dx = − cos x + C,
∫
cos x dx = sin x + C,
∫
∫
∫
dx
= tg x + C,
cos2 x
dx
= −ctg x + C,
sin2 x
∫
x
dx
1
= arctg + C, a ̸= 0,
2
2
x +a
a
a
∫
x − a
dx
1
+ C, a ̸= 0,
=
ln x2 − a2
2a x + a √
dx
√
= ln |x + x2 + a| + C, a ̸= 0,
x2 + a
∫
dx
x
√
= arcsin
+ C, a ̸= 0.
2
2
|a|
a −x
102
Некоторые интегралы, зависящие от параметра
1. Интеграл Дирихле
+∞
∫
sin (αx)
π
dx = sign α, α ∈ R,
x
2
0
где


1, если α > 0,
0, если α = 0,
sign α =

−1, если α < 0.
2. Интеграл Пуассона
+∞
√
∫
π
−αx2
, α > 0.
e
dx =
α
−∞
103
СОДЕРЖАНИЕ
Предисловие . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .3
1. Мера Жордана . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
1.1. Сеть в Rn . Кубы и элементарные множества . . . . . . . . . . . . . . . . . . 4
1.2. Мера элементарных множеств . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
1.3. Мера Жордана произвольного ограниченного
множества . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
1.4. Критерий измеримости. Цилиндрические множества . . . . . . . . . 12
1.5. Основные свойства меры Жордана . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
2. Кратные интегралы . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
2.1. Разбиение измеримых множеств . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
2.2. Интегральные суммы. Определение кратного
интеграла. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .18
2.3. Верхние и нижние суммы Дарбу. Критерий
интегрируемости . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22
2.4. Свойства кратных интегралов . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 26
2.5. Сведение двойного интеграла к повторному . . . . . . . . . . . . . . . . . . 29
2.6. Сведение тройного и n-кратного интеграла к
повторным. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .38
2.7. Замена переменных в кратных интегралах . . . . . . . . . . . . . . . . . . . 47
2.8. Криволинейные координаты . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 51
2.9. Несобственные кратные интегралы . . . . . . . . . . . . . . . . . . . . . . . . . . . 61
3. Геометрические и физические приложения кратных
интегралов . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 69
3.1. Геометрическое приложение кратных интегралов . . . . . . . . . . . . 69
3.2. Приложение кратных интегралов к геометрии масс . . . . . . . . . . 77
3.3. Приложение кратных интегралов к физике . . . . . . . . . . . . . . . . . . 84
4. Примеры контрольных работ . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 89
4.1. Контрольная работа №1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 89
4.2. Контрольная работа №2 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 90
4.3. Контрольная работа №3 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 94
Список литературы . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 97
104
Приложение. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 98
Download