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

Ми10 крив и пов

.pdf
Скачиваний:
8
Добавлен:
09.05.2015
Размер:
2.67 Mб
Скачать

условимся брать направление движения по кривой L так, чтобы область, ограниченная этой кривой, оставалась слева (рис. 1.3.4). Такое направление интегрирования называют положительным. Перемещение в противоположном направлении называют отрицательным.

Замечание 3. Если Y(x, y) 0, то интеграл принимает вид

X (x, y)dx и называется криволинейным интегралом по координате x .

AB

Если X (x, y) 0, то Y(x, y)dy называют криволинейным интегралом

AB

по координате y.

Рис.1.3.4

Поскольку определяемый нами интеграл (1.16) можно записать в виде суммы двух криволинейных интегралов по координатам x и y

X (x, y)dx +Y(x, y)dy = X (x, y)dx + Y(x, y)dy, то его называют еще

AB

 

AB

AB

 

 

составным криволинейным интегралом по координатам.

 

 

Подынтегральное выражение в левой части последнего равенства есть

скалярное произведение вектора

 

 

JG

G

JG

=G(X (x, y),Y(x, y)) и

 

 

F(x, y) = X (x, y) i

+Y(x,Gy) j

G

 

дифференциалаdr = dx i + dy j = (dx, dy) радиус-вектора r переменной

точки M кривой AB . Поэтому криволинейный интеграл от вектор-

G

 

JG

 

 

JG

функции F по кривой AB можно записать в векторной форме: F

dr.

AB

Если кривая L – замкнутая, то криволинейный интеграл

JG G G

v∫F dr называют циркуляцией векторного поля F по замкнутому

L

контуруL.

Возвращаясь к задаче о работе, можно сказать, что работа

JG

поля F при перемещении материальной точки по кривой AB криволинейным интегралом второго рода

силового

выражается

31

Α =

X (x, y)dx +Y(x, y)dy =

G

G

 

F

dr.

(1.17)

AB

JG

AB

 

 

 

Таким образом, если вектор F(x, y) задает силовое поле, то криволинейный интеграл второго рода выражает работу этого поля вдоль линии AB .

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

JG

G

JG

 

 

G

вектор-функции F(x, y, z) = X (x, y, z) i

+Y(x, y, z) j

+ Z(x, y, z) k,заданной

вдоль пространственной кривой AB , который аналитически записывают

так:

 

 

G

G

 

X (x, y, z)dx +Y(x, y, z)dy + Z(x, y, z)dz = F

dr.

(1.18)

AB

 

 

AB

 

 

1.2.3. Связь между криволинейными интегралами первого и второго рода

Пусть на кривой L установлено направление движения. Взяв точку M на этой кривой, проведем касательную к линии L в точке M и установим на касательной направление, соответствующее направлению движения по кривой.

Рис.1.3.5

Отложим по касательной в установленном направлении дифференциал

дуги ds , получим вектор dτ (рис. 1.3.5), проекциями которого служат дифференциалы dx, dy, dz . Обозначим через α(M ), β(M ),γ(M ) углы,

оразованные касательным вектором dτ с осями координат Ox,Oy,Oz . НапомнимGформулы, по которым определяются проекции вектора на оси: dx =| dτ | cosα(M ) = cosα(M ) ds, dy = cos β(M ) ds, dz = cosγ(M ) ds.

Поэтому составной криволинейный интеграл второго рода может

быть заменен интегралом первого рода по формуле

Xdx +Ydy + Zdz =(X cosα(M ) +Y cos β(M ) + Z cosγ(M ))ds.

L L

32

1.2.4. Существование и вычисление криволинейного интеграла второго рода

Теорема 1.2. Пусть кривая AB задана параметрическими уравнениями

x = x(t);

t

[α; β], причем функции x = x(t), y = y(t)

непрерывны, вместе

 

 

 

y = y(t),

 

 

 

 

 

 

со своими

 

производными первого порядка

на отрезке t [α; β] , где

α = t

A

, β = t

B

значения параметра t , соответствующие точкам A и B .

 

 

 

всякой вектор-функции

G

G

JG

Тогда

для

F(x, y)

= X (x, y) i +Y(x, y) j,

непрерывной вдоль кривой AB , существует криволинейный интеграл и

имеет место равенство

 

 

 

 

 

 

 

 

β

 

 

 

X (x, y)dx +Y (x, y)dy = ( X (x(t), y(t))x' (t) +Y (x(t), y(t)) y' (t))dt.

(1.19)

AB

 

 

 

 

α

 

 

 

Доказательство. Остановимся на втором утверждении теоремы. Покажем, что вычисление криволинейного интеграла второго рода сводится к вычислению определенного интеграла по формуле (1.19). Выберем произвольное разбиение кривой AB точками A0 = A, A1 ,..., An = B на

элементарные дуги Ai1 Ai длиной si . Пусть ti ,i = 0,..., n, – значения параметра кривой, соответствующие точкам Ai . Для этих значений верны

соотношения

t0 =α < t1 < t2 <... < tn1 < tn = β .

На каждой элементарной

дугеAi1 Ai

 

 

 

 

 

 

 

 

 

 

 

 

 

 

kk

 

произвольным

 

образом

выберем точку Mi (xi , yi

) .

По

определению для криволинейного интеграла по координате x имеем

 

 

 

 

 

 

 

n

kk

 

 

n

kk

 

 

 

 

 

 

X (x, y)dx = λ=maxlims 0 X (xi

, yi )xi

=maxlimx 0 X (xi , yi )xi

,

 

 

(1.20)

AB

 

i

i

i=1

 

 

 

i

i

i=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где xi = xi xi1 = x(ti ) x(ti1 ).

 

Последнюю

разность

преобразуем

по

формуле

Лагранжа

x

 

= x(t

) x(t

) = x' (τ

)t .Здесь

t

i

= t

i

t

и

 

 

 

i

 

i

 

i1

 

i

i

 

 

 

 

i1

 

ti1 <τi < ti , i =1,..., n..

 

 

 

 

 

 

kk

дуге Ai1 Ai выбрали

Так

 

как точку Mi (xi

, yi ) на

произвольным

образом,

 

то

выберем

ее так, чтобы ее координаты

 

 

 

 

 

 

 

 

 

k

 

k

 

 

 

 

 

 

 

соответствовали значению параметра τi : xi

= x(τi ),

yi = y(τi ) . Подставляя

найденные значения xi , yi , xi

в формулу (1.20)

и учитывая, что если

длина дуги si 0, то xi

0, следовательно ti 0, получим

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

X (x, y)dx = maxlimt 0 X (x(τi ), y(τi ))x' (τi )ti .

 

 

 

 

 

 

 

 

AB

 

 

 

i

i

i=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

33

t : X (x(t), y(t))x' (t), заданной

на

отрезке [α, β]. Ее предел равен

определенному интегралу. Таким образом, имеем

 

 

β

 

 

X (x, y)dx = X (x(t), y(t))x' (t)dt.

 

AB

α

 

 

 

 

β

Аналогично рассуждая, найдем

Y(x, y)dx = Y(x(t), y(t)) y' (t)dt.

 

 

AB

α

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

 

β

 

 

X (x, y)dx +Y(x, y)dy = (X (x(t), y(t))x' (t) +Y(x(t), y(t)) y' (t))dt.

AB

α

 

 

Доказательство первого утверждения теоремы смотрите, например [3]. Согласно доказанной теореме, для вычисления криволинейного

интеграла второго рода в подынтегральном выражении необходимо от переменных x и y перейти к параметру t кривой, для чего через t следует

выразить подынтегральные функции X (x, y), Y(x, y) и дифференциалы dx, dy . В доказательстве теоремы предполагалось, что начальной точке A кривой AB соответствует левый конец отрезка [α, β], являющегося областью изменения параметра кривой. Если же параметризация кривой не согласована с направлением ее обхода, то в формуле (1.19) определенный интеграл справа соответствует противоположному направлению обхода, т.е. левая и правая части формулы отличаются знаками. Для восстановления равенства можно в определенном интеграле поменять местами пределы интегрирования. Таким образом, действует следующее правило: при переходе от криволинейного интеграла второго рода к определенному интегралу нижний предел интегрирования должен соответствовать начальной точке кривой, а верхний предел интегрирования – конечной точке.

Рассмотрим частные случаи задания плоской кривой AB .

Если кривая AB задана уравнением, y = y(x), x [a,b], то в качестве

параметра кривой можно взять абсциссу x точки на кривой. В этом случае в соответствии с формулой (1.19) получаем

X (x, y)dx +Y(x, y)dy = b (X (x, y(x)) +Y(x, y(x)) y' (x))dx,

(1.21)

AB

a

 

где a и b – абсциссы точек A и B этой кривой.

 

Если кривая

AB задана уравнением, x = x( y), y [c, d], то в качестве

параметра кривой обычно берется ордината y точки кривой. Тогда

 

X (x, y)dx +Y(x, y)dy = d (X (x( y), y)x' ( y) +Y(x( y), y))dy.

(1.22)

AB

c

 

34

Здесь c и d — ординаты точек A и B кривой.

Теорема 1.2 без каких-либо затруднений переносится на случай пространственной кривой. Пусть пространственная кривая AB задана параметрическими уравнениями

x = x(t);

y = y(t); t [α, β],

z = z(t),

где x(t), y(t), z(t) – функции, непрерывные на отрезке [α, β] вместе со

своими производными, причем начальной точке A кривой соответствует значение параметра t =α, конечной точке B кривой — значение параметра t = β . Тогда криволинейный интеграл второго рода общего вида

(1.18) от

непрерывных

функций X (x, y, z),Y(x, y, z), Z(x, y, z)

вдоль

пространственной кривой

AB существует, и для него верно равенство,

аналогичное (1.19):

 

 

 

X (x, y)dx +Y(x, y)dy + Z(x, y, z)dz =

 

 

AB

 

 

β

 

 

 

=(X(x(t), y(t), z(t))x' (t) +Y(x(t), y(t), z(t))y' (t) +Z(x(t), y(t), z(t))z' (t))dt.

(1.23)

α

 

 

 

Пример

1.9. Вычислить криволинейный интеграл второго

рода

I = (x + y)dx + 2zdy + xydz вдоль пространственной кривой AB , заданной

AB

 

 

 

 

x = t;

параметрическими уравнениями

 

y = t2 ; t [1;2] .

 

 

 

 

 

z = 3 t,

Решение. Так как dx = dt, dy = 2tdt, dz = −dt, то, используя формулу (1.23) , находим

I = 2

((t +t2 ) + 2(3 t) 2t + t t2 (1))dt =

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

2

 

3

 

13

 

2

 

3

 

t4

 

2

35

 

 

 

 

 

 

 

 

 

 

 

 

= (13t 3t

 

t

 

)dt =

 

 

t

 

t

 

 

 

 

=

 

.

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

2

 

 

 

 

 

4

 

1

4

 

 

 

 

 

 

 

 

 

 

 

G

 

 

 

 

G

 

JG

Пример 1.10. Найти работу силового поля F(x, y) = (4x y) i

+5x2 y j при

перемещении материальной точки вдоль куска параболы y = 3x2

между ее

точками A(0;0) и B(1;3).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Решение. В соответствии с формулами (1.17), (1.21) имеем:

 

 

35

Α = X (x, y)dx +Y(x, y)dy = (4x y)dx +5x2 ydy =

AB AB

= 1

((4x 3x2 ) +5x2 3x2 6x)dx = 1

(4x 3x2 +90x5 )dx = (2x2 x3 +15x6 )

 

1

=16.

 

0

0

0

 

1.2.5. Свойства криволинейного интеграла второго рода

Рассмотрим основные свойства криволинейного интеграла второго рода общего вида.

Свойства линейности:

1) k X (x, y)dx + k Y(x, y)dy = k X (x, y)dx +Y(x, y)dy;

AB

 

AB

2) (X1 (x, y) ± X2 (x, y))dx = X1 (x, y)dx ± X2 (x, y)dx,

AB

AB

AB

(Y1 (x, y) ±Y2 (x, y))dy = Y1 (x, y)dy ± Y2 (x, y)dy.

AB

AB

AB

3) Свойство аддитивности: если на кривой ABвзята произвольная точка C , то

X (x, y)dx +Y(x, y)dy = X (x, y)dx +Y(x, y)dy + X (x, y)dx +Y(x, y)dy.

AB

AC

CB

4) При изменении направления обхода кривой AB криволинейный

интеграл второго рода вдоль этой кривой меняет знак:

X (x, y)dx +Y(x, y)dy = −X (x, y)dx +Y(x, y)dy.

AB

 

BA

Свойства 1) – 4) несложно доказать, используя определение криволинейного интеграла второго рода как предела интегральных сумм и известные свойства предела. Эти доказательства повторяют доказательства соответствующих свойств определенного интеграла.

5) Если кривая AB представляет собой отрезок прямой, параллельной оси Oy , то

X (x, y)dx = 0.

AB

Если же AB – это отрезок прямой, параллельной оси Ox , то

Y(x, y)dy = 0.

AB

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

36

В частности, когда AB – отрезок вертикальной прямой x = x0 , y [c, d], то вдоль него x' ( y) 0 и

X (x, y)dx +Y(x, y)dy = d (X (x( y), y)x' ( y) +Y(x( y), y))dy = d Y(x0 , y))dy.

AB

c

 

c

 

 

 

Аналогично, если кривая AB – отрезок горизонтальной прямой

 

 

y = y0 , x [a,b], то

 

 

 

 

 

X (x, y)dx +Y(x, y)dy = b (X (x, y(x)) +Y(x, y(x)) y' (x))dx = b X (x, y0 )dx.

AB

a

G

a

G

JG

Пример 1.11.

Найти циркуляцию векторного поля

 

 

F(x, y)

= (x y) i

+5x j

вдоль замкнутого контура L (направление положительное). Контур L

контур треугольника с вершинами O(0;0), A(1;0), B(0;1)

(рис. 1.3.6).

 

 

 

 

 

Рис.1.G3.6

G

JG

Решение. Циркуляция векторного поля F(x, y) = (x y) i

+5x j вдоль

замкнутого контура L определяется криволинейным интегралом второго

рода

JG

G

G

 

G

v∫

 

F

dr

= >(x y)dx +5xdy, ( dr = dx i

+ dy j ).

L

 

 

L

 

 

Так как контур

 

L :OABO состоит из трех отрезков OA, AB и BO (рис.

1.3.4) , то в силу аддитивности криволинейного интеграла имеем:

>(x y)dx +5xdy = (x y)dx +5xdy + (x y)dx +5xdy + (x y)dx +5xdy.

L OA AB BO

1) Вычислим интеграл (x y)dx +5xdy .

OA

Так как OA: y = 0, x [0,1] , то dy 0 на отрезке OA. Следовательно

(x y)dx +5xdy = 1

(x 0)dx + 5x 0 =

x2

 

 

1

=

1

.

 

 

 

OA

0

2

 

0

 

2

 

 

 

 

 

37

2)

Для того,

чтобы вычислить

интеграл

(x y)dx +5xdy ,

 

найдем

 

 

 

 

 

 

 

 

AB

 

 

 

 

 

x 1

 

 

 

 

y 0

 

уравнение прямой, проходящей через точки

A(1;0) и B(0;1) :

=

 

 

 

 

 

 

 

 

1 0

 

y = −x +1. Тогда

AB : y = −x +1, x [1,0]

 

 

 

 

 

0 1

или

учетом

 

заданного

направления

движения

по

контуру

 

L),

dy = −dx.

 

 

 

 

Отсюда

(x y)dx +5xdy = 0

(x + x 1)dx 5xdx = −3 0

xdx 0 dx = −3

x2

 

 

0

x

 

10

=

5

.

 

 

 

 

 

 

 

 

 

 

 

 

 

AB

 

1

 

 

 

 

1

 

1

2

 

1

 

 

 

 

 

 

2

 

 

 

 

(x y)dx +5xdy.

 

 

 

 

 

 

 

 

 

 

 

 

3) Вычислим

интеграл

Так

как

BO :x = 0,

 

y [1,0],

 

 

 

 

BO

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx 0 на отрезке BO , то получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(x y)dx +5xdy = 0

(0 y) 0 +5 0 dy = 0.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

BO

 

 

1

 

 

 

G

 

 

 

 

 

 

 

 

 

 

 

G

 

JG

Следовательно, циркуляция векторного поля

 

 

 

 

 

 

 

 

 

 

 

 

F(x, y) = (x y) i +5x j

вдоль замкнутого контура L определится равенством

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

>(x y)dx +5xdy =.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L

 

 

 

 

 

 

1

 

 

 

 

5

 

 

 

 

 

 

 

 

 

= (x y)dx +5xdy + (x y)dx +5xdy + (x y)dx + 5xdy =

+

+ 0 = 3.

 

 

 

 

 

 

 

 

OA

 

AB

 

 

BO

 

 

 

2

 

 

 

 

2

 

 

 

 

 

 

 

 

 

1.2.6. Формула Грина

Напомним некоторые понятия.

Множество G называется областью, если это множество удовлетворяет двум условиям:

1)каждая точка множества является внутренней точкой множества, то есть множество является открытым множеством;

2)любые две точки множества можно соединить ломаной, каждая точка которой принадлежит данному множеству, то есть множество связно.

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

Вдальнейшем, через D будем обозначать замкнутую область G. Рассмотрим случай криволинейного интеграла второго рода вдоль

простого замкнутого контура L в плоскости xOy . В этом случае контур L является границей некоторой плоской замкнутой области D .

38

Теорема 1.3. (Формула Грина) Пусть односвязная область D на плоскости xOy ограничена простым кусочно гладким контуром L , а

функции X (x, y) и Y(x, y) непрерывны в D вместе со своими частными производными. Тогда справедлива формула, называемая формулой Грина:

>∫

 

Y

 

X

 

(1.24)

X (x, y)dx +Y(x, y)dy = ∫∫

dxdy,

L

D

x

 

y

 

 

где контур L обходится в положительном направлении, т.е. при движении по контуру область D остается слева.

Формула Грина позволяет заменять анализ одного явления, происходящего в области (двойной интеграл), анализом другого явления, происходящего только на границе области (криволинейный интеграл).

Доказательство. Формула Грина фактически распадается на две независимые формулы:

1)>LX (x, y)dx = −∫∫D Xy dxdy;

2)>LY(x, y)dx = ∫∫D Yx dxdy.

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

Oy .

 

 

 

Рис.1.3.7

 

 

 

 

y = y1 (x) и

Это значит, что она ограничена

снизу и

сверху

кривыми

y = y2 (x), где

функции y1 (x), y2 (x) непрерывны

на

отрезке

[a,b] и

удовлетворяют

неравенству

y1 (x) y2 (x), x [a,b],

а

слева

и справа –

вертикальными

отрезками

прямых

x = a и

x = b

(рис.

1.3.7).

Положительное

направление

обхода

контура

 

L

 

зададим

последовательностью точек

ABFEA

этого

конура. (В

частном

случае

39

каждый из вертикальных отрезков

AE и BF может выродиться в точку).

Докажем, что в этом случае верно равенство

 

 

 

 

 

 

>X (x, y)dx = −∫∫X dxdy.

 

(1.25)

 

 

 

 

L

 

 

 

D

y

 

 

Вычислим двойной интеграл:

 

 

 

 

 

 

 

X dxdy =b dx

y (x)

X dy =

b dx X(x, y)

 

 

=b (X(x, y(y2 (x)) X(x, y(y1(x)))dx.

∫∫

2

 

yy==yy12((xx))

 

y

 

D

y

a

y (x)

a

 

 

 

a

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

На основании свойства аддитивности криволинейного интеграла

второго рода имеем:

 

 

 

 

 

 

 

>X (x, y)dx = X (x, y)dx + X (x, y)dx + X (x, y)dx + X (x, y)dx.

(1.26)

L

 

AB

 

BF

 

 

 

FE

EA

 

Интегралы

X (x, y)dx и X (x, y)dx раны нулю на вертикальных

 

 

 

BF

 

EA

 

 

 

 

 

 

отрезках EA:x = a,

и BF : x = b, так как на этих отрезках dx = 0.

 

Используя правило вычисления криволинейного интеграла второго рода, заключаем, что

 

X (x, y)dx = b

X (x, y( y1 (x))dx

(AB : y = y1 (x), x [a;b]),

 

 

AB

 

a

 

 

 

 

 

 

 

 

X (x, y)dx = a X (x, y( y2 (x))dx (FE : y = y2 (x), x [b;a]).

 

 

FE

 

b

 

 

 

 

 

 

 

Подставляя эти соотношения в (1.26), находим

 

>∫

 

b

1

 

a

2

 

b

1

2

X(x, y)dx =

 

 

 

 

X(x, y(y (x))dx+

 

X(x, y(y (x))dx =

 

(X(x, y(y (x)) X(x, y(y (x)))dx.

L

 

a

 

 

b

 

 

a

 

 

Сравнивая выражения для двойного и криволинейного интегралов, видим, что различие только в знаке. Что и доказывает равенство (1.25).

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

кривыми на конечное число замкнутых областей, правильных относительно оси Oy . Записывая формулу (1.25) для каждой частичной

области и суммируя полученные равенства, мы приходим к формуле (1.25) для всей области D.

Формулу Грина можно распространить на случай многосвязной замкнутой области D, ограниченной внешним контуром L0 и внутренними

контурами L1 , L2 ,..., Ln . Здесь положительным направлением обхода внешнего L0 контура следует считать движение против часовой стрелки, а

положительным направлением обхода внутренних контуров – движение по часовой стрелке, поскольку при таком движении область D, ограниченная этими контурами, остается все время слева (рис. 1.3.8).

40