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

Основы высшей математики для инженеров 2009

.pdf
Скачиваний:
226
Добавлен:
12.03.2016
Размер:
10.04 Mб
Скачать

5.9. Формула Грина – Остроградского…

271

Рассмотрим сначала

J1 DD

HP

 

b

>(x)

HP

dy.

 

 

 

 

 

 

 

dN

Ddx D

 

 

 

 

 

 

 

 

 

 

 

 

Hy

a

((x)

Hy

 

 

 

 

N

 

 

 

 

 

Так как при постоянном x

 

P(x, y) есть одна из пер

 

 

 

 

 

 

 

вообразных для

HP

 

, то

 

 

 

 

 

Рис. 5.21

 

 

 

 

 

Hy

 

 

 

 

 

 

 

>(x)

HP

dx P(x, y)|>(x)

 

 

 

 

 

D

P(x, >(x)) P(x, ((x)),

 

((x) Hy

 

((x)

 

 

 

 

 

 

 

 

 

 

 

 

т. е.

bb

J1 DP(x, >(x))dx DP(x, ((x))dx.

aa

Но по способу вычисления криволинейного интеграла

 

b

 

 

 

 

 

DP(x, >(x))dx

DP(x, y)dx

 

DP(x, y)dx,

 

M62

 

 

a

 

 

 

 

 

AmB

BmA

b

DP(x, ((x))dx DP(x, y)dx.

a

 

AnB

 

Дуги AnB и BmA дают в совокупности границу L области N, проходи мую в положительном направлении, поэтому

J1

DP(x, y)dx DD

HP

dN.

(1)

 

 

 

Hy

 

LN

Аналогично показывается, что

DD

HQ

dN DQ(x, y)dy.

(2)

Hx

N

L

 

 

 

Вычитая почленно равенства (1) и (2), получим

HQ

HP

DD

 

 

 

dN DQ(x, y)dy P(x, y)dx.

 

 

 

Hx

Hy

 

N

L

Доказанная формула может быть применена для вычисления пло щади плоской области при помощи криволинейного интеграла. Поло жим P(x, y) J 0, Q(x, y) J x. Тогда формула примет вид

DDdN D xdy.

NL

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

Так как слева стоит площадь области N, то получаем формулу

N D xdy.

L

Аналогично, полагая Q(x, y) J 0, P(x, y) y, получим

N D ydx.

L

П р и м е р. Вычислить площадь фигуры, ограниченной эллипсом

x a cost,

y b sin t.

Решение. Если обход эллипса совершается против часовой стрелки, так, чтобы область, ограниченная им, оставалась слева, то 0 ( 2'.

2'

2'

2'

1 cos2t

 

N Da cost b cost dt ab Dcos2t ab D

dt 'ab.

2

0

0

0

 

 

 

5.9.2. Независимость криволинейногоM62интеграла от пути интегрирования

Пусть функции P(x, y) и Q(x, y) непрерывны вместе с частными

производными HP и HQ в некоторой области G плоскости XOY. Рас

Hy Hx

смотрим в G две произвольные точки А и В (рис. 5.22). Эти точки можно соединить различными линиями, лежащими в G. Нетрудно проверить, что интегралы

DPdx Qdy, и

DPdx Qdy,

L1

L2

Рис. 5.22

вообще говоря, различны.

 

Ставится вопрос: когда интеграл DPdx Qdy принимает одинако

L

вое значение на любой кривой L, соединяющей точки А и В, т. е. когда интеграл не зависит от пути интегрирования.

Теорема 1. Для того, чтобы криволинейный интеграл DPdx Qdy в

L

некоторой области G не зависел от пути интегрирования, необходимо

5.9. Формула Грина – Остроградского…

273

и достаточно, чтобы интеграл по любому замкнутому контуру, лежаще2 му в G, был равен 0.

Доказательство. Достаточность. Пусть интеграл по замкнутому контуру в области G равен 0. Пусть А и В — любые точки G. Соеди ним эти точки в G произвольными путями AnB и AmB.

По условию

 

 

DPdx Qdy

D D 0.

 

 

AmBnA

AmBBnA

 

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

 

Рис. 5.23

DPdx Qdy DPdx Qdy DPdx Qdy.

 

AmB

BnA

AnB

Необходимость. Пусть интеграл не зависит от пути интегрирова ния, тогда

DPdx Qdy DPdx Qdy,

AnB

M62

AnB

 

AmB

т. е.

 

 

DPdx Qdy

DPdx Qdy 0,

 

AmB

но

 

 

 

D D,

 

AmB

BmA

т. е.

 

 

 

D D 0,

 

AnB BmA

но

 

 

0 D D

DPdx Qdy.

 

AnB BnA AnBmA

Теорема доказана.

Теорема 2. Для того, чтобы криволинейный интеграл DPdx Qdy не

L

зависел от пути интегрирования в односвязной области G, необходимо и

достаточно, чтобы в каждой точке этой области HP HQ .

Hy Hx

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

Достаточность. Пусть в области G HP HQ . По формуле Остро

Hy Hx

градского – Грина

HQ

HP

(3)

DD

 

 

 

dN DPdx Qdy,

 

 

 

Hx

Hy

 

 

 

 

 

N

L

 

 

 

 

где N — площадка, ограниченная контуром L. Эта площадка целиком

принадлежит G в силу односвязанности. Так как везде в G

HP

 

HQ

,

 

 

 

 

 

 

 

Hy Hx

то интеграл слева в (3) равен 0, следовательно, и интеграл по контуру L равен 0. Таким образом, интеграл по любому замкнутому контуру равен 0 и, следовательно, по теореме 1 интеграл не зависит от пути интегрирования.

Необходимость. Пусть интеграл DPdx Qdy не зависит от пути интег

ка (x0, y0), что

 

 

 

 

L

 

 

 

Hy

M62Hx . Пусть, например,

 

Hy

 

рирования. Покажем, что

HP

HQ в G. Пусть это не так. Тогда есть точ

 

 

 

Hy

Hx

 

 

 

 

HP(x0, y0)

 

 

HQ(x0, y0)

 

 

HP(x0, y0)

 

HQ(x0, y0) 0. В силу непрерывности HP и HQ существует окрестность

Hx Hy Hx

точки (x0, y0), в которой имеет место неравенство HP HQ 0. Рассмот

Hy Hx

рим круг N столь малого радиуса, что он целиком содержится в этой ок рестности. За контур L примем окружность, ограничивающую этот круг. В силу свойства (3) двойного интеграла и по формуле Остроградско го – Грина

 

HQ

 

HP

DPdx Qdy DD

 

 

 

dN 0.

Hx

Hy

L

 

 

 

N

 

 

Таким образом, интеграл по замкнутому контуру отличен от 0, а это противоречит утверждению теоремы 1. Теорема доказана.

Вспомним, что условие HP HQ является необходимым и доста

Hy Hx

точным условием того, чтобы выражение P(x, y)dx Q(x, y)dy было дифференциалом du некоторой функции u(x, y). Если функция u(x, y)

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

275

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

(x2, y2)

DP(x, y)dx Q(x, y)dy u(x2, y2) u(x1, y1),

(x1, y1)

где (x1, y1) и (x2, y2) — начальная конечная точка пути интегрирования.

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

Криволинейный интеграл по длине дуги. Вычисление массы кривой. Ин теграл по площади поверхности (поверхностный интеграл первого рода). Вычисление поверхностного интеграла первого рода.

5.10.1. Криволинейный интеграл по длине дуги. Вычисление массы кривой

Пусть в плоскости XOYM62задана дуга l, соединяющая точки А и В. Пусть в точках дуги определена непрерывная функция f (x, y). Разо бьем дугу АВ на n частей точками A0 A, A1, …, An B. Длину дуги

Ai Ai 1 (i 0, 1, …, n 1) обозначим через -li. На каждой из дуг Ai Ai 1 выберем точку Mi(xi, yi) и составим сумму

n 1

n 1

 

 

F f (Mi)-li F f (xi, yi)-li.

 

 

i 0

i 0

 

 

Если существует предел этой суммы при

 

условии, что все -li 0, который не зависит

 

ни от способа разбиения дуги на части, ни от

 

выбора точек Mi(xi, yi), то он называется кри

 

волинейным интегралом по длине дуги, или

 

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

Рис. 5.24

 

 

n 1

 

D f (x, y)dl lim

F f (xi, yi)-li.

 

AB

n

i 0

 

-li 0

 

 

 

 

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

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

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

Если дуга АВ задана уравнением y ((x), x [a, b], то

b

D f (x, y)dl D f [x,((x)]1 (2(x) dx.

AB a

Физический смысл: если f (x, y) — плотность точек кривой, то

D f (x, y)dl

AB

выражает массу дуги.

П р и м е р. Найти массу дуги кривой y ln x между точками с абс циссами x 1, x 2, если плотность + x2 .

 

 

2

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

1

 

 

m D+(x, y)dl

D x2 dl D x2

1 (ln x) 2 dx D x2 1

 

dx

x2

AB

AB

1

 

 

 

 

1

 

 

 

 

2

M62

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D x 1 x dx

1

(5 5 2 2).

3

1

 

 

 

5.10.2. Интеграл по площади поверхности (поверхностный интеграл первого рода)

Рассмотрение функций нескольких переменных естественно при водит к дальнейшему обобщению понятия интеграла — к поверхност ным интегралам, которые конструктивно вводятся так же, как криво линейные интегралы.

Пусть на поверхности S, ограниченной замкнутой кривой L, зада на функция f (M) f (x, y, z).

1.Разобьем S на n частей -S1, -S2 , …, -Sn.

2.На площадке разбиения -Si выберем произвольно точку Mi(xi, yi, zi) и вычислим

f (Mi)-Si f (xi, yi, zi)-Si

(i 1, 2, , n).

 

3. Составим сумму

 

 

 

n

n

 

 

F f (Mi)-Si F f (xi, yi, zi)-Si.

(1)

i 1

i 1

 

 

(1) — интегральная сумма.

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

277

ОПРЕДЕЛЕНИЕ. Если существует предел суммы (1) при n , -Si 0, который не зависит ни от способа разбиения, ни от выбора точек Mi, то он называется интегралом по площади поверхности:

 

 

 

n

 

DD f (x, y, z)dS

lim

 

F f (xi, yi, zi)-Si.

(2)

S

n

 

i 1

 

max -S i

0

 

 

 

 

Если +(x, y, z) — поверхностная плотность, то масса поверхно сти S есть

m DD+(x, y, z)dS.

S

В этом состоит физический смысл интеграла.

Интеграл (2) обладает теми же свойствами, что и двойной.

1. Если S S1 S2 , то

M62DD DD DD.

S S1 S

2. Если f (x, y, z) 0, то

DD f (x, y, z)dS 0.

S

Пусть поверхность площади S задана уравнением z ((x, y), где ((x, y) — непрерывная вместе с частными производными 1 го порядка функция.

Теорема существования. Для всякой непрерывной функции f (x, y, z) су2 ществует поверхностный интеграл первого рода.

5.10.3. Вычисление поверхностного интеграла первого рода

Пусть N — проекция поверхности S на плоскость XOY, а u f (x, y, z) — непрерывная функция, определенная во всех точках по верхности S. Разобьем N на части -N1, -N2 , …, -Nn. Тогда S разобьется на части -S1, -S2, , -Sn, при этом -Si — часть поверхности, которая проектируется в -Ni.

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

 

 

-Si DD 1 (x2(x, y) (y2(x, y) dN.

 

 

 

 

(-Ni )

 

 

 

 

 

 

 

По теореме о среднем

 

 

 

 

 

 

-S

i

1 (2

(x , y ) (2

(x , y ) -N

,

 

 

 

 

x

i i

y

i i

i

 

 

 

где xi, yi

 

координаты

некоторой точки

 

Рис. 5.25

Pi(xi, yi) -Ni. Возьмем в каждой части -Si

 

точку Mi(xi, yi, zi), где zi ((xi, yi).

 

 

 

 

 

 

Для данного разбиения поверхности S на части и для данного выбора

точек Mi составим интегральную сумму

 

 

 

 

 

n

 

n

 

 

1 (x2(xi, yi) (y2(xi, yi) -Ni. (3)

F f (xi, yi, zi)-Si F f (xi, yi, ((xi, yi))

i 1

 

i 1

 

 

 

 

 

 

 

 

Перейдем к пределу при n , считая, что -Ni 0. Сумма (3) являет ся интегральной суммойM62по области N для непрерывной функции

f (x, y, ((x, y)) 1 (x2 (y2 . Поэтому ее предел есть двойной интеграл по области N:

DD f (x, y, ((x, y))1 (x2(x, y) (y2(x, y) dN.

N

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

 

n

 

lim

F f (xi, yi, zi)-Si DD f (x, y, z)dS.

n

i 1

S

-S i 0

 

 

Таким образом, окончательно,

DD f (x, y, z)dS DD f (x, y, ((x, y))1 (x2(x, y) (y2(x, y) dN.

SN

 

П р и м е р. Вычислить массу

m части

S параболоида вращения

1

2

 

2

 

2

2

 

z

 

(x

y

 

), вырезанную цилиндром x y

 

4, если плотность в ка

 

 

 

2 ждой точке пропорциональна квадрату расстояния от этой точки

до оси OZ.

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

279

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

Задача о потоке жидкости. Определение поверхностного интеграла вто рого рода. Вычисление поверхностного интеграла второго рода.

5.11.1. Задача о потоке жидкости через поверхность

Рассмотрим движение жидкости в пространстве. Пусть скорость v частицы, протекающей через данную точку М пространства, зависит только от этой точки, т. е.

v v(x, y, z)

или

v vx(x, y, z)i vy(x, y, z)j vz(x, y, z)k.

Вычислим поток жидкости через поверхность S за единицу времени,

считая, что плотность жидкости + 1.

 

1. Пусть v const,

а

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

 

плоская площадка. За единицу времени час

 

тицы жидкости, находящиеся на площадке

 

S, переместятся в направлении v на расстоя

 

ние, равное по длине

|

v

|,

будут на пло

 

 

 

M62т. е.

 

щадке S1. Количество жидкости, которое

 

пройдет за единицу времени через S, чис

Рис. 5.26

ленно равно объему цилиндра: Q S h, где h — высота.

Пусть n — единичный вектор нормали к S. Угол между векторами

n

и

v

обозначим через (. Тогда

 

 

 

 

h |

 

| cos ( |

 

| |

 

 

| cos ( (

 

,

 

).

 

 

 

 

v

v

n

v

n

 

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

(1)

 

 

 

 

 

Q (

 

,

 

)S.

 

 

 

 

 

v

n

2. Пусть v v(x, y, z), а S — произвольная поверхность. Пусть в произвольной точке М этой поверхности определен единичный нор мальный вектор n cos) i cos* j cos L k , направляющие косинусы которого являются непрерывными функциями координат точек по верхности.

Найдем количество жидкости Q, протекающей в единицу времени через эту поверхность. Поскольку v меняется от точки к точке по величи не и направлению, то непосредственно применить формулу (1) нельзя. Поэтому разобьем поверхность S на n частей -S1, -S2 , …, -Sn и выберем

За количество жидкости, протекающей через поверхность S в еди ницу времени, принимается предел, к которому стремится правая часть последнего равенства при n и max -Si 0:
n
i 1
i 1

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

в каждой части точку Mi(xi, yi, zi). Пусть ni i cos)i j cos*i k cos Li — единичный вектор нормали к поверхности S в точке Mi.

Будем считать, что в пределах -Si скорость v const и равна ее значению vi в точке Mi:

vi vx(xi, yi, zi)i vy(xi, yi, zi)j vz(xi, yi, zi)k.

Будем также предполагать, что область -Si практически плоская. Тогда

-Qi = (vi, ni)-Si

nn

Q F-Qi = F(v, ni)-Si.

Q lim

F(

v

,

ni)-Si.

n

 

 

 

 

v

nM62v v

max -S i

0 i 1

 

Расписывая скалярное произведение (v, ni) в координатной форме, по лучим

Q lim F[ x cos)i y cos*i z cos L z]-Si.

(2)

n

max -S i 0 i 1

В равенстве (2) справа записана интегральная сумма для непрерывной функции vx(x, y, z) cos) vy(x, y, z) cos* vz cos L. Поэтому

Q DD(

v

,

n

)dS DD(vx cos) vy cos* vz cos L)dS.

(3)

S

 

 

 

S

 

Итак, количество жидкости, протекающей в единицу времени че рез данную поверхность S, или, как говорят, поток жидкости через эту поверхность есть интеграл по поверхности S, определенный в соответ ствии с формулой (3).

5.11.2. Определение поверхностного интеграла второго рода

Пусть F P(x, y, z)i Q(x, y, z)j R(x, y, z)k

некоторая

вектор

функция в пространстве, S — некоторая поверхность. Предположим,

что в каждой точке поверхности S определен

единичный

вектор