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

высшая математика

.pdf
Скачиваний:
73
Добавлен:
18.02.2016
Размер:
5.88 Mб
Скачать

 

 

 

 

 

∫∫ f (x, y)dx dy =

∫∫ f (ρ

cos θ ,ρ sin θ )ρ dρ dθ .

 

(5)

 

 

 

 

 

D

D

 

 

 

 

 

 

Если область интегрирования D ограничена двумя лучами θ

= α ,

θ

= β

(α

< β ),

выходящими из полюса, и двумя кривыми ρ = ρ

1(θ ) и

ρ

= ρ

2(θ

),

где ρ 1(θ ) и ρ 2(θ )

– однозначные функции при α ≤ θ ≤ β

 

и

ρ

1 (θ) ≤ ρ

2 (θ) ,

то двойной интеграл справа в формуле (5) вычисляется

так:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

â

ρ

2 (θ )

 

 

 

 

 

 

 

∫∫ F(ρ , θ )ρ dρ dθ = dθ

F(ρ , θ )ρ dρ ,

 

 

 

 

 

 

 

D

á

ρ

1 (θ )

 

 

где F(ρ, θ

) = f (ρ cos θ , ρ sin θ ), причем

сначала вычисляется интеграл

ρ2 (θ)

F(ρ , θ )ρ dρ , в котором θ считается постоянным.

ρ1 (θ)

Пример 3. Перейдя к полярным координатам, вычислить двой-

ной интеграл

∫∫

x 2 + y 2 dx dy ,

если

 

D

первая

четверть круга

x2 + y2 a2.

 

D

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Решение.

Полагая x = ρ

cos θ ,

y = ρ

 

sin θ

по формуле (5), имеем:

∫∫

x 2 + y 2 dx dy = ∫∫

ρ 2 cos2 θ +

 

ρ 2 sin 2 θ ρ dρ

dθ =

D

 

 

 

 

 

D

 

 

 

 

 

 

 

 

 

 

 

 

 

π 2

a

ρ 2 dρ = 1

π

2

 

0 = a

3

θ 0π 2

= π a

3

 

 

=

dθ

 

dθ ρ 3

 

 

.

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

0

0

3 0

 

3

 

 

 

6

 

 

 

 

 

 

 

 

 

 

 

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

z = 1 – x2 y2,

y = x,

y = x 3 ,

z = 0 и расположенного в первом октанте.

Решение. Данное тело ограничено сверху параболоидом z = 1 – x2 y2. Область интегрирования D – круговой сектор, ограниченный дугой окружности x2 + y2 = 1, являющейся линией пересечения

параболоида с плоскостью z = 0, и прямыми y = x, y =

x

3 . Следова-

тельно, V = ∫∫ (1x 2 y 2 ) dx dy .

 

 

D

 

 

Перейдем к полярным координатам x = ρ cos θ ,

y = ρ

sin θ . Урав-

нение окружности x2 + y2 = 1 в этих координатах примет вид ρ = 1, подынтегральная функция равна 1 – ρ 2, а пределы интегрирования по θ определяем из уравнений прямых:

2 9 1

 

 

tg θ 1 = 1, т.е. θ 1 =

π

; tg θ 2 =

3 , т.е. θ 2 =

π

.

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

π

 

 

 

 

 

 

 

π

 

 

 

 

 

 

 

 

 

 

1

 

 

π

 

 

 

 

2

 

 

 

 

 

 

3

 

1

 

 

3

 

3

1

 

2

 

1

 

4

 

1

3

 

π

 

V = ∫∫ (1− ρ

 

)ρ dρ dθ = dθ

(ρ − ρ

 

) dρ =

 

 

ρ

 

 

ρ

 

 

 

 

dθ =

 

 

dθ =

 

.

 

 

2

 

4

 

 

 

4

48

D

 

 

 

 

 

 

 

 

π

 

0

 

 

 

 

 

π

 

 

 

 

 

 

 

 

0

 

π

 

 

 

 

 

 

 

π

 

4

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

Ответ: V =

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

48

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

! Задания для самостоятельной работы

 

 

 

 

 

 

1. Вычислить

∫∫ (3x 2

2xy +

y) dx dy ,

если область D ограниче-

 

 

 

 

 

 

 

 

D

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

на линиями x = 0, x = y2,

y = 2.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2. Вычислить

∫∫ y ln x dx dy ,

если область D ограничена линия-

ми xy = 1, y =

 

 

 

 

 

D

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x , x = 2.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3. Вычислить

∫∫ sin(x +

y) dx dy , если область D ограничена ли-

 

 

 

 

 

π

 

 

D

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ниями x = 0,

y =

,

y = x.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4. Вычислить ∫∫ ln(x 2 +

y 2 ) dx dy , если D – кольцо между окруж-

 

 

 

 

 

 

 

 

D

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ностями x2 + y2 = e2,

x2 + y2 = e4.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5. Вычислить

∫∫ x 2 +

y 2 dx dy , если область D ограничена ли-

 

 

 

 

 

 

 

 

D

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ниями x2 + y2 = a2,

x2 + y2 = 4a2.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6. Найти площадь фигуры,

ограниченной линией x3 + y3 = axy.

7. Вычислить

 

площадь

 

 

фигуры,

 

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

линией

x = y2 – 2y, x + y = 0.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

8. Найти

 

 

объем

тела,

 

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

 

 

 

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

x2 + y2 = a2, x2 + z2 = a2.

2 9 2

Лекция 47

Общее дифференциальное уравнение первого порядка. Составление дифференциальных уравнений

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

10. Общее дифференциальное уравнение (ДУ) пер-

вого порядка. ДУ первого порядка называется уравнение вида

F(x, y, y) = 0

или разрешенное относительно y:

y′ = f (x, y) ,

(1)

где х – аргумент, у (х) – искомая функция; F – заданная функция трех переменных x, y, y; f – заданная функция двух переменных х, у.

Уравнение

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

где х и у имеют тот же смысл, что и выше; dx, dy – дифференциалы; P(x, y) и Q(x, y) – заданные функции, называется уравнением, записанным в дифференциалах. Здесь у можно принять за аргумент, х(у) – за функцию; dy, dx = x (y) dy – соответственно за дифференциалы.

Решением ДУ называется такая дифференцируемая функция y = y(x), которая при подстановке в уравнение обращает его в тожде-

ство, например, F(x, y(x ), y(x )) 0 или y(x ) f (x, y(x )) .

Процесс нахождения всех решений ДУ называется интегрированием ДУ, а график решения у = у(х) ДУ называется интегральной кри-

вой этого уравнения.

 

 

 

 

 

Пример 1. Проинтегрировать простейшее

ДУ:

y′ =

f (x ) .

Имеем:

 

dy

= f (x )

dy =

f (x ) dx dy = f (x ) dx

y =

F (x ) + C ,

 

 

 

 

dx

 

 

 

 

 

где F (x) –

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

(неопределенный

интеграл)

функции

f (x), С – произвольная постоянная.

 

 

 

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

2 9 3

Задачей Коши, или начальной задачей, называют задачу нахождения решения у = у(х) уравнения (1), удовлетворяющего начальному условию

y(x0) = y0.

(2)

Геометрически это означает, что ищется интегральная кривая уравнения (1), проходящая через заданную точку M0(x0; y0) плоскости Oxy (рис. 1).

Для дифференциальных уравнений важное значение имеет вопрос о существовании решения задачи Коши и единственность этого решения.

Имеет место теорема Коши. Если правая часть уравнения (1) непрерывна и имеет

непрерывную производную fy

в области D, то решение ДУ (1) Рис. 1 с начальным условием (2), где (x0; y0) D, существует и един-

ственно, т.е. через точку (x0; y0) D проходит единственная интегральная кривая данного уравнения.

Если во всех точках решения y = ψ (x) уравнения (1) условие единственности не выполняется, то такое решение называется особым. Через каждую точку M0(x0; y0) особого решения, кроме этого решения, проходит и другое решение уравнения (1), которое не совпадает с y = ψ (x) в сколь угодно малой окрестности этой точки.

Особым решением является огибающая семейства интегральных кривых (если она существует), т.е. линия, которая в каждой своей точке касается хотя бы одной интегральной кривой.

Для существования особого решения уравнения (1) необходимо, чтобы не выполнялись условия теоремы Коши.

Пример 2. Рассмотрим уравнение

 

 

 

 

y = 3y 2 3 .

 

 

 

(3)

Правая часть уравнения (3) f (x, y) = 3y2/3 определена и непрерывна

во всех точках плоскости Oxy. Частная производная

f

=

2

обращается

y

3 y

2 9 4

Рис. 2

в бесконечность при у = 0, т.е. на оси Ox, так что при у = 0 нарушается одно из условий теоремы Коши.

Следовательно, в точках оси Ox возможно нарушение единственности. Непосредственной подстановкой убеждаемся, что функция y = (x + C)3 (семейство кубических парабол) есть решение уравнения (3). Кроме того, уравнение (3) имеет очевидное решение y 0 (рис. 2).

Таким образом, через каждую точку оси Ox проходит, по крайней мере, две интегральные кривые, и, следовательно, в точках этой оси нарушается eдинственность, т.е. y = 0

– особое решение уравнения (3).

Общим решением ДУ (1) называется функция

y = ϕ(x, C),

(4)

зависящая от одной произвольной постоянной С, и такая, что выполняются условия: 1) она удовлетворяет уравнению (1) при любых допустимых значениях постоянной С; 2) каково бы ни было начальное условие (2), можно подобрать значение C0 постоянной С, что решение y = ϕ (x, C0) будет удовлетворять заданному начальному условию (2). При этом предполагается, что точка (x0; y0) принадлежит области D, где выполняются условия существования и единственности решения задачи Коши (1), (2).

Частным решением ДУ (1) называется решение, получаемое из общего решения (4) при каком-либо определенном значении произвольной постоянной С.

Пример 3. Проверить, что функция y = Cex – общее решение уравнения y′ − y = 0 , и найти частное решение, удовлетворяющее начальному условию у(1) = –1.

Решение. Имеем

y = Cex и y′ = Ce x . Подставляя

в данное урав-

нение выражения y

и y, получаем Cex – Cex 0,

т.е. функция

y = Cex удовлетворяет данному уравнению при любых значениях постоянной С.

Зададим любое начальное условие y(x0) = y0. Подставляя в функцию

2 9 5

y = Cex, будем иметь y0 =

Cex0 , от-

куда

C = y0 ex0 .

Функция

y = y0 ex x0 удовлетворяет данному

начальному условию.

Поэтому

функция y = Cex – общее решение

данного уравнения.

 

При x0 = 1 и y0 = –1 получим

частное

решение y = – ex–1.

С геометрической точки зре-

 

ния общее решение определяет

 

семейство интегральных кривых,

 

которыми являются графики пока-

 

зательных функций; частное реше-

 

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

Рис. 3

ходящая через точку M0(1; –1) (рис.

 

3).

 

Иногда общее решение ДУ в виде (3) найти нелегко. В связи с этим введем следующее определение: соотношение вида Ф(х, у, С) = 0, неявно определяющее общее решение, называется общим интегралом ДУ первого порядка; соотношение, получаемое из общего интеграла при конкретном значении постоянной С, называется частным интегралом ДУ.

20. Составление дифференциальных уравнений. Если известно однопараметрическое семейство функций, то для него можно построить ДУ, решениями которого будут все функции этого семейства. Может оказаться, что семейство функций задает общее решение или общий интеграл ДУ.

Пример 4. Для семейства функций y =

C

найти ДУ.

x 1

 

чим

y′ +

y =

Решение. Продифференцировав исходное соотношение,

полу-

y′ = −

C

. Подставив сюда C = y(x – 1), будем

иметь

(x 1)2

y= 0 . Это и есть искомое уравнение; его общее решение

x 1

C

x 1 .

Пример 5. Для семейства функций x2 + y2 Cx = 0 найти ДУ.

2 9 6

Решение. Из условия задачи находим

 

x 2

+

y 2

C =

 

 

 

. Дифференци-

 

x

 

 

 

 

 

 

руем соотношение x2 + y2 Cx = 0 и в полученное уравнение подставляем значение С. Тогда 2x 2 + 2xy y′ − x 2 y 2 = 0 , т.е. 2xy y′ + x 2 y 2 = 0 – искомое ДУ; x2 + y2 – Cx = 0 – его общий интеграл.

Отметим, что к составлению и интегрированию ДУ приводят многие задачи математики и других наук – физики, биологии, химии, экономики и т.д.

В экономике ДУ часто может быть составлено исходя из экономического смысла производной. В качестве примера рассмотрим динамику рыночных цен (макромодель Домара), где независимой переменной служит время t. Допустим, что для конкретного продукта функции спроса Qd и предложения Qs имеют следующий вид:

 

 

 

 

 

 

 

 

Qd = α β P,

(5)

 

 

 

 

 

 

 

 

Qs = –γ+ δP,

 

где α ,

β, γ, δ– некоторые положительные постоянные, Р – цена продук-

та.

Цена равновесия

 

находится из условия Qd = Qs

и будет равна

 

 

 

P

 

 

=

α

+ γ

– конкретной положительной постоянной.

 

P

 

 

 

 

 

 

 

β + δ

 

 

 

 

 

 

 

 

Если случится, что начальная цена P (0) точно равна P , то, очевидно, что рынок будет в положении равновесия и не нужен динами- ческий анализ. В более интересном случае P (0) P , и тогда, если P будет достижимой, то только после соответствующего процесса приспособления, во время которого меняется не только цена, но и Qd , Qs как функции Р, причем все эти переменные можно трактовать как функции времени. Поставим вопрос, который касается динамики цен, так: пусть есть достаточно времени на то, чтобы произошел процесс приспособления; будет ли цена Р (t) со временем приведена до цены равновесия, т.е. будет ли Р (t) стремиться к P , когда t → ∞.

Чтобы ответить на этот вопрос, нужно вначале найти Р (t). Последнее, в свою очередь, требует описать структуру изменения цен. Вообще говоря, цены определяются через релятивное воздействие сил спроса и предложения на рынке. Для простоты допустим, что уровень (скорость) изменения цены (относительно времени) в каждый момент времени прямо пропорционален разности Qd – Qs в этот момент. Такая структура изменения цены может быть выражена так:

2 9 7

dP

=

j(Qd Qs ) ,

(6)

dt

 

 

 

где j – некоторый установленный положительный множитель (коэффициент) приспособления.

Подставляя (5) в (6), запишем (6) в виде:

 

dP

= j(α − β P + γ− δP ) = j(α + γ) j(β + δ)P ,

 

 

dt

 

 

 

 

 

 

или

 

 

 

 

 

 

 

dP

+ j(β + δ)P = j(α + γ) .

(7)

 

 

 

dt

 

 

 

 

 

Уравнение (7) и есть ДУ для цены продукта Р (t) в макромодели Домара. Решение уравнения (7) и его анализ будут приведены на следующей лекции после изложения методов интегрирования ДУ первого порядка.

! Задания для самостоятельной работы

1. Дано семейство функций y = x + C 1+ x 2 ,

где С – параметр. Найти ДУ с заданным семейством решений. 2. Для семейства функций х + у + С (1 – ху) = 0 найти ДУ.

3. Проверить, что функция y = х + С есть общее решение ДУ y= 1, и найти частное решение, удовлетворяющее начальному условию у(0) = 0. Дать геометрическое истолкование результата.

4.Показать, что для уравнения y′ = y12 в каждой точке оси Ox нарушается единственность решения.

5.Показать, что данные функции являются решениями указанных дифференциальных уравнений:

a) y =

sin x

, xy′ + y = cos x ;

x

 

 

б) y = 2 + C 1x 2 , (1x 2 )y′ + xy = 2x .

6. Составить ДУ кривых, у которых точка пересечения любой касательной с осью абсцисс имеет абсциссу, вдвое меньшую абсциссы точки касания.

2 9 8

Лекция 48

Интегрирование дифференциальных уравнений первого порядка

Рассмотрены основные ДУ первого порядка и методы их интегрирования.

Рассмотрим ДУ первого порядка

 

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

(1)

10. ДУ

с разделяющимися переменными.

Пусть

P(x, y) = P1(x), a

Q(x, y) = Q2(y), тогда уравнение (1) называют ДУ с

записанное в дифференциалах, где P(x, y), Q(x, y) – заданные функции

двух переменных х и у.

 

разделенными переменными.

Используя инвариантность формы дифференциала, получаем общее решение

P1(x )dx + Q2 (y) dy = C ,

т.е. ДУ с разделенными переменными решается простым интегрированием.

Пусть P(x, y) = P1(x)P2(y), a Q(x, y) = Q1(x)Q2(y), тогда уравнение (1) называют ДУ с разделяющимися переменными.

Поделив это уравнение на P2(y)Q1(x) ≠ 0, тем самым сведем его к уравнению с разделенными переменными:

P1(x )

Q2 (y)

 

 

P1(x )

Q2 (y)

 

 

 

dx +

 

 

dy = 0

 

 

dx +

 

 

dy = C .

Q

(x )

P

(y)

Q

(x )

P

(y)

1

 

 

2

 

 

 

1

 

 

2

 

 

При таком преобразовании можно потерять решения x = x0, где Q1(x0) = 0 и y = y0, где P2(y0) = 0. Эти случаи нужно рассматривать отдельно. А именно: найти xk и ym такие, что Q1(xk) = 0, P2(ym) = 0, и проверить, являются ли x = xk или y = ym решениями исходного уравнения и содержатся ли они в общем интеграле при каком-то значении Ck (соответственно Cm).

Пример 1. Решить уравнение x(1 + y2) dx + y(1+x2) dy = 0.

Решение. Разделим данное уравнение на (1 + x2)(1 + y2). Переменные разделяются:

2 9 9

2x dx

+

2y dy

 

=

0

ln( x 2 + 1) + ln( y 2 + 1) = ln C

1+ x 2

y 2 + 1

 

 

 

 

(x2 + 1)(y2 + 1) = C – общий интеграл данного уравнения.

Функции Q1(x) = 1 + x2 и Q2(y) = y2 + 1 при вещественных х и у в нуль не обращаются. Поэтому полученная формула содержит все решения рассматриваемого уравнения.

20. Однородные ДУ первого порядка. Функция Р(х, у) назы-

вается однородной степени m, если для любых t она удовлетворяет равенству

P(tx, ty) = tmP(x, y).

Пример 2. Проверить на однородность следующие функции:

 

 

 

P

(x, y) = 2x2 + 3y2, P

(x, y) =

 

x 3

3xy,

 

 

 

 

 

 

2y

 

 

 

 

 

1

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P3 (x, y) =

 

 

xy + y 2 e

 

y

, P4(x, y) = exy.

 

 

 

 

 

 

x 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Решение.

P

(tx, ty) = 2(t2x2) + 3(t2y2) = t2P

(x, y).

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

Аналогично, P

 

(tx, ty) = t2P

 

(x, y). P3 (tx

,ty ) =

tx ty + t 2 y 2 etx ty

=

P3 (x, y) ;

2

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t 2 x 2

 

 

P4 (tx ,ty ) = etxty =

(exy

)t 2 = (P4 (x, y)) t 2

t m P4 (x, y)

ни при каком m.

Значит, P1(x, y),

P2(x, y) – однородные с m = 2, P3(x, y)

однород-

ная с m = 0, P4(x, y)

– неоднородная функция.

 

 

 

 

 

ДУ (1) называют однородным, если Р(х, у) и Q(х, у) – однород-

ные функции одной степени.

 

 

 

 

 

 

 

 

 

 

 

 

Однородное уравнение сводится к уравнению с разделяющимися

переменными заменой у = ux, где u

новая искомая функция. Дей-

ствительно, положив t = x1 (x

 

y

=

1

 

P 1,

 

 

 

P(x, y)

 

x m

 

x

 

 

Получим уравнение

 

y

 

y

P 1,

 

dx +

Q 1,

 

dy = 0

 

 

 

x

 

x

0),

будем иметь

 

 

y

=

1

 

,

Q 1,

 

 

 

Q(x, y) .

 

x m

 

 

x

 

 

 

 

 

y

 

 

y

 

или

P1

 

 

dx +

Q1

 

 

dy =

0 .

 

 

 

 

 

x

 

 

x

 

3 0 0