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

Математика-2-й семестр (курс лекций)

..pdf
Скачиваний:
8
Добавлен:
05.02.2023
Размер:
11.03 Mб
Скачать

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

Перейдем вначале к повторному интегралу типа (5.3). Тогда a 0; b 0,5; y1 (x) x; y2 (x) 1 x . Поэтому

 

 

 

0,5

 

1 x

 

 

f (x, y)dxdy dx

f (x, y)d y .

 

D

 

 

0

 

 

x

 

 

Для перехода к интегралу типа (5.4)

требуется разбить область на две:

D1 c гра-

ницами c 0;

d

0,5; x1

( y) 0;

x1 (y) y и

 

1

1

1

 

 

 

 

2

D

с границами

c 0,5;

 

d

2

1;

x2 ( y) 0;

2

 

 

2

 

 

 

1

x2

(y) 1 y . Поэтому

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

0,5

y

 

1

 

 

1 y

 

f (x, y)dxdy dy f (x, y)dx

dy f (x, y)dx.

D

0

0

 

0,5

 

0

 

Пример 3. Пусть область D -

внутренность тре-

угольника с вершинами A(1,1) ,

B(5,5) ,

C(3, 2) . В двойном

интеграле f (x, y)dxdy перейти к повторным и расставить

D

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

Найдём уравнения прямых AB , BC , AC . Уравнение прямой AB можно записать в виде

 

 

 

 

 

x 1

 

y 1

 

или, что то же самое, в форме

 

 

 

 

 

1

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y x ;

прямой

AC в форме

x 1

 

y 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

1

 

 

 

 

 

или

y 1 x 1

; прямой CB

в

 

виде

 

 

 

 

 

 

 

 

2

 

2

 

 

 

 

 

 

x 3

 

y 2

или y

3 x 5

. Как для перехода к интегралу

 

 

 

2

 

3

 

 

 

2

2

 

 

 

 

 

 

 

вида (5.3), так и для перехода к интегралу вида (5.4) приходится разбивать область на две. Для интеграла вида: (5.3) соответствующие области задаются неравенствами:

311

D1

1 x 3;

0,5x 0,5 y x ,

D2

 

 

3 x 5;

1,5x 2,5 y x . Таким образом

 

 

 

 

 

 

3

x

 

5

 

x

 

 

f (x, y)dxdy dx

 

f (x, y)dy dx

 

f (x, y)dy.

 

D

1

0,5 x 0,5

 

3

1,5 x 2,5

 

Расставить пределы интегрирования, взяв внешний интеграл по y (то есть представить двойной интеграл в

виде повторного интеграла вида (5.4)), предлагается самостоятельно.

Пример 4. Пусть область D задана неравенствами y x2 , y x . Тогда

 

 

 

 

1

 

x

 

 

 

1

 

 

y

 

f (x, y)dx .

f (x, y)dxdy dx

f (x, y)dy dy

D

 

 

 

0

x2

 

 

 

0

 

y2

 

 

Пример 5. Изменить порядок интегрирования в ин-

теграле

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 x

 

 

 

1

 

x

 

 

 

 

2

 

 

f (x, y)dxdy dx

f (x, y)dy

dx

f (x, y)dx

 

D

0

0

 

 

 

 

1

0

 

 

 

 

 

 

Исходная область представлена в виде объединения

двух областей: D1

 

0 x 1; 0 y

 

 

и

D2

x

1 x 2; 0 y 2 x .

 

Таким

образом,

эта

 

область

ограничена

кривыми

y

 

 

 

y 2 x

 

 

x 0 .

 

 

 

x ,

и

 

Её

также

можно

 

 

задать

 

 

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

0 y 1,

y2 x 2 y . Поэтому

 

 

 

 

 

 

 

1

2 y

 

 

 

 

 

 

 

 

 

f (x, y)dxdy dy

f (x, y)dx.

 

 

D

 

 

 

 

 

0

y2

 

 

 

 

 

 

 

 

 

312

5.2.2.2. Вычисление тройных интегралов

Аналогично случаю двойного интеграла доказывается, что если V [a,b] [c, d] [e, f ] - параллелепипед, то

 

 

b d f

 

 

f (x, y, z)dxdydz

 

 

 

 

 

 

 

 

 

 

 

 

f (x, y, z)dz dy dx.

V

 

a c e

 

 

Пусть теперь V - область, расположенная между плоскостями x a, x b и для x [a,b] область V одно-

значно проектируется на плоскость YOZ и D - эта проекция. Тогда

 

b

 

 

f (x, y, z)dxdydz

f (x, y, z)dydz dx .

V

a

D

 

Если V - цилиндр с образующими, параллельными оси OZ , направляющей, лежащей в плоскости XOY и являющейся границей области D , ограниченный поверхно-

стями z z1 (x, y), z z2 (x, y) , то

 

 

 

 

z2

( x, y )

 

f (x, y, z)dxdydz dxdy

 

f (x, y, z) dz .

V

D

 

z1 ( x, y )

 

Пример 1. Пусть область V ограничена поверхно-

стями x 0 ,

y 0 , z 0, x 4, y 4,

z x2

y2 1. В трой-

ном интеграле f (x, y, z)dxdydz

перейти к повторным и

 

V

 

 

 

 

 

расставить пределы интегрирования.

 

 

 

Данная область есть цилиндр, ограниченный по-

верхностями

z 0, z x2

y2 1. Проекция этого цилиндра

на плоскость

XOY есть

квадрат

с границей x 0 , y 0 ,

x 4 , y 4 ,

которая одновременно является направляю-

щей цилиндра. Поэтому

 

 

 

 

 

 

 

4

4

x2 y2 1

 

f (x, y, z)dxdydz dx dy

 

f (x, y, z)dz.

 

V

0

0

 

0

 

313

Пример 2. Область V ограничена поверхностями

y 0,

z 0,

3x y 6,

3x 2y 12,

x y z 6.

В тройном

 

 

интеграле

f (x, y, z)dxdydz

перейти к

 

 

 

 

 

 

V

 

 

 

 

 

 

повторным и расставить пределы ин-

 

 

тегрирования.

 

 

 

 

 

 

 

 

Область однозначно проектиру-

 

 

ется

на

треугольник

y 0,

 

 

3x y 6, 3x 2y 12,

лежащий в плос-

 

 

кости

XOY , является цилиндром,

ограни-

 

 

ченным

 

поверхностями

z 0,

z 6 x y,

 

 

направляющая

которого

есть

указанный

 

 

выше треугольник. Поэтому

 

 

 

 

 

 

12 2 y

 

 

 

 

 

 

 

 

6

 

3

 

 

6 x y

 

 

 

 

f (x, y, z)dxdydz dy

dx

f (x, y, z)dz.

 

 

V

 

 

0

 

6 y

0

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

5.2.3. Замена переменных в кратных интегралах

5.2.3.1. Криволинейные системы координат

Положение точки на прямой, на плоскости, в R3 и в Rn можно определить различными способами. В частности, это можно сделать, задав её декартовы координаты. Иногда же бывает удобно фиксировать положение точки при помощи других величин, например, связанных с решаемой задачей. Выяснением этих вопросов для общего случая мы и займёмся.

Пусть D, D1 Rn - области, r : D1 D отображение

314

 

 

 

 

x (u ,u

 

,...,u

 

)

 

 

 

 

1

1

 

2

 

n

 

x r(u) r(u ,u

2

,...,u

n

) x2 (u1

,u2 ,...,un ) .

1

 

........................

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xn

(u1

,u2

,...,un )

Если r - биективное (взаимно однозначное) отображение, то будем говорить, что задана криволинейная система координат, так как в этом случае положение точки x D однозначно определяется точкой u D1 . Если вектор-

функция r дифференцируема, то криволинейную систему координат будем называть регулярной. Заметим, что в этом случае, по теореме о производной обратной функции [раздел 2.5], обратное отображение, осуществляемое век- тор-функцией r 1 , дифференцируемо.

Система вектор-функций x r(u) r(u1 ,u2 ,...,un ) при ul const образует, как и в случае декартовых координат,

систему координатных поверхностей. Пересечения координатных поверхностей образуют координатные поверхности меньшей размерности. В частности при n 2 , отображение r r(u,v) задаёт криволинейную систему координат на плоскости, а кривые

(x, y)T r(u,C2 ) x(u,C2 )i y(u,C2 )j, (x, y)T r(C1,v) x(C1,v)i y(C1,v)j ,

образуют координатные линии. Аналогично при n 3 отображение r r(u,v, w) задаёт криволинейную систему ко-

ординат в пространстве R3 , поверхности

(x, y, z)T r(u,v,C3 ) x(u,v,C3 )i y(u,v,C3 )j z(u,v,C3 )k ,

(x, y, z)T r(u,C2 , w) x(u,C2 , w)i y(u,C2 , w)j z(u,C2 , w)k ,

(x, y, z)T r(C1,v, w) x(C1,v, w)i y(C1,v, w)j z(C1,v, w)k ,

образуют координатные поверхности, а их пересечения, то есть кривые

315

(x, y, z)T r(u,C2 ,C3 ) x(u,C2 ,C3 )i y(u,C2 ,C3 )j z(u,C2 ,C3 )k,

(x, y, z)T r(C1,v,C3 ) x(C1,v,C3 )i y(C1,v,C3 )j z(C1,v,C3 )k,

(x, y, z)T r(C1,C2 , w) x(C1,C2 , w)i y(C1,C2 , w)j z(C1,C2 , w)k,

образуют систему координатных линий.

Длины векторов rul ,l 1, 2,...,n, то есть числа hi rul , l 1, 2,..., n, называются коэффициентами Ламе криволи-

нейной системы координат. Если векторы rul ,l 1, 2,...,n, по-

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

ли перпендикулярны векторы ru (u, v) ,

rv (u, v) . Аналогично

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

ной, если перпендикулярны векторы

ru (u,v, w) , rv (u,v, w) ,

 

плоскости равны

rw (u,v, w) . Коэффициенты Ламе на

hu

 

2

 

 

2

 

, hv

 

 

 

 

 

 

2

 

 

2

(xu )

 

 

( yu )

 

 

 

 

 

 

(xv )

 

( yv )

 

 

hu

 

 

2

 

 

2

 

 

 

2

,

 

(xu )

 

 

 

( yu )

 

 

(zu )

 

 

hw

 

 

2

 

 

 

2

 

 

2

 

.

 

 

 

 

 

(xw )

 

( yw )

 

 

(zw )

 

 

 

 

 

 

, а в R3 соответственно

hv

 

2

 

2

 

2

,

(xv )

 

( yv )

 

(zv )

 

Заметим, что для ортогональной криволинейной системы координат модуль определителя матрицы Якоби r (производной матрицы) [раздел 2.1] равен произведению коэффициентов Ламе.

5.2.3.2. Полярная система координат на плоскости

Наиболее часто используемой криволинейной системой координат на плоскости является полярная система координат. Положение точки в этой системе координат определяется длиной радиус-вектора точки и углом ме-

жду радиус-вектором точки и осью. Если в роли оси полярной системы взять ось OX , то в координатном виде пе-

316

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

x cos ,

ся по формулам . В векторной форме то же са-

y sin ,

мое записывается в виде

x

x( , )

 

cos

( cos ) i ( sin ) j .

 

r( , )

 

 

 

y

y( , )

 

sin

 

Угол при этом может быть выбран из любого по-

луинтервала длиной

2 . Чаще всего бе-

рут полуинтервалы

[0, 2 ) ,

[

 

,

3

) ,

2

2

 

 

 

 

 

[ , ) . Полярная система координат яв-

ляется ортогональной. Действительно, вычисляя скалярное произведение векторов

r (cos , sin )T , r ( sin , cos )T ,

получаем требуемое. Коэффициенты Ламе для полярной системы координат равны h 1 , h .

5.2.3.3. Сферическая и цилиндрическая системы

координат в R3

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

жду радиус-вектором точки и осью OZ , углом между проекцией радиус-вектора

точки на плоскость XOY и осью OX. Формулы перехода в координатной форме приобретают вид

317

x cos sin ,y sin sin ,

z cos .

При этом 0 , 0 2 , 0 . В векторной форме то же самое записывается в виде

xyz

x( , , )

 

cos sin

 

 

 

 

 

 

 

r( , , ) y( , , )

 

sin sin

 

 

 

cos

 

 

z( , , )

 

 

 

( cos sin )i ( sin sin ) j cos k .

Сферическая система координат является ортогональной. Действительно, вычисляя скалярное произведение векторов

r (cos sin , sin sin , cos )T , r ( sin sin , cos sin ,0)T ,

r ( cos cos , sin cos , sin )T ,

получаем требуемое. Коэффициенты Ламе для сферической системы координат равны h 1 , h sin , h .

Второе обобщение полярной системы координат называется цилиндрической системой координат. Положение точки в этой системе координат определяется длиной проекции

радиус-вектора точки на плоскость XOY ,

углом между этой проекцией и осью OX ,

координатой z . Формулы перехода в координатной форме приобретают вид

x cos ,y sin ,z z.

При этом 0 , 0 2 , z . В векторной форме то же самое записывается в виде

318

x

x( , , z)

 

cos

 

 

 

 

 

 

 

 

 

y

( , , z) y( , , z)

 

sin

 

 

 

 

 

z

 

 

z

z( , , z)

 

 

 

( cos )i ( sin )j zk .

Цилиндрическая система координат также ортогональна. Предлагается проверить это самим. Коэффициенты Ламе для цилиндрической системы координат равны h 1 , h , hz 1.

5.2.3.4. Замена переменных в интегралах Теорема 5.15. Пусть f (x) f (x1, x2 ,..., xn ) - функция,

заданная в области D Rn , r : D1 D - биективное (осуще-

ствляющее взаимно однозначное соответствие) дифференцируемое отображение,

 

 

x (u ,u

 

,...,u

 

)

 

 

1

1

 

2

 

n

 

x r(u) r(u1 ,u2

,...,un )

x2 (u1

,u2 ,...,un )

........................

 

 

 

 

 

 

 

 

 

 

 

 

 

xn

(u1

,u2

,...,un )

.

Тогда

f (x)dx f ((u)) r (u) du,

D

D1

где r (u) - модуль якобиана r (u) (определителя матри-

цы Якоби или, что то же самое, производной матрицы r (u) ).

Доказательство. Пусть n 2 . Тогда взаимно одно-

значное

дифференцируемое отображение

D1 в D можно

записать

в

виде

r r(u,v) (x(u,v), y(u,v))T .

Разобьём об-

ласть

D1

на

части прямыми u uk , k 1, 2,..., n,

319

v vl ,l 1, 2,..., m,

параллельны-

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

разбиению соответствует раз-

биение области

D кривыми

 

 

 

r r(u ,v) (x(u ,v), y(u ,v))T ,

 

 

 

 

k

k

k

k 1,2,...,n,

r r(u,v ) (x(u,v ), y(u,v ))T ,

l 1, 2,..., m . При

 

 

l

l

l

 

 

этом прямоугольник D1kl с

вершинами

(uk ,vl ),

(uk 1 ,vl ),

(uk ,vl 1 ),

(uk 1 ,vl 1 ) перейдёт

в

криволинейный

четырёх-

угольник

 

Dkl , ограниченный

линиями

r(uk , v),

r(u,vl ),

r(uk 1 ,v), r(u,vl 1 ) .

Пусть

(uk , vl ) -

точка прямоугольника D1kl ,

xkl x(uk , vl ) ,

ykl y(uk ,vl ) . Рассмотрим интегральную сум-

му

 

 

 

n

m

 

S f (xkl , ykl )) (Dkl )

 

k 1

l 1

для вычисления интеграла от функции f по области D , в

 

 

которой (Dkl ) - площадь четы-

 

 

рёхугольника Dkl . Из геометри-

 

 

ческого

смысла

производной

 

 

[раздел 2.10] следует, что вектор

 

 

ru (uk ,vl )

является касательным к

кривой r(u, vl )

в точке (uk , vl ) , а вектор rv (uk ,vl ) будет каса-

тельным вектором кривой r(uk ,v) в той же точке. Далее,

r(uk 1 ,vl ) r(uk ,vl ) ru (uk ,vl ) uk ru (uk ,vl ) uk 1 ( uk ) ,

r(uk ,vl 1 ) r(uk ,vl ) rv (uk ,vl ) vl rv (uk ,vl ) vl 2 ( vl ) ,

где 1 ( uk ) и

2 ( vl )

- бесконечно малые более высокого

порядка малости, чем

uk и

vl . Можно показать, что

площади криволинейного четырёхугольника

Dkl и парал-

лелограмма,

построенного

на

векторах

320