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

Основы математики. Неопределенный и определенный интеграл

.pdf
Скачиваний:
207
Добавлен:
14.03.2016
Размер:
2.74 Mб
Скачать

Примеры

Вычислить

1. (x2 2xz)dx + (z2 2xz)dz,

l

где l – дуга параболы z = x2, пробегаемая от точки A(1,1) до точки B(1,1)

(рис. 13).

Решение. На плоскости xoz задано силовое (векторное) поле

F = ((x2 2xz),(z2 2xz)).

z

x

Рис. 13

Вычислим в некоторых точках значение силового вектора (рис. 13)

F(A) = {3,3}, F(0) = {0,0}, F(B) = {1,1}.

По условию линия интегрирования l задана явно уравнением z = x2,

тогда z′ = 2x. При движении из точки А в точку В аргумент х изменяется от –1 до 1. По формуле (4.5) получаем

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(x

2

2xz)dx + (z

2

2xz)dz = ((x

2

2x x

2

 

 

2

)

2

2x x

2

 

 

 

 

 

 

 

 

 

) + (x

 

 

 

 

2x)dx =

 

l

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

4

 

 

 

6

 

4x

5

 

 

1

= −14.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= (x2 2x3 + 2x5 4x4 )dx =

 

x

 

x

 

+

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

3

 

2

3

5

 

 

 

1

15

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Механический смысл результата есть работа по перемещению мате-

риальной точки из положения А в положение В в силовом поле F.

 

 

 

2. Вычислить

xdx + ydy + (x + y 1)dz

по

отрезку

 

прямой,

соеди-

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

няющей точки A(1,1,1) и B(2,3,4).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Решение

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Уравнение заданной прямой, проходящей через две точки, имеет вид:

 

 

 

 

 

 

x 1

=

y 1

=

z 1

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

60

Представим уравнение прямой l в параметрическом виде:

x = t +1,

 

 

y = 2t +1,

 

= 3t +1,

z

 

t ≤1,

0

dx = dt;

тогда = 2 ;

dy dtdz = 3dt.

При передвижении от точки В к точке А параметр t меняется от 1 до 0. По формуле (4.6) и по свойству определенного интеграла получим

 

0

 

 

 

 

 

 

 

xdx + ydy + (x + y 1)dz = (t +1)dt + (2t +1)2dt + (t +1+ 2t +11)3dt =

 

 

l

1

 

 

 

 

 

 

 

 

1

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= − (14t + 6)dt = −(7t2 + 6t)

 

= −13.

 

0

 

 

 

 

0

 

 

 

 

 

 

3. Вычислить криволинейный инте-

 

 

 

 

 

у

 

 

 

 

 

 

 

 

 

 

 

грал y2dx + x2dy, где l – верхняя половина

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

эллипса

x = acost,, пробегаемая по часо-

 

 

 

 

 

 

 

y = bsint

 

 

 

 

 

 

вой стрелке (рис. 14).

 

 

 

 

 

 

 

 

 

 

 

х

 

 

 

 

 

 

 

 

Решение. По формуле (4.6) имеем

Рис. 14

 

 

 

 

y2dx + x2dy = 0 (b2 sin2 t(asint) + a2 cos2 tbcost)dt =

 

 

 

 

 

 

i

π

 

 

 

 

 

 

 

0

 

0

 

 

4ab2.

 

= −ab2 sin3 tdt + a

2bcos3 tdt =

 

π

π

3

 

 

 

 

4.3.Приложения криволинейных интегралов

4.3.1.Вычисление площади плоской области

Площадь области D, ограниченной замкнутым контуром l, находится по формуле

61

S =

1

xdy ydx,

(4.9)

 

2

l

 

 

 

 

где направление обхода контура l выбрано так, что область D остается все

время слева от пути интегрирования.

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

Пример

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Вычислить площадь области D,

 

 

 

 

 

 

 

 

 

 

x = acost,

 

 

 

 

 

 

 

 

 

 

ограниченной эллипсом y = bsint,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0 ≤ t ≤ 2π.

 

 

 

 

 

 

 

 

 

 

 

 

Решение

 

 

 

 

 

 

 

Рис. 15

 

 

 

 

 

 

 

 

x = acost,

 

 

 

 

 

 

 

 

 

 

 

тогда dx = −asintdt,

 

 

 

 

 

 

По условию имеем y = bsint,

 

 

 

 

 

 

 

 

 

 

dy = bcostdt.

 

 

 

 

 

 

 

 

 

 

0 ≤ t ≤ 2π,

 

 

 

 

 

 

 

 

 

 

 

 

Используя формулу (4.9) получаем

 

 

 

 

 

 

 

1

xdy ydx =

1

2π

 

 

 

 

 

 

 

 

 

 

 

S =

(acostbcost + bsintasint)dt =

 

 

 

 

 

 

 

2

l

2

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ab

2π

 

ab

 

2π

ab

 

2π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

(cos2 t + sin2 t)dt =

 

dt =

t

= πab.

 

 

 

 

 

 

 

2

 

 

 

 

2

0

2

 

0

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4.3.2. Вычисление массы неоднородной кривой

Пусть средняя линия неоднородного стержня описывается уравнением линии l. Плотность стержня задается функцией (M ) , тогда масса m

стержня l вычисляется по формуле

m = µ(M )dl.

(4.10)

l

Пример

x = cost,

Найти массу четверти эллипса , расположенной в первом

z = 2sint

квадранте плоскости xoz, если линейная плотность массы (M ) = z.

62

Решение. Из формулы (4.10) сле-

дует, что m = zdl . Из уравнения кри-

l

вой l находим dl = sin2 t + 4cos2 tdt.

Очевидно, что параметр t меняется от 0 до π/2, тогда масса дуги

z

x

Рис. 16

 

π/2

 

π/2

 

 

m = zdl = 2sint

sin2 t + 4cos2 tdt = −2 1+ 3cos2 td(cost).

l

0

0

 

 

Используя метод замены переменной cost = u и определив новые пределы интегрирования, получим

0

 

 

1

 

 

 

 

 

 

 

 

3

 

ln(

 

+ 2).

 

m = −2

 

1/ 3 + u2 du = 2

 

1/ 3 + u

2 du = 2 +

 

 

 

3

3

3

 

 

1

 

0

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4.3.3. Вычисление площади цилиндрической поверхности

Пример

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Найти площадь цилиндрической

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

 

 

поверхности x = cost, , расположенной в

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z = 2sint

 

 

 

 

 

 

 

 

 

 

 

 

 

первом октанте и заключенной между

 

 

 

 

 

 

 

 

 

 

 

 

 

плоскостями y = 0,

z = 0 (рис. 17).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Из формулы

(4.10) следует, что

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

S = zdl . Из уравнения кривой l находим

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 17

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dl = sin2 t + 4cos2 tdt.

Очевидно, что параметр t меняется от 0 до π/2, тогда

 

π/2

 

π/2

 

 

S = zdl = 2sint

sin2 t + 4cos2 tdt = −2 1+ 3cos2 td(cost).

l

0

0

 

 

63

Рис. 18

Используя метод замены переменной cost = u и оп ределив новые пределы интегрирования, получим

0

 

1

 

 

 

3

 

ln(

 

+ 2).

S = −2

 

1/ 3 + u2 du = 2

 

1/ 3+ u2 du = 2 +

 

 

3

3

3

 

 

1

0

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

4.3.4.Вычисление статических моментов

ицентра тяжести кривой

Пусть масса распределения вдоль кривой АВ с плотностью µ(x, y). Тогда статические моменты относительно координатных осей определяются формулами

Mx = yµ(x, y)dl,

M y = xµ(x, y)dl.

(4.11)

l

l

 

Координаты центрра тяжести кривой l вычисляются по формулам

xc = My , m

Пример Найти координаты центра

однородной дуги циклоиды

x = a(t sint),

= (1cos ),

y a t

0 t 2π.

yc = Mx . m

тяжести у

(4.12)

х

Решение. Пусть точка (xc , yc ) – центр

тяжести (ЦТ) дуги. В силу симметрии кривой относительно прямой x = πa

получаем

 

xc

= πa . Найдем массу m. Из

уравнения циклоиды находим:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dl = a 2

1costdt, тогда масса

 

 

 

 

 

 

2π

 

 

 

 

 

 

 

 

2π

 

 

2π

t

 

t

 

2π

 

 

 

 

 

 

 

 

 

 

m = a

2

 

 

1costdt =

2

a

1costdt = 2a sin

dt = −4acos

 

= 8a.

 

 

2

 

0

 

 

 

 

0

 

 

0

2

 

0

 

 

 

 

 

 

 

 

 

 

Вычислим статичееский момент дуги относительно осии ох

64

 

 

 

2π

 

 

 

 

 

 

 

 

 

 

 

 

2π

 

t

 

32a2.

M =

ydx =

a(1cost)a

 

 

2

 

 

1costdt = 4a2

sin3

dt =

 

 

 

 

2

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

 

π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

Определяем вторую координату yc

по формуле (4.12)

 

 

 

 

 

y =

Mx

=

32

a2

:8a =

4

a.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c

m

3

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Итак, xc

= πa,

 

yc = 4a / 3.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4.3.5. Вычисление моментов инерции кривой

 

Моменты инерции Ix , Iy

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

осей ох, оу и начала координат определяются соответствующими формулами

Ix = y2µ(x, y)dl,

Iy = x2µ(x, y)dl,

I0 = (x2 + y2 )µ(x, y)dl.

(4.13)

l

l

l

 

Нетрудно заметить, что момент инерции относительно начала координат равен сумме моментов относительно осей координат I0 = Ix + Iy.

Пример Найти моменты инерции относительно координатных осей и начала

y = 2cost,

координат четверти однородной окружности , лежащей в пер-

z = 2sint

вом квадранте плоскости yoz.

Решение. В силу одинакового расположения однородной кривой по отношению к координатным осям имеем Iy = Iz. По формулам (4.13) получаем

π/2

Iy = Iz = z2dl = 4 sin2 t4sin2 t + 4cos2 tdt =

l

0

= 8

π/2

1cos2t

dt = 4

 

sin2t

 

π/2

 

2

t

2

 

 

= 2π

 

 

 

 

 

 

 

0

 

 

 

 

 

 

0

 

 

 

 

 

 

 

65

I0 = (y2 + z2 )ds = 4π/2(cos2 t + sin2 t)4sin2 t + 4cos2 tdt =

l

0

π/2

 

π/2

 

 

 

= 8 dt = 8t

 

= 4π.

0

 

0

 

Итак, Iy = Iz = 2π; I0 = 4π.

4.3.6. Вычисление работы силового поля

Пусть F = P(x, y)i + Q(x, y) j есть переменная сила. Функции P(x, y) и Q(x, y) непрерывны на кривой l. Работа, совершенная силой F вдоль пути l, определяется по формуле

A = P(x, y)dx + Q(x, y)dy.

(4.14)

l

 

Пример

 

Вычислить работу силы F = yzi + (y + z) j

при перемещении точки

массы m из точки O(0,0) в точку A(1,1) по прямой z = y, лежащей в плоскости yoz.

Решение. Из формулы (4.14) следует, что A = yzdy + (y + z)dz. Так

l

как мы интегрируем по прямой z = y и при перемещении из точки О в точку А у меняется от 0 до 1, получаем

1

1

 

y

3

 

 

 

1

 

4.

 

 

 

A = y y dy + (y + y)dy = (y2

+ 2y)dy =

 

+ y2

 

 

 

=

 

 

0

0

 

3

 

 

 

0

 

3

 

 

 

 

 

 

66

4.3.7. Условия независимости криволинейного интеграла по координатам от пути интегрирования

В некоторой односвязной плоской области G следующие утверждения равносильны:

1. Интеграл P(x, y)dx + Q(x, y)dy , взятый по любому замкнутому

l

контуру l, лежащему в G, равен нулю.

2. Интеграл P(x, y)dx + Q(x, y)dy не зависит от пути интегрирова-

lAB

ния, а зависит только от положения начальной А и конечной В точек.

3. Выражение Pdx + Qdy представляет собой полный дифференциал

некоторой функции U(x, y), однозначно определенной в области G.

4. Всюду в области G выполняется равенство P = Q .

y x

Условие независимости криволинейного интеграла второго рода от пути интегрирования можно сформулировать с помощью теоремы.

Теорема. Пусть функции P(x, y) и Q(x, y) непрерывны вместе со своими первыми частными производными в области G, ограниченной контуром l. Тогда для того, чтобы криволинейный интеграл

P(x, y)dx + Q(x, y)dy

l

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

P Q

всех точках области G выполнялось равенство y = x .

Эту теорему можно сформулировать иначе:

Криволинейный интеграл P(x, y)dx + Q(x, y)dy не зависит от пути

l

интегрирования тогда и только тогда, когда его подынтегральное выражение есть полный дифференциал некоторой функции U(x, y).

67

4.3.8. Нахождение функции по ее полному дифференциалу

Если в области G известен полный дифференциал некоторой функции U(x, y), т. е. P(x, y)dx + Q(x, y)dy = dU(x, y), то функцию U(x, y) можно найти, вычисляя соответствующий криволинейный интеграл

(x,y)

P(x, y)dx + Q(x, y)dy,

(x0 ,y0 )

где (x0 , y0 ) – координаты началь-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ной точки А; (x, y) – координаты

 

у

 

 

 

 

 

 

 

 

 

 

 

 

N(x0y)

 

B(x,y)

 

 

 

 

 

 

 

конечной точки В линии интегри-

 

 

 

 

 

 

 

 

 

 

рования, причем от точки А до точ-

 

 

 

 

 

 

 

 

 

 

ки В интегрирование можно прово-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

дить, например,

вдоль ломаной

 

 

 

 

A(x0y0)

 

 

K(x,y0)

 

 

 

 

 

 

 

 

 

 

 

 

 

АКВ (рис. 19), звенья которой АК и

 

0

 

 

 

 

 

 

 

х

 

 

 

 

 

 

 

 

 

 

КВ параллельны

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

 

 

 

 

 

 

Рис. 19

Точка К имеет координаты х и у0.

 

 

 

 

 

 

 

 

 

 

Тогда имеем: вдоль АК y = y0 ; dy = 0 , вдоль КВ x = const; dx = 0 .

Подставляя в исходный интеграл, получим:

 

 

 

 

(x,y)

 

 

 

 

x

y

U(x, y) = P(x, y)dx + Q(x, y)dy =P(x, y0 )dx +Q(x, y)dy +C .

(x0 ,y0 )

 

 

 

 

x0

y0

Аналогично можно начальную точку A(x0 , y0 ) и точку с текущими координатами B(x, y) соединить ломаной ANB (криволинейный интеграл от этого не изменится, так как подынтегральное выражение является полным дифференциалом dU).

При этом получим:

(x,y)

y

x

U(x, y) = P(x, y)dx + Q(x, y)dy = Q(x0, y)dy +P(x, y)dx +C.

(x0 ,y0 )

y0

x0

68

Пример. Найти первообразную функцию U(x, y), если известно,

что выражение (ey + x)dx + (xey 2y)dy является полным дифференциа-

лом искомой функции, т. е. dU = (ey + x)dx + (xey 2y)dy .

Решение. Выберем в качестве начальной точки начало координат, тогда получим:

(x,y)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

y

 

 

 

U(x, y) =

(ey + x)dx + (xey 2y)dy = (1+ x)dx + (xey 2y)dy + C =

 

(0,0)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

0

 

 

 

 

(1+ x)

2

 

x

 

 

 

 

 

 

 

 

 

 

 

(1+ x)

2

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

+ (xey y2 )

 

y

+ C =

 

 

 

 

 

2 + xey x y2 + C =

 

 

2

 

 

 

0

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

1+ 2x + x2

 

 

 

1

 

 

 

 

 

 

2 x + C =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

+ xey y

 

 

 

 

 

 

 

 

2

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

1

+ x +

x

2

1

+ xey

y2

x + C =

x2

+ xey y2

+ C .

 

 

 

 

 

 

2

2

 

2

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Окончательно получаем: U(x, y) =

x2

 

+ xey

y2 + C.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

Задания для самостоятельного решения

Задача 1.

Вычислить криволинейные интегралы I рода:

1.1. y2dl, l : x = a(t sint), y = a(1cost), 0 t 2π.

l

1.2. (x + y)dl, где l – контур треугольника с вершинами O(0,0),

l

A(1,0), В (0,1).

1.3. (x2 + y2 + z2 )dl, l : x = acost, y = asint, z = bt, 0 t 2π.

l

1.4. zdl, l : x = tcost, y = tsint, z = t, 0 t ≤ π.

l

69