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

Математический анализ II Учебное Пособие

.pdf
Скачиваний:
87
Добавлен:
14.04.2015
Размер:
7.4 Mб
Скачать

Вычисляя скалярные произведения ортов ρ0 , φ0 и z0 нетрудно убедиться, что ρ0 φ0 =ρ0 z0 = φ0 z0 = 0 , т.е. можно сделать вывод, что цилиндри-

ческая система координат ортогональна.

Если положить z = 0 , то мы придём к известным уже полярным координатам ρ, ϕ :

x= ρcosϕ

y= ρsinϕ

Нетрудно убедиться, что для полярных и цилиндрических координат выражения для Якобиана преобразования совпадают, т.е.

J (ρ,ϕ,z ) =J (ρ,ϕ) = ρ .

Замечание (о геометрическом смысле якобиана)

Рассмотрим ортогональную криволинейную систему координат Oξη, и пусть декартовы координаты x и y связаны с криволинейными ξ и η со-

отношениями:

x =x (ξ,η), y =y (ξ,η) ,

причём предполагается, что Якобиан преобразования

J (ξ,η) =xξyη′ −xηy ξ′ ≠ 0 .

Рассмотрим в плоскости xO y

простую область D , ограниченную кон-

туром K (рис. 16). В плоскости O1ξη ей соответствует некоторая область

, ограниченная контуром L .

 

η

 

 

y

 

L

 

D

 

 

(ξ,η)

 

 

 

 

0

K

x

01

ξ

 

рис. 16

Если функции x (ξ,η) и y (ξ,η) непрерывны и имеют непрерывные част-

ные производные первого порядка, а также непрерывные смешанные производные, то в этом случае можно доказать, что площади SD и S связаны

соотношением

SD = J (ξ ,η) S ,

где (ξ,η) -некоторая "средняя" точка области . Из этого соотношения вытекает очевидный геометрический смысл якобиана преобразования. А именно: мы можем сказать, что модуль якобиана преобразования J (ξ,η)

при переходе к криволинейным координатам даёт нам коэффициент искажения площадей отображаемых областей.

74

Можно доказать также, что для пространственного случая, когда осуществляется переход от декартовых координат x , y и z к криволинейным ξ, η, ζ при выполнении соответствующих предположений (непрерывность функций, их частных производных и т.п.) имеет место соотношение

vT

=

J (

ξ

,

η

,

ζ

)

vT ` , где vT - объём области T в системе координат O xy z ,

vT `

- объём отображенной области T в системе координат O1ξηζ . В дан-

ном случае, с геометрической точки зрения, на модуль Якобиана преобразования J (ξ,η,ζ ) можно смотреть как на коэффициент искажения объёмов.

7. Замена переменных в тройных интегралах

Пусть осуществляется переход от декартовых координат x и y к криволинейным ξ и η: x =x (ξ,η), y =y (ξ,η).

Предположим, что выполнены условия предыдущего параграфа относи-

тельно функций x (ξ,η), y (ξ,η), ξ =ξ(x ,y )

и η =η(x ,y )

 

и их частных про-

изводных, и имеет силу все сказанное об областях D и

 

. Пусть, кроме то-

го, в области D задана непрерывная функция f (x ,y ) .

 

 

 

 

 

 

 

 

 

 

Рассмотрим двойной интеграл I = ∫∫f (x ,y )dxdy .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Составим интегральную сумму для этого двойного интеграла.

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

σn = f (xk ,y k )

 

 

 

SDk .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k =1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Произвольному разбиению области D на ячейки D1, D 2 ,...,Dn

 

соответству-

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

 

на ячейки 1, 2 ,...,

 

 

n , а точкам (xk ,y k ) Dk соот-

ветствуют точки (ξk ,ηk ) k . Следовательно, можем написать:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

SDk =

J (

ξ

,

η

)

 

S k

,

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а тогда σn = f (x (ξk ,ηk ),y k (ξk ,ηk ))

J (ξk ,ηk )

 

 

 

 

S

k .

 

 

 

 

 

 

 

 

 

 

 

k =1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Измельчая дробление и устремляя ранг дробления к нулю, получим

limσ

n

=

∫∫

f

(

x (ξ,η),y (ξ,η)

)

 

 

J

(

ξ,η

)

 

dξdη .

 

 

 

 

 

 

 

 

 

 

 

 

n →∞

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

λ0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n →∞

n

 

∫∫

 

(

 

)

 

С другой стороны, принимая во внимание,

 

 

=

f

x ,y

dxdy ,

что limσ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

λ0

 

 

D

 

 

 

 

 

окончательно получим

∫∫f (x ,y )dxdy = ∫∫f x (ξ,η),y (ξ,η) J (ξ,η)dξdη.

D

75

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

I = ∫∫(x 2 +xy )dS , где D есть круг x 2 +y 2 R 2 .

D

Решение. Перейдём к полярным координатам x =r cosϕ, y =r sinϕ . В

полярных координатах уравнение окружности r =R

 

при любом ϕ (т.е. ϕ

не меняется от 0 до 2π ), а r является постоянным,

 

J (r ,ϕ)

 

=r , тогда по-

 

 

лучим

 

 

 

 

 

I = ∫∫(x 2 +xy )dxdy = ∫∫(r 2 cos2 ϕ +r cosϕ r sinϕ) r drdϕ =

 

D

 

 

 

 

 

= ∫∫(cos2 ϕ + cosϕsinϕ) r 3drdϕ = π(cos3 ϕ +cosϕsinϕ)dϕRr 3dr

= πR 4 .

0

0

4

Заметим, что для того, чтобы расставить пределы интегрирования, достаточно выяснить, как проходят (возрастают) через область D координатные линии r и ϕ .

Предположим далее, что нужно вычислить тройной интеграл I = ∫∫∫f (x ,y ,z )dxdy dz и пусть осуществляется переход к криволинейным

T

координатам ξ, η, ζ по формулам:

x =x (ξ,η,ζ ), y =y (ξ,η,ζ ), z =z (ξ,η,ζ ) .

Предположим, что относительно функций x (ξ,η,ζ ), y (ξ,η,ζ ) , z (ξ,η,ζ ) , ξ(x ,y ,z ),η(x ,y ,z ),ζ (x ,y ,z ) и их частных производных имеет место сказанное замечании п. 6. Допустим также, что в области T функция f (x ,y ,z )

непрерывна. Тогда аналогично вышесказанному относительно двойного интеграла имеет место соотношение

∫∫∫f (x ,y ,z )dxdydz = ∫∫∫f x (ξ,η,ζ ),y (ξ,η,ζ ),z (ξ,η,ζ ) J (ξ,η,ζ )dξdηdζ ,

T T

где J (ξ,η,ζ ) - якобиан преобразования.

Для расстановки пределов интегрирования по переменным ξ,η,ζ следует выяснить, как проходят координатные линии ξ, η, ζ через область T .

Рассмотрим теперь конкретные примеры.

76

Пример. Вычислить площадь поверхности сферы x 2 +y 2 +z 2 =R 2

(рис. 17).

 

 

z

 

 

 

 

 

 

 

 

 

Решение.

 

Вычислим восьмую

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R

z =

 

R 2 x 2 y 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

часть площади поверхности сферы,

 

 

 

 

 

 

 

 

 

 

 

 

 

лежащую в первом октанте:

R

 

 

 

R

 

 

 

 

 

 

1S = ∫∫ p 2 (x ,y ) +q 2 (x ,y ) +1dxdy

 

 

 

 

y

 

 

 

 

8

D

 

 

 

 

 

 

 

 

 

 

x

 

ϕ

 

 

 

 

 

 

 

 

,где

z (x,y )

 

 

 

 

 

 

z (x,y )

 

 

 

 

 

 

 

 

 

 

 

 

p(x,y ) =

, q(x,y ) =

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

рис. 17

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

В нашем случае z =

 

R 2 x2 y 2 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p (x ,y ) =

 

 

x

 

 

, q(x ,y ) =

 

 

 

y

 

 

 

.

 

 

R 2 x 2 y 2

 

R 2 x 2 y 2

Итак, имеем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

S

= ∫∫

 

 

 

x 2 +y 2

 

+1dxdy = ∫∫

R dxdy

 

 

.

 

 

8

 

R

2

x

2

y

2

R

2

x

2

y

2

 

 

 

D

 

 

 

 

 

 

D

 

 

 

 

 

 

 

Перейдём к полярным координатам

x =r cosϕ, y =r sinϕ, J (r ,ϕ) =r .

Ясно, что

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

π

 

 

 

 

 

 

 

 

 

 

 

 

R dxdy

 

 

 

 

 

 

 

R rdrdϕ

 

2

R

 

 

rdr

 

 

∫∫

 

 

 

 

= ∫∫

 

=R dϕ

 

 

 

.

R

2

x

2

y

2

 

 

R

2

r

2

 

 

R

2

r

2

D

 

 

 

 

D

 

 

 

 

 

 

0

0

 

 

 

 

 

Вычислим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R

 

rdr

 

 

 

 

 

1 R d (R 2 r 2 )

 

 

2

 

2

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

I Вн =

 

 

 

 

 

 

 

 

 

= −

2

 

 

 

 

 

 

 

= − R

 

 

r

 

 

 

=R ,

 

R

2

r

2

 

 

R

2

r

2

 

 

 

 

0

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

тогда

π

S =8R 2 2 dϕ = 4πR 2 кв. ед.

0

77

Пример. Найти объём тела, лежащего в первом октанте и ограниченного снизу конической поверхностью x 2 +y 2 =z 2 , сверху шаровой поверхностью x 2 +y 2 +z 2 =R 2 , а с боков координатными плоскостями x = 0 и y = 0 (рис. 18).

Решение. Искомый объём v = ∫∫∫dxdy dz . Перейдём к сферическим ко-

ординатам

 

T

к.л.ρ

 

z

x =r sinθ cosϕ

 

к.л.ϕ

y =r sinθ sinϕ

 

 

θ

 

 

 

к.л.

z =r cosθ

 

 

y

 

 

 

J (r ,θ,ϕ) =r 2 sinθ

x

Найдём уравнения конуса x 2 +y 2 =z 2

рис. 18

в сферических координатах:

 

(r sinθ cosϕ)2 +(r sinθ sinϕ)2 = (r cosθ )2 ,

откуда следует θ = π и θ =

3π

. Заметим, верхняя чаша конуса имеет урав-

 

 

4

 

 

 

 

4

 

 

 

 

 

 

 

 

 

нение θ = π , а нижняя

θ =

3π

 

. Нетрудно убедиться, что уравнение шаровой

4

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

поверхности x 2 +y 2 +z 2 =R 2

 

в сферических координатах r =R . Итак, ис-

комый объём

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

π

π

 

 

 

 

 

v = ∫∫∫r 2 sinθdrdϕdθ = 2 dϕ4 sinθdθRr 2dr .

 

 

 

T

 

 

 

 

 

 

 

 

0

0

0

 

 

 

 

Вычислим

 

 

 

 

 

 

 

 

 

π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R

2

=

 

R 3

; I 2Вн =

R 3 4

 

2 2

R

3

.

I 1Вн = r dr

 

 

3

 

3

sinθdθ =

6

 

0

 

 

 

 

 

 

0

 

 

 

 

Окончательно

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

π

 

 

 

 

 

 

 

v =

2

2 R

2 dϕ = (2

2 )πR 3

куб. ед.

 

 

 

 

6

 

 

 

 

 

0

 

12

 

 

 

 

 

78

Пример. Вычислить объём

тела, ограниченного поверхностью

x 2 +y 2 = (z 2)2 и плоскостью xO y

(рис. 19).

Решение. Тело симметрично, поэтому вычислим объём его четвертой части, лежащей в первом октанте:

 

z

 

 

1

v = ∫∫∫dxdydz .

 

 

 

 

4

 

(0,0,2)

 

 

T

 

 

Перейдём к цилиндрическим коор-

к.л.z

 

 

 

динатам:

 

 

0

1

2

y

 

 

 

 

 

 

 

 

 

x =r cosϕ

x

2 к.л.ϕ

1

к.л.r

y =r sinϕ

 

 

 

рис. 19

 

z =z

 

 

 

 

 

 

Уравнение конуса в цилиндрических координатах r 2 = (z 2)2 . Часть конуса, ограничивающая тело, имеет уравнение z = 2 r . Тогда будет

 

 

π

 

1 v = ∫∫∫rdr dϕdz = 2 dϕ2 dr

2r rdz .

4

T

0 0

0

Вычислим внутренние интегралы

2r

I 1вн = rdz =r z 02r = 2r r 2 ,

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

2

)dr

 

2

r 3

 

2

=

4

 

 

 

 

 

 

 

 

 

 

 

 

I 2 вн = (2r r

 

= r

 

3

 

 

 

3

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

π

 

 

 

 

 

 

 

 

 

 

 

 

И окончательно:

1 v =

4

2 dϕ =

2π .

 

 

 

 

 

 

 

 

 

 

 

 

 

4

3

0

 

3

 

 

 

 

 

 

 

 

 

 

Ответ:

v =

8π

куб. ед.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

79

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

§ 1. Криволинейные интегралы I рода

1. Определение криволинейного интеграла первого рода

Пусть в плоскости xO y лежит кривая A B , у которой существует длина дуги, и пусть в каждой точке этой кривой определена некоторая функция f (x ,y ) (рис. 1).

Разобьём кривую A B произвольным образом точками, следующими друг за другом от A к B , на n частей. Пусть

Sk - длина дуги M k M k +1 . Наибольшую из длин дуг Sk назовём рангом дробле-

ния λ . На каждой дуге M k M k +1 возьмём произвольную точку Pk (ξk ,ηk ) и

вычислим в ней значение функции f (ξk ,ηk ).

y

 

B =M n

 

 

 

 

M k +1

 

 

Pk (ξk ,ηk )

 

A =M 0

M k

0

x

 

рис. 1

Составим произведение f (ξk ,ηk ) Sk и просуммируем все такие произведения, т.е. образуем интегральную сумму (сумму Римана):

n1

σn = f (ξk ,ηk )ΔSk .

k=1

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

λ 0 : I = limσn .

λ0

Если этот предел существует, то он называется криволинейным ин-

тегралом от функции f (x ,y ) по кривой A B

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

I = f (x ,y )d S ,

 

A B

 

80

 

 

n 1

f (x ,y )dS = nlim→∞

f (ξk ,ηk ) Sk .

A B

λ0

k =0

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

f (x ,y )dS = f (x ,y )dS .

A B B A

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

I = f (x ,y ,z )dS .

A B

Теорема существования криволинейного интеграла I рода

Пусть кривая A B задана параметрическими уравнениями x =ϕ(t ), y =ψ(t ), где функции ϕ(t ) и ψ(t ) определены и непрерывны на промежутке [p ,q ] (p <q ) . Пусть в каждой точке кривой A B определена непрерывная функция f (x ,y ) . Тогда криволинейный интеграл первого рода от функ-

ции f (x ,y ) по кривой A B

существует и выражается через определённый

интеграл так:

q

 

f (x ,y )dS = f [ϕ(t ),ψ(t )] ϕt2 +ψt2dt .

A B

p

Заметим, что для пространственной кривой, заданной параметрическими уравнениями x =ϕ(t ) , y =ψ(t ), z =η(t ) на [p ,q ] , при соблюдении теоремы

мы получим

 

q

f (x ,y ,z )dS = f [ϕ(t ),ψ(t ),η(t )] ϕt2 +ψt2 +ηt2dt

A B

p

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

81

2.Применения криволинейных интегралов первого рода

1.Вычисление длины дуги кривой.

Из определения криволинейного интеграла первого рода следует, что длина дуги кривой A B , по которой ведётся интегрирование, равна

SA B = dS .

AB

2.Вычисление массы материальной кривой.

Найдём массу кривой A B , заданной уравнением y = f (x ) , где f (x ) - непрерывная функция на промежутке [a ,b], если плотность её в каждой точке кривой определена непрерывной функцией ρ = ρ(x ,y ) на том же

промежутке.

Разобьём кривую A B на n частей произвольным образом. Очевидно, что её масса

n 1

M ρ(ξk ,ηk ) Sk ,

k =0

где (ξk ,ηk ) - произвольная точка, принадлежащая каждому отрезку кривой A B , а Sk - длина дуги этого участка. В правой части выражения для M

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

M = ρ(x ,y )dS .

A B

3. Вычисление работы

Найдём работу силы по перемещению материальной точки вдоль кривой A B (рис. 2). Будем рассматривать кривую A B как направленную, тогда и на касательной к ней можно задать направление, совпадающее с

82

y

 

 

направлением

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

пусть

 

 

при движении по кривой от A к B точка M

 

 

 

 

 

 

 

 

предшествует точке N , тогда будем считать,

 

 

 

что секущая, проходящая от

M

к

N ,

имеет

 

N

 

направление, совпадающее с направлением

M

 

кривой A B . Касательная

τ

есть предельное

 

 

положение секущей при условии,

что N M .

A

 

 

 

x

Тогда направление с секущей переходит на ка-

0

 

 

рис. 2

 

сательную, и про такую касательную говорят,

 

 

что она имеет направление, совпадающее с на-

 

 

 

 

 

 

правлением кривой A B .

 

 

 

 

 

 

 

 

 

Найдём работу силы F по перемещению материальной точки M из точ-

ки A в точку B вдоль кривой (рис. 3).

 

 

 

 

 

 

 

 

 

 

y

 

 

Разобьём кривую A B

произвольным обра-

 

B зом на n частей и будем считать, что каждый

 

Fk

 

k -й участок кривой - прямолинейный отрезок,

 

θk

 

длина которого

Sk , сила F на каждом таком

 

 

участке постоянная по величине и направле-

 

 

τ

A

M

нию, и угол между силой F и участком кривой

x

равен θk . Тогда элементарная работа силы F на

0

 

рассматриваемом участке

A

k

=F

S

k

cosθ

k

.

 

рис. 3

 

Отсюда вся работа

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A = F cosθdS .

AB

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

1.Определение криволинейного интеграла второго рода

Пусть в плоскости xO y

лежит кривая A B , в каждой точке которой оп-

ределена функция f (x ,y ) (рис. 1).

Разобьём кривую A B

произвольным образом точками M 0 =A ,

M 1,M 2 ,... , M n , следующими друг за другом, на n частей, пусть d1 , d 2 , …,

dn - диаметры дуг M 0M 1 , M 1M 2 ,...,M n 1M n , наибольший из диаметров dk , равный λ , есть ранг дробления.

На каждой дуге M kM k +1 возьмём произвольную точку Pk (ξk ,ηk ) и вы-

числим в

ней значение функции f (ξk ,ηk ), т.е. составим произведение

f (ξk ,ηk )

xk .

Просуммируем все такие произведения, т.е. составим интегральную сумму (сумму Римана):

83