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

Ilin_Poznyak_-_Matanaliz

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

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

ординат x 1 , х2 , хЗ скалярное поле и будет, очевидно, функцией

переменных x 1 , х2 , хЗ :

и = и(х\ х2, хЗ).

Эта функция может рассматриваться как результат с;:перпози­

ции функции и(х, у, z) переменных х, у, z и функций (6.58). По­

 

 

 

 

 

 

 

 

 

ди

 

 

 

 

 

 

этому для вычисления производных - . мы можем применить

 

 

 

 

 

 

 

 

 

дх'

 

 

 

 

 

 

правило дифференцирования сложной функции. Обозначая

ди

 

через Ui, получим

 

 

 

 

+ ди .!!JL + ди ~.

 

 

дх'

 

и. =

ди дх

 

 

 

 

 

 

 

 

 

 

(6.63)

 

 

 

 

Z

дх axi

ду axi

az дх;

 

 

 

 

 

 

ди

ди

ди

 

 

 

 

 

d

и в

б

.. k

,

Так как - , - , -

 

-координаты вектора gra

 

 

азисе~,),

 

дх

ду

az

~

О

 

дх

ду

az

 

 

 

 

 

 

 

 

 

 

xyz,

 

 

 

 

 

 

связанном с системои

 

а - . , - . , - . -

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

 

 

 

 

 

 

 

дх'

дх'

дх'

 

 

 

 

 

 

ра Ti, то, очевидно, соотношение (6.63) может быть переписано

в следующей форме:

 

Ui = Ti grad и.

(6.64)

Используя формулы Гиббса (см. формулы (6.6)) для вектора grad и и формулы (6.64), получим

gradu = (TigradU)Ti = UiTi.

Итак, градиент скалярного поля и в криволинейных координа­

тах выражается следующим образом:

grad и = Ui Ti (Ui = ::J. (6.65)

На практике часто встречается случай ортогональной криволи­

нейной системы координат. В предыдущем пункте мы получили

(см. п. 30 предыдущего пункта) выражение для векторов T i вза­

имного базиса для ортогональной системы. Используя эти вы­

ражения и формулу (6.65), найдем для grad и в ортогональных

координатах следующую формулу:

1 ди

1 ди

1 ди

(6.66)

grad и = --- Tl + ---Т2

+ ---ТЗ.

Hf ax 1

H~ дх2

Hl дхЗ

 

Наряду с ортогональным базисом Ti рассматривают ортонор­

мированный базис ei = Ti/ Hi. Легко

видеть, что в

базисе ei

выражение для grad и имеет вид

 

 

 

1 ди

1 ди

1 ди

 

(6.67)

grad и = --- el +

---е2 + ---ез.

 

H 1 ax 1

Н2 дх2

Нз дхЗ

 

 

20. выlаJIcениеe nрои360дно'Й С'/ИЛЛРНО20 nолл и(М)

по на­

nра6леНU10 е 6 'кри60лине'йныlx 'Координатах. Пусть ei

-

контра­

вариантные координаты единичного вектора е в базисе Ti, так

что

172

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

ГЛ.6

Для производной дu мы получили в п. 2 § 2 этой главы следую­

де

щую формулу:

-дu = е· gгаdи

де

(см. формулу (6.40)). Подставляя в эту формулу выражение

для е в базисе ri и формулу (6.65) для grad и, получим

-ддеu = k rk)(Иiri ) = еk Иi(rkri ) = еk Иiдki = Иiеi .

Таким образом, производная скалярного поля и по направле­

нию е выражается в криволинейных координатах следующим

образом:

-дu = Иiеi .

(6.68)

де

 

3. Выражение дивергенции, ротора и производной по

направлению для векторного поля в криволинейных ко­

ординатах. Пусть в области П, в которой введены криволи­ нейные координаты, задано дифференцируемое векторное поле

р(М). При этих условиях дивергенция и ротор поля р опреде­

лены в каждой точке области П, и в каждой точке П по любому

направлению е может быть вычислена производная др. Дивер-

де генцию, ротор и производную по направлению в данной точ-

ке М мы будем относить к базису ri, r i в этой точке.

1О. Въtра;JfCение дивергенции векторного поля в криволинеiJ.­

ныlx координатах. После введения в области П криволинейных

координат х1 , х2 , х3 векторное поле р будет, очевидно, функцией переменных х1 , х2 , х3 :

р = р(хl, х2 , х3).

Эта функция может рассматриваться как результат суперпози­

ции функции р(х, У, z) и функций (6.58). Поэтому для вычисле-

ния производных др мы можем применить правило дифферен-

дх'

 

 

 

 

цирования сложной функции. Обозначая др через Pi'

получим

 

 

 

 

дх'

 

Р

=

др дх

+ др .!!.JL + др ~.

(6.69)

 

Z

дх axi

ду axi

az axi

 

Так как др = Ai,

др

= Aj,

др = Ak, где А-линейный опе-

дх

ду

 

az

АЬ..т + o(lb..rl)

(см. п. 3

ратор, определенный равенством Ь..Р =

§ 2 этой главы), то из соотношений (6.69) и свойств линейного

оператора получим

Pi =

дх.

ду.

az k)

=

А

ri·

(6.70)

А( - 1,

+ -) + -.

 

 

дх'

дх'

дх'

 

 

 

 

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

По определению divp = div А = riAri. Поэтому, согласно фор­

муле (6.70), в криволинейной системе координат дивергенция векторного поля р(М) может быть вычислена по формуле

divp = rip;

i

-

дР)

.

(6.71)

"

 

-

дх;

 

Найдем выражение для дивергенции для случая ортогональной криволинейной системы координат. Используя выражение для

векторов r i взаимного базиса для ортогональных криволиней­

ных координат и формулу (6.71), получим

divp = ~~PlTl + ~~P2T2 + ~lрзrз (Pi = :~). (6.72)

Формула (6.72) записывается также и другим способом. Обозна­

чим через pi координаты поля р в ортонормированном базисе

ei = ~ 1). Тогда после ряда преобразований выражение (6.72)

Н;

для div р примет следующий вид:

divp=

1

[д(р1Н2НЗ) +д(р2НзН1) +д(РЗН1Н2)].

(6.73)

 

Н1Н2Нз

дх1

 

 

дх2

 

дхЗ

 

20. выlаJIcениеe ротора векторного поля в криволинеuныx ко­

ординатах. По определению rot р =

rot А =

[r iAri]. Поэтому,

согласно формуле (6.70),

получим

 

 

дР)

 

 

 

 

rotp =

i

Pi]

(

Pi

(6.74)

 

 

 

= дхi

Найдем выражение ротора в ортогональной криволинейной си­

стеме координат. Используя выражение для векторов r i взаим­

ного базиса для ортогональной системы и формулу (6.74), полу-

чим

1

[ТlРl] +

1

[Т2Р2] +

1

[rзрз]

)

(6.75)

rotp =

- 2

- 2

- 2

Pi =

дx~ .

 

H 1

 

Н2

 

Нз

 

 

 

 

 

В ортонормированном базисе ei =

т·

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

--.::..

 

 

 

 

 

 

Н;

 

 

 

 

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

 

 

 

 

 

 

 

 

{ _1_ [д(РЗНЗ) _

д(р2Н2)]

_1_ [д(Р1Нl) _

д(РЗНЗ)]

Н2Нз

дх2

 

дхЗ '

НЗН1

дхЗ

 

дх1

'

 

 

_1_ [д(Р2Н2)

_ д(Р1Нl)]}.

 

 

(6.76)

 

 

H 1H 2

дх1

 

дх2

 

 

 

30. выаJIcениеe для nроизводноu векторного поля по направ­ лению в криволинеuныx координатах. Воспользуемся формулой

др = Ае

(6.46)

де

'

 

1) в правой части этой формулы

суммирование по индексу i

не произво­

дится.

174

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

ГЛ.6

полученной нами в п. 3 § 2 этой главы. Пусть е = eiri. Тогда из

формулы (6.46) и свойств линейного оператора получим

др .

-= eLAri·

де

Так как Ari = Pi' где Pi = дР, то для производной векторного

дх'

поля Р по направлению е получим следующее выражение:

др

.

(6.77)

-

= eZPi.

де

4. Выражение оператора Лапласа в криволинейных

ортогональных координатах. Мы определили оператор Ла­

пласа tlu как повторную операцию div grad и. Используя выра­

жения (6.67) и (6.73) для градиента и дивергенции в криволиней­

ных ортогональных координатах, мы получим выражение для

оператора Лапласа.

В рассматриваемом случае векторным полем Р, дивергенцию которого нужно вычислить, является поле grad и. Подставляя

(6.67) в (6.73), получим

tlu = 1 [~(Н2НЗ~)+~(НЗНl~)+~(НlН2~)].

Н1 Н2НЗ ax 1 H 1 ax 1 дх2 Н2 дх2 дхЗ Нз дхЗ

(6.78)

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

линдрической и сферической системах координат.

1О. Цuлuндрu'Чес'Ка.я система 'Координат. В силу результатов в 10 п. 1 § 3 параметры Ламе для цилиндрических координат

имеют вид

Н1 = 1, Н2 = р, Нз = 1.

В таком случае из формул (6.67), (6.73), (6.76) и (6.78) вытекают

следующие равенства:

ди

1 ди

ди

grad u = -ер + --еер + - ez ,

др

р д<р

az

d ·

-

1 д (

Р

р) + 1 дР,!,

+ дР"

 

lVP -

--

 

Р

---

--,

 

rot Р = (~apz _ дР'!') е

р др

 

 

 

Р д<р

az

 

+ (дР

 

aPz) е

+ (~д(PP'!') _~ дР ) e z ,

 

 

 

 

р

_

 

 

 

 

р

р д<р az

Р

 

az

 

др

ер

р

др

р д<р

tlu =

~~ (рди) + ~ д2и + д2и.

 

 

р др

др

 

р2 д<р2 az 2

 

20. Сферu'Чес'Ка.я система 'Координат. В

этом

случае пара­

метры Ламе имеют вид

 

 

 

 

 

 

 

 

 

Н1 =

1,

Н2 =

psine,

НЗ =

р.

 

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

Следовательно,

 

 

ди

1 ди

1 ди

gгаdи = -ер + ---еер + --ее,

др

р sin е дер

р де

 

 

·

1

д (

Р

2р)

 

 

1

дРср

1 д (. В?)

 

 

d lVP = --

 

р

+ -.--- + -.-- Sln

е,

 

 

 

р2 др

 

 

 

 

рsше

дер

рsше де

 

 

го

t

-

1 (д(Siперср)

-

дРВ)

ер

+

 

 

 

 

Р - psine

де

 

 

дер

 

 

 

 

 

 

 

 

 

 

 

 

р _

~ д(РРср))ее + (~д(PPB) _

 

 

 

 

+(_l_дР

 

~ дРР)е

 

 

 

 

psine

дер

 

р

 

др

р

др

р де ер,

 

 

!:::.И =

~~ (р2 ди) + _l_~ (siпвдu) +

1

д2 и.

 

 

 

р2 др

др

 

р2 sin е де

де

р2 sin 2 е дер2

В заключение этой главы приведем сводку формул, связы­

вающих операции взятия градиента, дивергенции и ротора с ал­

гебраическими операциями:

10.gгаd(и ± v) = grad и ± grad v.

20.grad(и . v) = и grad v + v grad

и)

v grad u - u grad v

grad ( -

= ---"------"'---

v

v 2

div(p ± q) = div р ± div q. div( ир) = р grad и + и div р. div[pq] = q rotp - prot q. rot[p ± q] = rot р ± rot q.

гоt(ир) = игоtр - [рgгаdи].

и.

(v#O).

в справедливости этих формул читатель легко убедится са­

мостоятельно.

3 а к л ю ч и т е л ь н ы е з а м е ч а н и я . В этой главе мы по­

знакомились с основными операциями теории поля. При этом мы не опирались на какие-либо физические представления, по­

скольку нашей целью являлось построение математического ап­

парата теории. В следующей главе мы получим ряд важных

интегральных соотношений, связывающих некоторые операции

теории поля. Эти соотношения позволят нам указать физиче­

скую интерпретацию понятий и операций, введенных в настоя­ щей главе.

ГЛАВА 7

ФОРМУЛЫ ГРИНА, СТОКСА И

ОСТРОГРАДСКОГО

В этой главе мы получим важные формулы, играющие боль­

шую роль в различных приложениях и, в частности, в теории

поля. Эти формулы в определенном смысле представляют собой

обобщения на многомерный случай формулы Ньютона--Лейбница

для одномерных интегралов.

§1. Формула Грина 1)

1.Формулировка основной теоремы. Пусть D - конеч­

ная, вообще говоря, многосвязная область на плоскости Оху с

кусочно-гладкой границей L 2). Область D с присоединенной

границей L мы будем обозначать D. Справедлива следующая

ос'Нов'Нля теорема.

Теоре,м,а 7.1. Пусть фу'Н'Кции Р(х, у) и Q(x, у) 'Неnрерыв­

'НЫ в D и имеют 'Неnрерыв'Ные 'Част'Ные nроизвод'Ные первого nоряд'Ка в D. Если существуют 'Несобстве'Н'Ные и'Нтегралы по

области D от 'Ка;ждой из 'Част'Ных nроизвод'Ных фу'Н'Кций Р(х, у)

и Q(x, у) 3),

то справедливо соот'Ноше'Ние

 

JJ(~; -~:)dxdy= !PdX+QdY,

(7.1)

 

D

L

 

'Называемое

фор м у л о й

Гр и 'Н а. При этом стоящий в пра­

вой 'Части (7.1) и'Нтеграл представляет собой сумму и'Нтегра­

лов по связ'ныlM 'Комnо'Не'Нтам гра'Ницъ! L, 'На 'которыlx у'Каза'Но

та'Кое 'Наnравле'Ние обхода, при 'Котором область D остается

слева.

1)Дж. Гринанглийский математик (1793-1841).

2)Граница L называется кусочно-гладкой, если она составлена ИЗ конечно­

го числа гладких кривых. Если граница L состоит ИЗ конечного числа замк­ нутых кусочно-гладких кривых Li, ТО связную область D обычно называют М'Н,О20св.яз'Н,оU, а кривые L i называют св.яз'Н,ымu ",омnо'Н,е'Н,mамu границы.

3) Так как частные производные функций Р(х, у) и Q(x, у) существуют

лишь в открытой области D, то упомянутые интегралы являются несобст­ венными. При дополнительном предположении о непрерывности указанных

частных производных в D упомянутые интегралы переходят в собственные.

§ 1

ФОРМУЛА ГРИНА

177

Мы докажем сначала формулу Грина для специального, но

достаточно широкого класса областей. Затем мы установим ряд

вспомогательных утверждений, которые понадобятся для дока­ зательства сформулированной теоремы.

2. Доказательство формулы Грина для специального

класса областей. Пусть D - односвязная конечная область с кусочно-гладкой границей L. Будем считать, что каждая пря­

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

По предположению, существуют несобственные интегралы

от частных производных функций Р(х, у) и Q(x, у). Это озна­ чает, что для любой системы областей {Dn }, монотонно исчер­

пывающих область D, справедливо, например, соотношение

n---+оо J!

дх

J! дх

lim

aQ dxdy =

aQ dxdy

D n

 

D

(аналогичные соотношения справедливы и для других частных производных функций Р(х, у) и Q(x, у)).

Опишем построение специальной системы областей {Dn},

монотонно исчерпывающих область типа К. Эта система пона­

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

формулы Грина для областей ука-

У

У2

занного типа.

 

 

 

Пусть сегмент [а, Ь] оси ОХ

У2 - Еn

 

представляет собой проекцию на

 

А

эту ось области D (рис. 7.1). Про­

 

Yl + Еn

-1-

ведем через точки а и Ь прямые,

 

Уl 1

параллельные оси Оу. Каждая из

 

 

1

этих двух прямых пересекается с

 

1

границей L лишь в одной точке.

О

а йn Х

Эти две точки А и В пересече-

ния указанных прямых с грани-

 

Рис. 7.1

цей L ,разделяют L на две кривые

 

 

L' и L ,которые, очевидно, представляют собой графики непре­ рывных и кусочно-дифференцируемых на сегменте [а, Ь] функ­

ций Уl(Х) и У2(Х) соответственно. Отметим (см. рис. 7.1), что Уl (х) :::;; У2 (х) (равенство имеет место лишь при х = а и х = Ь).

Рассмотрим далее последовательность сегментов [аn, Ьn] та­

ких, что а < аn < ЬN < Ь, аn --+ а, ЬN --+ Ь при n --+ 00. Пусть,

кроме того, при любом n сегмент [аn, Ьn] содержится в сегменте

n+l, Ьn+1]. Выберем число сп > О так, чтобы графики L~ и L~

функций Уl (х) + сп и У2 (х) - сп были расположены в области D

и не пересекались.

178 формулы ГРИНА, СТОКСА И ОСТРОГРАДСКОГО ГЛ.7

Границей области D n является кривая, составленная из ли­

ний L~ и L~ и отрезков вертикальных прямых, проходящих через точки аn и ЬN (см. рис. 7.1). Область Dn+l строится ана­ логичным образом, только вместо сегмента [аn, ьn] берется сег­

мент [аn+l, ьn+1], а число сn+l > О выбирается меньшим числа

сп. Очевидно, что если сп --+ О, то построенная система областей

{ D n} монотонно исчерпывает область D.

Докажем следующее утверждение.

Теоре,м,а 7.2. Пусть в области D типа К функции Р(х, у)

и Q(x, у) удовлетворяют условиям теоремы 7.1. Тогда для этои области и для функции Р(х, у) и Q(x, у) справедлива формула

Грина.

Д о к а з а т е л ь с т в о. Достаточно убедиться в справедливо­

сти равенств

f Q dy,

-!!~:dxdy =

f Pdx.

 

!! ~~dx dy =

(7.2)

 

 

 

 

D

L

D

L

 

Так как указанные равенства доказываются однотипно, мы про­ ведем доказательство второго из них. Рассмотрим двойной ин­

теграл

 

!!~:dxdy.

(7.3)

 

 

 

D n

 

Для области D n

и для подынтегральной функции дР в инте­

 

 

ду

грале (7.3) выполняются все условия, при которых действует

формула повторного интегрирования. По этой формуле имеем

 

Ьn

У2(Х)-Еn

!! ~:dx dy = ! dx

!

~:dy =

D n

аn

Yl(X)+E n

 

 

Ьn

 

Ьn

 

J Р(х, У2(Х) -

сп) dx - J Р(х, Уl(Х) + сп) dx (7.4)

Левая часть соотношений (7.4) при n --+ 00 имеет предел, равный

интегралу !! ~: dx dy. В силу равномерной непрерывности

D

функции Р(х, у) в замкнутой области D, каждое из слагаемых в правой части (7. 4) имеет при n --+ 00 предел, равный для пер-

ь

ь

вого слагаемого JР(х, У2(Х)) dx и для второго

JР(х, Уl(Х)) dx.

а

а

§ 1

ФОРМУЛА ГРИНА

179

Первый из этих двух интегралов представляет собой при ука­ занном на рис. 7.1 направлении обхода границы криволинейный

интеграл

- J Р(х, у) dx,

L"

а второй интеграл - криволинейный интеграл

J Р(х, у) dx.

L'

МЫ видим, что правая часть соотношений (7.4) при n --7 (х)

имеет предел, равный

- JР(х, у) dx.

L

Таким образом, вторая из формул (7.2) доказана. Справедли­ вость первой из формул (7.2) устанавливается аналогично (нуж­

но спроецировать D на ось Оу и повторить проведенные рассу­

ждения). Теорема доказана.

3. Инвариантная запись формулы Грина. Пусть функ­

ции Р(х, у) и Q(x, у) удовлетворяют условиям теоремы 7.1 в ко­

нечной связной области D с кусочно-гладкой границей L. Опре-

делим в области D = D + L векторное поле р, координаты ко­

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

равны Р(х, у) и Q(x, у). Очевидно, при условиях, наложенных на функции Р(х, у) и Q(x, у), поле р будет непрерывным в обла­

сти D и непрерывно-дифференцируемым в D. Найдем ротор

этого векторного поля. Используя выражение rot р в ортонор­

мированном базисе i, j, k, получим

rotp = (a Q _ aP)k.

дх ду

Из этого соотношения получим

aQ _

дР = krotp.

(7.5)

дх

ду

 

3 а м е ч а н и е 1. Перейдем в плоскости Оху к новому ор­

тонормированному базису i', j' и к новой декартовой системе

координат Ох'у', связанной с этим базисом. П1СТЬ векторное по­ ле р имеет в этом новом базисе координаты Р и Q'. Очевидно, в новой системе координат функции р' и Q' удовлетворяют усло­

виям теоремы 7.1. Кроме того, так как в новом базисе rotp =

= (a Q' _ aP')k то

дх'

ду'

'

 

 

 

 

 

aQ'

дР'

= krotp.

(7.6)

 

 

-

- -

дх' ду'

180 ФОРМУЛЫ ГРИНА, СТОКСА И ОСТРОГРАДСКОГО ГЛ.7

Так как скалярное произведение k rot р представляет собой ин­

вариант, то из (7.5) и (7.6) следует, что выражение aQ - дР не

дх ду

меняет ни значения, ни формы при переходе к новому ортонор­ мированному базису, т. е. также представляет собой инвариант.

С помощью этого замечания мы можем сделать следующий важный вывод: интеграл, находящиuся 6 ле60U 'Части формулы

Грина (7.1), имеет ИН6ариантныu хара'Ктер - его зна'Чение и

форма не меняются при переходе 'к НО60й. де'КартО60й. системе 'Координат. Действительно, при таком преобразовании коорди­

нат абсолютная величина якобиана преобразования равна еди­

нице. Согласно же замечанию, подынтегральное выражение не меняет ни значения, ни формы. Обратимся теперь к интегралу

§ Pdx + Qdy,

(7.7)

L

 

находящемуся в правой части формулы Грина. Убедимся, что

этот интеграл maKJlCe имеет ИН6ариантныu хара'Ктер - его

зна'Чение и форма не меняются при переходе 'к НО60й. де'Карто60й. системе 'Координат.

Пусть t - единичный вектор касательной в точках грани­ цЫ L, направление которого согласовано с направлением обхода на L, cos о: и sin о: - координаты вектора t. Выберем в качестве параметра на L длину дуги [, причем на каждой связной ком­ поненте границы возрастание параметра l согласовано с направ­ лением обхода на этой компоненте. При условиях, наложенных

на L, функция t(l) будет кусочно-непрерывноЙ. При сформули­

рованных выше условиях векторное поле р будет непрерывным

на L, а его координаты Р и Q представляют собой непрерывные

функции от [.

Заметим, что после выбора направления обхода и параметра

на кривой L криволинейный интеграл второго рода (7.7) преоб­

разуется в криволинейный интеграл первого рода. При этом Р и

Q вычисляются в точках L, а dx = cos о: dl, dy = sin о: dl. Таким образом,

§ Р dx + Q dy = §(Р cos о: + Q sino:) dl =

§ pt dl.

(7.8)

L

L

L

 

Соотношение (7.8)

показывает, что интеграл (7.7) действитель­

но имеет инвариантный характер: скалярное произведение pt -

инвариант, параметризация с помощью длины дуги не связана

с системой координат. Кроме того, в новой декартовой системе

координат Ох'у' имеем

pt dl = (Р' cos 0:' + Q' sin 0:') dl = р' dx' + Q' dy',

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