
Галкин С.В. Теория функций комплексного переменного и операционное исчисление- учебное пособие для вузов (2011) (1)
.pdfП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