metod1_KI_dlya_studentov
.pdfМосковский государственный технический университет имени Н.Э. Баумана
Факультет «Фундаментальные науки» Кафедра «Математическая физика и вычислительная математика»
Б.Т. Добрица, И.В. Дубограй, О.В. Скуднева.
Кратные интегралы и их приложения.
Москва
1
Модуль 1.
Занятие 1. Определение двойного интеграла. Геометрический смысл двойного интеграла. Теорема существования. Свойства
двойного интеграла. Сведение двойного интеграла к повторному и вычисление его в декартовых координатах. 1.Определение двойного интеграла. Теорема о его
существовании.
|
Пусть на координатной плоскости XOY дана |
||
|
замкнутая, ограниченная область D с кусочно- |
||
|
гладкой границей, имеющая площадь S(D) и |
||
|
конечный диаметр d, и пусть в этой области |
||
|
задана функция f(x,y). |
||
|
Разобьём область D на n элементарных |
||
|
областей, не имеющих общих внутренних |
||
|
точек, с площадями i и диаметрами di . |
||
Рис. 1.1 |
( di - наибольшее расстояние между точками |
||
|
элементарной области). В каждой |
||
элементарной области выберем произвольную точку Pi xi , yi , в которой |
|||
вычислим значение функции |
f Pi . |
|
|
|
n |
|
|
Составим сумму: Sn |
f Pi |
i . |
(1.1) |
i 1
Эта сумма называется интегральной .
Определение. Двойным интегралом f x, y d от функции f(x,y) по
D
области D называется предел интегральной суммы (1.1) при неограниченном возрастании числа разбиений n области D и стремлении наибольшего из диаметров разбиения к нулю (если этот предел существует).
|
|
|
|
n |
f x, y d f x, y dxdy |
lim Sn |
lim |
f xi , yi i (1.2) |
|
D |
D |
n |
n |
i 1 |
max di 0 |
max di 0 |
Область D называется областью интегрирования.
Теорема существования двойного интеграла. Если функция f(x,y)
непрерывна в замкнутой ограниченной области D (с кусочно-гладкой границей), имеющей конечную площадь и диаметр, то она интегрируема по этой области, то есть существует предел интегральной суммы (1.1) при
n и max di 0 , не зависящий от способа разбиения D на элементарные области и от выбора точек Pi .
2
2.Свойства двойного интеграла.
1.Линейность. Если функции f(x,y) и g(x,y) интегрируемы по области D, то их линейная комбинация тоже интегрируема по этой области и
f x, y g(x, y) dxdy f x, y dxdy g x, y dxdy
D D D
2.Аддитивность. Если область D разбита на две части D1 и D2 , не имеющие общих внутренних точек, то
f x, y dxdy f x, y dxdy f x, y dxdy
D |
D1 |
D2 |
3. Интеграл от единичной функции. |
Если функция f x, y 1 в |
области D, то двойной интеграл от единичной функции численно равен площади области интегрирования S(D).
1dxdy S D
D
4.Интегрирование неравенств. Если функции f(x,y) и g(x,y)
интегрируемы по области D, и f(x,y) <g(x,y) на D, то
f x, y dxdy g x, y dxdy .
D |
D |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Если |
f x, y g(x, y) |
на D, то f x, y dxdy g x, y dxdy . |
|||||||||||||
|
|
|
|
D |
|
D |
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
f x, y |
|
dxdy , |
||||
Замечание. Имеет место неравенство |
f x, y dxdy |
|
|
|
|||||||||||
|
|
||||||||||||||
|
|
|
|
|
|
|
D |
|
|
D |
|
|
|
|
|
которое следует из того, что |
|
f x, y |
|
f x, y |
|
f (x, y) |
|
и свойства 4. |
|||||||
|
|
|
|
||||||||||||
5. Теорема об оценке двойного интеграла. Если функция f(x,y) |
|||||||||||||||
интегрируема по области D, и в этой области |
|
m f x, y M , то |
m S D f x, y dxdy M S D .
D
6.Теорема о среднем. Если функция f(x,y) непрерывна в области D, то существует хотя бы одна точка P0 D , для которой выполняется
равенство f x, y dxdy f P0 S D .
D
3.Геометрический смысл двойного интеграла.
Если f(x,y)>0, то двойной интеграл f x, y dxdy
D
численно равен объёму тела, нижним основанием которого является область
3
Рисунок 1. 2
D , ограниченного цилиндрической поверхностью с образующей, параллельной оси OZ, и поверхностью z f x, y . См. рис. 1.2
V f x, y dxdy |
(1.3) |
D |
|
Рассмотрим пример непосредственного вычисления двойного интеграла. Пример 1.1. Вычислим двойной интеграл по области D от функции z xy ,
0 y 1
если D : .
0 x 1
Решение. Разобьём квадрат на n2 частей. Тогда i 1n 1n n12 .
Выберем точки в правой верхней вершине каждого квадрата. (См.рис.1.3). Значение функции z xy в этих точках будет равно
i |
|
j |
i 1, 2, |
, n; j 1, 2, |
, n . |
|
|
||||
n |
n |
По определению (см. (1.2)) имеем
Рис. 1.3
|
|
|
|
|
|
|
1 |
|
|
n |
n |
|
1 |
|
|||
x y dxdy lim |
i j lim |
1 2 3 ... n 1 2 3 ... n |
|||||||||||||||
2 |
2 |
||||||||||||||||
D |
|
|
|
n n |
|
|
i 1 |
j 1 |
n n |
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
1 |
n |
|
n 1 |
2 |
|
|
|
1 |
|
|
|
|
|
|||
lim |
|
|
|
|
|
|
. |
|
|
|
|
||||||
2 |
|
|
2 |
4 |
|
|
|
|
|||||||||
n n |
|
|
|
|
|
|
|
|
|
|
В подавляющем большинстве случаев вычисление пределов интегральных сумм представляет собой сложную и трудновыполнимую задачу. Для достижения результата эффективнее использовать повторный интеграл.
4. Повторный интеграл.
|
|
|
|
x2 |
y2 x |
f x, y dy (1.4а) |
Повторным называется двойной интеграл вида |
dx |
|
||||
|
|
|
|
x1 |
y1 x |
|
y |
|
x y |
|
|
|
|
или 2 |
dy |
2 |
f x, y dx (1.4б). |
|
|
|
y1 |
|
x1 y |
|
|
|
|
4
При вычислении первого из указанных повторных интегралов сначала
|
|
|
|
|
|
|
|
|
|
y2 x |
f x, y dy |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
берётся внутренний интеграл |
|
по переменной y , при этом |
x |
||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
y1 x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
играет роль параметра. Затем полученную функцию интегрируют по |
|
|
|||||||||||||||||||||||||||||||
переменной x. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
x2 |
|
y2 |
x |
|
x2 |
y2 |
x |
|
|
|
x2 |
|
|
|
|
|
2 |
|
|
|
|
|
|
|
1 |
|
|
, |
|
||||
|
|
|
f x, y dy |
|
|
|
|
x, y dy dx |
|
|
|
|
x |
F |
|
|
x |
|
|||||||||||||||
|
dx |
|
|
|
|
|
f |
|
|
F |
|
|
x, y |
|
|
|
|
x, y |
dx |
|
|
||||||||||||
|
|
y1 x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x1 |
|
|
x1 y1 x |
|
|
|
x1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
где |
F x, y - первообразная функции |
f x, y по y , т.е. |
Fy x, y f x, y . |
||||||||||||||||||||||||||||||
Аналогично, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
y2 |
|
x2 y |
|
|
y2 |
x2 |
y |
|
|
|
|
|
y2 |
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
1 |
|
|
, |
||
|
|
|
|
f x, y dx |
|
|
|
f x, y dx dy |
|
|
|
y |
|
F |
|
||||||||||||||||||
|
dy |
|
|
|
|
|
|
|
F |
x |
, y |
|
x y , y |
dy |
|
||||||||||||||||||
|
|
x1 y |
|
|
|
|
|
y |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
y1 |
|
|
|
y1 x1 |
|
|
|
|
|
y1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
где |
Fx x, y f x, y . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
x |
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dx |
|
x |
dy |
|
|
|
|
|
|
|||||
Пример1.2. Вычислим повторный интеграл : |
|
2 |
. |
|
|
|
|
||||||||||||||||||||||||||
1 |
|
1 |
|
|
y |
|
|
|
|
|
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x
Решение. Вычислим сначала внутренний интеграл, при этом множитель x2 , играющий роль параметра (не зависящего от переменной интегрирования), вынесем за знак этого интеграла.
2 |
x |
x2 |
2 |
x |
dy |
||||
dx |
|
|
dy x2dx |
|
|
||||
y |
2 |
y |
2 |
||||||
1 |
1 |
|
|
1 |
1 |
|
|
||
|
|
|
|
|
|
||||
|
x |
|
|
|
x |
|
|
|
x |
4 |
|
x |
2 |
|
|
2 |
9 |
|
|
|
|
||||||||
|
|
|
|
|
|
|
. |
|||
|
|
|
|
|
||||||
|
4 |
2 |
|
|
1 |
4 |
|
|||
|
|
|
|
|
|
|
|
|
|
21
x2 y
x |
|
2 |
|
|
dx x2 |
||
|
x |
||
1 |
|
|
|
|
|
1 |
|
x |
|
||
|
|
|
1 |
2 |
|
|
|
|
||
|
dx x3 |
x dx |
|
|
|||
x |
1 |
|
|
|
|
|
5.Переход от двойного интеграла
кповторному. Расстановка пределов
интегрирования.
Для вычисления значения двойного интеграла (1.1) необходимо перейти к повторному интегралу вида (1.4а) или (1.4б) . Способ перехода зависит от вида области D.
5
Рисунок 1.4
Область, ограниченная слева и справа прямыми
x x1 и x x2 , а сверху и снизу - кривыми y y1 x и y y2 x , такими, что любая прямая x x0 x1 x0 x 2 , проведённая через область D, пересекает каждую из кривых y y1 x и y y2 x в одной точке,
называется y-правильной. При этом точку А называют точкой входа в область, а точку В – точкой выхода из области.
(см.рис 1.4) То есть область D является y-
правильной, если
D x, y | x1 x x2 , y1 x y y2 x
Определив пределы изменения переменных по области D, расставляем пределы интегрирования в повторном интеграле:
|
|
|
|
|
x |
|
y x |
|
|
|
|
f x, y dxdy 2 |
dx |
2 |
f x, y dy . |
||
|
|
|
|
D |
x1 |
|
y1 x |
|
|
|
|
|
|
(1.4а) |
|
|
|
|
|
|
|
|
Рисунок 1.5 |
|
||
Аналогично вводится понятие x-правильной области (см. рис. 1.5). |
||||||||
В этом случае имеем D x, y | y1 |
y y2 , x1 y x x2 y и |
|
||||||
|
y |
|
x y |
|
|
|
|
|
f x, y dxdy 2 |
dy |
2 |
f x, y dx |
(1.4б) |
|
|||
D |
y1 |
|
x1 y |
|
|
|
|
|
Пример 1.3. Расставим двумя способами пределы интегирования в двойном
интеграле f x, y dxdy , где D – треугольник с вершинами O(0,0),
D
A(0,1),В(1,1).
Рисунок 1.5 Решение. Построим область D. Она одновременно является как
y-правильной, так и x-правильной. Поэтому первый способ расстановки пределов при переходе к повторному интегралу выглядит следующим образом.
Для всех точек данной области переменная x изменяется между x1 0 и x2 1, а переменная y для любого из этих x меняется
от y y1 x 0 до y y2 x x , т.е.
Рис. 1.6
6
D x, y | 0 x 1,0 y x .
Определив пределы изменения переменных по области D, расставляем пределы интегрирования в повторном интеграле:
|
x |
|
y x |
|
f x, y dxdy 2 |
dx |
2 |
f x, y dy |
|
D |
x1 |
|
y1 x |
|
1 |
x |
|
|
|
dx f x, y dy . |
|
(См.(1.4а)) |
||
0 |
0 |
|
|
|
Изменим порядок интегрирования, то есть расставим пределы интегрирования в повторном интеграле вторым способом (См.(1.4б)). Из рис. 1.6 следует, что для всех точек данной области переменная y
изменяется между y1 0 и y2 |
1, а переменная x для любого из этих y |
||||||
|
1 |
|
|
2 |
|
|
|
меняется от |
x x y y до x |
x y 1. Т.е. D |
x, y | 0 |
y |
1, y x 1 |
||
Теперь расставляем пределы интегрирования: |
|
|
|
||||
|
y |
|
x y |
1 |
1 |
|
|
f x, y dxdy 2 |
dy |
2 |
f x, y dx dy f x, y |
dx . |
|
||
D |
y1 |
|
x1 y |
0 |
y |
|
|
Пример 1.4. Изменим порядок интегрирования в повторном интеграле
2 |
2 |
x |
|
||
dx |
|
|
|
|
f x, y dy . |
0 |
2 x x2 |
Решение. Из данного интеграла следует, что область интегрирования ограничена прямыми x 0 , x 2 и кривыми y 2x x2 , y 2x .
Построим область D (См. рис. 1.7). Она является y-правильной, но не является x-правильной (См. рис .1.4). Для получения x-правильных областей придётся разбить область D на три части. Определим пределы, в которых меняются переменные по каждой из трёх областей для повторных интегралов
вида (1.4б).
0 y 1
D1 : y2 x 1 1 y2 ,4
|
0 y 1 |
|
|||
|
|
|
|
|
|
D2 |
: |
|
|
|
|
1 y2 |
x 2 |
||||
|
1 |
||||
|
|
|
|
|
7
Рис. 1.7
|
1 y 2 |
|
||
|
2 |
|||
|
|
|
|
|
D3 |
: y2 |
|||
|
|
|
x 2 |
|
|
|
|||
|
|
4 |
|
|
По свойству аддитивности имеем
f |
x, y dxdy f x, y |
dxdy f x, y dxdy f x, y dxdy . |
|
|||||||||||||||||
D |
|
|
|
|
D1 |
|
|
|
|
|
|
D2 |
|
|
|
D3 |
|
|
|
|
Переходя к повторным интегралам вида (1.4б), расставляем пределы |
|
|||||||||||||||||||
интегрирования в каждом из трёх слагаемых и получаем ответ: |
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
2 x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
1 |
1 1 y |
2 |
|
1 |
|
2 |
2 |
2 |
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|||||||||||
dx |
|
|
f x, y dy dy |
|
|
|
|
f x, y dx dy |
|
|
f x, y dx dy f |
x, y dx . |
||||||||
0 |
2 x x2 |
|
|
y2 |
|
|
|
|
|
|
|
|
|
y2 |
|
|
||||
0 |
|
|
|
0 |
|
1 1 y2 |
1 |
|
||||||||||||
|
|
|
|
|
|
4 |
|
|
|
|
|
|
|
|
|
4 |
|
|
||
Пример 1.5. Вычислим |
x2 |
y2 dxdy , |
где область D ограничена линий |
|||||||||||||||||
|
|
|
|
|
|
D |
|
|
|
|
|
|
|
|
|
|
|
y x2 и y 2x 3 .
Решение.
Рис. 1.8
1 x 3 и x2 y 2x 3.
Линии, задающие область D , пересекаются в точках с абсциссами x 1
и x 3 . Построим область D и выясним, является ли она правильной (См. рис.1.8). Сравнив данную область с областями, представленными на рис. 1.3 и рис. 1.4, видим, что наша область является y-правильной, но не является x –правильной.
Перейдём от двойного интеграла к повторному вида (1.4а) и вычислим его.
Учтём, что по области D переменные меняются в следующих пределах:
Расставим пределы интегрирования и вычислим интеграл.
|
3 2 x 3 |
3 |
|
|
|
3 |
|
|
2 x 3 |
|||
|
|
|
|
|||||||||
x2 y2 dxdy dx |
|
x2 y2 dy |
x2 |
y |
y |
|
|
|
|
dx |
||
|
|
|
|
|||||||||
D |
1 |
x |
2 |
1 |
|
3 |
|
|
x |
2 |
||
|
|
|
|
|
|
|
|
|
8
3 |
|
|
2x 3 |
3 |
|
x |
2 |
|
3 |
|
|
4 |
|
2x 3 |
4 |
|
|
5 |
|
|
7 |
|
3 |
|
||
|
x2 |
2x 3 |
|
x2 x2 |
|
|
dx |
x |
x3 |
|
|
x |
|
|
x |
|
|
|
|
|
6592 |
|||||
|
|
|
|
|
|
|
2 4 |
|
|
|
|
|
|
|
||||||||||||
|
|
|
3 |
|
|
3 |
|
|
|
|
2 |
|
|
5 |
21 |
|
|
|
35 |
|||||||
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Занятие 2
Криволинейные координаты на плоскости. Замена переменных в двойном интеграле. Геометрический смысл модуля и знака Якобиана. Двойной интеграл в полярных координатах. Вычисление площади плоской фигуры с помощью двойного интеграла.
Ауд.: ОЛ-6 №№ 2160, 2162, 2164, 2166, 2171, 2177, 2183 или ОЛ-5 №№ 8.43,
44, 46, 48, 50, 56, 62.
1.Замена переменных в двойном интеграле.
Вычисление двойного интеграла может оказаться затруднительным, и тогда ищут другую систему координат, в которой интегрирование упростилось бы. При этом происходит преобразование переменных, а также и области интегрирования.
Рисунок 2.1
Пусть взаимно однозначное отображение области D' плоскости UOV
x x u,v |
|
|
|
в область D плоскости XOY задано функциями y y u,v , |
(2.1) |
|
|
|
|
причём обе функции непрерывно дифференцируемы в области D' . (см. рис 2.1) При этом отображении площадь элементарной части
1i ui vi области D1 и площадь i xi yi соответствующей части
9
области D связаны так, что i |
|
|
J u,v |
|
1i , где |
J u,v - якобиан этого |
|||||
|
|
||||||||||
отображения, отличный от нуля в любой точке области D' , кроме, быть |
|||||||||||
может, конечного числа точек (линий), т.е. |
|
||||||||||
|
D x, y |
|
|
x |
x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|||||
J |
|
|
u |
v |
0 |
|
|||||
|
y |
y |
|
||||||||
|
D u, v |
|
|
|
|
|
|
(2.2) |
|||
|
|
|
|
u |
v |
|
|
|
|
|
|
Тогда справедлива формула замены переменных для двойного интеграла:
f x, y dxdy f x u,v , y u,v |
|
J |
|
dudv |
(2.3) |
|
|
|
|||||
|
|
|||||
D |
D ' |
|
Замечание.При переходе в новую систему координат площадь плоской
фигуры D (См. свойства двойного интеграла) |
вычисляется по формуле: |
|||||
S dxdy |
|
J |
|
dudv |
(2.4) |
|
|
|
|||||
|
|
|||||
D |
D ' |
|
||||
2.Геометрический смысл модуля и знака Якобиана. |
Модуль якобиана есть предел отношения площадей бесконечно малых площадок D и D΄.
Геометрический смысл знака якобиана состоит в том, что в случае положительного якобиана ориентации контуров, ограничивающих области D' и D при отображении совпадают, а при отрицательном — отличаются.
3.Двойной интеграл в полярной системе координат (как частный случай).
Связь между переменными в декартовой и полярной системах координат:
(см. рис. 2.2)
x r cos |
, tg |
y |
, r2 |
x2 y2 , |
|
|
|||||
|
|||||
y r sin |
|
x |
|
(2.5) |
|
r 0, , |
|
|
|
||
0 2 |
|
||||
Рисунок 2.2 |
|
|
|
|
Якобиан при переходе в полярную систему координат: J , r . (Убедитесь в этом, вычислив якобиан
самостоятельно по формуле (2.2)) |
|
|
Теперь формула (2.3) примет вид: |
|
|
f x, y dxdy f r cos , r cos rd dr |
(2.6) |
|
D |
D |
|
10
Рисунок 2.3