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

chirskii-lectures-4sem

.pdf
Скачиваний:
60
Добавлен:
08.04.2015
Размер:
3.67 Mб
Скачать

x(u,v + ∆v)x(u,v)= xv (u,v)v + o(v) y(u,v + ∆v)y(u,v)= yv (u,v)v + o(v) x(u + ∆u,v)x(u,v)= ux (u,v)u +o(u) y(u + ∆u,v)y(u,v)= uy (u,v)u + o(u)

x(u + ∆u,v + ∆v)x(u,v + ∆v)= ux (u,v + ∆v)u + o(u) y(u + ∆u,v + ∆v)y(u,v + ∆v)= uy (u,v + ∆v)u + o(u) x(u + ∆u,v + ∆v)x(u + ∆u,v)= xv (u,v + ∆v)v +o(v) y(u + ∆u,v + ∆v)y(u + ∆u,v)= yv (u,v + ∆v)v + o(v)

 

 

При малых u, v производные

x

,

x ,

y ,

y , вычисленные в точках

 

u

 

 

 

 

 

 

 

 

 

 

 

v

u

v

 

 

 

 

(u + ∆u,v),(u,v + ∆v),

мало

отличаются от

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

производных,

вычисленных в точке (u,v),

 

 

 

 

 

 

 

 

 

 

 

мало

отличаются

от

поэтому r u (u,v), r v (u,v)

 

 

и

 

 

 

 

 

 

 

 

 

 

 

 

соответственно,

и

r u (u + ∆u,v), r v (u + ∆u,v)

r u (u,v + ∆v), r v (u,v + ∆v),

рассматриваемый

четырёхугольник

представляет

собой

«почти

параллелограмм».

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

u ,

 

v

 

 

 

 

 

r

r

 

 

 

 

 

 

u

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

v

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

u

v = J

, т.е. равна

 

J

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u v

11

Поэтому

 

при

 

 

преобразовании

 

интегральная

сумма

n

 

 

 

 

 

 

 

 

 

 

 

 

f (x(ui ,vi ), y(ui ,vi

))

 

J (ui ,vi

)

 

uv

близка

к

интегральной

сумме

 

 

i=1

 

 

 

 

 

 

 

 

 

 

 

 

n

)пл.(Di

),

 

 

 

 

 

 

 

 

 

 

f (xi , yi

и т.к.

 

 

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

интегральные суммы для

i=1

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

Замечание. Утверждение теоремы сохранится, если условие взаимной

однозначности отображения x = x(u, v), y = y(u, v)

нарушится на

множестве

нулевой площади.

 

 

 

 

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

 

 

Пусть

требуется

вычислить

∫∫ f (x, y)dxdy

по области

D, которая

 

 

 

D

 

 

задаётся

в

полярных

координатах

условиями

α φ β

 

 

 

 

 

 

 

 

 

 

 

r r(φ)

Сделаем замену переменных

x = r cosφ

y = r sinφ

При этой замене нарушается взаимная однозначность отображения. Точке (0,0) соответствует целый отрезок [α, β] на оси φ . Однако точка и

отрезок имеют нулевую площадь, и теорема справедлива. Осталось

вычислить

 

 

J.

x

= cosφ,

y

= sinφ,

x

= −r sinφ,

y

= r cosφ.

 

r

r

φ

φ

 

 

 

 

 

 

 

 

 

 

J =

 

cosφ

r sinφ

 

= r cos2

φ + r sin 2 φ = r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

sinφ

r cosφ

 

 

 

 

 

 

 

 

 

 

Следовательно, ∫∫ f (x, y)dxdy = αβ dφ0r (φ ) f (r cosφ, r sinφ)rdr .

D

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

Рассмотрим пример. Найти I = 0ex2 dx .

12

Решение. I — это несобственный интеграл, и прежде всего

следует

установить его сходимость. По определению,

ex2 dx = 1 ex2 dx +

ex2 dx .

 

0

0

1

Первый из интегралов — собственный, второй — сходится по 1-й теореме о

сравнении,

так

 

как

 

 

при

 

x 1

 

 

справедливы

неравенства

x2 x x2

≤ −x,ex2 ex , а ex dx , очевидно, сходится.

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Обозначим

IR

= Rex2 dx

 

(очевидно,

IR I,R →∞). Тогда, поскольку

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

обозначение

переменной интегрирования можно

выбрать

произвольным,

IR2 = IR IR = Rex2 dx 0Rey2 dy = ∫∫ex2 y2 dxdy , где SR

— квадрат, а CR ,C 2R

 

 

0

 

 

 

 

 

 

SR

 

 

 

 

 

 

 

 

 

 

 

четверти круга, соответственно, радиусов R

и

 

2R . Так как ex2 y2 0 , то по

свойствам 2

и

3 двойного

 

интеграла

 

∫∫ex2 y2 dxdy IR2

∫∫ex2 y2 dxdy . В

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

CR

 

 

 

 

 

 

C 2R

интеграле ∫∫ex2 y2 dxdy перейдем к полярным координатам:

 

 

 

CR

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x2 y2

 

π/ 2

R

r2

 

πR

 

r2 1

2

 

π

 

r2

 

R

π

(1 e

R 2

).

 

 

 

 

 

 

∫∫e

 

dxdy = dφe

 

rdr =

e

 

 

dr =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

4

e

 

 

=

4

 

CR

 

 

0

0

 

 

2

0

 

 

 

 

 

 

0

 

 

 

Аналогично,

∫∫ex2 y2 dxdy =

π

(1 e2R 2 )и

π(1 eR 2 dx)IR2 π

(1 e2R 2 dx). При

 

 

 

C 2R

 

 

 

4

 

 

 

 

 

 

4

 

 

 

 

 

 

4

 

стремлении R к получаем, что limIR2 =

π

, то есть I2

=

π

,I =

π .

R→∞

4

 

 

4

 

2

 

 

 

 

 

 

4. Тройные интегралы

Рассмотрим кубируемую область в трехмерном пространстве V R3 .

Разбиение T на части Vi осуществляется непрерывными поверхностями.

Диаметр разбиения определяется аналогично двумерному случаю. Также, по аналогии, можно определить для функции f (x, y,z), разбиения T области V и

13

выбранных

точек

Mi (xi , yi ,zi ) Vi

интегральную

сумму

σ(f ,T ,{M i })= n

f (M i )µ(Vi

), где µ(Vi ) обозначает объем области Vi .

i=1

 

 

 

 

Определение. Пусть I R такое число, что ε > 0 δ > 0

T,d(T)< δ

{µi }σ(f ,t{µi })I <ε . Тогда мы говорим, что f интегрируема на V, число I

есть интеграл f по области V и обозначаем это так: I = ∫∫∫f (x, y,z)dxdydz .

V

Как и в случае двойного интеграла, выполняются аналогичные свойства 1-6. Можно доказать, что если f (x, y,z) непрерывна на V, то она интегрируема на V. Точно также можно убедиться в том, что если точки разрыва f лежат на конечном числе непрерывных поверхностей, и разбивающих V на кубируемые области, то f интегрируема на V.

Вычисление тройного интеграла производится по следующему правилу.

Теорема.

 

Пусть

V

задана

следующими

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

ψ1 (x, y)z ≤ ψ2 (x, y),(x, y) D . D — квадрируемая

замкнутая

область на

плоскости, ψ1 ,ψ2

— непрерывные функции на D. Тогда

 

 

 

 

 

ψ

(x, y )

 

 

 

 

 

∫∫∫ f (x, y, z)dxdydz = ∫∫dxdy 2

f (x, y, z)dz

 

 

 

 

V

 

D

ψ1 (x, y )

 

 

 

 

 

Замечание.

Если

область

D

задана

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

a x b,φ1 (x)y ≤ φ2 (x), где φ1 ,φ2

— непрерывные функции на [a,b], то

 

 

b

φ2 (x ) ψ2

(x, y )

 

 

 

 

 

∫∫∫ f (x, y, z)dxdydz = dx dy

 

f (x, y, z)dz.

 

 

 

 

V

 

a

φ1 (x )

ψ1 (x, y )

 

 

 

 

 

Сформулируем общую теорему о замене переменных.

 

 

Теорема.

Пусть отображение x = x(u,v, w),

y = y(u,v, w),

z = z(u,v, w)

устанавливает

взаимно

однозначное

соответствие между

областями

,(u,v, w) и

 

V ,(x, y, z) V ,

причем

функции

x, y, z

непрерывно

14

 

x

x

x

 

 

 

 

 

 

u

v

w

 

 

дифференцируемые и J =

y

y

y

 

0 ни в одной точке . Пусть

u

v

w

 

 

 

 

z

z

z

 

 

 

u

v

w

 

 

f (x, y, z) — непрерывная функция. Тогда

∫∫∫ f (x, y, z)dxdydz = ∫∫∫ f (x(u,,v, w), y(u,v, w), z(u,v, w)) J dudvdw

 

 

V

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Как и для двойного интеграла, теорема остается верной в случае

нарушения ее условий на множестве нулевого объема.

 

 

 

 

 

 

 

Пример 1. Переход к цилиндрическим координатам. Он

осуществляется с помощью функций: x = r cosφ, y = r sin φ,z = z .

 

 

 

 

При

 

 

 

 

этом

 

 

 

 

 

якобиан

 

равен

 

cosφ

 

r sin φ

0

 

 

cosφ

r sin φ

 

= r cos2 φ+ r sin2 φ= r .

 

 

 

 

 

 

 

 

 

 

 

 

 

J =

sin φ

 

r cosφ

0

 

=

 

 

 

 

0

0

 

1

 

 

sin φ

r cosφ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример 2. Переход к сферическим координатам осуществляется

функциями x = r cosφcosψ, y = rsin φcosψ,z = rsin ψ.

 

 

 

 

 

 

 

Якобиан

 

 

 

 

 

 

 

 

 

преобразования

 

 

 

 

 

равен

 

 

cosφcosψ

r sinφcosψ

r cosφsinψ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

J =

sinφcosψ

r cosφcosψ

r sinφsinψ

 

=

 

 

 

 

 

 

 

 

sinψ

 

 

0

 

 

 

 

 

r cosψ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(разложение

 

 

 

 

 

 

 

 

по

 

 

 

 

 

 

3-й

 

 

 

 

строке)

= r cosψ

 

cosφcosψ

 

r sinφcosψ

 

+sinψ

 

r sinφcosψ

r cosφsinψ

 

=

 

 

 

 

 

 

 

 

 

sinφcosψ

 

 

r cosφcosψ

 

 

 

r cosφcosψ

r sinφsinψ

 

 

 

 

(выделим общие множители у столбцов)

 

 

 

 

 

 

 

 

r2 cos3 ψ

 

cosφ

 

 

sin φ

 

+ r2 sin2 ψcosψ

 

sin φ

cosφ

 

= r2 cos3 ψ + r2 sin2 ψcosψ =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

sin φ

 

 

cosφ

 

 

 

 

 

 

 

 

 

 

 

 

cosφ

sin φ

 

 

 

 

= r2 cosψ(cos2 ψ + sin2 ψ)= r2 cosψ.

15

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

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

Рассмотрим спрямляемую (т.е. имеющую длину) кривую AB на плоскости (A, B – точки плоскости). Для простоты, считаем, что эта кривая задана параметрически x = x(t), y = y(t),t [T0 ;T1 ], причем x(t), y(t) – непрерывно дифференцируемые на отрезке функции такие, что каждому значению параметра соответствует единственная точка кривой.

Тогда длина кривой выражается формулой l = T1

(x(t))2

+(y(t))2 dt .

 

 

 

 

T0

 

 

Под

разбиением

T кривой AB будем понимать множество точек

A = A0 , A1 ,..., An = B ,

лежащих на этой кривой и занумерованных в

направлении от A к B. Пусть li - длина кривой Ai1 , Ai .

 

Диаметр d(T) определим как d(T )= max li .

 

 

 

 

 

 

i=1,...,n

 

 

Пусть функция

f (M )= f (x, y) определена на кривой AB. Выберем на

каждом

участке

Ai1 Ai

кривой точку M i

и

образуем сумму

σ(f ,T ,{M i })= n

f (M i )li

, называемую интегральной.

 

 

i1

 

 

 

 

 

Определение.

Пусть

I R . Если ε > 0 δ > 0 T : d(T )<δ {M i }

σ(f ,T,{M i })I < ε , то величина I называется криволинейным интегралом

первого типа по кривой AB и обозначается так: I = f (x, y)dl .

AB

Важное замечание. Если бы мы совершали движение по кривой не от A к

B,

16

а от B к A, то в разбиении T с

выбранными

точками

{M i }

изменилась

бы только

нумерация

отрезков и

точек

M i ,

а сама

интегральная сумма

не

изменилась

бы, поскольку в ее определении

фигурирует лишь длина li участка,

которая не зависит от того, в каком

направлении

проходится

участок.

Это

означает,

 

что

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

 

 

AB

BA

 

 

 

В этом важнейшее отличие от обычного определенного интеграла, который менял бы знак при изменении направления обхода кривой.

Сформулируем теорему, сводящую новый пока объект - криволинейный

интеграл к обычному определенному интегралу.

 

Теорема.

 

Пусть f (M ) - непрерывная

на

кривой AB функция (т.е.

ε > 0 δ > 0 M1 , M 2 - точек кривой таких,

что расстояние между M1 , M 2

меньше δ,

 

f (M1 )f (M 2 )

 

< ε ). Пусть кривая

AB параметризована так:

 

 

 

 

 

 

- непрерывные на [T0 ;T1 ]

x = x(t), y = y(t),t [T0 ;T1 ], где x(t), y(t), x (t), y (t)

функции, причем каждому значению параметра соответствует единственная

 

 

 

 

T1

 

2

2

dt .

 

 

 

 

 

 

 

 

точка кривой. Тогда f (M )ds = f (x(t), y(t)) (x (t))

 

+ (y (t))

 

 

 

 

 

AB

T0

 

 

 

 

 

 

 

 

Доказательство.

 

Интегральная

 

 

сумма

n

f (M i )li

= n

f (M i )ti

(x(t))2 +(y(t))2 dt для

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

интеграла

i=1

 

i=1

ti 1

 

 

 

 

 

 

 

 

первого

типа

отличается

от

интегральной

суммы

f (x(ti ), y(ti ))

(x(ti ))2 + (y(ti ))2 ti

 

для

 

 

интеграла

17

T

 

 

 

t

 

 

1

f (x(t), y(t))

(x(t))2 + (y(t))2 dt лишь тем, что

i

(x(t))2

+ (y(t))2 dt несколько

T0

 

 

 

ti 1

 

 

отличается

от

(x(ti ))2 + (y(ti ))2 ti , именно

,

этот

интеграл равен

(x(ξi ))2 + (y(ξi ))2 ti .

Нетрудно доказать, что при d(T )0 пределы этих сумм равны (строгое

доказательство опущено). Это означает, что теорема доказана.

Отметим, что изменение направления обхода кривой означает одновременную смену пределов интегрирования и знака величины dt, что не изменяет величину интеграла в правой части этого равенства.

Из свойств криволинейного интеграла отметим следующие 2 основных:

1.

(λ1 f1 (M )+ λ2 f2 (M ))dl = λ1 f1 (M )dl + λ2

f2 (M )dl

при условии,

что

 

AB

 

AB

AB

 

 

 

существуют f1 (M )dl и f2 (M )dl .

 

 

 

 

 

AB

AB

 

 

 

2.

Если

AB, BC -

кривые, удовлетворяющие

условиям теоремы,

то

 

f (M )dl = f (M )dl + f (M )dl .

 

 

 

 

AC

AB

BC

 

 

 

Свойство 2 позволяет определить криволинейные интегралы 1-го типа для кусочно-гладких кривых (т.е. кривых, состоящих из конечного числа частей, каждая из которых удовлетворяет условиям теоремы). В частности, можно определить криволинейный интеграл и для замкнутых кривых.

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

Рассмотрим, как и в параграфе 1, кривую AB, которую пока считаем незамкнутой.

18

Пусть проекция этой кривой на ось x представляет собой отрезок [a; b].

Пусть точки A0 ,..., An дают разбиение кривой AB. Рассмотрим их проекции x0 ,..., xn , лежащие на отрезке [a; b] и обозначим xi = xi xi1 , i =1,..., n .

Пусть P(x, y) - определена на AB. Пусть M i -

точка, лежащая на кривой

 

 

 

 

 

n

)xi .

 

 

 

между Ai1 и Ai . Положим σ = P(M i

 

 

 

 

 

 

 

 

i=1

 

 

 

 

Определение. Пусть

I R . Если

ε > 0

δ > 0 T : d(T )<δ

и {Mi }

выполняется

 

σ I

 

<ε ,

то говорят,

что I

- это

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

интеграл

 

 

второго типа P(x, y)dx .

AB

Точно также, рассматривая проекции на ось y, определим Q(x, y)dy .

AB

Интеграл общего вида P(x, y)dx +Q(x, y)dy определяется, как сумма этих

AB

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

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

соответствии со следующей теоремой.

 

 

 

Теорема. Пусть L – кривая, заданная уравнениями:

x = x(t)

, t [T ,T ],

 

y = y(t)

0

1

 

 

P(x, y) и Q(x, y) – непрерывные на L функции. Тогда:

Pdx + Qdy = T1 (P(x(t), y(t))x(t)+ Q(x(t), y(t))y(t))dt .

L

T0

Теорема без доказательства.

Примечание 1.

19

a) Если кривая L задана явным уравнением

y =φ(x), a x b ,

где φ(x) -

непрерывно дифференцируемая функция,

то предыдущая

формула

принимает вид: Pdx +Qdy = b P(x,φ(x))dx + b Q(x,φ(x))φ(x)dx .

L a a

b)

Если L

задана уравнением x =ψ(y), c y d , то Pdx +Qdy =

 

 

 

L

 

d

d

 

 

 

 

= P(ψ(y), y)ψ (y)dy + Q(ψ(y), y)dy .

 

c

 

c

c)

Если L -

отрезок прямой x = x0 , то Pdx 0 для любой функции P,

 

 

 

L

если L - отрезок прямой y = y0 , то Qdy 0 для любой функции Q.

L

Примечание 2.

Пусть

α - угол, составляемый

вектором касательной к кривой и

положительным направлением оси x.

Тогда dx = ds cosα, dy = ds sin α . Поэтому

Pdx +Qdy = (P cosα +Q sin α)ds .

 

L

L

 

Заметим, что при изменении направления обхода угол α изменяется на π +α . При этом cos(π +α)= −cosα, sin(π +α)= −sin α , и интеграл в правой части написанного выше равенства меняет свой знак.

20

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]