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

матан

.pdf
Скачиваний:
41
Добавлен:
27.05.2015
Размер:
3.81 Mб
Скачать

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

571

так как Qa, Qn — непрерывная функция, то интеграл в правой части существует.

Справедливо равенство:

 

Qa, Qn dS =

(P cos(Qn, x)+Q cos(Qn, y)+R cos(Qn, z)) dS =

S

 

 

 

S

+ (u, v) + (u, v) |r u×r v| du dv =

=

P κ(u, v)

 

 

 

 

(y, z)

 

 

 

 

(z, x)

 

(x, y)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P (u, v)

+ Q(u, v)

+ R(u, v)

du dv. (5)

 

 

=

 

 

 

 

 

(y, z)

 

 

 

(z, x)

 

 

(x, y)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Здесь в правой части равенства стоит обычный двойной интеграл по области G, в которой функции P, Q и R зависят от x = x(u, v), y = y(u, v), z = z(u, v).

Рассмотрим частный случай потока. Предположим, что гладкий кусок S поверхности S взаимнооднозначно проектируется на каждую из координатных плоскостей x = 0, y = 0, z = 0, соответственно в области Sx, Sy, Sz, причем каждая из этих

областей измерима по Жордану. Таким образом, гладкий кусок S описывается любой из трех функций:

x = f1(y, z), (y, z) Sx,

y = f2(z, x),

(z, x) Sy,

z = f3(x, y),

(x, y) Sz.

Обозначим через Sx , Sy , Sz , соответствующим образом ори-

ентированные проекции, ориентированной поверхности S . Нормаль Qn к куску поверхности S образует угол с осью Oz, косинус которого равен:

 

cos(Qn, z) = ±

 

 

1

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

где

1 + p2 + q2

 

p =

∂f3

, q =

∂f3

 

 

∂x

 

∂y

 

 

 

 

 

 

 

 

 

и знак + или выбирается в зависимости от ориентации S . Тогда

R cos(Qn, z) ds =

S

572 Глава 20. Теория интегрирования в Rn

= R(x, y, f3(x, y))

±

 

 

 

 

 

 

 

 

1 + p2

+ q2

dx dy =

 

 

1 + p2 + q2

Sz

 

1

 

 

 

 

 

 

 

 

 

 

 

 

= ± R(x, y, f3(x, y)) dx dy = R(x, y, f3(x, y)) dx dy.

Sz

 

 

Sz

 

 

 

 

Здесь последний интеграл берется по ориентированной области.

ОПРЕДЕЛЕНИЕ 12.5. Введем обозначение

 

R(x, y, z) dx dy = R(x, y, f3(x, y)) dx dy.

S

Sz

Интеграл в левой части последнего равенства называется поверхностным интегралом II-го рода.

Аналогичные рассуждения справедливы и для остальных двух интегралов:

 

P cos(Qn, x) ds =

P (f1(y, z), y, z) dy dz =

S

 

Sx

=P (x, y, z) dy dz,

 

S

 

 

Q cos(Qn, y) ds =

Q(x, f2(x, z), z) dz dx =

S

 

Sy

 

=Q(x, y, z) dz dx.

S

Таким образом, поток вектора Qa через ориентированную поверхность S может быть вычислен по формуле:

Qa, Qn ds = P (x, y, z) dy dz + Q(x, y, z) dz dx+

S

S

+R(x, y, z) dx dy.

(6)

ЗАМЕЧАНИЕ. Если поверхность S может быть разрезана на конечное число частей S = kSk, каждая из которых проектируется однозначно на все три координатные плоскости, то для того, чтобы вычислить поток вектора Qa через поверхность S , нужно вычислить потоки Qa через все Sk и их сложить.

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

573

ЗАМЕЧАНИЕ. Связь поверхностных интегралов второго рода с поверхностными интегралами первого рода дана уже в определении.

(P cos α + Q cos β + R cos γ)ds =

S

=P (x, y, z) dy dz + Q(x, y, z) dx dz + R(x, y, z) dx dy,

S

где cos α, cos β, cos γ — направляющие косинусы нормали, направленной в соответствии с выбранной стороной поверхности (т.е. ориентацией).

12.5.Формула Гаусса-Остроградского

Пусть G R3 — ограниченная область с кусочно-гладкой

границей S, и на G определено векторное поле

Q Q Q

Qa(x, y, z) = P i + Qj + Rk.

Будем предполагать, что функции

P (x, y, z), Q(x, y, z), R(x, y, z), ∂P∂x (x, y, z), ∂Q∂y (x, y, z), ∂R∂z (x, y, z)

непрерывны в G.

ОПРЕДЕЛЕНИЕ 12.6. Дивергенцией векторного поля Qa называется скалярная функция

div Qa = ∂P∂x + ∂Q∂y + ∂R∂z .

Далее будем считать, что поверхность S ориентирована с помощью единичной внешней нормали (внешней относительно к G).

ОПРЕДЕЛЕНИЕ 12.7. Область G будем называть элементарной областью относительно Oz, или z-областью, если существуют такие функции ϕ(x, y) и ψ(x, y), определенные в области Λz Oxy, что

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

Заметим, что граница z-области может быть представлена следующим образом

∂G = S = S1 S2 S0,

574

Глава 20. Теория интегрирования в Rn

где S1 — представима в виде z = ψ(x, y), S2 — представима в виде z = ϕ(x, y), а S0 — поверхность, являющаяся частью цилиндра с основанием Λz и образующими, параллельными оси Oz.

Точно так же определяются элементарные области относительно осей Ox и Oy.

ОПРЕДЕЛЕНИЕ 12.8. Области, элементарные относительно всех координатных осей, будем называть элементарными.

ТЕОРЕМА 12.1 (Гаусса7-Остроградского). Пусть

G – элементарная область и на G заданы функции P (x, y, z), Q(x, y, z), R(x, y, z), непрерывные вместе со своими частными производными

 

∂P

(x, y, z),

∂Q

(x, y, z),

∂R

(x, y, z).

 

 

 

∂y

 

∂z

 

 

∂x

 

 

 

Тогда

div Qa dx dy dz = Qa, Qn ds,

 

 

(7)

 

G

 

 

S

 

 

 

где интеграл в правой части формулы берется по внешней стороне границы ∂G = S, (т.е. Qn – внешняя нормаль к S).

ЗАМЕЧАНИЕ. Формула (1) называется формулой Гаусса— Остроградского, или, в векторном анализе, теоремой о дивергенции.

ЗАМЕЧАНИЕ. Формулу Гаусса-Остроградского иначе можно записать в следующей форме

 

∂P

 

∂Q

 

∂R

 

 

(

 

+

 

 

+

 

)dxdydz =

 

P dydz + Qdzdx + Rdxdy. (8)

 

∂x

 

∂y

 

∂z

 

 

G

 

 

 

 

 

 

 

S

+

 

 

 

 

 

 

 

 

Доказательство. a) Покажем справедливость формулы (7) для случая, когда P = 0, Q = 0 и G является z-областью, т.е. существуют функции λ1(x, y), λ2(x, y), (x, y) Λz, такие, что

G = {(x, y, z) : (x, y) Λz, λ1(x, y) ≤ z ≤ λ2(x, y)}.

Пусть S+ — граница области G ориентированная с помощью внешней нормали. Тогда для области G имеет место равен-

7Гаусс Карл Фридрих (30.4.1777 – 23.2.1855) – немецкий математик, астроном, физик и геодезист. Род. в Брауншвейге (Германия). В его творчестве органично сочетались исследования по теоретической и прикладной математике. Работы Гаусса оказали большое влияние на дальнейшее развитие высшей алгебры, теории чисел, дифференциальной геометрии, теории притяжения, классической теории электричества и магнетизма, геодезии, теоретической астрономии.

∂R∂z dx dy dz =

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

575

ство:

 

 

∂z dx dy dz =

dx dy

λ2(x,y)

∂z (x, y, z) dz =

 

 

 

∂R

 

 

∂R

G

 

 

Λz

λ1(x,y)

 

 

=

R(x, y, λ2(x, y)) dx dy − R(x, y, λ1(x, y)) dx dy.

Λz

 

 

 

Λz

 

 

 

Введем обозначения: σ1+ — ориентированная поверхность графика функции z = λ1(x, y), σ2+ — ориентированная поверхность графика функции z = λ2(x, y), σ+ — ориентированная боковая поверхность, т.е. S+ = σ1+ σ2+ σ+. Обозначим через σ1+,z, σ2+,z — проекции кусков σ1+ и σ2+ на плоскость z = 0

ориентированные соответствующим образом (σ+ — ориенти-

рована отрицательно, σ+

 

 

 

 

1,z

— положительно). Пользуясь фор-

 

 

2,z

 

 

 

 

 

мулой предыдущего параграфа, получаем

R(x, y, λ2(x, y)) dx dy =

 

R(x, y, λ2(x, y)) dx dy =

 

 

 

 

 

 

+

 

 

 

 

 

 

 

Λz

 

 

 

 

σ

2,z

 

 

 

 

 

 

 

 

=

R(x, y, z) dx dy =

 

 

R(x, y, z) cos(Qn, z) ds.

 

+

 

 

 

 

 

 

 

 

 

 

σ2

 

 

 

σ2

 

Аналогично получаем

 

 

 

 

 

R(x, y, λ1(x, y)) cos(Qn, z) ds = σ

R(x, y, z) cos(Qn, z) ds.

Λz

 

 

 

 

 

1

 

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

 

∂R

 

 

 

 

 

 

dx dy dz =

 

 

 

 

 

∂z

 

 

G

 

 

 

 

 

= R(x, y, z) cos(Qn, z) ds + R(x, y, z) cos(Qn, z) ds.

σ1

 

 

 

 

 

σ2

 

Замечая, что cos(Qn, z) = 0 на боковой поверхности σ, имеем:

R(x, y, z) cos(Qn, z) ds+

G

σ1

576

 

 

Глава 20. Теория интегрирования в Rn

+

R(x, y, z) cos(Qn, z) ds + R(x, y, z) cos(Qn, z) ds =

σ2

 

 

σ

 

=

R(x, y, z) cos(Qn, z) ds.

 

 

S

b) Аналогичным образом можно для x-области и y-области

доказать формулы:

 

 

∂x dx dy dz =

P (x, y, z) cos(Qn, x) ds,

 

∂P

 

G

 

S

Q(x, y, z) cos(Qn, y) ds.

 

∂y dx dy dz =

 

∂Q

 

G

 

S

 

Если область G является одновременно x-областью, y-областью и z-областью, то складывая полученные формулы, получаем

формулу Гаусса-Остроградского:

 

 

∂x +

∂y

+ ∂z

dx dy dz =

 

 

∂P

∂Q

∂R

 

G

=(P cos(Qn, x) + Q cos(Qn, y) + R cos(Qn, z)) ds.

S

Можно доказать, что формула Гаусса-Остроградского справедлива для любой ограниченной области, граница которой состоит из конечного числа кусочно-гладких поверхностей. В данном случае мы ограничимся лишь формулировкой.

ТЕОРЕМА 12.2. Пусть

граница

∂G

ограниченной

области G состоит из конечного числа кусочно-

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

на

G

 

заданы функции

P (x, y, z), Q(x, y, z), R(x, y, z),

непрерывные

вместе со

своими частными производными

 

 

 

 

∂P

(x, y, z),

∂Q

(x, y, z),

∂R

(x, y, z).

 

 

 

∂z

 

∂x

∂y

 

 

Тогда

div Qa dx dy dz =

Qa, Qn ds,

 

 

G

 

 

∂G

 

 

где n – внешняя нормаль к ∂G.

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

577

ЗАМЕЧАНИЕ. Формула Гаусса-Остроградского позволяет найти выражение для объема области через соответствующий поверхностный интеграл. Положим в формуле (8) P =

x, Q = y, R = z. Тогда

3|G| = 3 dx dy dz =

G

(x cos(Qn, x) + y cos(Qn, y) + z cos(Qn, z)) ds =

=xdydz + ydzdx + zdxdy.

S+

ПРИМЕР 3. Пусть G = {x2 + y2 + z2 ≤ R2}. Тогда

4

 

1

 

x2

y2

z2

 

R

 

 

πR3 =

 

 

 

 

+

 

+

 

ds =

 

 

|S|.

3

3

R

R

R

3

S

Следовательно, |S| = 4πR2.

12.6.Геометрический смысл дивергенции

ТЕОРЕМА 12.3. Пусть в области G определено непрерывно дифференцируемое вектороное поле Qa = Qa(M). Пусть M0 G и D — область с кусочно-гладкой границей

∂D = S такая, что M0 D, D G. Обозначим через S+

— поверхность S, ориентированную с помощью внешней нормали, а через d(D) – диаметр D. Тогда

div Qa(M0) = lim

S

Q

.

(9)

 

 

Qa, ds

 

 

 

 

+

 

 

 

 

d(D)0

 

 

 

 

 

 

|D|

 

 

 

Доказательство. По формуле Гаусса–Остроградского полу-

чаем

 

 

 

 

 

div Qa dx dy dz =

Qa, Qn ds =

 

Q

 

 

Qa, ds .

 

D

S

S

+

 

 

С другой стороны, по интегральной теореме о среднем получаем

div Qa dx dy dz = div Qa(M)|D|,

D

578 Глава 20. Теория интегрирования в Rn

где M – некоторая точка, принадлежащая D. Таким образом

div Qa(M) =

S

.

 

 

Qa, Qn

dS

 

+

 

 

 

 

|D|

 

 

Переходя к пределу при d(D) 0, получаем требуемое.

ЗАМЕЧАНИЕ. В правой части формулы (9) стоят величины, не зависящие от выбора системы координат. В результате получаем, что дивергенция векторного поля не зависит от выбора системы координат.

ЗАМЕЧАНИЕ. Точки векторного поля Qa в которых div Qa = 0 называются источниками векторного поля. Как видно из последней формулы, поток такого векторного поля через любую достаточно малую поверхность, окружающую источник, не равен нулю.

ПРИМЕР 4. Пусть u = f(x, y, z) — дважды непрерывно

дифференцируемая функция в области G, где S = ∂G — кусочно-гладкая поверхность. Пусть u = (fx, fy, fz). Тогда по формуле Гаусса—Остроградского поток u через S вычис-

ляется по формуле

 

 

 

 

 

 

 

 

 

 

 

 

 

∂n ds = u, Qn ds = div u dx dy dz =

 

∂u

 

 

 

 

 

 

 

 

 

 

 

 

S

 

S

∂x2

+ ∂y2

 

 

G

 

=

 

+ ∂z2 dx dy dz.

 

 

 

2f

 

2f

 

 

2f

 

 

G

 

 

 

 

 

 

 

 

 

 

 

 

Отметим, что выражение

 

 

 

 

 

 

 

 

 

 

 

2f

2f

2f

 

 

 

 

+

 

+

 

 

= ∆f

 

 

 

∂x2

∂y2

∂z2

называется оператором Лапласа8 от функции f(x, y, z).

12.7.Формула Стокса

Пусть S R3 — дважды непрерывно дифференцируемая поверхность, а Qr = Qr(u, v) — ее представление. Будем считать, что (u, v) D, где D — плоская ограниченная область

8Лаплас Пьер Симон (23.3.1749 – 5.3.1827) – французский математик, физик и астроном. Род.

вНормандии (Франция). Ему принадлежат многочисленные фундаментальные работы по математике, экспериментальной и математической физике и небесной механике. В области математики Лаплас получил значительные результаты по теории дифференциальных уравнений, теории вероятностей, теории потенциала и др.

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

579

для которой справедлива формула Грина. Пусть ∂D — один простой кусочно-гладкий контур, Γ0 — положительно ориен-

тированный контур ∂D, и

u = u(t);

v = v(t); a ≤ t ≤ b

— представление Γ0. Зададим через

Qn = Qr u × Qr v |Qr u × Qr v|

— ориентацию на S. Таким образом Qn = (cos α, cos β, cos γ). Также обозначим S — поверхность S с выбранной на ней нормалью Qn, а Γ — контур с представлением Qr = Qr(u(t), v(t)),

где t [a, b].

ОПРЕДЕЛЕНИЕ 12.9. При введенных выше обозначениях говорят, что контур Γ ограничивает поверхность S, или, что поверхность S натянута на Γ.

Пусть в области G R3 такой, что S G, задано векторное поле

 

 

 

 

Q

 

 

 

 

 

Q

 

 

 

Q

 

Qa(x, y, z) = P (x, y, z)i + Q(x, y, z)j + R(x, y, z)k.

 

Будем предполагать, что функции P, Q, R непрерывно диф-

ференцируемы. Напомним, что

 

 

 

 

 

 

 

 

 

 

 

∂R

∂Q

∂P

 

 

∂R

 

 

∂Q

∂P

 

rotQa =

 

 

Qi +

 

 

 

 

 

Qj +

 

 

 

 

Qk =

∂y

∂z

∂z

∂x

∂x

∂y

 

 

 

 

 

 

Q

 

Q

 

Q

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

j

 

k

 

 

 

 

 

 

 

 

 

 

 

P

Q

 

R

 

 

 

 

 

 

 

 

 

= det

 

 

 

.

 

 

 

 

 

 

 

 

 

 

∂x

∂y

 

∂z

 

 

 

 

 

 

Справедливо следующее утверждение.

ТЕОРЕМА 12.4 (формула Стокса9). Пусть функции

P, Q, R непрерывно дифференцируемы в G. Тогда справедлива формула Стокса:

rota,Qn ds =

Q

Qa, dr ,

S

Γ

т.е. поток вектора rotQa через ориентированную поверхность S равен циркуляции вектора Qa вдоль границы Γ этой поверхности, ориентированной соответственно ориентации S .

9Стокс Джорж Габриэль (13.8.1819 – 1.2.1903). Род. в Скрине (Ирландия). Работал преимущественно в области оптики, гидродинамики и математической физики. Особо следует отметить вывод Стоксом одной из важнейших формул в векторном анализе.

580 Глава 20. Теория интегрирования в Rn

ЗАМЕЧАНИЕ. Другими словами, справедливо

 

 

 

 

 

∂y

∂z

dy dz + ∂z

∂x

dz dx+

 

 

 

 

 

 

 

 

∂R

 

∂Q

 

∂P

∂R

 

 

 

S

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

∂Q

∂P

dx dy = P dx + Q dy + R dz,

 

 

 

 

 

 

 

 

 

 

 

 

∂x

 

∂y

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Γ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂R

 

 

 

∂Q

 

 

∂P

 

 

∂R

 

 

 

 

 

 

 

 

 

 

cos(Qn, x) +

 

 

cos(Qn, y)+

 

 

∂y

∂z

 

 

∂z

∂x

 

S

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂Q

 

 

∂P

 

cos(Qn, z) ds =

 

 

 

 

 

 

 

 

 

 

 

 

 

P dx + Q dy + R dz.

 

∂x

∂y

 

Γ

Отметим, что в случае, если S — область в плоскости (x, y) и Qn = (0, 0, 1), тогда формула Стокса принимает вид:

 

∂x

∂y

dx dy = P dx + Q dy,

 

 

∂Q

∂P

 

S

 

 

 

 

Γ

т.е. получаем формулу Грина. Таким образом, формула Стокса это обобщение формулы Грина на случай кривых поверхностей, лежащих в пространстве R3.

Доказательство. Справедливы следующие равенства

P (x, y, z) dx =

Γ

b

=P [x(u(t), v(t)), y(u(t), v(t)), z(u(t), v(t))]xt(u(t), v(t))dt =

a

=

P [x(u, v), y(u, v), z(u, v)]

∂x

(u, v)du +

∂x

(u, v)dv =

 

 

∂u

∂v

Γ0

=

 

 

 

 

 

 

 

 

 

 

 

P

∂x

(u, v)du + P

∂x

(u, v)dv.

 

 

 

 

 

 

 

∂u

∂v

 

 

Γ0

 

 

 

 

 

 

 

 

 

 

Применяя к правой части данного равенства формулу Грина, получаем

P ∂u(u, v)du + P

∂v (u, v)dv =

∂u

P ∂v

∂v

P ∂u

dudv =

 

∂x

∂x

 

 

∂x

 

∂x

 

Γ0

D