9388

Используем независимость предела интегральных сумм для интегрируемой функции от способа разбиения исходной области на

Поэтому для двойного интеграла в декартовых координатах используется ещё одно обозначение

∫∫ f ( x, y)dxdy .

D

Достаточным условием существования двойного интеграла является непрерывность функции f (x, y) в области D , граница которой

представляет собой кусочно-непрерывную кривую.

Пользуясь рассуждениями предыдущего параграфа, мы приходим к геометрическому смыслу двойного интеграла. Он состоит в том, что для функции z = f (x, y) ³ 0 двойной интеграл равен объёму V цилиндрического тела, ограниченного поверхностью z = f ( x, y) . Итак,

V = ∫∫ f ( x, y)ds .

D

Если же функция z = f (x, y) £ 0 , то V = − ∫∫ f ( x, y)ds

D

Рис. 49.3

Полученный результат можно использовать для вычисления объёма

xOу является область D (рис. 49.3).

62

Объём такого тела равен разности объёмов двух цилиндрических тел

V = ∫∫ f1 (x, y)dxdy - ∫∫ f2 (x, y)dxdy .

49.3. Свойства двойного интеграла. Поскольку двойной интеграл вводится по той же схеме, по которой вводился определённый интеграл функции одной переменной на отрезке, то свойства этих интегралов аналогичны. Они вытекают из рассмотрения соответствующих интегральных сумм. В частности, двойной интеграл суммы двух (или любого конечного числа) функций по области D равен сумме двойных интегралов по этой области каждой из функций в отдельности

∫∫( f (x, y)+ g(x, y))dxdy = ∫∫ f (x, y)dxdy + ∫∫ g(x, y)dxdy .

D D D

Постоянный множитель можно выносить за знак двойного интеграла, т.е. при c = const верно равенство

∫∫cf (x, y)dxdy = c∫∫ f (x, y)dxdy .

D D

Если область D разбита на конечное число областей D1 , D 2 , ..., D n без общих внутренних точек, то

аналогичное неравенство

∫∫ f (x, y)dxdy ³ ∫∫ g(x, y)dxdy .

D D

Если в качестве подынтегральной функции рассмотреть функцию, тождественно равную единице ( f (x, y) ≡ 1) , то в интегральной сумме все

слагаемые будут представлять площади частичных областей

интегральных сумм. Тем самым, получаем способ вычисления площади плоской области D с помощью двойного интеграла

∫∫1d x d y = S .

D

Из последних двух свойств вытекает, что если m ≤ f ( x, y) ≤ M в области D , то

m S ≤ ∫∫ f ( x , y ) d x d y ≤ M S .

D

Это даёт оценку отношения двойного интеграла к площади области

m ≤ S1 ∫∫D f ( x, y )d xd y ≤ M .

Введём понятие среднего значения функции f (x, y) в области D

f ср = S1 ∫∫D f ( x , y ) d x d y .

Его геометрическая интерпретация состоит в следующем. Если дана

которого равен объему исходного цилиндрического тела.

Поскольку функция, непрерывная в замкнутой и ограниченной области D , принимает все свои промежуточные значения между наименьшим и наибольшим, то в этой области найдётся точка (ξ , η ) , в которой достигается среднее значение функции, т.е.

∫∫ f ( x, y)dxdy = f (ξ, η) S .

D

Последнее свойство является аналогом теоремы о среднем значении для определённого интеграла функции одной переменной.

64

Лекция 50. Вычисление двойного интеграла в декартовых

иполярных координатах

50.1.Вычисление двойного интеграла в прямоугольных декартовых координатах. Чтобы вычислить двойной интеграл, следуя введенному выше определению, мы должны были бы находить предел интегральных сумм, доказывая при этом его независимость от разбиения области и выбора точек в каждой подобласти. Такой способ при необходимости применяется в приближенных вычислениях двойных интегралов. Однако, подобно тому, как определенный интеграл по отрезку, являющийся пределом соответствующих интегральных сумм, находится с помощью формулы Ньютона-Лейбница, двойной интеграл может быть вычислен путем двух операций последовательного интегрирования.

Для этого нужно, прежде всего, точно описать область, по которой ведётся интегрирование. Область D назовем правильной в направлении

выполняются условия c ≤ y ≤ d , ψ1 ( y ) ≤ x ≤ ψ2 ( y) .

Это означает, что для правильной в направлении оси Oy (Ox) области всякая прямая, параллельная оси Oy (Ox) и проходящая через

внутреннюю точку области, пересекает границу области ровно в двух точках.

На рис. 50.1 изображены области, правильные в направлении обеих осей.

65

Рис. 50.1

Пусть z = f ( x, y ) ³ 0 – непрерывная на этой области функция. Напомним, что объём тела можно вычислить как интеграл от площадей его плоских сечений. А именно, если для каждого x [a,b] известна площадь S (x) сечения тела плоскостью, перпендикулярной оси Ox , то его объём

b

выражаетcя формулой V = ∫ S ( x )dx .

a

Используем эту формулу для вычисления двойного интеграла, помня, что он равен объёму цилиндрического тела. Сечение цилиндрического тела плоскостью x = x0 представляет собой криволинейную трапецию ABCE

(рис. 50.2).

Рис. 50.2

Если эту трапецию спроецировать параллельно оси Ox на плоскость

Рис. 50.3

66

Интеграл в правой части формулы (50.1) носит название повторного интеграла. Для его нахождения сначала вычисляется внутренний интеграл

вычисляется внешний интеграл. Для таких областей используется порядок интегрирования, при котором внешним является интеграл по переменной y . Рассуждая так же, как для областей, правильных в направлении оси Oy

Если область D не является правильной в выбранном направлении, то её разбивают на правильные области. Например, область на рис. 50.4 разбита прямой y = y0 на области D1 и D2 , являющиеся правильными

в направлении оси Oy .

Рис. 50.4

Если область D является прямоугольником со сторонами, параллельными осям координат, для координат точек которого

Заметим, что внешние пределы в повторном интеграле всегда постоянны, а во внутреннем интегррале пределы бывают постоянными числами только в том случае, если D – прямоугольник.

Рассмотрим для примера вычисление с помощью двоойного интеграла объема тела, ограничеенного гиперболическим параболоиддом z = xy

цилиндром, y = x2 и плоскостью y = x (рис. 50.5).

Рис. 50.5

Как видим, область D является правильной в направлении обеих осей, тем самым можно использовать любой порядок интегрирования.

Заменим двойн ой интеграл повторным по формуле (50.1)

Найдём сначала внут ренний интеграл по переменной y , считая

переменную x константой, затем проинтегрируем по x полученную функцию, зависящую уже только от переменной x

пределы останутся прежними, так как в области D переменная y меняется от 0 до 1, а для нахождения функций ψ1( y) и ψ2 ( y) следует выразить x через переменную y из уравнений линий, ограничивающих D . По формуле (50.2) получается, что

68

50.2. Вычисление двойного интеграла в полярных координатах.

Все проведённые до сих пор вычисления выполнялись в прямоугольной декартовой системе координат. Использование полярной системы координат даёт возможность упростить вычисление двойного интеграла во многих случаях, например, если область интегриррования является сектором круга или кольца.

Введём на плоскости полярную систему координат, полюс которой совпадает с началом декартовой системы координат, а полярная ось – с осью Ox . В этом случае полярные координаты связаны с декартовыми

Эти соотношения позволят перейти в записи уравнений границы области от декартовых координат к полярным.

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

69

В правой части последней формулы мы видим предел интегральной суммы для функции , поэтому получаем формулу вычисления двойного интеграла в полярных координатах

Коротко можно сказать, что при переходе в дво йном интеграле к

производится повторным интегрированием, только теперь по переменным r и ϕ . Область интегрирования D должна быть в этом случае правильной по отношению к полярной системе координат (рис. 50.7). Это означает, что

она будет ограничена лучами ϕ = ϕ1 и ϕ = ϕ2 , а также кривыми r = r1 (ϕ) и

r = r2 (ϕ) , то есть для точек области выполнены условия ϕ1 ≤ϕ≤ ϕ2 , r1(ϕ) ≤ r ≤ r2 (ϕ) .

70

Рис. 50.7

Внешний интегр ал берётся по переменной ϕ , и формула замены двойного интеграла повторным в полярных координатах выглядит следующим образом:

Рис. 50.8

Если полюс принадлежит области интегрирования (рис. 50.8), то в этой

Рис. 50.9

71