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

Методы теории функций в геометрии. Учебное пособие

.pdf
Скачиваний:
0
Добавлен:
12.08.2026
Размер:
989 Кб
Скачать

§ 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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

s0

s

s

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

rr

 

r

 

 

rr

 

 

r

rr

r drr

 

 

 

 

 

 

 

 

 

 

= lim

A

 

+ lim ±α

 

 

 

 

 

 

= A lim

 

= A

 

 

,

 

 

 

 

 

 

 

 

 

 

s0

 

s

s0

 

 

 

 

 

s

 

 

s0

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 , имеетвид

 

 

 

 

axxayx

azx

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

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]