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

Галкин С.В. Теория функций комплексного переменного и операционное исчисление- учебное пособие для вузов (2011) (1)

.pdf
Скачиваний:
340
Добавлен:
23.03.2016
Размер:
1.65 Mб
Скачать

П2. СКАЛЯРНЫЕ И ВЕКТОРНЫЕ ПОЛЯ

П2.1. Характеристики скалярных и векторных полей

Говорят, что в области (плоской или пространственной) задано скалярное поле (M), если в этой области задана скалярная функция (M). Говорят, что в области (плоской или пространственной) задано векторное поле a (M), если в этой области задана векторная функция a (M). Например, масса или температура частиц в комнате – скалярные поля, скорость или силы взаимодействия частиц – векторные поля.

В интегралах первого рода (двойных, криволинейных, поверхностных) имеют дело со скалярным полем – распределением масс точек кривой или поверхности в пространстве. В интегралах второго рода вычисляются характеристики векторных полей: работа векторного поля (силового поля) в криволинейном интеграле, поток векторного поля в поверхностном интеграле. Рассмотрим подробнее основные характеристики скалярных и векторных полей.

П2.1.1. Характеристики скалярных полей

Линии уровня плоского поля (x, y) – кривые, на которых значения функции постоянны: (x, y) = С = const, например линии равной высоты, нанесенные на географическую карту (h (x, y) = 0 – уровень моря, h = 7000 м – немногие горные вершины, h = – 10 000 м – самые глубокие океанские впадины).

Поверхности уровня пространственного поля (x, y, z) – поверхно-

сти, на которых значения функции постоянны: (x, y, z) = С = const, например поверхности с равной температурой воздуха в атмосфере. Любая линия на поверхности уровня поля – это линия уровня этого поля.

Пример. Задано поле x2 y2 z2 C. При положительной

постоянной С поверхности уровня – однополостные гиперболоиды, при С = 0 поверхность уровня – конус, при С < 0 поверхности уровня – двуполостные гиперболоиды.

Линии или поверхности различных уровней поля не пересекаются, так как каждая линия или поверхность уровня поля соответ-

180

ствует своей константе. Чем чаще (гуще) поверхности или линии уровня поля, тем интенсивнее его изменение.

 

 

 

,

 

,

 

Градиентом поля называется – вектор grad

x

y

.

 

 

 

 

z

Теорема. Градиент скалярного поля ортогонален его поверхности уровня.

Доказательство. Пусть точка (x, y, z) остается на поверхности уровня g(x, y, z) = 0 при вариациях переменных. Тогда равенство превращается в тождество, а тождество можно дифференцировать:

dg(x, y, z)

g dx

g dy

 

 

0.

g dz gradg dr

 

x

y

z

 

 

Вектор dr (x, y, z) – это вектор, касательный в точке (x, y, z) к любой кривой, лежащей на поверхности уровня, проходящей через эту точку. Поэтому в точке (x, y, z) вектор градиента скалярного поля ортогонален всем касательным к линии уровня, проходящим через эту точку. Следовательно, он ортогонален касательной плоскости к поверхности уровня и направлен по нормали к поверхности уровня.

Производная скалярного поля по направлению l определяется как

g

 

limt 0

g M tl g(M )

.

 

l

 

t

 

M

 

 

 

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

g

 

 

 

l

 

 

 

 

 

 

 

 

gradg

 

 

.

l

 

M

 

| l

|

 

 

 

 

Пример. Найдем производную скалярного поля g(x, y, z) = = x2 + y2 + z3 по направлению {1, 3, 2} в точке (1, 0, 4). Имеем

181

 

2 y,

 

3z2 ,

g

 

 

 

 

 

 

 

 

 

 

gradg 2x,

 

 

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

1

 

3

 

2

1,0,4

98

 

 

 

 

 

 

 

 

2,0, 48

 

,

 

 

,

 

 

 

 

 

.

14

14

14

 

14

 

 

 

 

 

 

 

П2.1.2. Характеристики векторных полей

Векторная линия – линия, в каждой точке которой вектор поля направлен по касательной к ней. Уравнения векторной линии легко получить из условия коллинеарности векторов поля

a(M ) P(x, y, z)i Q(x, y, z) j R(x, y, z)k

и касательной

 

dyj dzk.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dr dxi

 

 

Эти уравнения выглядят так:

 

 

 

 

 

 

 

 

dx

dy

dz

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

P(x, y, z)

Q(x,

y, z)

R(x, y, z)

Пример. Уравнение векторных линий векторного поля

 

 

 

 

 

 

 

a(M ) yi xj

 

 

можно записать в виде

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

 

 

dy

.

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

Интегрируем дифференциальное уравнение:

 

dx

 

dy

,

xdx ydy, xdx ydy 0,

d(x2 y2 ) 0.

 

y

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

Положив С > 0, получим линии уровня

x2 y2 C, т. е. ок-

ружности.

Векторной трубкой называется поверхность, образованная векторными линиями.

182

П2.2. Формула Грина

Теорема (формула) Грина. Пусть G – плоская односвязная область с кусочно-гладкой границей L. Пусть функции P(x, y), Q(x, y) непрерывны и имеют непрерывные частные производные по своим переменным в области G и на ее границе L. Тогда справедлива

формула Грина

 

 

Q

 

P

P(x, y)dx Q(x, y)dy

x

dxdy.

L

G

 

y

Доказательство. Назовем плоскую область D (в плоскости OXY) правильной, если любая прямая, параллельная координатной оси (OX или OY), пересекает область не более чем в двух точках.

Покажем, что область G можно представить как объединение

конечного числа правильных областей G n Dk .

k 1

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

Поэтому доказательство может быть проведено для правильной области D (рис. П2.1). Пусть D – правильная область. Так как P, Q могут быть произвольными функциями, то формула Грина сводится к двум формулам

L

и

L

P(x, y)dx

Q(x, y)dy

 

 

 

P

 

dxdy

D

 

 

y

Q dxdy.

D x

183

Рис. П2.1

Эти формулы надо доказать. Докажем первую формулу, представив контур L в виде объединения дуг L1, L2 (вторая формула доказывается аналогично). Сведем двойной интеграл к повторным интегралам, а затем к двум определенным интегралам, вычисляя внутренний интеграл. Далее два определенных интеграла сведем к криволинейным интегралам, которые объединим в криволинейный интеграл по контуру:

 

 

P

b

( x)

 

P

 

 

b

y) | ((xx)) dx

 

dxdy

 

 

 

P(x,

y

 

dy dx

D

 

 

( x)

 

y

 

 

a

 

 

 

a

 

 

 

 

 

 

b

 

 

b

 

 

 

 

 

 

P(x, (x))dx P(x, (x))dx

 

 

 

a

 

 

a

 

 

 

 

 

 

P(x, y)dx P(x, y)dx

 

 

 

 

L2

 

L1

 

 

 

 

 

 

P(x, y)dx P(x, y)dx

 

P(x, y)dx P(x, y)dx.

 

L2

 

L1

 

 

 

L1 L2

 

 

L

184

П2.3. Формула Остроградского – Гаусса

Теорема (формула) Остроградского Гаусса. Пусть компо-

ненты векторного поля

a(M ) P(x, y, z)i Q(x, y, z) j R(x, y, z)k

непрерывны и имеют непрерывные частные производные в пространственно односвязной замкнутой области V и на ее кусочногладкой границе . Тогда справедлива формула Остроградского –

Гаусса

P(x, y, z)dydz Q(x, y, z)dxdz R(x, y, z)dxdy

 

 

 

 

 

 

 

 

 

P

 

Q

 

R

 

x

y

dv.

V

 

 

 

z

Отметим, что левая часть формулы представляет собой поток векторного поля a(M ) через поверхность – «количество» вектор-

ного поля, протекающее через поверхность в единицу времени. Доказательство. Формула Остроградского – Гаусса в си-

лу произвольности функций P, Q, R состоит из трех частей, в каждую из которых входит одна из компонент векторного поля P, Q, R. В самом деле, можно взять P = 0, Q = 0 и доказывать отдельно часть формулы, в которую входит только R. Остальные части формулы (при P = 0, R = 0, Q = 0, R = 0) доказываются аналогично.

Будем доказывать часть формулы Остроградского – Гаусса:

 

R(x, y, z)dxdy

R

z dv.

 

V

 

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

185

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

Итак, будем доказывать соотношение

 

R(x, y, z)dxdy

R

z dv

 

V

 

для цилиндрического тела V, проектирующегося в область D на плоскости OXY (рис. П2.2). Пусть «верхняя» граница цилиндрического

тела – поверхность 1 – описывается уравнением z z1(x, y), «нижняя» граница – поверхность 2 – описывается уравнением z z2 (x, y). Боковую поверхность цилиндрического тела, параллельную оси OZ, обозначим 3.

Рис. П2.2

186

Сразу отметим, что поток векторного поля через боковую поверхность равен нулю. Действительно,

R(x, y, z)dxdy R x, y, z cos d 0,

3

3

так как нормаль на боковой поверхности ортогональна оси OZ и cos 0 .

Отметим также, что на «верхней» поверхности циллиндрического тела 1 cos 0, а на «нижней» поверхности 2 cos 0.

Поэтому при переходе от поверхностного интеграла по поверхности 2 к двойному интегралу по области D и обратно надо менять

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

надо. Для доказательства соотношения

 

 

 

 

 

 

R

 

 

 

 

 

R(x, y, z)dxdy z dv

 

 

 

 

 

 

 

V

 

 

 

выполним следующие выкладки:

 

 

 

 

R

 

z1

( x, y)

R

 

 

 

 

 

 

 

 

dxdy R(x,

y,

z1(x,

y) dxdy

z dxdydz

z dz

 

 

 

 

 

 

 

 

 

 

V

D z2 ( x, y)

 

 

D

 

 

 

R(x,

y, z2 (x,

y) dxdy

R(x,

y, z)dxdy R(x,

y, z)dxdy

D

 

 

 

1

 

 

2

 

 

R(x,

y, z)dxdy R(x, y,

z)dxdy R(x,

y,

z)dxdy

1

 

 

2

 

 

3

 

 

 

R(x,

y, z)dxdy.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В результате выкладок получено искомое соотношение. Таким образом, часть формулы Остроградского – Гаусса доказана.

187

Замечание. Формулу Остроградского – Гаусса можно записать в «полевом» виде:

П (a) a n d div a dv

V

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

П2.4. Дивергенция и ее свойства

Дивергенция векторного поля (расходимость) есть

P Q R .x y z

Дивергенция – это характеристика векторного поля, инвариантная относительно системы координат. Покажем это. Дадим ин-

вариантное определение дивергенции.

Рассмотрим произвольную точку M в пространственной области V. Выберем ее окрестность VM – шар радиусом r с центром в

точке M. Обозначим M ее границу – сферу радиусом r. По теореме о среднем для тройного интеграла

diva dv diva(M )VM , M VM ,

VM

где VM – объем шара.

По формуле Остроградского – Гаусса

 

 

1

 

 

1

 

 

 

diva

M

 

divadv

 

 

a

nd .

V

V

 

 

M VM

 

M M

 

 

Стягиваем окрестность VM к точке M, получаем дивергенцию векторного поля в точке M инвариантное определение дивергенции:

188

diva(M ) lim

a n d

 

M

.

 

VM M

VM

 

 

 

Отсюда следует, что дивергенция векторного поля в точке M имеет смысл объемной плотности потока векторного поля через окрестность этой точки и характеризует мощность источника или стока векторного поля в точке M. Если diva(M ) > 0, то точка M

источник векторного поля; если div a(M ) < 0, то точка M – сток

векторного поля. Если в некоторой области дивергенция равна нулю, то в этой области нет ни источников, ни стоков, поток векторного поля через границу такой области равен нулю – «сколько поля втекает в область, столько и вытекает из нее».

Пример. Определим расположение источников и стоков векторного поля

a x2 yi xzyj xz2k.

Выясним, является ли точка M(1, 2, 3) источником или стоком. Вычислим diva 2xy xz 2xz 2xy xz. Все точки, для которых

выполнено неравенство 2xy + xz > 0, – источники; все точки, для которых выполнено неравенство 2xy + xz < 0, – стоки. На поверхности 2xy + xz = 0 нет ни источников, ни стоков. Точка M – источ-

ник, так как diva(M ) 4 3 7 0.

Сформулируем и докажем свойства дивергенции.

1. Дивергенция, как всякий линейный оператор, обладает свой-

ством линейности:

div 1a1 2a2 1diva1 2diva2 .

Докажем это:

div 1a1 2a2 div ( 1P1 + 2 P2 )i ( 1Q1 + 2Q2 )j ( 1R1 + 2 R2 )k

x ( 1P1 2 P2 ) + y ( 1Q1 2Q2 ) z ( 1R1 2 R2 )

1

 

P1

 

Q1

 

R1

 

2

 

P2

 

Q2

 

R2

 

 

 

 

x

y

z

 

 

x

y

z

 

1diva1

2diva2 .

 

 

 

 

 

 

 

 

 

 

 

 

189