Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Математика / Шилкин Лекции.PDF
Скачиваний:
33
Добавлен:
11.05.2015
Размер:
7.25 Mб
Скачать

 

ЛЕКЦИЯ 32

 

 

 

x + y + z = 2

.

Пример 1. Вычислить ∫∫∫ xdxdydz, T :

 

T

z = 1, x = 0, y = 0

 

Решение. Область T представляет собой пирамиду, поднятую на высоту z = 1 :

 

 

 

 

1 1x

 

2xy

 

1

 

1x

 

 

 

 

 

 

 

 

 

 

1

 

1

 

 

3

 

 

 

 

 

 

2

 

 

1

 

 

 

1

.

∫∫∫ xdxdydz = dx

dy

 

xdz = xdx

(1 x y)dy =

 

 

 

x

 

x

 

+

 

 

 

x dx =

 

2

 

 

 

2

 

24

T

 

 

 

0

0

 

 

1

 

 

 

0

 

0

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

1

, x y 2 x

 

 

 

 

 

 

Пример 2. Вычислить ∫∫∫zdxdydz, T :

0

2

.

 

 

 

 

 

 

 

 

 

 

1

x

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0 z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

2

 

 

 

 

 

 

1

Решение.

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 x

1x2 y2

 

 

 

 

1

 

2 x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

2 x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

I = 2dx dy

zdz =

2dx (1 x2

y2 )dy =

 

1

 

 

2( y yx2

 

 

 

y3 )

 

 

dx =

 

 

2

2

 

 

 

3

 

 

x

 

 

0

 

x

0

 

 

 

0

x

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

12

 

 

 

8

 

 

 

 

 

 

1

 

 

 

 

 

1

1

 

 

 

 

5

 

 

 

 

 

 

 

 

 

7

 

 

 

 

 

 

 

3

 

 

 

3

 

3

 

 

3

 

 

 

 

2

 

 

 

 

4

 

 

2

 

 

 

 

.

 

 

 

=

 

 

(2x 2 x

 

 

 

 

x

 

x + x

 

+

 

x

 

)dx =

 

 

 

 

 

 

 

x

 

 

x

 

 

 

 

 

=

 

 

 

 

 

 

 

 

2

 

3

 

 

 

3

 

2

 

2

 

6

 

 

 

192

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Двойной и тройной интегралы в полярной и цилиндрической системе координат

Пусть двойной интеграл преобразуется от прямоугольных координат x, y к криволинейным координатам u, v , связанным с прямоугольными соотноше-

ниями

x = x(u, v), y = y(u, v) , где функции

x(u, v), y(u, v) имеют непрерывные ча-

стные производные в области D

плоскости uO v

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

области D

не обращается в нуль

Ι =

u

 

v

 

0 . При этом устанавливается вза-

 

y

 

y

 

 

 

 

 

 

 

 

 

 

 

 

u

 

 

v

 

 

 

имнооднозначное соответствие между точками области D плоскости XOY и D

плоскости

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

uO v .

 

 

 

 

 

 

 

 

 

 

 

 

 

Формула преобразования двойного интеграла в этом случае имеет вид

∫∫ f (dx, y)dxdy = ∫∫ f ( x(u, v), y(u, v)) | I | dudv .

 

 

 

 

 

D

 

 

 

 

D

 

 

 

 

 

 

 

 

 

 

 

 

 

В

 

 

случае

полярных

 

 

координат

x = r cosϕ, y = r sin ϕ

 

x

x

 

 

cosϕ

r sinϕ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ι =

r

ϕ

 

=

= r cos

2

ϕ + r sin

2

ϕ = r и тогда преобразование двойно-

y

y

 

sinϕ

r cosϕ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

ϕ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

прямоугольных к полярным

осуществляется по формуле

∫∫ f (dx, y)dxdy = ∫∫ f (r cosϕ, r sinϕ)rdrdϕ = ∫∫F(r,ϕ)rdrdϕ .

D

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Если область ограничена двумя лучами ϕ = α, ϕ = β (α < β ) , выходящи-

ми из полюса и двумя кривыми r = r1 (ϕ), r = r2 (ϕ) , где r1 (ϕ), r2 (ϕ) однозначные функции при α ϕ β , то двойной интеграл сводится к повторному по сле-

β r2 (ϕ )

дующей формуле: ∫∫F(r,ϕ)rdrdϕ = dϕ F (r,ϕ)rdr .

α r1 (ϕ )

Пример.

 

Найти

площадь

 

фигуры,

 

 

ограниченной

линиями

x2 6 x + y2 = 0, x2 10x + y2 = 0, y = x

3

, y = 3 x .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Решение. Вычислим в полярной системе координат. Уравнения линий,

ограничивающие область D , следующие: tgϕ

= 1

,

tgϕ =

3 ϕ1

=

π ,

ϕ2 = π ,

 

 

 

 

 

 

 

 

3

 

 

 

 

 

6

3

 

 

 

 

 

 

 

 

 

 

 

 

 

π

 

 

 

 

 

 

 

 

 

S = ∫∫dxdy = ∫∫rdrdϕ = 3

 

10 cosϕ

r1 = 5cosϕ, r2

= 10cosϕ . Таким образом,

dϕ

rdr =

 

 

 

 

 

 

 

 

D

 

 

π

 

6 cosϕ

 

 

 

 

 

 

 

 

 

 

 

 

 

6

 

 

π 3

2

 

π 3

 

 

 

1

 

 

 

π 3

π 3 .

 

 

 

 

 

 

 

 

 

 

 

= 32 cos

ϕdϕ = 16 (1 + cos 2ϕ)dϕ =

16 ϕ +

 

sin 2ϕ

 

= 8

 

 

 

2

 

 

 

π 6

 

 

π 6

 

 

 

 

 

 

π 6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Вообще полярную систему координат удобно применять там, где x, y в уравнениях границ области D присутствуют в четных степенях.

Отметим, что модуль якобиана равен коэффициенту изменения бесконечно малой площади при переходе от плоскости XOY к XOY , то есть

dxdy = rdϕdr .

Тройной интеграл в цилиндрической системе координат

На плоскости XOY проекции точки M ставятся в соответствие полярные координаты r,ϕ . Связь цилиндрической и декартовой систем координат сле-

x = r cosϕ, 0 ϕ 2π

дующая:

y = r sinϕ, 0 r ≤ ∞ . (Если r = const , а ϕ, z произвольные, то получаем

 

z = z, − ∞ < z < ∞

 

 

 

 

 

цилиндр радиуса r .)

 

 

 

 

 

 

Для

 

 

цилиндрической

системы

координат

 

( x, y, z)

 

cosϕ

r sinϕ

0

 

 

неотрицателен ( r равен 0

только на

 

 

 

I(r,ϕ, z) =

=

sinϕ

r cosϕ

0

 

= r

 

(r,ϕ, z)

 

0

0

1

 

 

 

 

 

 

 

 

 

 

 

границе), поэтому ∫∫∫ f ( x, y, z)dxdydz =∫∫∫ f (r cosϕ, r sinϕ, z)rdϕdtdz .

 

 

 

 

V

 

 

 

T

 

 

Существуют и другие системы координат в пространстве, например, сфе-

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

Пример 1.

Вычислить объем тела V , ограниченного поверхностями

z = x2 + y2 , z = 1 в

цилиндрической системе координат.

Решение. r

2 = z, z = 1, ϕ = 0, ϕ = 2π V = ∫∫∫dxdydz = ∫∫∫rdϕdrdz =

V T

2π

1

1

1

 

 

r

2

 

3

4

 

 

1

 

π

.

 

 

 

 

= dϕrdr dz = 2π r(1 r

2

 

 

 

 

 

 

=

 

 

 

 

 

 

 

)dr = 2π

2

4

 

 

 

2

0

0

r 2

0

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример

2.

 

Вычислить

 

массу тела

с

плотностью

µ = 5z

V : { x2 + y2

+ z 2

= 16, x2

+ y2

= 9z 2 , x = y = 0, x 0, y 0, z 0} .

 

 

 

 

 

Решение.

 

M = ∫∫∫µ( x, y, z)dxdydz .

Проекция

тела V

на плоскость

 

 

 

 

 

 

 

 

 

 

V

 

 

 

 

 

 

 

 

 

 

 

 

 

XOY дает четверть окружности

x2 + y2

=

 

144

, x 0, y 0 . Перейдем к цилиндри-

10

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ческой

системе координат:

0 ϕ

π , 0 r 12

, 1 r z

16 r 2 .

Тогда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

10

3

 

 

 

π 2

12

10

16r 2

 

 

 

12

 

 

 

 

 

2

 

 

 

 

 

 

 

 

5π

10

 

 

 

r

 

 

 

 

 

M = dϕ

rdr 5zdz =

 

 

 

16 r

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

dr = 72π .

 

 

 

 

0

0

 

 

1

r

 

4

0

 

 

 

 

9

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ЛЕКЦИЯ 33

Дифференциальные уравнения

При решении многих задач математики физики, экономики, естествознания часто приходится отыскивать неизвестную функцию по данному соотношению между ней, ее производными и независимыми переменными.

Так, например, в процессе распада радиоактивного вещества установлено, что скорость распада пропорциональна с коэффициентом пропорциональности к наличному количеству нераспавшегося радиоактивного вещества. Обозначая через x0 массу радиоактивного вещества в начальный момент времени t = 0 , а

через x его массу в момент времени t , получаем dxdt = −kx, k > 0 . Знак “-” берется в уравнении потому, что масса убывает. Если это уравнение переписать в

виде

dx

= −kdt

,

то можно проинтегрировать обе части, получив соотношение

x

 

 

 

 

ln x = −kt + ln c .

 

Учитывая, что при t = 0 x = x0 , получаем, что c = x0 и тогда

x = x0 ekt . Эта функция является решением дифференциального уравнения, и

мы можем определить количество радиоактивного вещества в любой момент времени t . Значение этого позволяет производить датировку событий, имевших место тысячи и миллиарды лет тому назад.

К такому же уравнению приходим при определении численности населения данного региона в зависимости от времени t . Необходимо только знать коэффициенты рождаемости и смертности для данного региона.

Обыкновенным дифференциальным уравнением n-го порядка назы-

вается соотношение вида F( x, y, y,K, y( n) ) = 0 , где F - функция, определенная в некоторой области D Rn+2 , x - независимая переменная, y - искомая функция переменной x ( y = y( x)) , а y, y′′,K, y( n) - ее производные.

Порядком ДУ называется порядок старшей входящей в нее производной.

Если уравнение

F( x, y, y,K, y( n) ) = 0 можно записать в виде

y( n) = f ( x, y, y,L, y( n1) ) , где

f определена в некоторой области D Rn+1 , то го-

ворят, что ДУ разрешено относительно старшей производной y( n) и называется

ДУ в нормальной форме.

ДУ, в котором неизвестная функция y зависит от одной переменной x ,

называется обыкновенным.

Если же ДУ содержит неизвестную функцию многих переменных и ее частные производные, то оно называется уравнением в частных производных.

Мы ограничимся изучением ОДУ и слово “обыкновенные” далее будем опускать.

Так, например, уже рассмотренное нами уравнение радиоактивного распада является ОДУ 1-го порядка в нормальной форме.

Уравнение y′ − (2 xy)2 ln y′ = 0 - ОДУ 1-го порядка;

Уравнение y′′′ = 1 y x - ОДУ 3-го порядка в нормальной форме:

Уравнение

2 u

+

2 u

= 0 - уравнение 2-го порядка в частных производных.

 

x2

 

y2

 

Решением (или интегралом) ДУ F( x, y, y,K, y( n) ) = 0 называется всякая действительная функция y = y( x) , определенная на интервале (a, b) такая, что

1)

y( x) n раз непрерывно

дифференцируема

на (a, b) ,

то есть

y( x) C ( n+1) (a, b)

 

 

 

 

2)

Точка ( x, y( x), y( x),K, y( n) ( x)) D Rn+2

для всех

x (a, b) , где

D - об-

ласть определения функции F

 

 

 

 

3) F( x, y( x), y( x),K, y( n) ( x)) 0 x (a, b) .

 

 

 

Всякому решению ДУ F( x, y, y,K, y( n) ) = 0

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

рая кривая y = y( x) x (a, b) , которая

называется интегральной кривой ДУ

F( x, y, y,K, y( n) ) = 0 .

Процесс нахождения решения ДУ называется интегрированием этого

уравнения.

 

 

 

 

Пример

1. Решением уравнения y′′′ = 1 ( y′′)2

является функция

y = −sin x + 2 x + 1 , так как y′ = −cos x + 2, y′′ = sin x, y′′′ = cos x,

1 sin 2 x = cos x,

 

π

,

π

 

x

2

.

 

 

 

2

 

Пример 2. Решением уравнения xyy′ + x2 + 2 xy = 0 являются функ-

ции ln | x + y | + x +x y = c .

Упр. Проверить это.

ДУ, как правило, имеет бесконечное множество решений. Чтобы выделить одно какое-то решение ДУ необходимо задание некоторых дополнительных условий. Так, в первоначальном примере таким условием явилось задание первоначального количества вещества.

ДУ 1-го порядка. Задача Коши.

Рассмотрим ДУ 1-го порядка. Оно получается из F( x, y, y,K, y( n) ) = 0 при

n = 1 : F ( x, y, y) = 0

(3)

или

y′ = f ( x, y)

(3a)

где

f ( x, y) - известная функция ( x, y) D R2 .

Уравнение (3а) всегда можно записать в так называемой дифференциаль-

ной форме:

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

 

(4) где P( x, y), Q( x, y)

заданные функции в плоскости области D . Легко получить из (4) (3а) и наобо-

рот. Так, если

y′ =

dy

подставить в (3а), то получится

f ( x, y)dx 1 dy = 0

аналог

dx

(4), а если в (4) разделить обе части на Q( x, y)dx , получим

y

 

P( x, y)

 

уравне-

= − Q( x, y)

 

ние типа (3а).

ления функции

Далее будем рассматривать уравнения типа (3а) или (4).

Решением (или интегралом) уравнения (3) называется действительная функция y = y( x) , определенная на интервале (a, b) такая, что

1)y( x) непрерывно дифференцируема на (a, b)

2)Точка ( x, y( x), y( x)) D R3 для всех x (a, b) , где D - область опреде-

F

3) F ( x, y( x), y( x)) 0 x (a, b) .

Упр. Сформулировать понятие решения для (3а).

Уравнение y′ = f ( x, y) имеет бесчисленное множество решений (почему?). Чтобы из множества этих решений выделить одно, так называемое частное решение, надо задать некоторые дополнительные условия. Таким условием

является условие Коши y( x0 ) = y0 , ( x0 , y0 ) D

(5)

Задача отыскания частного решения уравнения (3а), удовлетворяющего начальному условию (5), называется задачей Коши для этого уравнения

С геометрической точки зрения задача Коши для ДУ (3а) означает следующее: из множества интегральных кривых выделить ту, которая проходит через заданную точку M0 = ( x0 , y0 ) D .

Теорема 1. Пусть функция f ( x, y) определена, непрерывна и имеет непрерывную частную производную fy в области D . Тогда найдется интервал

( x0 δ , x0 + δ ) , на котором существует единственное решение y = y( x) ДУ (3а), удовлетворяющего условию y( x0 ) = y0 .

Единственность здесь понимается в том смысле, что если есть 2 решения y = y1 ( x), y = y2 ( x) , то они совпадают в окрестности точки x0 . При нарушении

условий

теоремы

через

точку ( x0 , y0 ) D могут проходить несколько инте-

гральных кривых.

 

f

1

 

 

Пример. y

 

 

 

= 2 y .

Функции f ( x, y) = 2 y , y =

y

определены и непре-

 

рывны для y > 0 ,

то есть условия теоремы Коши выполнены. Пусть

y(1) = 0 .

Поверим, что функция

y = ( x 1)2 является решением задачи Коши.

Действи-

тельно,

y(1) = 0, y′ = 2( x 1), 2( x 1) = 2 ( x 1)2 = 2( x 1)

( x 1 > 0) , но если рас-

смотреть это уравнение в области y 0 (а именно так определена правая часть

уравнения), то добавляется еще решение y = 0 и к

семейству парабол

y = ( x c)2 , являющимся общим решением уравнения y′ = 2

y , добавляется еще

прямая

y = 0 , то есть через каждую точку оси OX проходит 2

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

кривые. Здесь все дело в том, что f =

1 → ∞ .

 

 

y

y y0

 

Общим решением ДУ (3а)в некоторой области существования и единст-

венности решения задачи Коши называется семейство функций

y = ϕ( x, c) , за-

висящее от параметра c и удовлетворяющее условиям:

 

1.

Функции y = ϕ( x, c) являются решением уравнения (3а) при лю-

бом допустимом значении параметра c

2. Рои любом начальном условии y( x0 ) = y0 , ( x0 , y0 ) D c0 такое, что функция y = ϕ( x, c0 ) удовлетворяет условию ϕ( x0 , c0 ) = y0 .

Решения, получающиеся из общего при определенном значении параметра c , называются частными.

Для выделения из общего решения частного достаточно использовать начальные условия.

Часто общее решение ДУ y′ = f ( x, y) находится в виде Φ( x, y, c) = 0 и на-

зывается общим интегралом ДУ.

При фиксированном значении c = c0 общий интеграл Φ( x, y, c) = 0 превра-

щается в частный Φ( x, y, c0 ) = 0 .

При интегрировании ДУ действие выделения неопределенного интеграла называют “квадратурой”, чтобы не было путаницы с термином “интегрирование уравнения”.

Из уравнения y′ = f ( x, y) следует, что угловой коэффициент yкасательной к интегральной кривой в каждой ее точке ( x, y) D равен значению функции f ( x, y) в этой точке. Таким образом, в каждой точке ( x, y) D можно указать направление касательной к интегральной кривой, проходящей через точку ( x, y) . Поэтому, если через каждую точку ( x, y) провести отрезок с угловым коэффициентом k = f ( x, y) , то получим так называемое поле направлений в области D . То есть ДУ (3а) определяет в области D поле направлений, а решению ДУ соответствует кривая, направление касательной которой в каждой ее точке совпадает с направлением поля в этой точке.

Найти решение приближенно можно с помощью метода изоклин. Изоклиной ДУ (3а) называется кривая, в каждой точке которой поле на-

правлений имеет один и тот же наклон.

 

Таким образом,

семейство

изоклин определяется равенством

f ( x, y) = k = tgα , где k -

параметр, α -

угол наклона поля направлений к оси

OX . Придавая параметру k близкие числовые значения можно получить сеть изоклин, с помощью которых приближенно строятся интегральные кривые ДУ

y′ = f ( x, y) .

Пример. y′ = x + y. f ( x, y) = x + y определена во всей плоскости, следовательно, поле направлений может быть построено во всей плоскости. Обозначим через α угол наклона к оси OX направление поля изоклина, в точках которой

α =0,

то есть

y′ = tgα = 0 имеет уравнение x + y = 0 или y = −x . Изоклина, в ко-

торой

α = π ,

имеет уравнение x + y =

3 и чтобы по данному полю направле-

 

6

 

3

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

надо на плоскости XOY взять любую

точку ( x0 , y0 )

за начальную и вести через нее линию так, чтобы она в каждой

своей точке шла в направлении поля.

Практически с помощью метода изоклин можно приближенно построить семейство интегральных кривых уравнения y′ = f ( x, y) . При этом чем больше

построено изоклин, тем точнее можно провести интегральные кривые, представить взаимное расположение интегральных кривых ДУ.

Упр. Методом изоклин приближенно построить интегральные кривые ДУ

y′ = xx +yy .

ДУ 1-го порядка, интегрируемые в квадратурах

1. Уравнение вида y′ = f ( x) , где f ( x) непрерывная на (a, b) функция. Очевидно, все решения этого простейшего ДУ исчерпываются соотноше-

нием y = f ( x)dx + c , где c - произвольная константа. (Здесь f ( x)dx + c обозна-

чает одну из первообразных, а не совокупность их. Далее в теории ДУ это подразумевается.)

Геометрически: все интегральные кривые уравнения 1. в полосе

{a < x < b, − ∞ < y < +∞}

получаются из одной из них, например,

y = f ( x)dx

сдвигом параллельно оси OY . Задавая в этой полосе любую точку

M0 ( x0 , y0 )

можно единственным

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

y = f ( x)dx + c0 проходит через эту точку.

Пример. y′ = 3 x2 . f ( x) = 3 x2 непрерывна для − ∞ < x < +∞ . Общее решение во всей плоскости XOY имеет вид y = x3 + c . Возьмем начальное условие x0 = 1, y0 = 3 . Тогда 3 = 1 + c0 c0 = 2 , то есть y = x3 + 2 .