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

книги из ГПНТБ / Сахарников Н.А. Высшая математика учебник

.pdf
Скачиваний:
98
Добавлен:
27.10.2023
Размер:
24.22 Mб
Скачать

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

При этом условии имеют место формулы

 

V(p==grad(p,

ya = diva,

у Ха —rota.

 

(9)

Действительно,

 

 

 

 

 

 

 

 

 

 

 

_3_

 

Зф

Зф

Зф

grad ф,

 

 

 

 

д у

 

 

д х

д у

d z

 

 

 

 

 

 

 

 

Va

.

3

 

. 3

ах - з ау-'- k аг) ~~ div а.

 

1

----- 3

J

-т—

 

 

 

д х

'

д у

d z

 

 

 

 

По правилу векторного умножения векторов имеем

 

 

i

j

к

 

у

Ха =

3

3

3

rot а.

д х

ду

d z

 

 

 

 

 

ах

ау

az

 

В результате выполнения дифференциальных операций вто­ рого порядка получим непосредственно из (9) следующие фор­ мулы :

уу Ф= у g

r a dф =

d i v g фr

a= dДф,

V X УФ = у Xg r a фd = r o t g r aф,d

y y a

= y d

i v a

= g r a d d i v a ( ,1 0 )

Ѵ УX а = у r o tа = d

i v r а,o t

У Xу Xа = у Xr o аt

= r o t r oа.t

203. Выражение градиента в ортогональных криволинейных координатах. Рассмотрим две системы координат — деікартову прямоугольную систему (х , у, z) с ортами i, j, к и какую-либо криволинейную ортогональную систему (u, v, w) и пусть е„, еѵ,

ew — соответствующие орты этой системы (см. п. 178). Формулы связи между координатами предполагаются известными:

х — х(и,

v,

w),

у = у(и,

v,

w),

z = z(u,

v,

w),

и = и(х,

у,

z),

ѵ = ѵ(х,

у,

z),

w = w(x,

у,

z).

Следовательно, имеет место соотношение

“Ь МуУи ф UzZu — 1,

(11)

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

и(х(и, v, w), у (и, v, w), z(u, v, w)) и.

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

По правилу дифференцирования сложной функции имеем

f x ~ і ѵ М х “Г I r P x ~f“ f w W X l

f y - ~ f u U y~ \ - f v v y ~ c‘ i w W y ,

f z ~ f u U z -t~ f x P z "T~ î w W z -

Умножим эти равенства соответственно на i, j, к и, сложив,

получим согласно

равенству

(7)

 

 

grad / (М) =- /ûgrad и t

/;grad v -f/^grad w.

(12)

Известно (см.

п. 184),

что

градиент скалярной

функции

и (X, у, z) есть вектор, направленный перпендикулярно поверх­

ности уровня и = и0. Поэтому он коллинеарен вектору еи.

Сле­

довательно, имеем

 

grad и = и'х i-\-Uyi-r u'zk -- h еш

(13)

где h — скалярный множитель, подлежащий определению. Производная векторной функции г = xi -j- yj + zk по ска­

лярному аргументу, например и, согласно п. 185 есть вектор, направленный по касательной к годографу вектора г (и), т. е. по касательной к координатной линии и. Поэтому

 

 

ru = x'u\ + y u \ + z uk = Hueu,

 

 

(14)

где Ни — соответствующий

коэффициент

Ламе

(см.

п.

178) си­

стемы (и, v, w).

 

 

 

 

 

 

 

 

 

 

В результате скалярного умножения равенств (13) и (14) полу­

чим равенство

-ЕуЦх~!г УѵЦу~т

 

= HJi,

 

 

(15)

 

 

 

 

 

левая часть

которого

равна

единице по

формуле (11). Поэтому

h = 1/Ни и

согласно

(13)

имеем

grad

 

j

 

 

и — -гг еи. Аналогично

 

 

1

 

 

 

 

1

п и

 

 

получим grad

ѵ =

 

 

 

 

 

 

 

 

-гг£ѵ, grad w =

-==- ew.

 

 

 

 

'

'

-L* V

 

 

* * W

 

 

 

 

Теперь из (12) следует выражение градиента скалярной функ­

ции точки в любых

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

grad /(М ) =

 

1 -ËL е -

1 df

1 - Ê L

е

(16)

На ди

и Ііѵ

дѵ

н„

dw

е“)-

204. Выражение дивергенции в ортогональных криволинейных координатах. Это выражение установлено ниже с помощью инва­ риантного определения дивергенции (см. и. 197). Для получения дивергенции вектора а (М) в точке М 0 (и0, v0) w0) рассмотрим криволинейный параллелепипед, ограниченный координатными поверхностями (рис. 147) и = и0, и = и0 du, v = ѵ0, v = v0 +

+ dv, w — w0, w WQ -f dw. Элементы длин дуг координатных линий и элемент объема найдены в п. 178. Поток вектора через каждую грань параллелепипеда представим в виде произведения площади соответствующей грани на среднее значение величины ап. Поэтому поток вектора а через противоположные грани M 0M 3N sN n и М ХМ 2N 2N х (площади которых равны dsudsw) приблизительно равен

JJ

anda æ (a nds0dsw)

\ — (ands0dsw) | ==

S1+S2

Kohdu

 

щ

(auHvHw) J

(auHvHw)

dv dw æ

d

(auHvHw) du dv dw.

Uoi-du

 

du

 

 

 

 

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

Аналогичные формулы имеют ме­ сто для потоков вектора через остав­ шиеся пары граней. Поток вектора через всю поверхность параллеле­ пипеда получится равным

^ a nda « [^ -(а „ Я „ Я J + ~ { a aHuHw) + ^

(ашЯ иЯ„)] dudvdw.

s

 

 

 

 

Теперь, положив в формуле (23)

п. 197

Ѵв = HuH vHwdu dv dw,

в результате сокращения дроби на du dv dw

после

предельного

перехода при стремлении к нулю

величин

du, dv и

dw получим

выражение дивергенции в ортогональных криволинейных коорди­ натах

^iv а = H UHVHW\_~~ди (a“HvHw) + ЛдГ (avHuHw) + -^Г (awHuHv) j . (17)

205.Выражение оператора Лапласа в ортогональных криволи­

нейных

координатах. Оператор Лапласа сікалярной функции

/ (и, V,

w) согласно формуле (6) есть дивергенция от ее градиента.

Поэтому, применяя формулу (17) к вектору а = grad / в форме (16), т. е. полагая в (17)

___1

df

a,

J __ dj__

1

df

a u ~~ Л и

du

H u дѵ

а“> Hw

dw

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

А/

I

д

(

IIуПg, _9/_\

I

(

HUHW d f \ ,

UuH-Jlw

ди

\

Пи

ди

/

' дѵ

\

IIѵ д ѵ J '

 

 

 

.

д

!

1IUIIV

df

 

 

(18)

 

 

~

d w V IIw

d w J J '

 

 

 

 

 

Выражение оператора Лапласа в цилиндрических координатах выведем из (18), положив согласно п. 178 Нр — 1, 7/ф = р, Нг — 1. Таким образом, получим

А/ =

д

1

02/ ,

02/

(19)

~Эр

р 2

9 ф 2

dz*

Выражение оператора Лапласа в сферических координатах получается из (18) при Hr ~ 1, 7/ф = г sin 0, Н% = г (см. п. 178). В результате получим

А/ = 7-2

' д

Л І

1 .

д Л

1

д

df

( 20)

~д7

дг

S ill2 0 "Г-

ÖCp2

S i n ö

sin О 90~

Иногда

в сферической

системе

координат

угол Ѳ отсчитывается

не от оси Oz,

а от плоскости Оху. Тогда новый угол Ѳ' =

— Ѳ

и в (20) вместо sin Ѳвсюду войдет cos Ѳ'.

Многие задачи естествознания приводят к уравнению А/ = 0,

называемому уравнением Лапласа. Оно

имеет решение / = —

в области, где г Ф 0. Действительно,

f r’ — ----—2 и А

что прямо следует из (20). Функция / = — называется фунда­

ментальным решением уравнения Лапласа.

206. Уравнение диффузии. Пусть в данной неподвижной среде неравномерно распределено некоторое вещество а. Рассмотрим процесс изменения концентрации этого вещества и (717, і) вслед­ ствие диффузии.

Выделим и закрепим в рассматриваемой среде область В, ограниченную поверхностью S. Подсчитаем двумя способами

количество

Q вещества а, переместившееся вследствие диффузии

в область В

извне за промежуток времени от tx до t2.

Согласно экспериментальному закону Фике, через элемент

поверхности da в сторону нормали

п за время dt переместится

вследствие диффузии количество вещества, равное

 

dQ — D

da dt,

( 21)

где D — положительная величина, называемая коэффициентом диффузии. Среда предполагается изотропной в том смысле, что D не зависит от направления п.

Из (21) следует, что

И ’Д ИD — do.

ОП

ди

Согласно формуле (8) п. 184 - grad u n; поэтому в соот­

ди

ветствии с формулой Остроградского имеем

1

2

 

Q J

dt П j div (D grad и) dr.

(22)

Это перемещение вещества повлечет за собой изменение кон­ центрации внутри В и может быть подсчитано иначе. Увеличение

концентрации и на du = dt за промежуток времени dt потре­

бует сообщения элементу dr количества вещества, равного

cdudx

■с 4 ^- dt dx,

 

 

 

dt

 

где с — положительный

коэффициент пористости среды. Вся

область В за время от tx до t2поглотит количество вещества

І *

Ш ( С -5Г

(23)

 

где / — интенсивность источников вещества а

в области В. При­

равнивая выражения (22) и (23), придем к равенству

^2

 

 

 

1 ^ 11 і [_С

^— diV (Z) grad и)J dx О,

которое выполняется для

л ю б о й области

В. Отсюда следует

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

с

ди

(24)

— = div (D grad и) 4- / ,

называемое уравнением диффузии в неоднородной среде (или урав­ нением теплопроводности, так как этому уравнению удовлетво­ ряет температура среды, • изменяющаяся вследствие теплопро­ водности).

Если среда однородна (с и D постоянны), то при отсутствіи! источников (/ = 0) концентрация и удовлетворяет уравнению теплопроводности

ди

- =

а2Ди.

(25)

dt

1

П 5

как сііѵ grad и = Аи. В случае установив­

где а = у

— , так

шегося процесса

—- 0 и из (25) следует,

что концентрация и

удовлетворяет уравнению Лапласа

 

 

 

Au::-0.

(26)

Лемма. Если / (М) непрерывна в области С и 1))] / (М) ско —

в

= 0, где область В есть произвольная частъ области С, то / (М) :

II « С.

Д о к а з а т е л ь с т в о .

Рассуждая от

противного,

пред­

положим,

что

в точке М 0 области С функция / (М) принимает

неравное

нулю

значение.

Пусть для определенности / (М0) > 0.

Неравенство / (М) О 0

будет

выполняться

в некоторой

окре­

стности В х точки М 0, так как / (М) непрерывна.

Поэтому имеем

И

і і (М)/ dcû > 0,

что противоречит

условию

леммы. Отсюда

“ в,

 

 

 

 

следует ее заключение.

 

 

 

к

Д о б а в л е н и е.

Следуя

логической

схеме, которая привела нас

уравнению диффузии, выведем

так называемое уравнение неразрывности.

 

Рассмотрим для этого некоторую массу жидкости (или газа) и обозначим

через V (М, t) и р (М , t) соответственно скорость течения жидкости и ее плот­ ность в точке М в момент времени t. Выделим мысленно в потоке жидкости произвольную область В , ограниченную поверхностью 5. Пусть Q — коли­

чество жидкости,

вытекающее из В за промежуток времени от до t 2. Через

элемент площади

do за время dt переместится в сторону внешней нормали п

количество жидкости, равное dQ = pxndodt. Поэтому

t,

12

 

Q = J dt J j

рѵл da = J dt J J [ div (pv) dx.

U S

t ,

B~

Это перемещение жидкости влечет за собой изменение плотности жидкости внутри В. Рассуждая так же, как п выше, получим

і, В

Знак минус поставлен потому, что если плотность уменьшается внутри В,

то

< 0, а Ç > Ü. Если же жидкость втекает в область В, то ~ > 0,

аQ < U.

Всилу закона сохранения вещества оба выражения Q равны между

собой. Отсюда с помощью леммы следует уравнение неразрывности

l - f d iv ( p v ) = 0,

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

Глава X II

РЯДЫ

§ 35. ЧИСЛОВЫЕ РЯДЫ

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

и

Лейбница. Существенный вклад в теорию рядов внесен Эйлером

и

Гауссом, позднее Коши и Больцано.

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

(см.

и. 2).

Задать

числовую

последовательность — это значит

задать

ее

общий

член как

функцию натурального аргумента

ап =

/

(«)•

 

 

 

Составленное из элементов данной последовательности выра­ жение

СО

 

аг + й2+ а3-г . . . -і- ап-f . . ., или £ ап,

(1)

п = 1

 

принято называть числовым рядом, или просто рядом. Отдельные элементы ап, из которых образовано выражение (1), называют

членами данного ряда.

Между членами ряда в выражении (1) стоит знак плюс; это указывает как будто на то, что все члены ряда надо сложить. Однако их бесконечно много, а действие сложения чисел опре­ делено лишь для конечного числа слагаемых. Возникают во­ просы — что же такое сумма ряда, всякий ли ряд имеет сумму? Ответы на эти вопросы даны ниже.

Сумма п первых членов ряда называется п-й частичной суммой ряда', обозначим ее sn. Имеем

s n — а і + а 2 + • • • а п - (2)

Составим последовательность частичных сумм ряда (1)

О п р е д е л е н и е . Если последовательность частичных сумм ряда (1) имеет конечный предел

 

 

 

lim sn = s,

 

(4)

 

 

 

п-*-СО

 

 

то ряд (1)

называется

с х о д я щ и м с я , а число s называется

с у м м о й

р я д а

(1). В этом случае можно написать

 

 

 

 

 

СО

 

 

s = ах+

Л-2+

• • • + йп т~ • • ., или

s = ^ ап.

(5)

 

 

 

 

п = 1

 

Ряд (1)

называется

р а с х о д я щ и м с я ,

если последователь­

ность его частичных сумм не стремится к конечному пределу (например, sn -> при п со).

Перефразируя определение, можно сказать, что бесконечный ряд называется сходящимся, если сумма его первых п членов при неограниченном возрастании числа п стремится к конечному пределу: и этот предел называется суммой ряда.*

О сумме бесконечного ряда можно говорить только тогда,

когда он сходится, и тогда частичная сумма sn является

при­

ближенным выражением для суммы ряда.

 

П р и м е р 1. Рассмотрим так называемый геометрический ряд

 

а+ + а?2+ • • .+ ад”+ . •

(6)

члены которого образуют геометрическую прогрессию со знаменателем q,

причем а ф 0. При q ф 1

сумма его п первых членов равна

 

 

S n

 

а aqn

a

 

aqn

 

 

 

------------------------ = 2 — -------------

—------------ .

 

 

 

 

1— q

1—g

1q

1)

Если \q\ <; 1,

ТО

qn

0 при

n oo

и,

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

 

lim

sn =

lim

aqn

\

a

 

 

 

 

b 11? '

 

n-yœ

 

 

 

 

 

В случае | g| < 1 ряд

(6) сходится и его сумма s =

------ .

2)

Если ] q I > 1, то

qn -► оэ при

оо и тогда переменная sn предела

не имеет: sn —» оо при п -*■ оо. Поэтому в случае | g | >> 1 ряд (6) расходится. 3) Если g = 1, то ряд (6) имеет вид а + а + ... + а + ••• В этом случае

sn = ап -> оо при п-у оо, т. е. ряд расходится.

* При рассмотрении пределов последовательностей переменная п всегда стремится к бесконечности. Поэтому символ п -*■ оо мы иногда будем опускать в записи предела последовательности. Например, равенство (4) можно запи­ сать так: lim s„ = s.

Точно так же ряд с общим членом ап можно представить символом ^ ап, если известно с какого номера начинается суммирование.

4)

Если q — —1, то ряд (6) имеет вид а а

а — я -f-

. . . В

этом слу­

чае sn = 0,

если

п — четное, и

sn — а, если п — нечетное. Следовательно,

sn предела не имеет, а ряд расходится.

 

1:

 

Итак, геометрический ряд (6) сходится только в случае | q\

 

 

 

 

 

оо

 

 

 

 

 

 

 

 

2

ч ” - 7 ^ '

 

 

«

 

 

 

 

п=0

 

 

 

П р и м е р

2.

Дап ряд

Его общий член можно представить

 

 

 

 

1

 

 

 

 

в виде ап =

 

 

 

1

 

 

 

 

п(/г+1)

п

п-\-1 .

Поэтому

■••+»" = ( Ч

) +

 

 

 

 

1

1

1 n-f-1

сумма ряда

5 =

+

+ - +(^ ~ 4 д )

= lim ^ 1 —

 

= 1.

Этот пример иллюстрирует

возможность

точного

суммирования ряда в отдельных случаях.

208. Основные свойства рядов. 1°. Если все члены сходящегося ряда (1), сумма которого равна s, умножитъ на одно и то же число Ь, то получим сходящийся ряд, и его сумма равна bs:

 

2 ban = b ^

а„.

(8)

 

П=1

 

 

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

 

<т„ Ьаг

4- Ъап Ъ(%

. . . -j- ап) = bsn

(9)

имеет конечный предел, равный bs:

lim ап--lim bsn b lim snbs.

( 10)

2°. Сходящиеся ряды можно почленноскладывать и вычитать

с о СО ОО

2

вя ± 2

ьп~ Ъ (ап ± ь п).

(il)

71=1

71= 1

71= 1

 

Для доказательства обозначим через Ап и В„ частичные суммы

соответственно рядов ап и У, Ъп, и через А и В — их суммы. Частичная сумма ряда из правой части равенства (11) ап = Ап + + В„ имеет предел, равный

lim сг„ = lim А „ lim В„ = А ~г В.

(12)

Для .выяснения других свойств рядов введем понятие под­ последовательности. Пусть {хп} — некоторая последовательность. Рассмотрим произвольную возрастающую последовательность на­ туральных чисел ftlt к 2, . . ., кп, . . . Выберем из последователь­

ности {хп} элементы с номерами к и

к2, . . ., кп, . . . и расположим

их в таком же порядке, как и числа кп:

З'йп

■• • ,

ап) ■• • •

Полученную последовательность

назовем подпоследователь­

ностью последовательности {хп}.

 

Лемма. Если числовая последовательность {.х„}

сходится к не­

которому числу Ь, т. е. lim хп =

Ь, то и любая

подпоследователь­

ность этой последовательности сходится к тому же

числу Ъ.

Действительно, так

как

lim хп -= Ъ,

то для

любого е > 0

существует число N такое, что при п ^

N выполняется неравен­

ство

\хп 6 I < е. Пусть

{х^} — некоторая

подпоследователь­

ность

последовательности {хп}.

Так как kN ^

N, то,

начиная

с номера kN, элементы

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

{ ^ я} удовлетво­

ряют неравенству ] xkn Ъ| <; е. Поэтому подпоследовательность {xk } сходится и имеет пределом число Ъ.

3°. Сходящиеся ряды обладают сочетательным свойством, т. е. можно объединять в группы любые стоящие рядом члены сходя­ щегося ряда, а затем производить суммирование по группам, причем эта операция не нарушит сходимости ряда и не изменит его суммы.

Действительно, пусть дан сходящийся ряд (1) и его сумма равна s. Построим ряд а, который получается из (1) путем объеди­ нения в группы некоторых рядом стоящих членов. Выполнить это объединение можно различно, например, так:

К + a-ù + («3 4- «4 + аъ) + «в - Г («7 + я8) -г . . •

(13)

Последовательность частичных сумм ряда а образует подпоследо­ вательность {skn} последовательности частичных сумм ряда (1).

В частности, для ряда (13) это будет s2, s5, se, s8, . . . Согласно лемме существует и равен s предел lim s* . Поэтому ряд а схо­

дится и его сумма равна сумме ряда (1).

Заметим, что правила действий с рядами не всегда совпадают с правилами действий с конечными суммами. Например, если дан сходящийся ряд вида (13), то, опустив скобки, мы получим новый ряд, который может оказаться и расходящимся. Так, ряд (1 1) + (1 — 1) + (1 1) + ..., очевидно, сходится, между тем как полученный из него опусканием скобок ряд 1 1 4 - -f-l — 1 + ... будет расходящимся.

Если в данном ряде (1) отбросить первые к членов, то получим новый ряд, который называется остатком ряда (1) после к-ѵо члена:

СО

 

Ъ ап =■■ak+1+ ak+2+ . . .

(14)

71=Ml

 

4°. Ряд и любой его остаток либо оба сходятся, либо оба рас

ходятся.

Соседние файлы в папке книги из ГПНТБ