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

Ilin_Poznyak_-_Matanaliz

.pdf
Скачиваний:
733
Добавлен:
02.06.2015
Размер:
14.78 Mб
Скачать

§ 2

ОСНОВНЫЕ ПОНЯТИЯ И ОПЕРАЦИИ

161

 

Сформулируем следующее определение.

 

 

Оnределе'Н,ие 3.

Век;торное поле р(М) называется

д и ф­

Ф е р е н Ц и р у е м ы м

в т О "l к; е М области [2, если прира­

щение поля ь..р в mO"lne М MO;JfCem быть представлено в сле­

дующей форме:

 

ь..р = АЬ..т + o(lb..rl),

(6.42)

где А - линейный оператор, не зависящий от ь..т (не зависящий

от выбора то"lК;И М').

Соотношение (6.42) мы будем называть условием дифферен­ цируемости поля р(М) в mO"lne М.

Докажем, что если векторное поле р(М) дифференцируемо в точке М, то представление (6.42) для приращения ь..р этого

поля в точке М единственно. Пусть

ь..р = АЬ..т + 01 (Ib..rl) и ь..р = ВЬ..т + 02(Ib..rl)

(6.43)

- два представления приращения ь..р в точке М. Из формул

(6.43) при ь..т # о получаем соотношение

(А - В)е = o(lt.rl)

(6.44)

It.rl

'

в котором е = t.r - единичный вектор, o(lb..rl) = 01 (Ib..rl) -

It.rl

-02(Ib..rl). Так как o(lt.rl) - бесконечно малый вектор при ь..т -+

It.rl

-+ О, а е-произвольный единичный вектор, то из (6.44) следу­ ет, что (А - В)е = О для любого е, т. е. А = В. Единственность представления (6.42) доказана.

Будем говорить, что век;торное поле р(М), заданное в обла­

сти [2, дифференцируемо в этой области, если оно дифферен­

цируемо в na;JfCaou mO"lne области [2.

Введем nонятие производной по направлению для век;торно­

го поля р(М).

Пусть поле р(М) задано в области [2, М -некоторая точка

[2, е - единичный вектор, указывающий направление в точке М.

Пусть далее М' -любая точка из [2, отличная от М и такая, что

--,

вектор ММ коллинеарен вектору е. Расстояние между точка-

ми М и М' обозначим через р.

Если существует предел

1. t.p

lШ­

р--+О Р

(ь..р = р(М') - р(М)), то этот предел называется производ­ ной поля р(М) в mO"lne М по направлению е и обозна"lается

символом др.

де

6 В. А. Ильин и э. г. Позняк, часть II

162

ОСНОВНЫЕ ОПЕРАЦИИ ТЕОРИИ поля

ГЛ.6

 

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

lim 6..р.

 

 

др =

(6.45)

 

де

р-+О р

 

Справедливо следующее утверждение. Пусть поле р(м) диффе­

ренцируемо в тО'Чо'Х:е М области S1. Тогда nроизводная др поляр

де

в этой тО'Чо'Х:е по любому направлению е существует и может

быть найдена по формуле

др = Ае

(6.46)

де

'

где А - линейный оператор, определенный соотношением (6.42).

Докажем это утверждение. Пусть е - любое фиксированное

направление и пусть точка м' берется так, что вектор tlr = ре и Itlrl = р. Подставляя это значение tlr в соотношение (6.42) и

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

tlp =

рАе + о(р).

 

Отсюда получаем формулу

Ае + о(р).

 

6..р =

(6.47)

рр

Из соотношений (6.45) и (6.47) вытекает формула (6.46). Утвер­

ждение доказано.

Пусть р(М) - дифференцируемое в точке М области s1 поле.

Тогда

tlp = Atlr + o(ltlrl).

Найдем матрицу линейного оператора А для случая орто­ нормированного базиса i, j, k. Мы будем считать, что с этим ба­

зисом связана декартова прямоугольная система коорди­

нат Oxyz.

Обозначим через Р, Q и R координаты векторного поля р(М) в базисе i, j, k. Очевидно, согласно формуле (6.46),

др =

др = Ai

др =

др =

А'

др =

др = Ak.

ai

дх

'aj

ду

З,

ak

az

Из этих формул и из соотношений (6.26) для матрицы коэффици­

ентов линейного оператора в ортонормированном базисе i, j, k

о

следует, что матрица А рассматриваемого оператора А имеет

вид

дР

дР

дР

 

 

 

 

дх

ду

az

 

о

aQ

aQ

aQ

 

А= дх

ду

az

(6.48)

 

aR

aR

aR

 

 

дх

ду

az

 

Введем понятие дивергенции и ротора дифференцируемого

в области s1 векторного поля р(М), т. е. такого поля, прира-

§ 2

ОСНОВНЫЕ ПОНЯТИЯ И ОПЕРАЦИИ

163

щение 6.р которого в каждой точке М области s1 может быть

представлено в виде

6.р = А6.т + o(l6.rl),

причем оператор А, вообще говоря, меняется при переходе от одной точки области s1 к другой. Иными словами, оператор за­ висит от точки М и не зависит, конечно, от 6.Т.

Назовем диверге'Н'Циеu и ротором поля р(М) в то'Чке М

области s1 диверге'Н'Цшо и ротор ли'Неu'Ного оператора А. Таким образом, по оnределе'Ншо

div р = div А, rot р = rot А.

(6.49)

3 а м е ч а н и е. При наших предположениях о дифференци­

руемости поля р(М) в области s1 дивергенция div р и ротор rot р

определены в каждой точке S1. Поскольку эти объекты являют­

ся инвариантами (не зависят от выбора базиса), то, очевидно, div р представляет собой скалярное поле, а rot р - векторное по­

ле в области S1.

Найдем выражения дивергенции, ротора и производной по

направлению для дифференцируемого векторного поля р(М),

считая, что в пространстве выбран ортонормированный базис i, j, k, с которым связана декартова прямоугольная система координат Oxyz. При этом, как и выше, будем считать, что поле

р(М) имеет координаты Р, Q, R в базисе i, j, k.

о

Так как матрица А линейного оператора А определяется в

рассматриваемом случае соотношением (6.48) и по определению divp = div А, rotp = rotA (см. (6.49)), то, согласно формулам

(6.27) и (6.29), получим

 

divp =

дР + дQ + дR

 

(6.50)

 

 

дх

ду

дz'

 

 

rotp = (дR _

дQ)i + (дР _

дR)j + Q _

дР)k.

(6.51)

ду

дz

дz

дх

дх

ду

 

Для вычисления производной векторного поля р(М) по направ­ лению е воспользуемся формулой (6.46) и свойствами линейного

оператора.

Пусть е = i cos СУ+j cos (3 +k cos r

1). Тогда, согласно (6.46),

получим

 

 

 

 

 

 

др = Ае =

cos СУ Ai + cos (3 Aj + cos r Ak =

 

 

де

 

 

cos СУ др + cos (3 др + cos r

 

= cos СУ др + cos (3 др + cos r др =

др.

m

~

~

 

&

~

&

1) Так как

е - единичный

вектор,

то

его

координаты имеют

вид

{COSct, cos{3, cos,}, где а, {3 и ,-углы, которые составляет этот вектор

с осями Ох, Оу и Oz соответственно.

6*

164 ОСНОВНЫЕ ОПЕРАЦИИ ТЕОРИИ поля ГЛ.6

Таким образом, производная др может быть вычислена либо по

де

формуле

 

 

др =

cos СУ др + cos (3 др + cos, дР,

(6.52)

 

 

де

 

дх

 

ду

az

 

 

либо, учитывая, что Р, Q, R-координаты р(М), по формуле

др

дР

дР

cos (3

 

дР).

+

 

 

 

(-

cos СУ + -

 

+ -

cos,

1,

 

 

де

дх

ду

 

 

az

 

 

 

 

 

 

 

+ (aQcos СУ + aQ cos (3 + aQ cos ,) j

+

 

 

 

дх

 

 

 

ду

 

 

az

 

 

 

 

 

 

aR

 

aR

 

aR)

(6.53)

 

 

 

+ (-

COSCY + -

cos(3 + -

cos, k.

 

 

 

 

дх

 

 

ду

 

az

 

 

4. Повторные операции теории поля. Будем считать,

что в области n евклидова пространства Е3 заданы скалярное

поле u (М) класса с2 1) и векторное поле р(М) класса с2 .

При этих предположениях grad u представляет собой диффе­

ренцируемое векторное поле в п, div р -

дифференцируемое ска­

лярное поле,

а rot р - дифференцируемое векторное поле. По­

этому возможны следующие повторные операции:

 

rot grad и,

div grad и,

grad div р,

div rot р,

rot rot р.

Докажем, что

 

 

 

 

 

rot grad u =

о, div rot р = о.

(6.54)

Для доказательства вычислим rot grad u и div rot р в декарто­

вой прямоугольной

системе координат. Так как в этом случае

координаты grad u равны

дu

дu

дu

 

- , - , - , то на основании формулы

(6.51) получим

 

дх

ду

az

 

 

 

 

 

 

rotgradu = ( д2u

д2u

)i + ( д2u

д2u )j +

ayaz -

azay

 

 

azax -

axaz

+ 2u _ д2u)k = о.

дхду дудх

Таким образом, первое из равенств (6.54) справедливо для декар­

товой системы координат. В силу инвариантности выражения rot grad и, первое из равенств (6.54) доказано. Перейдем к до­ казательству второго равенства (6.54). Обратимся опять к де­ картовой системе координат. В этой системе, согласно (6.51),

1) Функция принадлежит классу C k в области !1, если все ее частные про­

изводные порядка k непрерывны.

§ 3 ВЫРАЖЕНИЕ ОПЕРАЦИЙ В КРИВОЛИНЕЙНЫХ КООРДИНАТАХ 165

векторное поле rot р имеет координаты (aR _ aQ ),

(дР _ aR),

( aQ _

ду

az

az

дх

дР), где Р, Q, R-координаты векторар. Согласно (6.50)

дх

ду

 

 

 

дивергенция векторного поля rot р в декартовой прямоугольной

системе координат равна сумме производных компонент этого

поля по одноименным координатам. Таким образом,

divrotp = .!!...- (aR _ aQ ) + .!!...- (дР _aR) + ~ (aQ _

дР) = о.

дх ду

az

ду

az

дх

az

дх

ду

Таким образом, второе из равенств (6.54) справедливо для де­

картовой системы координат. В силу инвариантности выраже­

ния divrotp, второе из равенств (6.54) справедливо в любой си­

стеме координат.

Одной из основных повторных операций теории поля явля­ ется операция div grad и. Кратко эту операцию обозначают flu,

причем символ fl обычно называют оператором Лапласа 1) . Та­

ким образом,

flu = div grad и.

(6.55)

Вычислим оператор Лапласа в декартовой прямоугольной системе координат. В такой системе векторное поле grad u име-

ет координаты

дu

дu

дu

О

бращаясь к выражению

(

6.50

)

для

- , - , - .

 

 

 

 

дх

ду

az

 

 

 

 

 

 

 

 

дивергенции векторного поля, получим

 

 

 

 

 

 

flu =

д2 u

д2 u

д2u

 

(6.56)

 

 

-

 

+ -

+ -.

 

 

 

 

 

дх

2

ду2

az2

 

 

 

 

Повторные операции grad div р и rot rot р связаны соотношением

rotrotp = graddivp - flp,

(6.57)

где flp представляет собой вектор, координаты которого в бази­ се i, j, k равны flP, flQ, flR (Р, Q, R-координаты векторного

поля р в базисе i, j, k). В справедливости соотношения (6.57)

читатель легко убедится самостоятельно.

§3. Выражение основных операций теории поля

вкриволинейных координатах

1.Криволинейные координаты. Пусть Sl- область ев­

клидова пространства Е3 ; Х,-У, z - декартовы координаты в этом

пространстве. Пусть далее, Sl- область евклидова пространства

Е3; х1 , х2 , х3 - декартовы координаты в jj;з.

1) п. с. Лаплас - выдающийся французский астроном, математик и фи­

зик (1749-1827).

166

ОСНОВНЫЕ ОПЕРАЦИИ ТЕОРИИ поля

ГЛ.6

Рассмотрим взаим~о однозначное и взаимно непрерывное

отображение области Q на область Q, которое осуществляется посредством функций

с помощью указанного отображения в области Q вводятся кри­

волинейные координаты х1 , х2 , хз . Смысл этого наименования

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

точке М(х, у, z) области Q сопоставляются три числа х1 , х2 , хз . Более точно, точка М определяется тройкой чисел х1 , х2 , хз .

Этим объясняется наименование «координаты» точки М для

чисел х1 , х2 , хз . Во-вторых, если в правых частях соотношений (6.58) фиксированы две какие-либо координаты, например х2

и хз , то при переменном х1 эти соотношения определяют в обла­

сти Q некоторую линию, отличную, вообще говоря, от прямой. Эту линию естественно называть nоор­

дин.аmн.оЙ лин.иеЙ х1 , подчеркивая тем

самым, что в точках указанной линии

меняется лишь координата х1 . В полной

 

аналогии определяются

координатные

 

линии х2 и хз . Вообще говоря, коорди­

 

натные линии х1 , х2 и хз

не будут пря­

 

мыми. Этим И объясняется термин «кри­

 

волинейные координаты».

 

 

Мы выяснили, что через каждую

 

точку М области Q проходят три коор-

Рис. 6.2

динатные линии хl, х2 , хз

(рис. 6.2). По-

строим в точке М базис Ti, T i , естественным образом связанный

с координатными линиями, проходящими через эту точку. При

этом мы воспользуемся соотношениями (6.58). Очевидно, про-

изводные

дх

ду

az

вычисленные в точке М, представляют

- , - , - ,

 

ax1

ax1

ax1

 

собой координаты вектора касательной к линии х1 в этой точке.

Мы обозначим, этот вектор через Тl. Аналогичным способом мы

строим векторы Т2 и Тз касательных к линиям х2 и хз

соответ­

ственно. Таким образом,

 

 

 

 

 

{

дх

ду

az}

k =

1, 2, 3.

(6.59)

axk '

axk '

axk '

Tk =

 

 

 

Для того чтобы векторы Тl, Т2, Тз образовали базис, нуж­ но потребовать некомпланарности этих векторов. Достаточным условием для выполнения этого требования является, очевидно,

§ з ВЫРАЖЕНИЕ ОПЕРАЦИЙ В КРИВОЛИНЕЙНЫХ КООРДИНАТАХ 167

условие необращения в нуль якобиана

 

дх

ду

дz

 

дх1

дх1

дх1

Г(х, у, z)

дх

ду

дz

D(x 1 , х2 , х3)

дх2

дх2

дх2 '

 

дх

ду

дz

 

дх3

дх3

дх3

ибо этот якобиан равен смешанному произведению векторов Tl, Т2, тз· С помощью построенного базиса Tl, Т2, Тз стандартным

образом строится взаимный базис T 1 , т2 , тз.

Итак, если в области n введены криволинейные координаты x 1 ,

х2 , хЗ , то с каждой точкой М этой области естественным обра­ зом связываются базисные векторы Ti, T i . Рассмотрим примеры.

1О. Цuлu'Ндрu"tесnал система nоорди'Нат. Эта система коор­

динат вводится с помощью соотношений

х = р cos ср, у = р SlП ср,

z = z.

(6.60)

Таким образом, x 1 = р, х2 = ср, хз = z.

Известно, что указанные

координаты р, ср, z (или, что то же,

x 1 , х2 , хз ) изменяются в

следующих пределах 1):

 

О ~ р ~ +00, О ~ t.p < 27Г,

- 00 < z < 00.

Эти неравенства определяют в евклидовом пространстве коор­

динатами р, ср, z (или x 1 , х2, хз) бесконечную область S1, jj;з с

изображенную на рис. 6.3. Мы можем, следовательно, рассмат-

z

Линия <р

_--Г--.J.

O}-----j___

у

х

Линия

р

z

I

I

I

Q I 2n

/1- I

Рис. 6.3

ривать введение цилиндрических координат в евклидовом про­

странстве Ез как результат отображения области s1 простран­ ства jj;з в пространство Ез с помощью формул (6.60).

1) См. вып. 3 «Аналитическая геометрия» настоящего курса.

168

ОСНОВНЫЕ ОПЕРАЦИИ ТЕОРИИ поля

ГЛ.6

 

Очевидно, координатные линии р (или линии x 1 )

представ­

ляют собой прямые, проходящие через ось Oz, перпендикулярно

этой оси, координатные линии t.p (линии х2) - окружности С цент­

рами на оси Oz, плоскости которых параллельны плоскости Оху.

Координатные линии z (линии хз ) - прямые, параллельные оси Oz (см. рис. 6.3). Найдем векторы Tl, Т2, Тз и T1, т2 , ТЗ. Имеем

 

Tl=

{ дх

ду

az}

 

{

.

}

 

- , - , - =

cost.p,Slnt.p,O,

 

 

др

др

др

 

{ .

}

 

Т2 =

{ дх

ду

az}

=

 

- , - , -

-рsшt.p,

pcost.p, О ,

 

 

arp

arp

arp

 

{

}

 

 

ТЗ =

{ дх

ду

az}

=

 

 

- , - , -

О,

0,1 .

 

 

 

az

az

az

 

 

 

 

Подчеркнем, что выражения в фигурных скобках представ­

ляют

собой декартовы

координаты базисных векторов Tl, Т2

и ТЗ.

Непосредственно

можно

убедиться,

что базис Tl, Т2, ТЗ

ортогональный. Для вычисления «взаимного базиса воспользу­

емся формулами, приведенными в п. 1 § 1 этой главы. Имеем

 

 

 

cos t.p

 

sшt.p

О

 

 

 

 

Tl Т2ТЗ =

Р sш t.p

р cos t.p

О =

р,

 

 

 

 

 

О

 

О

 

1

 

 

 

 

 

[Т2ТЗ] =

{pcost.p,

psint.p,

О},

 

 

 

 

[ТЗТ1]

=

{-sint.p, cost.p, О},

 

 

 

 

[TIT2] =

{О, О,

р}.

 

 

 

 

 

Поэтому

 

 

 

 

 

 

 

 

 

 

 

 

Т

1

[Т2ТЗ]

 

{

cos

ср,

ср, О

}

,

 

 

 

= --- =

 

Sln

 

 

 

 

 

ТIТ2ТЗ

{1.

 

 

 

 

 

 

 

2

[ТЗТ1]

 

 

~ cos ср,

О},

 

Т

 

= --- =

 

-- Slnt.p,

 

 

 

ТIТ2ТЗ

 

р

 

 

 

 

 

 

 

ТЗ

= [ТIТ2]

=

{О, О, 1}.

 

 

 

 

 

 

 

ТIТ2ТЗ

 

 

 

 

 

 

 

 

 

20. Сфери'Ческа.я систе.м.а координат. Эта система координат

вводится с помощью соотношений

 

 

 

 

 

 

х = psinf1cost.p,

у = psinf1sint.p,

z =

pcosf1.

(6.61)

Таким образом, x 1 = р, х2

= ср, хз

= е. Известно, что указанные

координаты р, ср, f1 (или, что то же, xl, х2 , хз ) изменяются в

следующих пределах:

О ~ Р < +00, О ~ t.p < 2п, О ~ f1 ~ п.

(6.62)

§ З ВЫРАЖЕНИЕ ОПЕРАЦИЙ В КРИВОЛИНЕЙНЫХ КООРДИНАТАХ 169

Неравенства (6.62) определяют в евклидовом пространстве jijЗ с координатами р, ср, () (или xl, х2 , хЗ) бесконечную область S1,

изображенную на рис. 6.4. Мы можем поэтому рассматривать

введение сферических координат в евклидовом пространстве ЕЗ

Рис. 6.4

как результат отображения области П пространства jijЗ в про­

странство ЕЗ с помощью формул (6.61).

Очевидно, координатные линии р (линии х1 ) представляют со­

бой лучи, выходящие из начала координат; координатные линии t.p

(линии х2 ) - окружности с центрами на оси ОZ, плоскости кото­

рых параллельны плоскости Оху, координатные линии () (линии

хЗ ) - полуокружности, центры которых находятся в начале коор­

динат и плоскости которых проходят через ось Oz (см. рис. 6.4).

Найдем векторы Тl, Т2, Тз и Tl, т2 , ТЗ. Имеем

Тl =

{

sin()cost.p,

sin()sint.p,

cos()

},

Т2

=

{

-psin()sint.p, psin() cos ср,

О

},

ТЗ

=

cos () cos ср,

р cos () sin ср,

- р sin ()

}.

Непосредственно можно убедиться, что базис Тl, Т2, ТЗ ортого­

нальный. Для вычисления взаимного базиса воспользуемся фор­

мулами, приведенными в п. 1 § 1 этой главы. Имеем

 

sin () cos t.p

sin () sin t.p

cos ()

 

 

ТI Т2 ТЗ

= - Р sin () sin t.p

р sin () cos t.p

О

= _р2 sin(),

 

 

р cos () cos t.p

р cos () sin t.p

sin ()

 

 

[Т2ТЗ]={

_р2 sin2 () cos ср,

- р2 sin2 () sin ср,

- р2 sin () cos ()

},

[ТЗТl]={

psint.p,

р cos ср,

О

},

[ТIТ2]={

- Р cos () sin () cos ср, - р cos () sin () sin ср,

р sin2 ()

}.

170

 

ОСНОВНЫЕ ОПЕРАЦИИ ТЕОРИИ поля

 

ГЛ.6

Поэтому

 

 

 

 

 

 

 

 

 

 

 

т1

=

[Т2ТЗ]

=

{

sin () cos ср,

sin () sin ср,

cos ()

},

 

 

ТIТ2ТЗ

 

 

 

 

 

 

 

 

т2

=

[ТЗТ1]

=

{

 

1 sin

 

1 cos

о}

ТЗ

 

ТIТ2ТЗ

 

 

 

р sin 8 '

-

Р sin 8 '

- Sln [7

,

=

[ТI

Т

2]

=

{

- cos [7 cos ср,

cos [7 Sln ср, -

 

 

 

 

 

 

1

IJ

1

IJ·

1 . 1J

}

 

 

ТIТ2ТЗ

 

 

Р

 

Р

 

Р

 

30. Ортого1-t.аЛ'Ь1-t.ал nрuволu1-t.еU1-t.ал система nоординат.

Криволинейную систему координат мы будем называть орто­

го1-t.ал'Ь1-t.оU, если в каждой точке области n базис Ti, определя­

емый равенством (6.59), является ортогональным. Только что

рассмотренные цилиндрическая и сферическая системы коорди­ нат представляют собой примеры ортогональных криволиней­

ных координат.

Получим выражение для векторов T i взаимного базиса для

случая ортогональной системы координат.

Введем следующие обозначения:

Н1

= I ll,

Н2

= I 21, НЗ = Iтзl·

 

T

 

T

Величины Н1 , Н2 , НЗ обычно называются nоэффu'Цuе1-t.тамu или nараметрами Ламе 1).

Так как система координат ортогональная и тройка векторов

Тl, Т2, Тз правая, то

Используя эти соотношения и формулы, выражающие векторы

взаимного базиса через векторы тi (см. п. 1 § 1 этой главы),

получим

1

1

2

1

з

1

т = Н2Т1, Т = Н2Т2, Т = Н2Т1.

1

2

3

2. Выражение градиента и производной по направ­ лению для скалярного поля в криволинейных коорди­

натах. Пусть в области П, в которой введены криволинейные

координаты х1 , х2 , хз , задано дифференцируемое скалярное по­

ле и(М). При этих условиях grad u определен в каждой точке n и в каждой точке n по любому направлению е может быть

дu

вычислена производная - . Как градиент grad и, так и произ-

де

водную по направлению в данной точке М мы будем относить

к базису Ti, T i в этой точке, построение которого описано в пре­

дыдущем пункте.

1О. выlаJIcе1-t.uеe градиента сnаллр1-t.ого nолл в nрuволu1-t.еu1-t.ыlx nоординатах. После введения в области n криволинейных ко-

1) г. Ламе - французский математик (1795-1870).

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