Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Краткий курс математического анализа. Том 2

.pdf
Скачиваний:
396
Добавлен:
10.05.2015
Размер:
4.49 Mб
Скачать

230Гл. 5. Интегральное исчисление функций многих переменных

Вчастности, взяв поочередно две функции из трех P , Q и R тождественно равными нулю, будем иметь

P dy dz =

P cos α dS,

 

S+

S

 

Q dz dx =

Q cos β dS,

(47.8)

 

 

S+

S

 

R dx dy =

R cos γ dS.

 

S+

S

 

Интуитивный смысл этих формул состоит в том, что элемент площади dS (см. (46.25)) данной поверхности, умноженный на косинус угла, который он составляет с некоторой координатной плоскостью, «приближенно» равен площади его проекции на рассматриваемую координатную плоскость, как если бы речь шла о площади плоской фигуры и ее проекции. На рис. 49 изоб-

ражен случай, когда указанное проектирование производится на плос-

кость переменных x и y.

В этом случае

dx dy ≈ dS cos γ,

где dS — элемент площади поверхности S, dx dy — элемент площади проекции этой поверхности на плоскость переменных x, y, а γ — угол

между dS и указанной плоскостью, очевидно, равный углу между нормалью n к поверхности S и единичным ортом k к оси z.

Если за непрерывную единичную нормаль на поверхности S взять вектор n и обозначить через Sповерхность S, ориентированную с помощью этой нормали, то согласно определению поверхностного интеграла второго рода получим

 

a dS = a(n) dS.

 

(47.9)

 

S

S

 

 

 

 

Очевидно, что

 

 

 

 

 

 

a dS =

a(

n) dS =

an dS =

a dS.

(47.9)

 

(47.5)

S+

S

S

 

 

S

 

Таким образом,

 

a dS =

 

 

 

 

 

a dS,

 

(47.10)

S

S+

§ 47. Поверхностные интегралы

231

т. е. при изменении ориентации поверхности интеграл второго рода меняет знак.

Если векторная функция a непрерывна на поверхности S, то ин-

тегралы a dS и

a dS существуют, так как в силу сказанного вы-

S+

S

ше существуют интегралы, стоящие в правых частях равенств (47.5) и (47.9).

З а м е ч а н и е. Если поверхность кусочно-гладкая (см. п. 46.5), то поверхностные интегралы первого и второго рода по ней определяются как суммы соответствующих интегралов по ее гладким частям.

47.2. Формулы для представления поверхностного интеграла второго рода в виде двойного интеграла. Пусть S

гладкая поверхность. Если она задана своим векторным представлением S = {r(u, v); (u, v) G}, то обычно через S+ обозначается поверхность S, ориентированная с помощью вектора

Вспомнив, что

a dS =

(47.5)

S+

S

n = ru ×rv .

|ru ×rv |

g g

2

=

r

u ×

r

v|

,

получим

11

22 − g12

(46.16)|

 

 

 

an dS = a ru ×rv |ru ×rv| du dv =

|ru ×rv |

G

= (a, ru, rv) du dv, (47.11)

 

 

 

 

 

 

G

 

где (a, ru, rv) — смешанное произведение векторов a, ru, rv.

 

Если

 

 

 

 

 

 

 

 

a = (P , Q, R),

r(u, v) = (x(u, v), y(u, v), z(u, v)),

 

ru = (xu, yu, zu),

rv = (xv, yv , zv),

 

то формула (47.11) в координатной форме имеет вид

 

S

 

G

 

P

Q

R

 

 

 

a dS =

 

 

 

zu

 

(47.12)

 

xu yu

du dv,

 

+

 

 

 

yv

 

 

 

 

 

xv

zv

 

где P = P (x(u, v), y(u, v), z(u, v))

,

Q = Q(x(u, v), y(u, v), z(u, v)), R =

= R(x(u, v), y(u, v), z(u, v)).

Если ориентированная поверхность S получена с помощью склейки из ориентированных поверхностей Si, i = 1, 2, ..., k (см. п. 46.5), то по определению

k

 

a dS =

a dSi.

Si=1 Si

232Гл. 5. Интегральное исчисление функций многих переменных

47.3.Некоторые специальные случаи поверхностных интегралов второго рода. Если поверхность S имеет явное представ-

ление z = f (x, y), (x, y)

G

,

то

= 0 1 = 1

 

 

 

 

 

 

xv

yv

 

 

 

 

 

 

 

 

 

 

yu

 

 

0

 

 

 

 

 

 

 

xu

 

1

 

 

 

(здесь x = u,

 

 

 

 

 

 

+

по верхней стороне

y = v), и потому

для

интеграла

S (совпадающей в данном случае с S ) поверхности S (см. пример

в п. 46.6) в случае P ≡ 0 и Q ≡ 0 на S будем иметь

 

 

 

 

R dx dy =

 

 

R(x, y, f (x, y)) dx dy,

(47.13)

 

 

 

 

(47.12)

 

 

 

 

 

 

 

 

 

 

 

 

G

 

 

 

 

 

 

S

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а для интеграла по нижней стороне S (совпадающей здесь с S)

поверхности S

 

 

 

 

 

 

 

 

 

 

 

 

 

R dx dy =

 

R(x, y, f (x, y)) dx dy.

(47.14)

 

 

 

S

 

 

 

G

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

При

использовании координатной записи поверхностного интегра-

 

 

 

 

 

 

 

 

 

R dx dy иногда пишут просто

ла второго рода (см. (47.6)) вместо

R dx dy,

 

 

 

 

 

 

 

 

S

 

R dy dx,

 

 

 

R dx dy соответственно

 

S

 

а вместо

 

 

 

 

 

 

 

S

 

т. е. в случае

 

 

 

 

S

 

 

 

 

 

 

 

нижней стороны S поверхности S дифференциалы dx и dy пишутся в обратном порядке.

Равенства (47.13) и (47.14) получены в предположении непрерывности функции R и гладкости поверхности S (т. е. непрерывной дифференцируемости на G функции f ). Однако интегралы в правой части этих равенств имеют смысл, например, лишь при предположении непрерывности функций R и f. Поэтому естественно, если функция f (x, y) непрерывна на замыкании G области G, а функция R = R(x, y, z) непрерывна на поверхности S = {z = f (x, y), (x, y) G},

определить интегралы R dx dy и R dx dy равенствами

S

S

S

S

R dx dy =

 

R(x, y, f (x, y)) dx dy,

 

G

(47.15)

 

 

R dx dy =

 

R(x, y, f (x, y)) dx dy.

 

G

 

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

аследовательно, ее график является гладкой поверхностью (в этом случае x = u, y = v, |rx ×ry|2 = 1 + fx2 + fy2 = 0), то, согласно доказан-

§ 47. Поверхностные интегралы

233

ному выше, определение (47.15) равносильно для рассматриваемых здесь поверхностей определению (47.5).

Подчеркнем, что в сделанных здесь определениях

определяются

не в отдельности понятия верхней S и нижней S стороны поверхно-

сти S, а лишь два интеграла, которые обозначаются соответственно

символами

и

. По аналогии со случаем гладкой поверхности

S

S

первый из них называется поверхностным интегралом по верхней, а второй — по нижней стороне поверхности S. В случае гладкой поверхности S эти интегралы совпадают с определенными раньше (см. (47.13) и (47.14)).

Нам потребуется еще понятие интеграла по определенного вида цилиндрическим поверхностям без предположения о гладкости этих поверхностей.

Если Γ0 = {x(u), y(u); a u b} — кривая на плоскости пере-

менных x, y, то поверхность S = {x(u), y(u), v; a u b, 0 v h} называется цилиндрической поверхностью с направляющей Γ0 и об-

разующей длины h, параллельной оси Oz (см. замечание 13 в п. 42.1). Если Γ0 — гладкая кривая, т. е. функции x(u), y(u) непрерывно дифференцируемы и x2u + yu2 > 0 на отрезке [a, b], а r(u, v) = x(u) i +

+ y(u) j + v k, то

 

 

ru = (xu, yu, 0),

rv = (0, 0, 1),

 

 

0

0

1

 

 

 

ru ×rv =

xiu

yju

k

 

 

 

 

0 = yu i xu j = 0, (ru ×rv) k = 0.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Следовательно, если

n =

ru ×rv

, то n k = 0, т. е. единичная нормаль

 

 

 

 

|ru ×rv |

 

 

 

n = (cos α, cos β, cos γ) к поверхности S перпендикулярна к оси z,

и поэтому cos γ = 0. Отсюда

 

 

 

R dx dy =

 

R cos γ dS = 0,

R dx dy =

R dx dy = 0.

S+

S

 

 

 

 

S

S+

В случае отсутствия предположения о гладкости кривой Γ0 эти

равенства примем за определение интегралов R dx dy и

R dx dy

S+

S

по цилиндрической поверхности S рассматриваемого вида или по поверхности S, являющейся частью такой цилиндрической поверхности:

R dx dy =

def

0.

(47.16)

R dx dy =

S+

S

 

 

234 Гл. 5. Интегральное исчисление функций многих переменных

Аналогичным образом рассматриваются интегралы вида

P dy dz и

Q dz dx.

S±

S±

В заключение отметим, что в теории поверхностных интегралов понятие гладкой поверхности бывает удобно понимать в более общем смысле, чем это было определено в п. 46.1 (это делается по аналогии с расширением понятия гладкой кривой; см. п. 45.4 ).

Поясним это на примере поверхности S, имеющей явное представ-

ление:

S = {z = f (x, y); (x, y) G},

G — квадрируемая область, функция R непрерывна на поверхности S, т. е. непрерывна функция R(x, y, f (x, y)), (x, y) G.

Пусть существует такое непрерывно дифференцируемое преобразование параметров

x = x(u, v), y = y(u, v),

(x, y)

= 0, (u, v) D,

(47.17)

(u, v)

D — квадрируемая область, что функция z = f (x(u, v), y(u, v)) является непрерывно дифференцируемой на замыкании D области D (функция f предполагается только непрерывной). Тогда поверхность S, задаваемая представлением

x = x(u, v), y = y(u, v), z = f (x(u, v), y(u, v)) (u, v) D,

является непрерывно дифференцируемой, а если у нее нет особых точек, то в каждой ее точке существует единичная нормаль n = (cos α, cos β, cos γ). Сделав в интеграле, стоящем в правой части равенства (47.15), замену переменной (47.17) и воспользовавшись формулой (47.12) (при P = Q = 0), получим

R(x, y, z) dx dy =

R(x, y, f (x, y)) dx dy =

 

 

 

 

 

 

S

G

 

 

 

 

 

 

= R(x(u, v), y(u, v), f (x(u, v), y(u, v)))

(x, y)

du dv

=

 

 

(u, v)

 

 

 

 

 

47.12)

D

 

 

 

 

((47.8)

 

 

=

R(x, y, z) cos γ dS.

 

(47.12)

 

 

 

 

 

(47.8)

S1

 

 

 

Таким образом, несмотря на то, что в этом случае исходное явное представление поверхности S не является, вообще говоря, дифференцируемым, для нее, какова бы ни была непрерывная на ней функция F , при надлежащем выборе параметров справедлива формула

F (x, y, z) dx dy = F (x, y, z) cos γ dS.

(47.18)

S

S1

§ 48. Скалярные и векторные поля

235

Обычно в правой части этого равенства вместо S1 пишут S, рассматривая S и S1 как «одну и ту же поверхность» с разными параметризациями.

В результате преобразования параметров (47.17) из поверхности S получилась гладкая поверхность. Один из простейших примеров по-

добной ситуации дает полусфера, заданная как график функции

z = 1 − x2 − y2 , x2 + y2 1. Если сделать преобразование параметров x = cos ϕ cos ψ, y = sin ϕ cos ψ, 0 ϕ 2π, 0 ψ π/2,

то та же полусфера будет задаваться непрерывно дифференцируемым представлением

x = cos ϕ cos ψ, y = sin ϕ cos ψ, z = sin ψ, 0 ϕ 2π, 0 ψ π/2,

и, следовательно, для нее справедлива формула (47.18). Оказывается целесообразным следующим образом расширить по-

нятие гладкой поверхности. Поверхность {M (u, v), (u, v) G} называется гладкой, если существует такое преобразование параметров

u = u(u1, v1)), v = v(u1, v1), (u1, v1) G1,

что поверхность {M (u(u1, v1), v(u1, v1)); (u1, v1) G1} является гладкой в старом смысле, т. е. в смысле определения, данного в п. 46.1.

Поверхность, являющаяся объединением гладких в новом смысле поверхностей, называется, как и раньше, кусочно-гладкой.

§48. Скалярные и векторные поля

48.1.Основные понятия. Действительную функцию u(x, y, z), заданную на некотором множестве E R3, называют также скалярным полем на этом множестве, в отличие от векторных функций, называемых векторными полями.

Всякой дифференцируемой на области G R3 функции (или, что то же самое, дифференцируемому на G скалярному полю) u(x, y, z) соответствует векторное поле ее градиентов

u =

∂u

, ∂u∂y ,

∂u

.

 

 

 

 

(48.1)

∂x

∂z

Если ввести символический вектор =

 

,

,

, то равенство

 

∂x

∂y

∂z

(48.1) можно рассматривать как «формальное» произведение этого

 

 

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

u = ∂x, ∂y, ∂zu.

Иначе говоря, — это оператор (векторнозначная функция), определенный на множестве дифференцируемых на области G функций.

236 Гл. 5. Интегральное исчисление функций многих переменных

Уравнение касательной плоскости в точке (x0, y0, z0) к поверхности уровня функции u, т. е. к поверхности, задаваемой неявно уравнением u = const, имеет вид (см. (46.13))

∂u ∂u ∂u 0,

(x − x0) ∂x + (y − y0) ∂y + (z − z0) ∂z =

где значения частных производных функции u берутся в точке

(x0, y0, z0).

Отсюда в силу геометрического смысла коэффициентов уравнения плоскости и из формулы (48.1) видно, что градиент u перпендикулярен поверхности уровня функции u (т. е. касательной плоскости к этой поверхности уровня).

Если G — плоская область, G R2xy, то градиент функции u =

=u(x, y), определенный в G, имеет вид

u = ∂u , ∂u∂y .∂x

Поскольку уравнение касательной прямой к линии уровня u = = const функции u имеет вид (см. замечание в п. 40.1)

(x − x0) ∂u∂x + (y − y0) ∂u∂y = 0,

то градиент u функции u перпендикулярен ее линии уровня.

Если в области G R3xyz задано векторное поле a и существует

скалярное в G поле u(x, y, z), для которого векторное поле a является полем градиентов, т. е. a = u, то функция u называется потенциальной функцией (или потенциалом) векторного поля a.

Векторное поле, для которого существует потенциальная функция, называется потенциальным полем.

Если координаты вектора a обозначить P , Q, R, т. е. a = (P , Q, R), и если a = u, то

P =

∂u

,

Q =

∂u

,

R =

∂u

.

∂x

∂y

 

 

 

 

 

 

∂z

Поэтому потенциальность непрерывного векторного поля a означает, что выражение P dx + Q dy + R dz является полным дифференциалом:

du = P dx + Q dy + R dz.

В теории векторных полей имеются важные понятия дивергенции div a и вихря (или, что то же, ротора) rot a дифференцируемого

в области G векторного поля a = (P , Q, R). Определим эти понятия формулами с помощью символического вектора :

def

· a

(48.2)

div a =

§ 48. Скалярные и векторные поля

237

— как «скалярное произведение» векторов и a;

 

×

∂x

∂y

 

∂z

 

 

 

 

i

j

 

k

 

 

rot a = a

 

P

Q

 

R

 

(48.3)

 

 

 

def

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

— как «векторное произведение» этих же векторов.

В координатном виде эти определения записываются следующим образом:

 

 

 

 

 

 

def

 

∂P

 

∂Q

 

∂R

 

 

 

 

 

 

 

 

 

 

div a =

 

 

 

+

 

+

 

 

,

 

 

 

 

(48.4)

 

 

 

 

 

∂x

∂y

∂z

 

 

 

 

def

 

∂R

 

∂Q

 

∂P

 

∂R

 

 

∂Q

 

∂P

 

 

rot a =

 

 

 

 

i +

 

 

∂x j +

 

 

 

k.

(48.5)

∂y

∂z

∂z

∂x

∂y

Если G — плоская область и, следовательно, P = P (x, y), Q = = Q(x, y), R ≡ 0 (рассматриваемые векторные поля лежат в тех же плоскостях или пространствах, в которых множество G является областью), то

div a =

∂P

+

 

∂Q

,

(48.6)

∂x

∂y

 

 

 

 

 

 

rot a =

∂Q

∂P

k.

(48.7)

∂x

∂y

Определим еще понятия циркуляции и потока векторного поля a. Если Γ — кусочно-гладкий замкнутый контур и векторное поле a

задано на Γ, то криволинейный интеграл второго рода

a dr

(48.8)

Γ

называется также циркуляцией векторного поля a по этому контуру. Если кусочно-гладкая поверхность S ориентирована с помощью единичной нормали n, то для векторного поля a, заданного на по-

верхности S, поверхностный интеграл второго рода

a dS = an dS

(48.9)

S+ S

называется также потоком векторного поля через поверхность S.

Если Γ — кусочно-гладкая плоская кривая, a n — единичная нормаль к Γ, то для векторного поля a, заданного на Γ, интеграл

an ds

(48.10)

Γ

называется потоком векторного поля через контур Γ.

238 Гл. 5. Интегральное исчисление функций многих переменных

48.2.Формула Гаусса–Остроградского. Пусть G — область

впространстве R3xyz . Предположим, что на плоскости Rxy2 существует такая квадрируемая область D, ограниченная спрямляемой кри-

вой, что граница ∂G области G состоит из двух поверхностей S1 и S2, задаваемых явными представлениями соответственно z = ϕ(x, y) и z = ψ(x, y), где функции ϕ и ψ непрерывны на замыкании D области

D, ϕ(x, y) < ψ(x, y), (x, y) D, и, быть может, поверхности S0, являющейся частью цилиндрической поверхности, основанием которой

является граница ∂D области D, а образующая

параллельна оси Oz (см. п. 44.1):

 

∂G = S1 S2 S0.

(48.11)

В этом случае область G называется элементарной относительно оси Oz (рис. 50). Она

имеет вид

G = {(x, y, z) : (x, y) D, ϕ(x, y) < z < ψ(x, y)}.

(48.12)

Мы уже встречались с областями такого типа при изучении вопроса о сведении кратного интеграла к повторному. Обозначим для краткости границу ∂G обла-

сти G через S, тогда (см. (48.11))

S = S1 S2 S0.

(48.13)

Пусть на S задана функция F = F (x, y, z). Поверхностные интегралы

второго рода

 

 

F (x, y, z) dx dy от функции F по нижней стороне

 

 

 

 

 

 

 

 

 

 

 

S1

 

 

 

поверхности S1,

 

F (x, y, z) dx dy по верхней стороне поверхности

 

 

,

 

,

 

S2

 

 

S2 и

F (x

 

 

dx dy

по поверхности S0 (см. п. 47.3) называются

 

y

 

z)

 

S0

поверхностными интегралами второго рода по внешним сторонам этих поверхностей, а их сумма — интегралом по внешней стороне

поверхности S и обозначается F (x, y, z) dx dy, т. е.

S+

F (x, y, z) dx dy =

S+

= F (x, y, z) dx dy + F (x, y, z) dx dy + F (x, y, z) dx dy. (48.14)

S1

2

S0

S

 

§ 48. Скалярные и векторные поля

239

Аналогично определяется поверхностный интеграл

F (x, y, z) dx dy

S

no внутренней стороне поверхности S:

F (x, y, z) dx dy =

S

= F (x, y, z) dx dy + F (x, y, z) dx dy + F (x, y, z) dx dy. (48.15)

S1

S2

S0

Напомним (п. 47.3), что F (x, y, z) dx dy = 0 и, следовательно,

S0

это слагаемое можно было бы и не писать. Оно пишется для того, чтобы формулы (48.14) и (48.15) формально соответствовали формуле (48.13). Как будет видно из дальнейшего, это оказывается очень удобным.

Если поверхности S1, S2 и

S0 кусочно-гладкие,

то интеграл

F (x, y, z) dx dy представляет

собой поверхностный

интеграл по

S+

поверхности S, ориентированный с помощью внешней единичной нормали n. В этом случае, если

n = (cos α, cos β, cos γ),

(48.16)

то

 

 

F (x, y, z) dx dy =

F (x, y, z) cos γ dS.

(48.17)

S+

S

 

Здесь термин «кусочно-гладкая поверхность» понимается в смысле общего определения в п. 47.3.

Аналогично областям, элементарным относительно оси Oz, определяются области G, элементарные относительно осей Ox,

Oy, и интегралы

 

 

F dy dz,

 

F dy dz,

S+

S

(48.18)

F dz dx,

 

F dz dx

S+

S

 

по внешней и внутренней сторонам поверхности S.

Области, элементарные относительно всех координатных осей, называются элементарными (рис. 51). Примерами элементарных обла-

стей являются пирамиды, параллелепипеды, вообще любые выпуклые многогранники, шары, эллипсоиды и т. п.