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

книги из ГПНТБ / Дьяченко, А. Н. Интегральное исчисление функций нескольких переменных. Теория вероятностей и элементы математической статистики учебное пособие

.pdf
Скачиваний:
26
Добавлен:
19.10.2023
Размер:
7.54 Mб
Скачать

функция, определенная в области У). Аналогично тому, как это было сделано для плоских областей (см. § 1.5), вводятся понятия статических моментов и моментов инерции области V относительно координатных осей и координатных плоскостей.

Введем следующие обозначения: Sx, Sy, Sz — статические мо­ менты тела V относительно координатных осей х, у, г; Sxy, Syz,

Sxz — статические моменты тела

V

относительно

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

плоскостей Оху, Oyz, Oxz\ Jx, Jy,

Jz — моменты инерции тела V

относительно координатных осей х,

у,

z; Jху, Jyz,

,Jхг —■моменты

инерции тела V относительно координатных плоскостей Оху, Oyz, Oxz.

Возьмем в области V произвольную точку с координатами (х, у, г). Расстояния этой точки до координатных осей и координат­

ных плоскостей

соответственно равны

(рис. 27):

 

 

 

•)

4 = V V

+ z2 ; dy=]/~x2 + z*; dz = х2-f у2 ;

 

d-xy —

dyz= x,

dxz = y.

Применяя рассуждения, аналогичные тем, которые привели нас к формулам (1.57), (1.58), получим следующие формулы:

= J JJ/ (х,

у,

г) ]/ y2 + z2dv,

 

 

 

v

 

 

 

 

 

 

 

 

 

Sy= f j f f ( x ,

у,

z) ]/rx2 + z2dv,

 

 

 

v

 

 

 

 

 

 

 

 

 

Sz = ^ j f ( x ,

у,

z) У x2, -\-y2dv,

 

 

 

V

 

 

 

 

 

yf (*■ У

dv;

s xy= Ш zf (*• У’

г) dv’

V

 

 

 

 

 

 

V

 

 

V

 

У’

z^dv’

 

 

 

(2.23)

 

 

 

 

 

 

 

 

Jx= y y y 2+ z 2)f(x,

y,

z) dv,

 

 

 

 

 

V

 

 

 

 

 

 

 

 

 

Jy = ГJ J (x2+ z2)/ (x,

y,

z) dv,

 

 

 

Jz = § y ( x 2 + y2)f(x,

y,

z) dv,

 

 

 

V

 

 

 

 

 

 

 

 

 

Jxy = J j J z2f (X,

y,

 

z) dv,

, Jxz=

y \ y2f (x, y,

Z) dv,

V

 

 

 

 

 

 

V .

 

 

Jyz= J ffx 2f(x,

y,

z)dv.

 

 

 

 

V

 

 

 

 

 

 

 

 

 

Пример 5. Вычислить статические моменты

относительно плоскости

Оху цилиндра х%+

у2 = R2; 0 < 2 < а

(рис.

28), если плотность рас­

50

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

сти Оху:

]5 (х, у, z) — kz2 (k — const).

Применяя формулы (2.23), получим:

= J [ j г/(х, у, z) dv =

k J ]" J zMxdydz.

V

V

 

Перейдем в этом интеграле к цилиндрическим координатам:

x = /'co s 9 ; y = r sin ср;

z = z ;

|/ |= г.

Применив формулу (2.20), имеем:

 

Sxy = k J J [ rz3d<pdrdz.

(2.24)

J v:

 

 

Легко видеть, что интервалы изменения цилиндрических коорди­ нат ср, г, z не зависят друг от друга:

0 < ф < 2л; 0 < г < R; 0 < z < а.

Расставляя пределы в тройном интеграле (2.24), имеем:

г

,

2? j

?

j

г oj

, п

Я2 а4

=

ЫЯга*

.

Sxy =

k

йф

о-

rdr

z3dz =

к2л

----------2 4

------------4

 

 

о 1

 

о

 

 

 

Вычисление координат центра масс тел. Определение центра масс, сформулированное ранее (см. § 1.5), обобщим на случай трех­ мерной области V. Центром масс трехмерной области будем назы­ вать точку, обладающую следующим свойством: если в эту точку поместить материальную точку с массой, равной массе всей области, то статический момент этой материальной точки относительно лю­ бой координатной плоскости равен статическому моменту самой области относительно этой же координатной плоскости. Пусть (хц, г/ц, 2Ц) — координаты центра масс области V. Тогда согласно определению имеем:

>yz

у*=- м

>ху

м

м

где М — масса области

V.

 

3*

51

Используя формулы (2.23); (2.3), получим:

J j J xf (х, у, z) dv

 

 

j \J у! (X, у, г) dv

_v_____________

Уц

v

 

J j j / (*. У, г) dv

1 J J j f (X,

у, г) dv

j j j z / ( x , у, z) dv

 

= —------------------ .

 

j j

j f (x,

У,

z) dv

 

V

 

 

 

 

Пример 6. Вычислить координаты

центра масс

цилиндра, заданного

в примере 3. Вычислим сначала массу этого цилиндра:

М __ J J | kz2 dxdydz. v

Перейдем, как и в примере 5, к цилиндрическим координатам:

R

a

D2

knR2a3.

М — k j'

dcp j

rdr J z2dz=k2n

------------2

=

o

o

o

3

3

Из соображений симметрии относительно оси z (и цилиндр, и плот­ ность распределения массы симметричны относительно оси z) очевидно, что центр масс находится на оси z, т. е. хц = уц = 0. Вычислим гц, воспользовавшись результатом примера 5:

с

knR2al

 

^ х у

:

 

Dxy _ / knRZa*

) :(

knR2a3

а.

М ~ { 4

 

 

Таким образом, центр масс рассматриваемого цилиндра имеет

координаты (о, 0, -5 - а |

Г л а в а 3

КРИВОЛИНЕЙНЫЕ ИНТЕГРАЛЫ

3.1. Криволинейные интегралы первого рода

Задача, приводящая к понятию криволинейного интеграла пер­ вого рода. Пусть требуется найти массу тяжелой нити. Под нитью будем понимать тело, двумя размерами которого можно пренебречь по сравнению с третьим размером (длиной). Под словами «тяжелая нить» будем понимать кривую линию с распределенной вдоль нее массой. Рассмотрим плоскую кривую АВ. Пусть в каждой точке кривой определена непрерывная положительная функция / (х, у), равная плотности распределения массы вдоль кривой. Разделим кривую АВ на п частей (частичных кривых) следующими друг за другом точками M lt М 2, ■■.,Мп- 1 (рис. 29). В каждой k-ой части кривой Mk-\Mk (к — 1, 2, . . . , п) (будем считать, что точка М0

52

совпадает с точкой А, а точка Мп совпадает с точкой В) выберем

по произвольной точке Pk (xk, yk). Будем считать, что вдоль частич­ ной кривой Mk-\Mk масса распределена равномерно с плотностью,

равной значению функции f (х, у) в выбранной точке, т. е. f (xk, yk). Обозначим длину k-ой части кривой через Д/ь ее масса приближенно

равна величине / (xk, yk) Alk. Суммируя массы всех частичных кривых, получим приближенное значение массы всей кривой:

М = 2 •/(**, yk) Alk. k=i

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

%n = max{Alk} (k=\, 2, . . . , п).

Очевидно, что при Кп 0 стремятся к нулю одновременно длины всех частичных кривых. За массу кривой следует принять величину, которая не зависит ни от способа

дробления кривой, ни от выбора

 

 

промежуточных точек, а опреде­

 

 

ляется только плотностью распре­

 

 

деления и формой кривой. Заме­

 

 

тим, что влияние произвола обус­

 

 

ловленного

выбором внутренних

 

 

точек и способом дробления кри­

 

 

вой, уменьшается при уменьшении

 

 

ранга дробления кривой,

так как

 

 

в достаточно малой части кривой

Р и с.

 

значения

непрерывной

функции

 

мало друг от друга отличаются.

 

М, опре­

Поэтому за массу кривой естественно принять величину

деляемую

равенством:

 

 

 

 

 

М=

Пт ^

f(xk, yk) A/fc.

 

(3.1)

Определение криволинейного

интеграла первого рода.

Вычисле­

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

Пусть в плоскости Оху задана кривая L в каждой точке которой определена функция f (х, у). Разобьем кривую L на п частичных кривых, занумерованных последовательно от начальной точки А до конечной точки В: L1; L 2, . . . , Ln, не имеющих общих внутрен­ них точек. Обозначим длины этих частичных кривых соответственно АА12, ■■■, А/„ (заметим, что А 4 > 0 для всех k). Наибольшую из этих длин назовем рангом дробления кривой и обозначим через Кп. В каждой частичной кривой Lk (k — 1, 2, . . . , п)-выберем по

53

произвольной точке Pk (xk, yk) и составим интегральную сумму:

Ук)А1к-

k=i

Если при стремлении к нулю ранга дробления кривой сущест-

П

вует конечный предел суммы 2 f(xk> Ук)

не зависящий ни от

k=i

способа дробления кривой, ни от выбора точек Pk (xk, yk), то этот предел называется криволинейным интегралом первого рода от функции / (х, у) по кривой L от точки Л до В и обозначается

J f (X, у) dl или | f (х, у) dl.

L

АВ

 

Таким образом, по определению

[f{x, у) dl =

lim 2 f(xk> Ук)АК-

(3-2)

L

k=l

 

Используя это определение, формулу (3.1) можно записать в

виде

(3.3)

M = ff(x , у) dl,

L

т. е. если масса распределена вдоль кривой L, то она может быть представлена криволинейным интегралом первого рода от плотно­ сти, вычисленным вдоль рассматриваемой кривой.

Сформулируем без - доказательства теорему, указывающую на класс функций, для которых существует криволинейный интеграл первого рода.

Теорема. Если функция f (х, у) непрерывна вдоль кривой L,

а множество значений функции, вычисленных в точках кривой ог­ раничено, то существует криволинейный интеграл первого рода от функции f (х, у) вдоль кривой L.

Свойства криволинейного интеграла первого рода. Рассмотрим простейшие свойства криволинейных интегралов первого рода, считая, что все указанные ниже интегралы существуют.

Свойство 1.

Ja/(x,

y)dl = a^f(x,

y)dl,

(3.4)

l

'l

 

 

где a — любая постоянная.

 

 

 

Свойство 2.

 

 

 

f t/i (*. У) + /* (*. У)\ dl = J7i (x, y)dl + \ ft (x, y) dl.

(3.5)

L

L

L

 

Свойства 1 и 2 непосредственно следуют

из определения

криво­

линейного интеграла I рода.

 

и L2, не

Свойство 3. Если кривую L разбить на две части L,

имеющие общих внутренних точек, то

 

J f(x, у) dl =

f f(x,

у) dl + J7 (x, у) dl.

(3.6)

L

L,

L 2

 

54

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

ком или кривой или кривой L 2. Это позволяет разбить сумму

П

2 f ( x k> Ук) Alk на две суммы: в одну сумму включить слагаемые, fc=i

соответствующие

кривой L b

а в другую — слагаемые, соответст­

вующие кривой Ь 2, т. е.

 

2 f{xk-

Ук) Alk = ^ f (xk, yk) A 4+2/(*/S:> Ук)

fc= l

L ,

L „

Нижние индексы, стоящие под суммами в правой части равен­ ства, указывают, к каким кривым относятся слагаемые данной суммы. Переходя в этом равенстве к пределу при кп -* 0, получим доказываемое равенство (3.6).

Свойство 4. Величина криволинейного интеграла первого рода не изменится, если изменить направление интегрирования. Если А и В — конечные точки кривой АВ, то

J /(*• y)dl= | f(x, y)dl.

АВ

BA

Д о к а з а т е л ь с т в о . Рассмотрим формулу (3.2). Все сла­ гаемые суммы, стоящей в правой части равенства, имеют вид

/ {хьУк) Aik- При изменении направления интегрирования не ме­ няется ни значение функции f (xk, yk), ни длина k-ой частичной кривой Alk. Отсюда следует, что величина интеграла (3.2) также не изменится при изменении направления интегрирования.

Сведение криволинейного интеграла первого рода к определен­ ному интегралу. Пусть кривая L, вдоль которой вычисляется кри­ волинейный интеграл первого рода от функции / (х, у), задана па­ раметрическими уравнениями:

x = x(t); у = у (t); ( а < Д < Р ) .

Будем считать, что функции х (t) и у (t) непрерывны в [а, р] и имеют в этом промежутке непрерывные производные х' (/) и у' (7). Разобьем кривую L на п частей, этому разбиению соответствует раз­ биение промежутка [а, р] на я частей:

a = n 0 < ^ < ^ 2 < . . . < t n= р.

Вычислим длину k-оц частичной кривой, воспользовавшись фор­ мулой для вычисления дуги плоской кривой:*

*k t_______________________

А/ * = J У W (t)]2+ [y' (t)]2dt (k= 1, 2, . . . , я).

1k—\

*См. [2], формула (7.12).

55

Применяя

теорему

о среднем для определенного

интеграла,*

получим:

 

 

 

 

А/, =

V w (Ik)?+ W (h)? А4 .

(3.7)

гДе t k - i ^

^ V ’ ^ k ~ ^ k h - v

 

Составим интегральную сумму для функции f (х , у) вдоль кри­

вой L, при этом в качестве произвольных точек Pk (xk, yk) выберем точки, координаты которых определяются равенствами:

xk = x (lky, ук = у(Ък),

(3.8)

получим:

ljf(Xh> Ук)Ык.

к= 1

Подставив в эту интегральную сумму формулы (3.7) и (3.8), получим:

f(xk, ук) А1к =

2 f[x (tk), y {lk)\V\x'{lk)? + W (lk)? Atk. (3.9)

k=\

k=i

 

 

В правой части этого равенства стоит интегральная сумма от

функции f (х (t), у {t))V lx' ( 0 12 + W (О)2 по промежутку (а,

0].

Перейдем в равенстве (3.9) к пределу при

-* 0, заметив,

что

при этом стремится к нулю и величина Кп, определяемая равенст­ вом:

А,я = шах{А^) (k— l, 2, . . . . п).

Получим:

Р __________________

J f [х, у) dl = J f(x(t), у (/)) V lx' W]2+ W (t)?dt.

(3.10)

L

a

 

 

Это равенство позволяет вычислять криволинейный интеграл I

рода с помощью определенного интеграла.

(3

а,

Заметим,

что равенство (3.10) имеет смысл только при

так как именно при этом условии справедлива формула для вычис­

ления

длины

кривой, которой мы воспользовались. В случае, если

р < а ,

то для

того, чтобы формула (3.10) осталась верна, перед

интегралом, стоящим в правой части равенства, следует поставить знак «—» (это следует из того, что определенный интеграл меняет знак при изменении направления интегрирования). В некоторых случаях бывает полезно, чтобы формула, связывающая криволи­

нейный

интеграл I рода с определенным интегралом, была верна

* См.

[2!, формула (6.10).

5R

как для случая а > р , так и для случая р > а . Тогда формула (ЗЛО) может быть записана в следующем более общем виде:

j f(x,

y.)dl = Sign ф — a) j f(x(l), у it)) У [х' (t)]2+ [y' (t)]2dt. (ЗЛ1)

L

а

Здесь использована функция Sign х, определяемая равенством

1,'если х--\> 0;

Signx =

— 1, если х < 0.

В частном случае, если кривая L задана в явной форме у = у (х) (a<ix<i6), то формула (ЗЛО) принимает вид:

 

ь

 

 

 

 

 

l f { x , y ) d l = l f (х, у (х)) У 1 +

[у! (x)]2dx.

(3Л2)

Эта формула справедлива лишь при b

Это ограничение мо­

жет быть снято, если записать

эту

формулу в виде, аналогичном

(ЗЛ1):

ь

 

 

 

 

 

 

 

 

 

j

/ (X, у) dl = Sign (b—a) J

f (x,

у (x)) У \+ [y' (x)]2 dx.

(3.13)

L

a

 

 

 

 

Криволинейный интеграл первого рода по пространственной кривой. Понятие криволинейного интеграла первого рода, введен­ ное для плоской кривой L, легко обобщается на случай, когда L — пространственная кривая, в каждой точке которой задана функция / (х, у, z). Формула (3.2) в этом случае принимает вид:

(7(х, у, z) dl =

П

_

(3.14)

lim V

f(xk, yk, zk) Mk.

L

V °

 

 

Легко видеть, что приведенные выше свойства 1—4 справед­ ливы и для криволинейного интеграла первого рода, вычисленного вдоль пространственной кривой.

Если кривая L задана параметрически

x = x(t), y — y(t), z = z(t) (as^t^P ),

то имеет место формула, аналогичная формуле (ЗЛО):

J / (х, у, z) dl =

L

f>

= j / [х(0, у(1), г Ц ) \ У [x'(t)]2+\y'(t)]2+[z'(t)]2di. (3.15)

а

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

57

имеет место формула, аналогичная формуле (3.11):

J7(x,

у, z) Л = Sign (Р— a) \f[x{t),

y(t), 2(01 X

 

 

X J / V

(0 1 *+IУ' (012+ [2 ' {t)?dt.

(3.16)

Пример 1.

Вычислить

интеграл J = J xydl,

где L — часть

параболы

 

 

L

 

 

у= х2 от точки (0, 0) до точки (2, 4). Воспользуемся формулой (3.12):

 

 

2

_______

2

______

 

 

J =

Jxydl = f хх2 у " 1 +

(2л:)2 d x =

I X *V 1 + 4х2 dx.

 

 

L

0

 

 

0

 

_________

Сделаем в этом интеграле замену

переменных:

У 1

+ 4х2 = г.

откуда

 

 

 

 

 

 

 

 

 

1 4- 4х2 = z2;

xdx =

zdz;

х2 = -------—.

 

 

 

 

 

 

4

 

4

 

 

Получаем:

 

 

 

 

 

 

 

 

У 17

 

 

 

 

 

 

 

= —

( (г2-

1)z2dz ■

 

 

 

 

 

16

J

 

16 V 5

 

 

 

 

289 V

17

1 7 1 7

1

 

391 У 17 — 1

16

5

 

3

 

 

120

 

Пример 2.

Вычислить интеграл

 

 

 

 

 

 

 

 

/ = J x (R 2- y 2)dl,

 

 

 

 

 

 

L

 

 

 

 

 

где L — четверть окружности х2 + У2 =

Я2

от точки

(0,

Я) до точки

(Я, 0).

 

 

 

 

 

 

 

 

Зададим кривую L в параметрической форме:

 

 

 

x = R c o s t ;

y==Rsi nt;

 

л

 

 

 

0 < t < — >

 

 

Вычислим х' (t) и у' (t):

 

 

 

 

 

 

х'

(i) =

— Я sin t\ у' (t)

= Я cos /.

 

 

Воспользуемся формулой (3.10):

 

 

 

 

я

 

 

_______________

 

 

2

 

 

 

 

У = j

Я cos t (Я2 — Я2 Sin2 0 К Я2 Sin21+

Я2 cos21 dt =

 

Я4 J co s <(lsin2 /)

= Я4 j* (1 — sin2 /) d (sin t) ■

 

 

 

sin31

Я4.

 

 

 

 

Я4 ( sin t

3

 

 

 

 

\

 

 

 

 

 

58

Пример 3. Вычислить

интеграл по

пространственной кривой

L:

 

 

/

= j г (л:2 + У2) dl,

 

 

 

где L — часть винтовой

L

 

 

 

 

линии,

 

 

Л

 

x =

3cos^;

г/ =

3 sin ^;

z = 4 t ;

0 < t <

 

— .

 

 

 

 

 

 

 

2

 

Вычислим производные от функций х (t), у (t) и z (t):

 

х'

(t) = — 3siri t\ у' (t)

= 3cos t;

z' (t) =

4.

 

Воспользуемся формулой (3.15):

 

 

 

я

 

 

 

________

 

 

2

 

 

 

 

 

J = j* 4^ (9 cos21+ 9 sin2 t ) V^9 sin2 t +

9 cos2 ^ +

16 dt =

 

о

я

 

я

 

 

 

 

 

 

 

 

 

 

=

180 ГШ = 1 8 0 -^ -| 2

= 180^ - = 22,5я2.

 

 

J

 

2 L

«

 

 

 

о

 

0

 

 

 

 

3.2.

Криволинейные интегралы второго рода

 

Определение

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

Пусть

в плоскости Оху задана кривая L (АВ), в каждой точке которой определена функция f (х, у). Разобьем эту кривую на п частей (ча­ стичных кривых) L 1; L 2, . . . , Ln. Обозначим, как и ранее, через %п— ранг дробления кривой, т. е. наибольшую из длин частичных кривых Lk (k = 1, 2, . . . , п). Пусть ДхА и Аyk — проекции k-ой частичной кривой L на оси координат. В каждой из частичных кри­ вых Lk выберем по точке Р (xk, yk).

Составим интегральную сумму:

П

И Н ч , ук) Ач -

к= 1

Если при стремлении к нулю ранга дробления кривой сущест­ вует конечный предел этой суммы, не зависящий ни от способа раз­ биения кривой L на части, ни от выбора точек Pk (xk, yk), то этот предел называется криволинейным интегралом второго рода по ко­

ординате х от функции / (х,

у) вдоль

кривой L и

обозначается

J / (х, у) dx\ таким образом,

 

 

 

 

 

1

 

 

 

 

 

Г / (х, у) dx =

lim

2

f

Ук) АЧ-

(3-17)

L

 

k=\

 

 

 

Аналогично определяется криволинейный интеграл второго

рода по координате у:

 

 

 

 

 

Г / (х, у) dy =

lim

2

/ (Ч> Ук) АУк-

(3-18)

L

 

Л=1

 

 

 

59

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