Методы теории функций в геометрии. Учебное пособие
.pdf
§ 2. Линейные операторы
В заданной системе координат между множеством линейных функций и матрицами вида (6.11) существует взаимно одноr-
значное соответствие. Суммой линейных операторов A и B
называют такой линейный оператор C , для которого выпол-
няется равенство
Cr(vr) = Ar(vr)+ B(vr), или C = Ar + B.
Аналогично определяется оператор разности линейных опера-
торов и умножение линейного оператора на скаляр λ :
Cr(vr)= λAr(vr), или C = λAr.
Произведением двух линейных операторов A и B называют
такой линейный оператор C , для которого выполняется ра-
венство
Cr(vr)= Ar(Br(vr)), или C = Ar B.
Легко проверяются свойства линейности (6.5) для опера-
тора C во всех случаях, как и то, что все рассмотренные операции над линейными операторами описываются аналогичными операциями над соответствующими матрицами.
Произведение линейного оператора A на самого себя назы-
вают квадратом (второй степенью) и обозначают A2 . Анало- r r
гично через A3 обозначают произведение A A2 и т. д. По определению полагают, что
A0 = E .
Сложение операторов ассоциативно, коммутативно, умножение операторов ассоциативно, но не коммутативно. Умножение и
сложениеоператоровобладают свойством дистрибутивности. |
||
r r r |
|
|
Пусть u,v, w — три произвольных вектора. Смешанное |
||
rrr r |
r r |
r |
произведение uvA(w) |
векторов u,v, A(w) является скалярной |
|
функцией трех векторных аргументов, линейной по всем трем аргументам. Для первых двух это очевидно (так как смешан-
101
Глава VI. Векторное поле
ное произведение обладает свойством линейности для каждого сомножителя), для третьего аргумента в силу соотноше-
ний (6.5) имеем |
|
r r |
|
r r |
rrr |
r rr |
rr |
||
uvA(λw)= uv |
(λA(w))= λuv |
(A(w)). |
||
Линейной по всем аргументам является и скалярная функция
трех векторных аргументов: |
|
rr r rr r r |
r rr r |
uvA(w)+vwA(u )+ wuA(v ).
Эта функция является кососимметричной по всем аргументам, rrr
поэтому она отличается от смешанного произведения uvw только на скалярный множитель. Обозначая этот скалярный
множитель через I1 , получаем |
|
|
|
rrr r rr r r |
r rr r |
rrr |
(6.12) |
uvA(w)+ vwA(u ) |
+ wuA(v ) = I1 |
uvw. |
Так как величина I1 однозначно определяется функцией, стоящей в левой части равенства (6.12), а эта функция, в свою очередь, определяется только линейным оператором A , то скаляр I1 однозначно определяется оператором A и называется пер-
вым инвариантом линейного оператора. Подставляя в равен- |
||||||||
|
|
|
|
|
r |
r |
|
|
ство (6.12) вместо векторов u,v, w орты декартовых координат |
||||||||
r r r |
|
|
|
rrr r |
rrr r |
rrr r |
|
|
i , j, k , |
получаем |
ijA(k )+ jkA(i )+ kiA( j )= I1 , |
или |
|||||
r |
r r r |
r r |
r |
r |
r |
|
|
|
i |
A(i )+ j |
A(j )+ k |
A(k )= I1. |
В этом равенстве слева стоят |
||||
координаты a11 , a22 , a33 векторов (6.7). Получаем выражение первого инварианта в координатах
I1 = a11 + a22 + a33 ,
называемого следом матрицы (6.11).
Смешанное произведение векторов ur, A(vr), Ar(wr) является скалярной трилинейной функцией векторов
Поэтому имеет место равенство |
|
|
|
rr r r r rr r r r r r |
r |
r r |
rrr |
uA(v )A(w)+ vA(w)A(u )+ wA(u )A(v ) = I2 |
uvw, |
||
также r r u,v, w .
102
§ 2. Линейные операторы
так как слева стоит трилинейная кососимметричная функция,
а I2 — скаляр, определяемый этой функцией. Следовательно, r
величина I2 определяется только линейным оператором A и
называется вторым инвариантом линейного оператора. На-
конец, справедливо равенство |
|
|
|
||
|
r r |
r r r r |
rrr |
|
(6.13) |
|
A(u )A(v )A(w) = I3 |
uvw, |
|
||
в |
котором смешанное |
произведение |
векторов |
r |
r r |
A(u ), A(v ), |
|||||
r |
r |
|
|
|
|
A(w), стоящее в левой части равенства (6.13), определяет три-
линейную кососимметричную скалярную функцию от векто- r r r
ров u,v, w . Это функция определяет скаляр I3 , а так как эта трилинейная rфункция однозначно определяется линейным
оператором A , то и скаляр I3 определяется линейным опера- r
тором A . Этот скаляр называют третьим инвариантом ли-
нейного оператора. Подставляя в равенство (6.13) вместо век- |
|||||||
r r |
r r |
|
|
|
r r r r r r |
|
|
торов u,v, w орты i , j, k , получаем |
I3 = A(i )A( j )A(k ). |
Сле- |
|||||
довательно, |
|
|
|
|
|
|
|
|
|
a11 |
a12 |
a13 |
|
|
|
|
|
|
|
||||
|
I3 = |
a21 |
a22 |
a23 |
|
. |
|
|
|
a31 |
a32 |
a33 |
|
|
|
§ 3. Дифференцируемость векторных полей |
|
||||||
Пусть ar(rr) |
— некоторое векторное поле. Обозначим че- |
||||||
рез a приращение вектора поля при переходе из точки rr0 в |
|
точку r : ar = ar(rr)− ar(rr0 ). Положим |
rr = rr − rr0 . |
Определение. Векторное поле a |
называется дифференци- |
руемым в точке rr0 , если его приращение a в этой точке может быть представлено в виде
103
Глава VI. Векторное поле
где |
|
r |
|
ar = A( |
rr)+αr |
|
rr |
|
, |
|
(6.14) |
||||
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
A — некоторый линейный оператор, определённый для |
||||||||||||||
данного поля и для данной точки rr0 , а вектор α зависит от |
rr |
||||||||||||||
и стремится к нулю при r → 0 . Приращение |
a векторного |
||||||||||||||
поля состоит из двух частей: главной части A( |
rr), линейной |
||||||||||||||
относительно приращения |
r независимого переменного, |
и |
|||||||||||||
бесконечно |
малой относительно |
приращения |
r |
величины |
|||||||||||
αr |
|
rr |
|
. Главная часть приращения |
|
a называется дифферен- |
|||||||||
|
|
|
|||||||||||||
циалом векторного поля и обозначается через |
da. |
Следова- |
|||||||||||||
|
|
|
|
r |
r |
r |
|
r |
|
на dr , |
|
|
|
||
тельно, da = A( |
r ), или, заменяя |
|
|
|
|
||||||||||
dar = A(drr).
Определение. Векторное поле называется дифференцируемым, если оно дифференцируемоr в каждой своей точке. Ли-
нейный оператор A называется производной от векторного поля a .
Ограничимся рассмотрением только дифференцируемых векторных полей. Определим для векторного поля производную по направлению. В фиксированной точке M зададим направлением l и проведем через эту точку, в направлении l кривую rr = rr(s). В точке s кривой вектор поля есть ar = ar(s),
где |
ar(s)= ar(rr(s)). |
Придавая параметру s |
приращение |
s , |
|||||||||||||||||
получим приращение вектора поля |
ar = ar(s + s)− ar(s). |
Ис- |
|||||||||||||||||||
пользуя равенство (6.14), |
вычислим предел отношения |
ar |
|||||||||||||||||||
при |
s → 0 . Имеем |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
s |
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
r |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
lim |
ar = lim |
(A( |
rr)+αr( |
|
|
rr |
|
)) |
= |
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
s |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
s→0 |
s |
s |
→0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
r |
rr |
|
r |
|
|
rr |
|
|
r |
rr |
r drr |
|
|
|||||
|
|
|
|
|
|
|
|||||||||||||||
|
= lim |
A |
|
+ lim ±α |
|
|
|
|
|
|
= A lim |
|
= A |
|
|
, |
|||||
|
|
|
|
|
|
|
|
||||||||||||||
|
|
s→0 |
|
s |
s→0 |
|
|
|
|
|
s |
|
|
s→0 |
s |
ds |
|
|
|||
104 |
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
§ 3. Дифференцируемость векторных полей
так |
|
как rr → 0 при s → 0 . Следовательно, αr |
|
rr |
|
→ 0 и |
|||||
|
|
||||||||||
|
rr |
|
→1. Таким образом, получаем |
|
|
|
|
|
|||
|
|
|
|
|
|
|
|||||
|
s |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dar |
r drr |
|
|
|
|
|
|
|
|
|
|
|
= A |
|
. |
(6.15) |
|||
|
|
|
|
ds |
|
||||||
|
|
|
|
ds |
|
|
|
|
|
||
Существование производной доказано и найдено её выражение, из которого видно, что эта производная не зависит от выбора кривой, а зависит только от точки M и от направления l в этой точке, поэтому ее называют производной векторного поля по направлению l . Её значение совпадает со значением
линейного оператора A от орта drr данного направления. ds
Используя формулу (6.15), найдём матрицу, соответст-
вующую линейному оператору A — производной векторного поля. Продифференцируем a в направлении оси Ox , вектор-
ное уравнение которой имеет вид rr = (s − x0 )i + rr0 . В коорди-
натах эта прямая имеет уравнение x = s, y = y0 , z = z0 . |
Следо- |
|||||||
вательно, для поля ar(x, y, z) |
вдоль |
данной прямой |
имеем |
|||||
ar(s)= ar(s, y0 , z0 ), и производная по |
s совпадает с частной |
|||||||
|
r |
(x, y, z): |
∂ar |
|
∂ar |
|
||
производной по x от a |
∂s = |
∂x , и в силу формулы |
||||||
(6.15) |
|
|
r |
|
|
|
|
|
|
|
|
r |
r |
|
|
|
|
|
|
∂∂ax = |
A(i ), |
|
|
(6.16) |
||
так как в данном случае |
∂rr |
r |
|
|
|
|
|
|
|
∂s |
= i . Аналогично |
|
|||||
r |
|
|
|
r |
|
|
|
|
|
r |
r |
|
r |
r |
|
||
∂∂ay = |
A |
(j ), |
∂∂az = |
A(k ). |
(6.17) |
|||
|
|
|
|
|
|
|
|
105 |
Глава VI. Векторное поле
Умножая формулу (6.16) на орт i , получаем слева ∂∂axx , а справа
r |
|
|
|
|
∂ax |
|
||
ir A(iv)= a11 . |
Следовательно, a11 = |
. Умножая же формулу |
||||||
|
||||||||
|
|
|
|
|
|
∂x |
|
|
|
|
r |
r |
r |
r |
|||
(6.16) |
на орт |
j , получим |
∂ay |
= j |
A(i )= a12 и т.д. Из формул |
|||
∂x |
||||||||
|
|
|
|
|
|
|
||
(6.16) |
и (6.17) получаем, что матрица, соответствующая линей- |
|||||||
|
|
r |
|
|
|
|
||
номуоператору A , имеетвид |
|
|
|
|
||||
∂ax∂x∂ay∂x
∂az∂x
∂ax ∂y
∂ay
∂y
∂az ∂y
∂ax |
|
||
∂z |
|
|
|
|
|
||
∂ay |
|
(6.18) |
|
. |
|||
∂z |
|||
|
|
||
∂az |
|
|
|
∂z |
|
|
|
|
|
||
Если через cosα, cos β, |
cosγ |
обозначить координаты ор- |
||||||||||||||||
та данного направления, то производная по направлению от ar |
||||||||||||||||||
в координатах имеет вид |
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
dar |
= ∂ar cosα + |
∂ar cos β + |
∂ar cosγ = |
|
|||||||||||||
|
|
|
||||||||||||||||
|
ds |
|
|
∂x |
|
∂y |
|
∂z |
|
|
|
|
||||||
|
r |
∂a |
x cosα |
+ |
∂a |
x cos β + |
∂a |
|
|
|
+ |
|||||||
= i |
|
|
|
|
x |
cosγ |
||||||||||||
|
|
∂x |
|
∂y |
|
|
|
∂z |
|
|
|
|||||||
|
r |
|
∂a |
y |
|
|
∂a |
y |
|
|
∂a |
y |
|
|
||||
+ j |
|
|
|
cosα + |
|
|
cos β + |
|
|
|
cosγ |
+ |
||||||
|
|
|
∂y |
|
|
|
|
|||||||||||
|
|
|
∂x |
|
|
|
∂z |
|
|
|||||||||
|
r |
|
∂a |
z cosα + |
∂a |
|
|
cos β + |
∂a |
|
|
|
|
|||||
+ k |
|
|
|
|
z |
|
|
z cosγ . |
||||||||||
|
|
|
∂x |
|
∂y |
|
|
∂z |
|
|
|
|||||||
§ 4. Дивергенция векторного поля
Определение. Дивергенцией векторного поля a называют первый инвариант I1 линейного оператора A , являющегося
106
§ 4. Дивергенция векторного поля
производной данного векторного поля. Дивергенция векторного поля ar обозначается diva и называется также расхожде-
нием векторного поля.
Из определения дивергенции следует ее выражение в декартовых координатах. Действительно, в координатах она равна следу матрицы (6.18). Таким образом,
r |
|
∂a |
|
∂ay |
|
∂a |
|
|
diva |
= |
|
x + |
|
+ |
|
z . |
(6.19) |
|
∂y |
|
||||||
|
|
∂x |
|
∂z |
|
|||
Из определения дивергенции следует, что значение diva , найденное по формуле (6.19), не зависит от выбора системы координат. Справедливы следующие свойства дивергенции:
divc = 0, cr = const;
div (ar+br)= divar+ divbr;
r |
r |
ruuuuur |
div (ϕa)=ϕdiva |
+ a gradϕ, |
|
где ar и br — некоторые векторные поля, а ϕ — скалярное по-
ле. Докажем, например, третью из этих формул. Имеем в силу формулы (6.19)
div(ϕar)= ∂(ϕ∂xax ) + ∂(ϕ∂yay )+
=ax ∂∂ϕ + ay ∂∂ϕ + az ∂∂ϕ +ϕ
xy z
∂(ϕaz ) = ∂z
∂∂axx +ϕ ∂∂ayy +ϕ ∂∂azz ,
что и доказывает формулу. Аналогично доказываются остальные равенства.
Дадим интерпретацию дивергенции на гидродинамической модели векторного поля. Рассмотрим течение идеальной жидкости, такое, что в каждой точке вектор a является вектором скорости движущейся частицы. Обозначим через S ориентированную поверхность, т. е. поверхность, в каждой точке которой задан орт нормали n . Для замкнутых поверхностей положительная ориентация выбирается таким образом, чтобы орт
107
Глава VI. Векторное поле
нормали nr был направлен во внешнюю (по отношению к замкнутой поверхности) часть пространства.
Зафиксируем на поверхности S произвольную точку M и рассмотрим в этой точке элемент поверхности dS . Количество жидкости, протекающей за единицу времени через элемент dS , равно (с точностью до бесконечно малых высших
порядков) объему цилиндра с основанием dS |
и образующей |
ar(M ) . Высота h этого цилиндра равна проекции образующей |
|
на нормаль к элементу dS . Следовательно, h = ±a nr , а объем |
|
цилиндра равен абсолютной величине выражения |
|
r |
(6.20) |
dQ = a ndS. |
|
Выражение dQ положительно, если угол между вектором ско-
рости а движущейся жидкости и ортом n — острый; и отрицательно в противном случае. Ориентированный элемент dS поверхности имеет две стороны: положительную, которая видна с конца вектора n , и отрицательную — противоположную. Таким образом, положительное значение величины (6.20) означает, что при своем движении жидкость проходит через поверхность с отрицательнойr стороны на положительную (в направлении орта n ). Отрицательное же значение величины (6.20) означает, что, наоборот, текущая жидкость проходит через поверхность с положительной стороны на отрицательную.
Суммируя все такие элементарные потоки dQ по всей поверхности S , приходим к интегралу
r |
r |
(6.21) |
Q = ∫∫a |
ndS, |
S
который даёт разность между количеством жидкости, протекающей через поверхность в положительном направлении (в направлении нормали) и количеством жидкости, протекающей в обратном направлении. Величина Q называется потоком
векторного поля ar через поверхность S . В координатах (если обозначить через cosα, cos β, cosγ компоненты вектора nr )
интеграл (6.21) имеет вид
108
§ 4. Дивергенция векторного поля
Q = ∫∫(ax cosα +ax cos β + az cosγ )dS.
Обозначим dSr = ndSr . Длина вектора dS равна площади эле-
мента поверхности. Направление вектора dS совпадает с направлением вектора n . При таком обозначении интеграл
(6.21) можно записать в виде Q = ∫∫ar dS.
S
Рассмотрим случай замкнутой ориентированной поверхности S . Если вектор n в каждой точке поверхности направлен во внешнюю часть пространства, то движение жидкости в положительном направлении — в сторону нормали — означает, что жидкость вытекает из области V , ограниченной поверхностью S . Напротив, движение в отрицательном направлении — в сторону, противоположную направлению нормали n , — означает, что жидкость втекает в область V через поверхность
S . Вместе с тем сам поток вектора Q = ∫∫ar dS через замкну-
S
тую поверхность дает разность между количеством жидкости, вытекающей из области V , и количеством жидкости, втекающей в эту область. Если этот интеграл равен нулю, то из области V вытекает ровно столько же жидкости, сколько в нее и втекает. Если же этот интеграл положителен, то это означает, что из области V вытекает жидкости больше, чем втекает в нее. Так как в рассматриваемом примере жидкость несжимаемая, то количество жидкости внутри поверхности S неизменно. Поэтому тот факт, что вытекает ее больше, чем втекает, означает наличие источников внутри области V. Числен-
ную величину Q называют мощностью всех источников
внутри V.
Отрицательное значение интеграла (6.21) также означает, что втекает внутрь области V больше жидкости, чем вытекает из нее, следовательно, внутри V должны быть стоки, поглощающие избыток жидкости.
109
Глава VI. Векторное поле
Отношение VQ , где V — объем, ограниченный поверхно-
стью S , дает среднюю мощность источников, отнесенную к единице объема, ее называют также плотностью источников (при отрицательном Q это будет средняя мощность стоков).
Зафиксируем некоторую точку M пространства и заключим ее внутрь замкнутой поверхности. Предел отношения VQ , ко-
гда поверхность S стягивается в точку M , называют истин-
ной плотностью источников в точке M , он равен диверген-
ции поля ar в точке M . Для доказательства используем теорему Остроградского, дающую следующую интегральную
формулу для любых функций P(x, y, z), |
Q(x, y, z), R(x, y, z), |
||||
имеющих непрерывные частные производные: |
|
|
|||
|
|
∂P + |
∂Q + |
∂R |
|
∫∫(Pdydz +Qdzdx + Rdxdy)=∫∫∫ |
dxdydz . |
||||
S |
V |
∂x |
∂y |
∂z |
|
Применяя формулу Остроградского к случаю, когда функции
P,Q, R совпадают с координатами ax , ay , az векторного поля ar , получим:
|
|
|
∂a |
|
|
∂ay |
|
∂a |
|
∫∫(ax dydz + ay dzdx + az dxdy)=∫∫∫ |
|
x |
+ |
|
+ |
|
z dxdydz . |
||
|
∂y |
|
|||||||
S |
V |
|
∂x |
|
|
∂z |
|||
Интеграл, стоящий в левой части этого равенства, взятый по замкнутой поверхности S , характеризует поток вектора ar через эту поверхность. Под знаком интеграла, стоящего в правой части, взятого по объёму V , содержится дивергенция векторного поля. Если обозначить через dV элемент объема dx dy dz , то правый интеграл может быть записан в виде
∫∫∫divadVr , а вся формула Остроградского переписывается
V
следующим образом:
110
