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

Письменный Д.Т. Конспект лекций по высшей математике - Полный курс

.pdf
Скачиваний:
26430
Добавлен:
10.03.2016
Размер:
20.05 Mб
Скачать

или, коротко,

d2y1

-- = F2(x;y1;Y2; ... ;yn)· dx2

Продифференцировав полученное равенство еще раз и заменив значе-

ния производных~'... , ~11;' из системы (52.1), получим

d3y1

- - =F3(х;у1;у2; ... ;уп)· dx3

Продолжая этот процесс (дифференцируем - подставляем - получа­

ем), находим:

(52.3)

dn

dxl(t =Fn(x;y1;y2; ... ;yn)·

Из первых (п-1) уравнений системы (52.3) выразим функции у2, уз, ...

... , Уп через х, функцию У1 и ее производные у'1 , у1/1 , ••• , у1(n-1) . Полу-

чим:

2

--

о/2 (х·' у1,. у'1.>".,.

1

'

 

 

 

 

 

.1,

 

 

 

y(n-1))

 

 

у

 

_ ф (х·у

·у'·

·у(п-1))

 

 

.

~-~-

..3.. .' .. ~'.. ~:::::

. ~... • . :

 

 

 

 

 

 

 

 

 

(52.4)

 

у

n

-

.1,

(х· у

. у'.

. y(n-1))

 

 

-

o/n

'

1, 1'

· · · '

1

·

Найденные значения yz, уз, ... , Уп подставим в последнее уравнение си-

стемы (52.3). Получим одно ДУ п-го порядка относительно искомой

функции У1= a;,'fP = Ф(х;у1;у~;.";у~п-l)). Пусть его общее решение

есть

У1 = ср1(х;с1;сz;".;сп)·

Продифференцировав его (п - 1) раз и подставив значения произ-

водных у~,у~', ... , Yin-l) в уравнения системы (52.4), найдем функции

У2, Уз,···, Yn=

У2 = 1Р2(х; с1; с2; ... ;Cn), ... , Yn = 1Рп(х; С1; с2; ... ; Сп)·

При.мер 52.1. Решить систему уравнений

{~ =4y-3z,

~~ = 2уЗz.

370

Q Решение: Продифференцируем первое уравнение: у" = 4у' - Зz'.

Подставляем z' = 2у- 3z в полученное равенство: у" = 4у' - 3(2уЗz), у" - 1 + = 9z. Составляем систему уравнений:

{у'= 4y-3z,

 

 

 

у" -

4у' + = 9z.

 

Из первого уравнения системы выражаем z через у и у':

 

 

4у-у'

 

(52.5)

 

 

z= --3-.

 

Подставляем значение z во второе уравнение последней системы:

11

4 1 + 6 - 9(- у')

'

у

- у

У-

3

т. е. у" - у' - = О. Получили одно ЛОДУ второго порядка. Решаем

его: k 2 - k- 6 = О, ki = -2, k2 = 3 и у = с1е-2х +с2е3х - общее решение

уравнения. Находим функцию z. Значения у и у' = 1е-2х + с2е3х)' =

= -2с1е-2х + Зс2е3х подставляем в выражение z через у и у' (форму­

ла (52.5)). Получим:

Таким образом, общее решение данной системы уравнений имеет

вид у = С1е-2х + с2е3х, z = 2с1е-2х + ic2e3x.

8

Заме-чание. Систему уравнений (52.1) можно решать методом ин­ тегрируемъtх комбинациft. Суть метода состоит в том, что посредством арифметических операций из уравнений данной системы образуют так называемые интегрируемые комбинации, т. е. легко интегрируемые уравнения относительно новой неизвестной функции.

Проиллюстрируем технику этого метода на следующем примере.

Пример 52.2. Решить систему уравнений:

{~~=у+ 1,

~х + dt = 1.

Q Решение: Сложим почленно данные уравнения: х'+у' = х+у+2, или

(х+у)' = (х+у)+2. Обозначим х+у = z. Тогда имеем z' = z+2. Решаем

полученное уравнение:

d+z 2

= dt, ln(z + 2) - lnc1

= t, z + 2 =

et

Z

 

 

С1

'

z + 2 = с1et, или х + у = с1et - 2.

Получили так называемый nepвuft интеграл системы. Из него

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

371

(52.6).

уменьшить на единицу число искомых функций. Например, у = с1et -

- 2 - х. Тогда первое уравнение системы примет вид

х' =

c1et - 2 - х + 1, т. е. х' + х = c1et - 1.

Найдя из него х

(например, с помощью подстановки х = uv), найдем

и у.

Заме'Чание. Данная система «позволяет» образовать еще одну ин­

тегрируемую комбинацию: х' - у' =

у - х, т. е. (х - у)' =

-(х - у).

Положив х - у = р, имеем: р' = -р,

или 1:Е = -dt, lnp -

ln с2 = -t,

 

р

 

р = c2 e-t, или х - у = c2 e-t. Имея два первых интеграо-~а системы, т. е.

х +у = с1et - 2 и х - у = c2 e-t, легко найти (складывая и вычитая пер-

вые интегралы), что х =

1

с1е

t

1

с2е

- t

-

1

- 1

t

-

1

- t

-

1

.

2

 

+ 2

 

 

, у - 2с1е

 

2с2е

 

 

 

52.3.Системы линейных ДУ с постоянными коэффициентами

Рассмотрим еще один метод интегрирования нормальной системы

уравнений (52.1) в случае, когда она представляет собой систему ли­

неtiн'Ых однороднЪtх ДУ с постоянными коэффициентами, т. е. систему

вида

d

~~·:.":_+ ~:'~~

+_ ~~~~~·

 

 

{ :

 

 

dx

= ап1У1 + ап2У2 + ···+ annYn·

 

Для простоты ограничимся рассмотрением системы трех уравне-

ний с тремя неизвестными функциями У1, У2 и у3:

 

 

 

~ = анУ1 + а12У2 + а1зУз,

 

 

11

 

 

 

 

 

п_

= az1Y1 + а22У2

+ аzзУз,

(52.6)

 

 

dx

 

 

~

= аз1У1 + аз2У2

+ аззУз,

 

 

 

dx

 

где все коэффициенты aij

( i, j = 1, 2, 3) -

постоянные.

 

Будем искать частное решение системы (52.6) в виде

 

 

У1 =а· ekx, У2 = f3 ekx,

Уз= [ ·ekx,

(52.7)

где а, (3, /, k-постоянные, которые надо подобрать (найти) так, чтобы функции (52.7) удовлетворяли системе

Подставив эти функции в систему (52.6) и сократив на множитель ekx =j:. О, получим:

ak = ана + ai2f3 + а~з/,

{fЗk = а21а + а22/З + а2з/, f'k = аз1а + аз2/3 + азз/,

372

или

{ (а11 - k)a + a12f3 + а13')' =О,

 

 

а21а: + (а22 - k)(З + а2з')' =О,

(52.8)

 

аз1а: + аз2f3 + (азз - k)f' =О.

 

Систему (52.8) можно рассматривать как однородную систему трех ал­

гебраических уравнений с тремя неизвестными а, (3, 1'· Чтобы эта си­

стема имела ненулевое решение, необходимо и достаточно, чтобы опре­ делитель системы был равен нулю:

а1з

 

 

а2з

=О.

(52.9)

азз - k

~Уравнение (52.9) называется харакmерисmическим уравнени­

ем системы (52.6). Раскрыв определитель, получим уравнение тре­

тьей степени относительно k. Рассмотрим возможные случаи.

Слу'Ча11 1. Корни характеристического уравнения действительны и

различны: ki, k2, k3 . Для каждого корня ki (i = 1, 2, З) напишем систе­ му (52.8) и определим коэффициенты a:i, fЗi, Т'i (один из коэффициентов

можно считать равным единице). Таким образом, получаем:

 

для корня ki

частное решение системы (52.6): ур) = a 1ek1x, у~1) =

_ (3

ek1x

у(1)

_

"'ek1x.

-

1

' 3

-

11

'

 

для корня k2

- Yi2) = a2ek2x' у~2) = f32ek2x' у~2) = ')'2ek2x;

для корня kз - ур) = а:зеkзх' у~з) = fЗзеkзх' у~з) = ')'зеkзх.

Можно показать, что эти функции образуют фундаментальную систе­

му, общее решение системы (52.6) записывается в виде

 

У1 = c1a1ek

1x + c2a:2ek2x + с3а3еkзх,

 

У2 = c1f31ek

1x + c2f32ek2x + сзfЗзеkзх,

(52.10)

Уз = С1')'1 ekix + C2')'2ek2x + С3')'3еkзх.

При.мер 52.3. Решить систему уравнений:

{~ =у1 -у2,

1:1п. - -4у1 + У2· dx -

Q Решение: Характеристическое уравнение (52.9) данной системы

имеет вид

 

 

 

1- k

-1

1 =

о,

1 -4

1- k

 

или 1-2k+k2 -4 =О, k2 -2k-3 =О, k1 = -1, k2 = 3. Частные решения

данной системы ищем в виде ур) = а1ekix, у~1) = (31 ek1x и у?) = a 2ek2 x,

у~2) = (32 ek2 x. Найдем a:i и f31 (i = 1, 2).

373

При k1 = -1 система (52.8) имеет вид

 

 

{ (1- (-l))a1 - /31

=О,

 

т. е.

 

 

-4а1 + (1 - (-l))/31 =О,

 

 

 

 

Эта система имеет бесчисленное множество решений. Положим а1

= 1,

тогда (31 = 2. Получаем частные решения

 

 

(1)

 

 

 

 

У1

= е

 

 

 

 

При k2 = 3 система (52.8) имеет вид

 

 

 

{-2а2 - /32 =О,

 

 

 

-4а2 -

2(32 = О.

 

 

Положим а2 = 1, тогда (32 = -2. Значит, корню k2

= 3 соответствуют

частные решения:

 

 

 

 

 

у~2)

= езх

и

у;2) = _2езх.

 

 

Общее решение исходной системы, согласно формуле (52.10), запи-

шется в виде: У1 = С1е-х + с2е3х, У2 = 2с1е-х - 2с2е3х.

8

Cлy"taii, 2. Корни характеристического уравнения различные, но

среди них есть комплексные: k1

= а + ib, k2 = а -

ib, k3 • Вид частных

решений в этой ситуации определяют так же, как и в случае 1.

За.ме"lание. Вместо полученных частных решений можно взять их линейные комбинации (п. 50.1, случай 3), применяя формулы Эйле­

ра; в результате получим два действительных решения, содержащих

функции вида еах · cos Ьх, еах · sin Ьх. Или, вьщеляя действительные и

мнимые части в найденных комплексных частных решениях, получим

два действительных частных решения (можно показать, что они тоже являются решениями уравнения). При этом понятно, что комплексно­ сопряженный корень k2 = а - ib не даст новых линейно независимых

действительных решений.

Пример 52.4. Найти частное решение системы

\~=У1+У2,

~ = -у1 + У2 -у3,

0:1& 3У2 +Уз,

dx =

удовлетворяющее начальным условиям: У1(О) = 7, У2(О) = 2, Уз(О) = 1.

374

О Решение: Составляем и решаем характеристическое уравнение:

1-k

1

 

 

о

 

 

 

-1

 

1-k

 

 

-1

=0,

 

о

 

3

 

1- k

 

 

 

(1 - k) ·11; k

 

1-=.\ \- 1·\~1

1-=.\ 1=о,

 

(1 - k)(k2 - 2k + 4) -

(k - 1) =о,

(1 - k)(k2 - 2k + 5) =о,

 

k1 = 1,

k2

= 1+ 2i,

= 1- 2i.

 

Для k1 = 1 получаем:

 

 

 

 

 

 

 

 

о. а1 + /31 +о. 'Yl =о,

 

{-а1 + О · /31

-

'У1 = О,

 

 

О · а1 + 3/31

+ О · "(1

=

О

 

(см. (52.8)). Отсюда находим: /31 =

О,

а1

=

1 (положили), 'У1

-1.

Частное решение системы: yi1) = ех, у~1) = О, у~1) = -ех.

 

Для k2 = 1+2i получаем (см. (52.8)):

-2ia2 + /32 =О,

{-а2 - ~i/32 - 'У2 = О,

3/32 - 2i"(2 = о.

Отсюда находим: а2 = 1 (положили), /32 = 2i, 2 = 3. Частное комп­

лексное решение системы:

Yi2) = e(l+2i)x, у~2) = 2ie(1+2i)x, у~2) = 3e(I+2i)x.

В найденных решениях выделим действительную (Re) и мнимую (Im) части:

yi2) = e(1+2i)x = ex(cos2x + isin2x),

Re yi2) = ех cos 2х, Im yi2 ) = ех sin 2х;

у~2) = 2ie(l+2i)x = ex(2icos2x - 2sin2x),

Rey~2) = -2exsin2x, Imy~2) = 2excos2x;

у~2) = 3e(l+2i)x = ex(3cos2x + i3sin2x),

Rey~2) = 3ех cos2x, Imy~2) = 3ех sin2x.

Как уже отмечено, корень k3 = 1 - 2i приведет к этим же самым реше­

ниям.

375

Таким образом, общее решение системы имеет вид

У1 = с1ех + с2ех cos2x + сзех sin2x,

У2 = с1 ·О - 2с2ех sin + 2сзех cos 2х,

Уз = -С1 ех + 3с2ех cos + 3сзех sin 2х.

Выделим частное решение системы. Пр~ заданных начальных услови~

ях получаем систему уравнений для определения постоянных с1 , с2, сз:

7 = С1+С2+0,

 

 

{ 2 - О+ 2сз,

==> С1 = 5, С2 = 2,

С3 = 1.

1 = -с1 + 3с2 + О,

 

 

Следовательно, искомое частное решение имеет вид

 

У1 = 5ех + 2ех cos + ех sin 2х, У2 = -4ех sin + 2ех cos 2х,

Уз = -5ех + 6ех cos + 3ех sin 2х.

Случаii, 3. Характеристическое уравнение имеет корень k кратно­

сти т (т = 2, 3). Решение системы, соответствующее кратному корню,

следует искать в виде:

 

а) если т = 2, то У1 = (А+ Bx)ekx, У2

= (С+ Dx)ekx, Уз =

= + Fx)ekx;

= (D + Ех + Fx2)ekx,

б) если т = 3, то у1 = (А+ Вх + Cx 2)ekx, у2

у3 = (G + Нх + Nx2 )ekx.

Это решение зависит от т произвольных постоянных. Постоянные А, В, С, ... , N определяются методом неопределенных коэффициентов. Выразив все коэффициенты через т из них, полагаем поочередно один из них равным единице, а остальные равными нулю. Получим т ли­

нейно независимых частных решений системы (52.6).

При.мер 52.5. Решить систему уравнений:

!

~-

У1 - У2 +Уз,

dx -

rJш. -

У1 + У2 - Уз,

dx -

 

01& -

2

 

dx -

-у2 + УЗ·

Q Решение: Составляем и решаем характеристическое уравнение

1-k -1 1

1 1-k -1 =0,

о -1 2-k

376

(1- k)(2 - 2k - k + k2-1) - 1(-2+k+1) =О, kl = 2, k2 = = 1. Корню

kl = 2 соответствует система (см. (52.8)):

-а:1 - /31

+ 'Yl =о,

{/31 =о,

{а:1 - /31 -

'Yl = о,

:==} а:1 - 'Yl =о.

-/31 =о,

 

 

Полагая 'Yl = 1, находим а:1 = 1. Получаем одно частное решение ис­

ходной системы: ур) = е2х, у~1) =О, у~1) = е2х.

Двукратному корню k = k2 = k3 = 1 = 2) соответствует реше­

ние вида у~2'3) = + Вх)ех, у~2'3) = + Dx )ех, у~2'3) = + Fx )ех.

Подставляем эти выражения (решения) в уравнения исходной системы:

{В· ех +(А+ Вх)ех =(А+ Вх)ех - (С+ Dx)ex + + Fx)ex, D · ех +(С+ Dx)ex =(А+ Вх)ех +(С+ Dx)ex - + Fx)ex, F · ех + + Fx)ex =-(С+ Dx)ex + 2(Е + Fx)ex,

или, после сокращения на ех "1 О и группировки,

(D - F)x +В+ С - Е = О,

{ - F)x - D - Е =О,

(D - F)x +С+ F - Е = О.

Эти равенства тождественно выполняются лишь в случае, когда

D-F=O,

B-F =О,

В+С-Е=О,

A-D-E=O,

C+F-E=O.

Выразим все коэффициенты через два из них (m = 2), например че­

рез А и В. Из второго уравнения имеем F = В. Тогда, с учетом пер­

вого уравнения, получаем D = В. Из четвертого уравнения находим

Е =А - D, т. е. Е =А - В. Из третьего уравнения: С= Е - В, т. е. С = А- В - В, или С = А- 2В. Коэффициенты А и В - произвольные.

Полагая А = 1, В = О, находим: С = 1, D = О, Е = 1, F = О.

Полагая А= О, В= 1, находим: С= -2, D = 1, Е = -1, F = 1.

Получаем два линейно независимых частных решения, соответ-

ствующих двукратному корню k = 1:

у(2) _ ех

у(2) _ ех

у(2) _ ех и

1 - ,

2 - ,

3 -

ур) =хех, у~з) = (-2+х)ех, у~з) = (-1+х)ех.

Записываем общее решение исходной системы:

 

У1 = с1е2х + с2ех + сзхех, У2 = с2ех + сз(х -

2)ех,

Уз= с1е2х + с2ех + сз(х - 1)ех.

8

Глава XI. ДВОЙНЫЕ И ТРОЙНЫЕ

ИНТЕГРАЛЫ

1 Лекции 44-46 1

§ 53. ДВОЙНОЙ ИНТЕГРАЛ

53.1. Основные понятия и определения

Обобщением определенного интеграла на случай функций двух пе­

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

у

Пусть в замкнутой области D плос­

кости Оху задана непрерывная функция

 

 

z = f(x;y). Разобьем область D на п

 

«элементарных областей» Di (i = 1, п),

 

площади которых обозначим через лsi,

 

а диаметры (наибольшее расстояние ме­

- 0--i---------x-

жду точками области) -

через di (см.

рис. 214).

 

 

 

В каждой области D i

выберем про­

Рис. 214

извольную точку

Mi(x,; Yi), умножим

значение f(xi; Yi)

функции в этой точке

 

на дS; и составим сумму всех таких произведений:

n

f(x1;Y1)ЛS1 + f(x2; У2)ЛS2 + ···+ f(хп; Yn)ЛSn = L f(xi;Yi)дSi.

i=l

(53.1)

Эта сумма называется интегра.л:ьноi1 сум.моii, функции f(x; у) в обла­

сти D.

Рассмотрим предел интегральной суммы (53.1), когда п стремит­

ся к бесконечности таким образом, что maxdi -+О. Если этот предел существует и не зависит ни от способа разбиения области D на ча­

сти, ни от выбора точек в них, то он называется iJвoii,нuм интегралом

от функции f(x;y) по области D и обозначается// f(x;y)dxdy (или

// f(x;y)ds). D

D

~Таким образом, iJвohнoii, интегра.л. определяется равенством

 

JL"Jo

n

 

// f(x;y)dxdy =

Lf(xi;y;). лsi.

(53.2)

D

(max ci;--+0)

i=l

 

~В этом случае функция f(x;y) называется интегрируемоit в

об,11,асти D; D - область интегрирования; х и у - пере­

менные интегрuровани..я; dx dy (или dS) - эл.емент n.лoщaiJu.

378

Для всякой ли функции f (х; у) существует двойной интеграл? На

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

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

Теорема 53.1 (достаточное условие интегрируемости функции). Если функция z = f(x;y) непрерывна в замкнутой области D, то она

интегрируема в этой области.

Заме'Чаnи.я.

1. Далее будем рассматривать толь­

ко функции, непрерывные в области инте­

грирования, хотя двойной интеграл может

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

функций.

2. Из определения двойного интегра­ ла следует, что для интегрируемой в обла­

сти D функции предел интегральных

у

Лу;

сумм существует и не зависит от спосо­

о

Лх;

х

 

 

ба разбиения области. Таким образом, мы

Рис. 215

можем разбивать область D на площадки

 

прямыми, параллельными координатным осям (см. рис. 215). При этом ЛSi = Лхi · Луi, равенство (53.2) можно записать в виде

jj f (х;у) dx dy =

 

n

·Ji..IIJo

L, f (xi; у;) · Лх; · Луi.

D

(max d;--+0)

i=l

53.2. Геометрический и физический смысл двойного

интеграла

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

Объем цилиндрического тела

Рассмотрим тело, ограниченное сверху поверхностью z = f(x; у) ~ ~О, снизу - замкнутой областью D плоскости Оху, с боков - цилин­ дрической поверхностью, образующая которой параллельна оси Oz, а

направляющей служит граница области D (см. рис. 216). Такое тело

называется цилиндрu'Ческим. Найдем его объем V. Для этого разобьем

область D (проекция поверхности z = f(x;y) на плоскость Оху) произ­ вольным образом на п областей D;, площади которых равны ЛSi (i = = 1,п). Рассмотрим цилиндрические столбики с основаниями Di, огра­ ниченные сверху кусками поверхности z = f(x; у) (на рис. 216 один из

379