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

kr_kriv_pov_int

.pdf
Скачиваний:
12
Добавлен:
09.02.2015
Размер:
675.9 Кб
Скачать

не зависит от способа параметризации (согласованного с `+) кривой `+. Действительно, пусть заданы две параметризации `+

~r = ~r1(t) = [x(t), y(t)]T , t [t(1),

t(2)],

t(1) < t(2);

(3.13)

r = ~r2(t) = [¯x(τ), y¯(τ)]T ,

τ [τ(1),

τ(2)],

τ(1) < τ(2),

(3.14)

согласованные с направлением

кривой, и

t

=

h(τ), t(1) =

h(τ(1)),

t(2) = h(τ(2)).

 

 

 

 

 

Пусть также h(τ) монотонная дифференцируемая функция.

Òàê êàê t(1) < t(2) è τ(1) <

τ(2), то в этом случае h0(τ) ≥ 0. Êàê

и в раз. 3.2, замена переменных

t = h(τ) в интеграле (3.11) приводит к

равенству:

 

t(2)

 

Z

[fx(x(t), y(t))x0(t) + fy(x(t), y(t))y0(t)]dt =

t(1)

τ(2)

Z

=[fx x(h(τ)), y(h(τ)) x0t(h(τ)) + fy x(h(τ)), y(h(τ)) yt0(h(τ))]h0(τ)dτ =

τ(1)

τ(2)

Z

=fx x¯(τ), y¯(τ) x¯0(τ) + fy x¯(τ), y¯(τ) y¯0(τ) dτ,

τ(1)

которое и означает независимость интеграла (3.12) от способа параметризации.

Если же параметризация (3.13) согласована с `+, а параметризация (3.14) c `, то точке A соответствуют значения параметров t(1) è τ(2), à точке B значения t(2) è τ(1). В этом случае для функции h(τ) справедливы

равенства: t(1) = h(τ(2)), t(2) = h(τ(1)) è

t(2)

Z

fx(x(t), y(t))x0(t) + fy(x(t), y(t))y0(t) dt =

t(1)

τ(2)

Z

=fx(¯x(τ), y¯(τ))¯x0(τ) + fy(¯x(τ), y¯(τ))¯y0(τ) dτ =

 

 

τ(1)

 

 

 

 

 

 

 

τ(2)

 

 

 

 

 

Z

 

Z

 

 

 

 

 

 

= −

 

x(¯( ) ¯( ))¯0

( ) +

y(¯( ) ¯( ))¯0( )

 

= −

~

→− )

(f(P ),

τ

 

 

 

 

`

d` ,

 

f x τ , y τ x τ

f x τ , y τ y τ dτ

 

 

 

(1)

 

 

 

 

 

 

 

что соответствует соотношению (3.9).

57

на две части: ограниченную

3.4. Теорема Грина

Кривая ` в определении криволинейных интегралов может быть замкнутой. В этом случае точки A и B (концы `) совпадают, и в качестве точки A можно вообще взять любую точку `. Криволинейный интеграл второго рода по замкнутой кривой `+ принято обозначать символом

I (f(P ), →− )

~

d` .

 

`+

 

Замкнутая кривая ` разбивает плоскость R2

часть D и неограниченную R2 \D. Будем далее считать, что направление на `+ согласовано с D так, что при движении по ` в выбранном направлении часть D плоскости остается слева.

Криволинейный интеграл второго рода (в R2) по замкнутой кривой ` тесно связан с двойным интегралом. Справедливо следующее утверждение.

Теорема 3.11. (Теорема Грина). Пусть D ограниченная замкну-

тая область в R2

с кусочно-гладкой границей `, и на ` выбрано направ-

ление, согласованное с D. Пусть функция

~

R

2

R

2

непрерывно

f(P ) :

 

 

 

 

 

 

 

 

дифференцируема на области D. Тогда справедливо равенство:

 

 

 

 

 

 

ZZ ∂x

∂y

 

 

= I

 

f(P ), →−

 

 

 

 

 

 

 

 

 

 

 

 

∂fy

 

∂fx

 

 

~

 

 

 

 

 

 

 

 

 

 

 

 

D

 

 

 

 

 

 

∂D

 

 

 

 

 

 

 

 

 

 

 

ãäå ~

 

 

 

 

 

 

 

dS

 

 

 

d` ,

 

 

2

 

(3.15)

~

~;

(x, y)

декартовы координаты в

,

~ ~

}

 

f = fx(x, y)i + fy(x, y)j

 

 

 

 

 

 

 

 

 

 

 

R

{i, j

 

соответствующий этим координатам базис, а символом ∂D обозначена ориентированная граница `+ области D.

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

Правильной будем называть область, являющуюся одновременно правильной относительно осей Ox и Oy (см. раз. 1.4). Правильной будет, на-

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

Докажем сначала теорему Грина для случая, когда D правильная

относительно оси

Ox

область, и

~

~(1)

 

 

 

 

(1)

(1)

~

ãäå

 

 

 

 

 

 

f = f

 

(P ) = fx (x, y) + fy

(x, y)i,

 

fx(1) = fx, fy(1) = 0. Итак, пусть

 

 

 

 

 

 

 

 

 

 

 

 

D

 

(

x, y

 

a

 

x

 

b, g

 

 

x

 

y

g x) .

 

 

 

 

=

 

) :

 

 

 

1(

 

) ≤

 

(1) 2( (2)

(3)

(4)

 

На рис. 3.1 изображена область

D и части `+ , `+ , `+ , `+ , íà êî-

58

торые в этом случае разбивается граница ` области D. Направления на `(1)+ ,

`(2)+ , `(3)+ , `(4)+ , согласованные с D, указаны стрелками. По теореме 1.12 справедливо равенство:

I1 ZZ

∂x

∂y

! dS =

 

∂fy(1)

∂fx(1)

 

y. .

. .

`(4)+

.

y = g2

(x)

 

.

 

.

 

`+(3)

 

.

D

 

. .`+(2)

. `+(1)

 

.

 

 

 

 

.

.

 

 

D

 

∂fx dy dx.

.

y = g1(x) .

 

 

 

.

. .

=

∂fx(1) dS =

.....

O a

b

x

 

b

 

g2(x)

 

 

.

 

 

 

 

 

ZZ

 

 

 

 

 

 

 

Z

 

 

 

 

Z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ðèñ. 3.1

 

 

 

 

 

 

∂y

 

 

 

 

 

 

 

 

 

 

∂y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D

 

 

 

 

 

 

 

 

a

 

 

g1(x)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Вычисляя внутренний интеграл, получим:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

"fx(x, y) y=g1

 

 

 

 

 

 

 

 

 

b

 

 

fx

x, g2(x)

− fx

x, g1(x)

 

dx.

I1 = − Z

(x)# dx = − Z

 

 

 

 

 

 

 

 

 

 

 

 

y=g2

(x)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~(1)

d` , который совпадает с суммой таких

Вычислим теперь интеграл

 

(f

, →− )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

`+

(2)

(3)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1)

 

 

(4)

 

 

 

 

 

 

 

 

 

(1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

H

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

же интегралов по кривым `+ , `+ , `+ , `+ . Для кривой `+

возможна такая

параметризация:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

`+(1) : ~r = ~r1(t) = [t, g1(t)]T ,

 

t [a,

b].

 

 

 

 

 

 

Поэтому по формуле (3.11)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Z

 

 

 

 

 

 

 

b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

 

 

 

 

,

→− ) =

 

 

Z

h

 

(1)

(

 

 

 

( ))

0 +

 

 

(1)

(

 

 

 

( )) 0

i

 

 

=

Z

(

 

 

( ))

`+

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

(f

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~(1)

 

 

d`

 

 

 

f

 

 

t, g

1

 

t t

 

 

f

 

 

 

t, g

1

t g

1

 

dt

 

 

 

f

 

t, g

1

t dt.

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

(1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Для кривой `(3)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(3.16)

в качестве параметризации возьмем:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

`(3) : ~r = ~r3(t) = [t, g2(t)]T ,

t [a, b].

 

 

 

 

 

 

Тогда аналогично (3.16) получим равенство:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Z

 

 

 

 

 

 

 

 

 

 

 

Z

 

 

 

 

 

 

 

 

 

 

 

b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~(1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Za

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(

f

 

, d`

 

 

 

 

 

 

 

 

 

d`

 

 

 

 

 

 

f

 

 

t, g

2

t

 

dt.

 

 

 

(3.17)

 

 

 

 

 

 

 

 

→− ) = −

 

(f , →− ) = −

x(

 

 

(

))

 

 

 

 

 

 

 

 

(3)

 

 

 

 

 

 

 

 

 

 

 

(3)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

`+

 

 

 

 

 

 

 

 

 

 

`

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

59

Для кривой `+(2)

выбираем параметризацию:

 

 

 

 

 

 

 

 

 

Тогда

 

 

 

 

`+(2) : ~r = ~r2(t) = [b, t]T ,

 

 

t [g1(b), g2(b)].

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

g2(b)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

`+

 

 

,

 

 

 

 

 

g1(b) h

x

 

( ) 0

+ y

 

( ) 0

= 0

 

 

 

 

 

(f

→− ) =

 

 

 

 

 

 

 

 

 

 

 

~(1)

 

d`

 

 

 

 

 

 

f

(1)

 

 

b, t b

 

 

f

(1)

b, t t

dt

 

 

.

(3.18)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Точно так же

 

 

 

 

 

 

 

 

 

 

Z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~(1)

d`

 

 

 

 

.

 

 

 

 

 

 

 

(3.19)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(f

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

, →− ) = 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(4)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

`+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Складывая равенства (3.16) (3.19), получим

 

 

 

 

 

 

 

 

Z (f

 

 

 

 

 

 

 

b

[fx(t, g1(t)) − fx(t, g2(t))] dt = I1,

 

 

 

 

, →− ) = Za

 

 

 

 

 

~(1)

 

d`

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(3.20)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

`+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

что доказывает теорему в рассматриваемом случае.

 

 

 

 

 

Пусть теперь

D

 

область,

 

 

правильная

относительно оси

Oy è

~ ~(2)

 

 

 

 

 

 

~

. Пример подобной

 

 

 

 

 

 

 

 

 

 

f = f

(P ) = fy(x, y)j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

области изображен на рис. 3.2. В этом слу-

 

 

y. .

 

`+(4).

 

 

 

÷àå

 

 

 

 

 

 

 

 

 

 

 

 

 

! dS =

 

 

 

 

d

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

. . .... ..

 

I2 ZZ

 

∂x

 

 

∂y

 

 

 

 

 

 

 

 

 

.

`+. .

 

 

 

 

 

 

 

 

 

 

 

x = h1(y)

 

 

 

 

 

 

∂fy(2)

 

∂fx(2)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

. .

 

 

 

(3)

 

 

 

D

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1)

D

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂fy

 

 

d

 

 

h2(y)

 

 

 

 

 

 

 

 

 

 

 

c

`+

.

 

x = h2(y)

ZZ

 

Z

 

Z

 

∂fy

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

... .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

=

 

 

dS =

 

 

 

 

 

 

 

 

 

dx

dy =

 

 

 

 

 

. . . ..... . .

 

 

 

 

∂x

 

 

 

 

 

 

 

∂x

 

 

 

 

 

 

 

`+(2)

 

 

.

 

D

 

 

 

 

c

 

 

h1(y)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

O

 

 

 

 

 

x

 

d

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= Z

 

fy(h2(y), y) − fy(h1(y), y) dy.

 

 

 

 

 

 

 

Ðèñ. 3.2

 

 

 

c

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Для кривых `+(2)

è `+(2), как и в предыдущем случае, получим равенства:

 

 

 

 

 

 

 

 

Z

(

 

 

→− ) = (

 

 

 

 

→− ) = 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~(2)

, d`

 

 

 

 

Z

 

~(2)

, d`

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f

 

 

 

 

 

f

 

 

 

 

 

 

 

 

 

 

 

 

 

(2)

 

 

 

 

 

 

 

 

 

 

(4)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

`+

 

 

 

 

 

 

 

 

 

`+

 

 

 

 

 

 

 

 

 

 

 

 

 

Äëÿ `(1)è `(4)+ параметризациями в рассматриваемом случае будут:

`(1): ~r = ~r1(t) = [h1(t), t]T , t [c, d];

60

Поэтому

 

 

 

`+(3) : ~r = ~r3(t) = [h2(t), t]T ,

t [c, d].

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Z

 

 

 

Z

 

 

d

 

 

 

 

 

 

 

 

 

 

 

,

→− ) =

 

~(2)

→− ) =

Z

 

 

(2)

( ( ) ) 0 +

 

(2)

( ( ) ) 0 =

(f

 

(

 

 

 

 

~(2)

 

 

`

 

c

 

 

 

 

 

 

 

 

 

`+

 

d`

f , d`

f

 

h

 

t , t h

 

f

 

 

 

 

 

 

 

 

 

 

x

 

1

 

1

 

y

1

 

(1)

 

 

 

(1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d

y

1( )

 

;

Z (f

, →− )

= − Z

f

h

t , t dt

 

~(2)

d`

 

 

d

Z

=fy h2(t), t dt.

c

(3)

c

 

`+

 

Складывая полученные равенства, получим и в этом случае соотношение:

Z (f , →− ) =

ZZ

 

∂x

∂y

!

 

 

~(2)

 

 

∂fy(2)

∂fx(2)

 

 

`+

d`

D

 

 

 

 

 

 

dS.

(3.21)

 

 

 

 

 

 

 

 

 

 

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

формулы (3.20) и (3.21). Так как ~

 

 

~(1)

 

 

~(2)

(P ), то складывая

 

 

f(P ) = f

(P ) + f

равенства (3.20) и (3.21) получим:

∂x

∂y

 

 

 

Z (f,

→− ) =

ZZ

 

 

 

~

d`

 

∂fy

 

∂fx

dS.

 

(3.22)

D

 

 

 

 

 

 

 

`+

 

 

 

 

 

 

 

 

 

 

Таким образом, теорема Грина для правильной области доказана. Наконец, пусть область D представляет собой объединение конеч-

ного числа правильных областей. По-

 

 

 

 

 

 

 

кажем, что теорема Грина справедлива

y. .

(1)

 

 

 

 

 

и для такой области D. Рассмотрим

 

.

.

 

 

 

 

. .

 

 

простейший случай,

когда D является

 

`+

 

 

 

 

 

 

 

 

`+

 

D2

 

объединением двух

правильных

îáëà-

 

 

 

 

 

 

 

D1

 

 

.

 

 

 

(1)

 

 

 

 

 

 

 

 

 

 

 

 

 

.

ñòåé D1 è D2 (рис. 3.3). Обозначим `

 

 

 

. .

 

 

.

ту часть границы D, которая входит

 

 

 

 

 

(2)

 

и в границу D1, à

`

(2) оставшуюся

 

 

 

 

 

`+

.

 

 

 

 

 

 

 

 

.

 

 

 

 

...

часть границы D (являющуюся частью

O

 

 

 

 

 

x

границы D2). Общую границу

D1 è

 

 

 

 

 

 

 

D2 обозначим `. Введем на границе D

 

 

Ðèñ. 3.3

 

 

направ-

 

 

 

 

 

 

 

 

 

 

ление, согласованное с D, а на кривой ` произвольное направление.

61

Тогда из (3.22) для правильных областей D1 è D2 получим:

ZZ

∂xy

∂yx

dS = ZZ

∂xy

∂yx

dS + ZZ

∂xy

∂yx dS =

 

∂f

 

∂f

 

∂f

 

∂f

 

 

∂f

 

∂f

D

 

 

 

 

D1

 

 

 

 

D2

 

 

 

 

 

~→− ~ →−

=(f, d`) + (f, d`).I I

 

 

∂D1

 

∂D2

 

 

С другой стороны,

 

Z

 

 

Z

 

I

 

 

 

 

 

~

 

~

 

~

d`

(

f, d`

f, d`

 

 

→− ) = (

 

→− ) + (f,

→− );

∂D1

 

`+

 

 

(1)

 

 

 

 

 

 

`+

 

(3.23)

(3.24)

(

 

→− ) = (

 

→− ) + (

 

→− ) = (

 

→− ) − (f, →− )

 

Z

 

Z

 

Z

 

Z

 

 

Z

 

 

 

~

 

~

 

~

 

~

 

~

d` .

(3.25)

 

f, d`

f, d`

f, d`

f, d`

 

∂D2

 

(2)

 

(1)

 

(2)

 

 

(1)

 

 

 

 

`+

 

`

 

`+

 

 

`+

 

 

Складывая равенства (3.24) и (3.25), получим:

(

 

→− ) +

(

 

→− ) = (

 

→− ) + (

 

→− ) = (f, →− )

I

 

 

I

 

Z

 

Z

 

 

I

 

 

~

 

 

~

 

~

 

~

 

~

d` ,

 

f, d`

 

f, d`

f, d`

f, d`

 

∂D1

 

 

∂D2

 

(1)

 

(2)

 

 

∂D

 

 

 

 

 

 

`+

 

`+

 

 

 

 

и теперь из равенства (3.23) следует теорема Грина для области D. Аналогично рассматривается случай, когда D является объединением трех или

большего числа правильных областей. Замечание 3.4. Теорема Грина и формула Грина (3.15) справедливы

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

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

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

4.1. Поверхностный интеграл первого рода

Пусть в пространстве R3 задана некоторая гладкая поверхность Σ. Напомним, что под гладкой поверхностью Σ мы понимаем множество тех

62

(и только тех) точек P в R3, радиус-векторы ~r которых удовлетворяют соотношению:

~r = ~r(u, v) = [gx(u, v), gy(u, v), gz(u, v)]T , (u, v) D,

(4.1)

где D ограниченная замкнутая область в R2, (u, v) координаты точек из D, à gx, gy, gz непрерывно дифференцируемые функции. Дополнительно предполагается, что для всех точек D

 

 

 

 

 

 

 

~r 0

, ~r

0

~

,

 

 

 

 

 

 

 

 

(4.2)

ãäå

 

 

 

 

 

 

[

u

 

v

] 6= 0

 

 

 

 

 

 

 

 

 

 

∂ux

,

 

∂uy

,

∂uz

 

, ~r v0 =

∂vx

,

∂vy

,

 

∂vz

 

 

~r u0 =

 

 

 

 

 

∂g

 

∂g

 

∂g

 

T

 

 

 

 

∂g

 

∂g

 

∂g

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

частные производные функции ~r(u, v) по u и v соответственно.

Для гладкой поверхности Σ определено понятие площади, которую сейчас будем обозначать |Σ|. Согласно формуле (ref1.78)

|Σ| = ZZ k[~r u0 , ~r v0 ]k dS.

(4.3)

D

 

Пусть дана также некоторая функция f(P ), определенная, по край-

ней мере, для всех точек Σ. Разобьем поверхность Σ на части

σi, i =

= 1, 2, . . . , n, таким образом, что для каждой из этих частей определена

площадь | σi|. Будем считать, что каждой части

 

σi соответствует такая

часть Di области D, что

Pi σi и вычислим

 

 

σi = (x, y, z) S : (u, v)

Di .

Выберем произвольно точки

 

 

сумму

 

 

 

 

n

 

 

 

Xi

 

 

 

f(Pi )| σi|,

 

(4.4)

 

=1

 

 

называемую интегральной суммой для функции f(P ), поверхности Σ, выбранного разбиения поверхности Σ и выбранных точек Pi .

Рангом заданного разбиения назовем число λ = max diam(Δσi).

1≤i≤n

Определение 4.1. Число I называется поверхностным интегралом первого рода функции f(P ) по поверхности Σ, если для любого ε > 0

существует такое δ > 0, что для любого разбиения с рангом λ, удовлетворяющим неравенству λ < δ, и любого выбора точек Pi справедливо

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

n

X

i=1

f(Pi )| σi| − I < ε.

63

Функция f(P ) называется в этом случае интегрируемой по Σ. Для поверхностного интеграла первого рода используется обозначение:

ZZ

I = f(P )dσ.

Σ

Замечание 4.1. Поверхностный интеграл первого рода обозначается так же, как и двойной интеграл. Обычно из контекста бывает ясно, о каком из них идет речь. Использование же для этих интегралов одного обозначе- ния обусловлено тем, что поверхностный интеграл первого рода является обобщением понятия двойного интеграла и сводится к нему в том частном случае, когда gx(u, v) = u, gy(u, v) = v, gz(u, v) ≡ 0; в этом случае x = u,

y = v и Σ часть плоскости z = 0, совпадающая с D.

Свойства поверхностного интеграла первого рода аналогичны свойствам как криволинейного интеграла первого рода, так и двойного интеграла. В частности, для поверхностного интеграла первого рода справедливы теоремы, аналогичные теоремам 1.1 1.9.

Для вычисления поверхностного интеграла отметим, что в соответствии с (4.3)

| σi| = ZZ

 

~r u0

, ~r v0

 

dS.

 

 

 

 

 

 

 

 

 

 

Di

Так как функции ~r 0u è ~r 0v непрерывны, то по теореме о среднем для двойного интеграла

 

 

σ

 

 

 

~r 0

u , v , ~r 0

u , v

 

 

D ,

значения пара-

 

(˜u , v˜ ) некоторая

 

 

 

D

 

v

i ˜i)

|

u ,

v

 

|

 

i| =

 

ui ˜i)

 

 

i|

 

 

i i

 

 

 

 

 

 

i

 

 

 

i

i

 

ãäå

 

 

 

точка из

 

 

. Обозначим

 

 

 

метров u и v, соответствующие точке Pi , т. е. пусть f(Pi ) = f(xi , yi , zi ), xi = gx(ui , vi ), yi = gy(ui , vi ), zi = gz(ui , vi ).

Таким образом, интегральная сумма (4.4) может быть записана в виде:

 

n

 

 

 

 

 

 

 

 

 

 

 

Xi

 

 

 

~r 0 u , v , ~r 0 u , v

 

 

 

D ,

 

f x , y , z

 

 

 

 

 

( i i

i )

 

 

˜i)

|

i|

 

=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ò.

е. является обобщенной

интегральной

 

суммой для функции

 

λ

 

 

 

 

 

 

Σ

 

f

gx(u, v), gy(u, v), gz(u, v)

 

~r u0

(u, v), ~r v0 (u, v)

и области D. Можно так-

è ðàíã µ = max diam(ΔDi) соответствующего разбиения области D будет

же доказать, что если ранг

 

разбиения поверхности

 

достаточно мал, то

1≤i≤n

также малым. Поэтому справедлива следующая теорема.

Теорема 4.1. Если Σ гладкая поверхность, задаваемая соотношениями (4.1), а функция f(P ) непрерывна в замкнутой области Ω, содер-

64

жащей Σ, то

ZZ f(P )dσ =

ZZ f gx(u, v), gy(u, v), gz(u, v)

~r

 

 

 

0 (u, v), ~r 0 (u, v) dS.

u v

Σ

D

(4.5)

Остановимся подробнее на важном частном случае, когда поверхность Σ является частью графика дифференцируемой функции g(x, y) декарто-

вых координат x, y. В этом случае в качестве параметрического задания Σ можно принять

~r = ~r(x, y) = [x, y, g(x, y)]T , (x, y) D.

Тогда

 

 

 

 

 

 

 

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂g

 

 

 

 

 

 

 

 

 

 

 

 

∂g

 

 

 

 

 

 

 

 

~r u0 = ~r x0

= 1, 0,

 

 

 

; ~r v0 = ~r y0 = 0, 1,

 

 

 

 

 

;

 

 

 

 

 

 

 

 

∂x

 

∂y

 

 

 

 

 

 

 

~r u0

, ~r v0

 

= −∂x, −∂y

, 1 ;

 

 

~r u0

, ~r v0

 

= s

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 +

 

 

∂x

 

 

+

∂y .

 

 

 

 

∂g

 

T

 

 

 

 

 

 

 

 

 

2

 

2

 

 

 

 

 

∂g

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂g

 

 

 

 

∂g

 

 

Поэтому в таком случае

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ZZ f(P )dσ = ZZ f x, y, g(x, y) s

 

 

 

 

 

 

 

 

 

 

(4.6)

1 +

∂x

+

∂y

 

dx dy.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂g

2

 

∂g

 

2

 

 

 

 

 

 

Σ

 

 

 

D

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример. Вычислим площадь поверхности сферы SR радиуса R, применяя формулу (4.5). Возьмем параметрические уравнения сферы в виде:

~r = ~r(ϕ, ϑ) = [R cos ϕ sin ϑ, R sin ϕ sin ϑ, R cos ϑ]T ,

(ϕ, ϑ) D = {(ϕ, ϑ) : 0 ≤ ϕ ≤ 2π, 0 ≤ ϑ ≤ π} .

Тогда

 

~r 0

 

 

R

 

ϕ

 

ϑ, R

 

 

ϕ

 

ϑ,

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ϕ

= [− sin

 

sin

 

 

cos

 

sin

 

 

0] ;

 

 

 

 

 

 

~r ϑ0 = [R cos ϕ cos ϑ, R sin ϕ cos ϑ, −R sin ϑ]T ;

 

T

 

 

~r 0

, ~r 0

R2

 

ϕ sin2 ϑ, R2 sin ϕ sin2 ϑ,

 

 

R2 sin ϑ cos ϑ

;

 

ϕ

ϑ = −

 

cos

~r ϕ0

, ~r ϑ0

 

= R2 sin ϑ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

π

 

 

 

 

|SR| = ZZ

dσ = ZZ

2

 

 

 

 

 

2

Z0

Z0

 

 

 

2

 

R

sin ϑdϕdϑ = R

 

sin ϑdϑ = 4πR

.

 

SR

 

 

D

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

65

4.2.Двухсторонние и односторонние поверхности

Пусть по-прежнему в R3 задана гладкая поверхность Σ, описываемая параметрическими уравнениями:

~r = ~r(u, v) = [gx(u, v), gy(u, v), gz(u, v)]T , (u, v) D.

 

 

 

Вектор N~ = ~r u0

, ~r v0

 

является вектором, ортогональным к поверхности

Σ

(ñì. ðàçä. 1.8). Â

 

 

N~ = ~0. Пусть ~n+ = ~n(u, v) =

1

 

N~

 

силу уравнения (4.2),

6

 

 

 

 

 

 

 

~

 

 

 

 

 

 

 

kNk

 

вектор единичной нормали к поверхности Σ в точке с декартовыми ко-

íûì

 

 

Σ

 

 

~n

= −~n+ .

ординатами

gx(u, v), gy(u, v), gz(u, v)

. Отметим, что наряду с ~n+ , единич-

 

вектором, ортогональным к

 

, является также и вектор

 

 

Таким образом, в каждой точке Σ определены два взаимно противоположных единичных нормальных к Σ вектора. При этом каждая из функций ~n+(u, v) è ~n(u, v) непрерывна в области D.

Пусть `D гладкая замкнутая кривая, лежащая в области D, параметрическими уравнениями которой являются:

`D :

(v = h2(t);

t [t(1), t(2)].

 

u = h1(t);

 

Òàê êàê `D замкнута, то h1(t(1)) = h1(t(2)) è h2(t(1)) = h2(t(2)). Кривой `D соответствует кривая `S с параметрическими уравнениями:

`S : ~r = ~r(t) = [gx(h1(t), h2(t)), gy(h1(t), h2(t)), gz(h1(t), h2(t))]T , t [t(1), t(2)],

лежащая на поверхности Σ. Так как ~r(t(1)) = ~r(t(2)), то кривая `S также замкнута.

Для точек кривой `S функции ~n+(u, v) è ~n(u, v) являются непре-

рывными функциями ~n+(h1(t), h2(t)) è ~n(h1(t), h2(t)) параметра t при t [t(1), t(2)]. Это означает, в частности, что определены векторы

~n(1)

=

~n

 

h

t(1)

)

, h

2(

t(1)

)

lim

~n

 

h

t

, h t

,

+

 

+

1(

 

 

 

 

 

 

= t→t(1)+0

 

+

1( )

2( )

~n(2)

= ~n

+

h

 

t(2)

 

, h (t(2))

=

lim

~n

h (t), h (t)

 

+

 

 

 

1(

(1))

 

2 (2)

 

 

t→t(2)−0

 

 

+

1

2

и аналогичные векторы ~n

 

è ~n

для функции ~n(h1(t), h2(t)).

Оказывается, что, несмотря на то, что для замкнутой кривой `D âû-

полняются равенства

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

hi(t(1)) =

 

lim

 

hi(t) = hi(t(2)) =

lim

 

hi(t),

i = 1, 2,

 

 

t→t(1)+0

 

 

 

 

 

 

 

 

t→t(2)−0

 

 

 

66

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]