книги из ГПНТБ / Дьяченко, А. Н. Интегральное исчисление функций нескольких переменных. Теория вероятностей и элементы математической статистики учебное пособие
.pdfфункция, определенная в области У). Аналогично тому, как это было сделано для плоских областей (см. § 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, к цилиндрическим координатам:
2я |
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, не имеющих общих внутрен них точек. Обозначим длины этих частичных кривых соответственно А1Ъ А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± или кривой 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