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

pismennyy_d_t_konspekt_lekciy_po_vysshey_matematike_polnyy_k

.pdf
Скачиваний:
86
Добавлен:
12.05.2015
Размер:
11.99 Mб
Скачать

51.2. Метод вариации произвольных постоянных

Рассмотрим ЛНДУ (51.1). Его общим решением является функ­

ция (51.3), т. е.

у= у* +у.

Частное решение у* уравнения (51.1) можно найти, если известно об­

щее решение у соответствующего однородного уравнения (51.2), мето­ дом вариации nроизвол.ьнwх постоянных (метод Лагранжа), состоя­ щим в следующем. Пусть у= с1у1 (х)+с2у2(х) - общее решение уравне­

ния (51.2). Заменим в общем решении постоянные с1

и с2 неизвестными

функциями с1 (х) и с2(х) и подберем их так, чтобы функция

у*= с1(х) · У1(х) + с2(х) · у2(х)

(51.6)

была решением уравнения (51.1). Найдем производную

(у*)' =с~(т)у1(т) +ci(x)y~(x) +с~(х)у2(х) +с2(х)у~(х).

Подберем функции с1 (т) и с2(т) так, чтобы

 

с~ (х) · У1 (х) + с~(х) · У2(х) =О.

(51.7)

Тогда

(у*)'= с1(х) · у~(х) + с2(х) · у~(х), (у*)"= с~ (х) ·у~ (х) + с1 (х) · у~'(х) + с;(х) · у~(х) + с2(х) · у~(х).

Подставляя выражение для у*, (у*)' и (у*)" в уравнение (51.1), полу­

чим:

с~(х) ·у~ (х) + с1 (х) · у~(х) + ~(х) · у~(х) + с2(х) · у~(х)+

+а1 (х) (с1 (х)у~(х) +с2(х)у~(х)] +а2(х) [с1(х)у1 (х) +с2(х)у2(х)] = f (х),

или

с1(х) · [у~'(х) +а1(х) ·у~(х) +а2(х) ·yi(x)]+

+с2(х) [у~(х) +а1(х)у~(х) +а2(х)у2(х)] +с~(х)у~ (х)+с~(х)у~(х) = f (х).

Поскольку у1 (х) и у2(х) - решения уравнения (51.2), то выражения в

квадратных скобках равны нулю, а потому

с~ (х) ·у~ (х) + с~(т) · у~(х) = J(x).

(51.8)

Таким образом, функция (51.6) будет частным решением у* уравне­ ния (51.1), если функции с1 (х) и с2(х) удовлетворяют системе уравне­

ний (51.7) и (51.8):

{с~~х) · Yi(x) + с~(х) · у2(х) =О,

(51.9)

с~ (х) ·у~ (х) + ~(х) · у~(х) = J(x).

360

Определитель системы 1~~~~~ ~~~~~ 1-::J. О, так как это определитель

Вронского для фундаментальной системы частных решений у1( х) и У2(х) уравнения (51.2). Поэтому система (51.9) имеет единственное ре­ шение: ci(x) = rp1(x) и с~(х) = rp2(x), где rp1(x) и rp2(x) - некоторые функции от х. Интегрируя эти функции, находим с1 (х) и с2(х), а затем по формуле (51.6) составляем частное решение уравнения (51.1).

Пример 51.1. Найти общее решение уравнения у" + у = _1_.

COSX

Q Решение: Найдем общее решение у соответствующего однородного

уравнения у" +у= О. Имеем: k2 +1 =О, k1 = i, k2 = -i. Следовательно,

у = с1 · cos х + с2 · sin х. Найдем теперь частное решение у* исходного уравнения. Оно ищется в виде (51.6): у* = с1 ( х) ·cos х + с2( х) ·sin х. Для нахождения с1 (х) и с2(х) составляем систему уравнений вида (51.9):

сНх) · cosx + с~(х) ·sinx =О,

 

 

 

 

 

{ci(x) · (-sinx) + с~(х) ·cosx =

-

 

1

-.

 

 

 

 

 

 

 

cosx

 

Решаf'м ее:

 

 

 

 

 

 

 

 

 

 

Д =

cosx

sinxl

=

2

 

2

Х = 1,

.

cosx

COS

Х +

 

l

-sшх

 

 

 

 

 

 

 

 

Л1 = 1 01

sinx 1= -tgx,

 

Л2=1 со~х

 

_01 1=1,·

--

COSX

 

 

 

-SШХ

 

 

 

 

 

 

 

 

 

с;(х) = ~1 = - tgx,

с1(х) = J(- tgx) dx = ln 1cosxl;

~(х) = ЛЛ2 = 1, с2(х) = ! 1 · dx = х.

Запишем частное решение данного уравнения: у*

= ln 1cos xl · cos х +

+х · sin х. Следовательно, общее решение данного уравнения имеет вид

у= (у+ у*)= с1 ·cosx + с2 ·sinx + cosx · ln 1cosxl + х ·sinx.

е

При нахождении частных решений ЛНДУ может быть полезной

следующая теорема.

 

 

Теорема 51.2 (о наложении решений). Если правая часть уравне­

ния (511} представляет собой сумму двух функций: /(х) = fi(x)

+

+ /2(х). а Yi и У2

- частные решения уравнений у"+ а1 (х) ·у'+

+ а2(х) ·у= fi(x)

и у"+ al(x) ·у'+ а2(х)у = /2(х) соответственно,

то функция у* = Yi + У2 является решением данного уравнения.

361

~ ,Действительно,

(yi +у;)"+ ai(x) · (yi +у;)'+ az(x) · (yi +у;)=

= ((у;)"+ ai (х) ·(yi)' + а2(х) · Yi) + ((у;)"+ ai (х) ·(у;)'+ а2(х) ·у;) =

= fi(x) + f2(x) = f(x). 8

51.З. Интегрирование ЛНДУ второго порядка

спостоянными коэффициентами

иправой частью специального вида

Рассмотрим ЛНДУ второго порядка с nостоянн'ЫМи коэффициен­

тами, т. е. уравнение

j у"+Р. у'+ q. у= f(x), J

(51.10)

где р и q - некоторые числа.

Согласно теореме 51.1, общее решение уравнения (51.10) предста­ вляет собой сумму общего решения у соответствующего однородного уравнения и частного решения у* неоднородного уравнения. Частное

решение уравнения (51.10) может быть найдено методом вариации про­

извольных постоянных (п. 51.2).

Для уравнений: с постоянными коэффициентами (51.10) существу­

ет более простой способ нахождения у*, если правая часть f(x) урав­

нения (51.10) имеет так называемый «специальный вид»:

1. f(x) = Pn(x) · e0 z

или

11. f(x) = e0 z · (Pn(x) · cosj3x + Qm(x) · sin/Зx).

Суть метода, называемого методом неоnреде.л.енн'ЫХ коэффициен­

тов, состоит в следующем: по виду правой части f (х) уравнения (51.1 О)

записывают ожидаемую форму частного решения с неопределенными коэффициентами, затем подставляют ее в уравнение (51.10) и из полу­ ченного тождества находят значения коэффици~нтов.

Слу'ЧаiJ. 1. Правая часть (51.10) имеет вид J(т) = Pn(x) · е°'х, где о: Е 1R, Pn(x) - многочлен степени n. "Уравнение (51.10) запишется в

виде

у"+р. у'+ q. У= Pn(x). eaz. J

(51.11)

1

В этом случае частное решение у* ищем в виде:

 

 

ly* = Хт. Qn(x). еах, 1

(51.12)

где r - число, равное кратности о: как корня характеристического

уравнения k2 + pk + q = О (т. е. r - число, показывающее, сколько раз а является корнем уравнения k2 + pk + q =О), а Qn(x) = A0 xn + + А1х"'-1 + ... + An -.-- многочлен степени п, записанный с неопреде­

ленными коэффициентами Ai (z = 1, 2, ... , п).

362

Q а) Пусть о: не является корнем характеристического уравнения

k2 + pk + q =о,

т. е. о:-:/:- ki,2· Следовательно,

r =О, у*= Qп(х) · ео:х, (у*)'= Q~(x) · ео:х + Qп(х) · ео:х ·о:,

(у*)"= Q~(x) · ео:х + 2Q~(x) · ео:х ·о:+ Qп(х) · ео:х · о:2

После подстановки функции у* и ее производных в уравнение (51.11),

сокращения на ео:х, получим:

Q~(x) + (2о: + p)Q~(x) + (о:2 + ро: + q) · Qп(х) = Pn(x).

(51.13)

Слева - многочлен степени п с неопределенными коэффициентами,

справа - многочлен степени п, но с известными коэффициентами. При­

равнивая коэффициенгы при одинаковых степенях х, получим сисге­ му (п + 1) алгебраических уравнений для определения коэффициентов

Ао,А1, ".,Ап.

б) Пусть о: являегся однократным (простым) корнем характери­

стического уравнения k2 + pk + q =О, т. е. о:= ki -:/:- k2.

В этом случае искать решение в форм(' у• = Qп(х)ео:х нельзя, т. к.

о:2 + ро: + q =О, и уравнение (51.13) принимает вид

Q~(x) + (2о: + р) · Q~(x) = Pn(x).

В левой части - многочлен степени (п - 1), в правой части - много­ член степени п. Чтобы получить тождество многочленов в решении у*,

нужно иметь многочлен степени (п + 1). Поэтому частное решение у* следует искать в виде у* = х · Qп(х)ео:х (в равенстве (51.12) положить

r = 1).

в) Пусть о: является двукратным корнем характеристического

уравнения k2+pk+q =О, т. е. о:= k1 = k2. В этом случае о:2 +pa+q

и 2о:+р =О, а поэтому уравнение (51.13) принимает вид Q~(x) = Pn(x).

Слева стоит многочлен степени п- 2. Понятно, чтобы иметь слева

многочлен степени п, частное решение у* следует искать в виде

у•= x 2 Qn(x)eo:x

(в равенстве (51.12) положить r = 2).

 

C.tty'Ч,aii, 2. Правая часть (51.10) имеет вид

 

f(x) = ео:х · (Pn(x) · соs/Зх + Qт(х) sin/Зx),

где Pn(x) и Qт(х) - многочлены степени пит соответственно, о: и

/3 - действительные числа. Ураnнснис (51.10) запишется в виде

1 у"+ру' + qy = ео:х · (Pn(x) · соs/Зх + Qт(х) · sin/Зx)· 1

(51.14)

363

Можно показать, что в этом случае частное решение у* уравнения

(51.14) следует искать в виде

 

1у*= xr · е0ж · (М1(х) · соs/Зх + N1(x) · sin/Зx), j

(51.15)

где r - число, равное кратности а+ /Зi как корня характеристического

уравнения k 2 +pk + q =О, М1(х) и N1(xy- многочлены степени l с не­

определенными коэффициентами, l - наивысшая степень многочленов

Рп(х) и Qт(х), т. е. l = max(n, m).

Замечания.

1. После подстановки функции (51.15) в (51.14) приравнивают мно­

гочлены, стоящие перед одноименными тригонометрическими функци­

ями в левой и правой частях уравнения.

2.Форма (51.15) сохраняется и в случаях, когда Рп(х) = О или

Qт(х) :=О.

3.Если правая часть уравнения (51.10) есть сумма функций вида

1 или 11, то для нахождения у* следует использовать теорему 51.2 о

наложении решений.

Пример 51.2. Найти общее решение уравнения у"-2у' +у= х-4.

0 Решение: Найдем общее решение у ЛОДУ у" - 2у' +у = О. Характе­

ристическое уравнение k 2 - 2k + 1 = О имеет корень k 1 = 1 кратности 2.

Значит, у = с1 еж + с2 х · еж. Находим частное решение исходного

уравнения. В нем правая часть х - 4 = - 4) · е0 ж есть формула ви­ да Р1 (х) ·е0·ж, причем а = О, не является корнем характеристического

уравнения: а =/:- k1 Поэтому, согласно формуле (51.12), частное реше­

ние у* ищем в виде у*= Q 1 (x) · е0 ж, т. е. у*= Ах+ В, где А и В -

неопределенные коэффициенты. Тогда (у*)'= А, (у"')"= О. Подставив у*, (у*)', (у*)" в исходное уравнение, получим -2А+Ах+В = х-4, или Ах+ (-2А +В) =х - 4. Приравнивая коэффициенты при одинаковых

степенях х, получаем систему уравнений:

{А= 1,

-2А+В = -4.

Огсюда А= 1, В= -2. Поэтому частное решение данного уравнения имеет вид у*= х-2. Следовательно, у= (у+у*) = с1еж +с2хеж +х-2-

искомое общее решение уравнения.

8

При.мер 51.3. Решить уравнение у" - 4у' + 13у = 40 · соsЗх.

0 Решение: Общее решение ЛНДУ имеет вид у = у + у*. Находим

решение однородного уравнения у: у" - 4у' + 13у = О. Характеристиче­

ское уравнение k 2 - 4Jc + 13 = О имеет корни k1 = 2 + Зi, k2 = 2 - Зi.

Следовательно, у = е2ж · (С1 cos Зх + с2 · sin Зх).

364

Находим частное решение у*. Правая часть ЛНДУ в нашем слу­

чае имеет вид f(x) = е0 :r (40соsЗх +О· sinЗx). Так как а= О, {3 = З,

а + (Зi =

Зi не совпадает с корнем характеристического уравнения,

то r = О.

Согласно формуле (51.15), частное решение ищем в виде

у* = А cos Зх + В sin Зх. Подставляем у* в исходное уравнение. Имеем:

(у*)' = -ЗА sin Зх + ЗВ cos Зх, (у*)" = -9А cos Зх - sin Зх. Получаем:

-9АсоsЗх - 9BsinЗx - 4(-ЗAsinЗx + ЗВсоsЗх)+

+lЗ(АсоsЗх + BsinЗx) = 40соsЗх,

или

(-9А-12В+ lЗА) соsЗх+ (-9В + 12А+ lЗВ) sinЗx = 40 cosЗx+O·sinЗx.

Отсюда имеем:

- 12В = 40,

 

 

 

 

{12А + 4В =О.

 

Следовательно, А= 1, В= -З. Поэтому у*= соsЗх -

ЗsinЗx. И нако­

нец, у = е2х ( с1 · cos Зх + с2 · sin Зх) + cos Зх - З sin Зх -

общее решение

уравнения.

 

Пример 51.4. (Д.л.я самосто.ятельного решения.) Для предложен-

ных дифференциальных уравнений получить вид частного решения:

а) у" - Зу' + = 5 + ех;

6) у" - 2у' + у = 2;

в) у"+ 4у = sin + cos 7х;

г) у" +у = 5 cos - х sin 2х;

д) у" - Зу' = х2 - 1 + cos х.

Ответы: а) А + х · В · ех; 6) А; в) х(А cos + В sin 2х) + С cos 7х + + D sin 7х; г) (Ах+ В) cos + (Сх + D) sin 2х; д) х(Ах2 + Вх +С)+

+ D cos х + Е sin х.

51.4.Интегрирование ЛНДУ n-го порядка (n > 2)

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

Рассмотрим линейное неоднородное ДУ п-го (п > 2) порядка y(n) + а1(х) · y(n-l) + а2(х) · y(n- 2) + ... + ап(х) ·У= f(x),

где а1 (х), а2(х), ... , an(x), f(x) - заданные непрерывные функции на (а;Ь).

Соответствующее ему однородное уравнение имеет вид y(n) + а1 (х) ·y(n-l) + ... + ап(х) ·У= О.

365

Теорема 51.З (о структуре общеrо решения ЛНДУ n-ro поряАка).

Общее решение у ЛНДУ п-го порядка равно сумме частного решения

у" неоднородного уравнения и общего решения у соответствующего ему однородного уравнения, т. е. у = у* + fj.

Частное решение у* ЛНДУ п-го порядка может быть найдено, если

известно общее решение у однородного уравнения, методом вариации

произвольных постоянных. Оно ищется в виде

у*= с1(х) · У1(х) + с2(х) · У2(х) + ... + Сп(х) · Уп(х),

где Yi(x), i = 1,п, - частные решения, образующие фундаментальную

систему, однородного уравнения.

Система уравнений для нахождения неизвестных с,(х) имеет вид

C~YI + с~у2 + сЗуз + ... + C~Yn = О, с~у~ + с~у~ + с3у~ + ... + с~у~ =О,

с~у~' + с~у~ + с3У~ + ... + с~у~ =О,

c~yin-l) + c~y~n-l) + cЗy~n-l) + ... + c~y~n-1) = j(x).

Однако для ЛНДУ п-го порядка с постоянными коэффициентами,

правая часть которого имеет специальный вид, частное решение у* мо­

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

Метод подбора частного решения у* уравнения

у<п) +P1Y(n-l) + ···+ PnY = f(x),

где Pi - числа, а правая часть f (х) имеет специальный вид, описанный

в п. 51.3 для случая n = 2, переносится без всякого изменения и на случай уравнения, имеющего порядок п > 2.

Пример 51.5. Решить уравнение y1 v -у'==; 2х.

О Решение: Находим fj:

k4 - k =о, k(k - 1). (k2 + k + 1) =о,

'4 = -~2 ± v'2Зi,

fj = с1 + с2ех + е

-lx(2

J3

J3 )

.

 

сз cos Тх + С4 sin Тх

Находим у*: f(x) = ( = е0·х·Р1(х)), r = 1, у*= х(Ах+В) = Ах2 +Вх,

отсюда

(у*)'= 2Ах +В, (у*)"= 2А, (у*)111 = 0, (y*) 1 v =О.

366

Тогда -(2Ах +В) = 2х. Отсюда А= -1, В= О и получаем у* = 2 .

Следовательно, функция

является общим решением уравнения.

 

§ 52. СИСТЕМЫ fJ.ИФФЕРЕНЦИАЛЬНЫХ

УРАВНЕНИИ 52.1. Основные понятия

Для решения многих задач математики, физики, техники (задач

динамики криволинейного движения; задач электротехники для не­ скольких электрических цепей; определения состава системы, в ко­

торой протекают несколько последовательных химических реакций 1

порядка; отыскания векторных линий поля и других) нередко требу­

ется несколько функций. Нахождение этих функций может привести к

нескольким ДУ, образующим систему.

Cucmeмoil ДУ называется совокупность ДУ, каждое из которых содержит независимую переменную, искомые функции и их производ­

ные.

Общий вид системы ДУ первого порядка, содержащей п искомых

функций у1,У2, ... ,yn, следующий:

~~~~~~~~~.2.;.· ·""; .~~;. ~~;. ~~~: : : ~~~?.~-~·

{ F . у . у .

 

 

. у . у' . у' .

. у' ) - о

·

n '

1,

2" "'

n•

1•

2""'

n -

Система ДУ первого порядка, разрешенных относительно произ­

водной, т. е. система вида

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

rJ1n.

- f

1

(х· у

 

. у .

. у

п

)

,

 

 

dx

-

 

 

,

1, 2" ·"

 

 

 

~ - f

2

(х·У ·У ·

·У

n

)

,

 

 

dx

-

 

 

,

1, 2" ·"

 

(52.1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

OJm.

-

f

 

(х·У

·У

·

· Уп )

 

 

 

dx

 

n

 

,

 

1,

2" ",

 

 

 

,

 

f§j называется нормаJtьноii cucmeмoi1 ДУ. При этом предполага­

е1'ся, что число уравнений равно числу искомых функций.

Заме-чание. Во многих случаях системы уравнений и уравнения

высших порядков можно привести к нормальной системе вида (52.1).

367

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

d2x

-

F (x·y·z·t·x'·y'·z')

"'di,2-

1

,,,, '''

2

= Р.2(х·у·х'·у'·z')

d y

-;jtI"

 

 

' '''. ' '

2

 

3

,,,, '''

di'1"-

z

-

F

(x·y·z·t·x'·y'·z')

описьшающая движение точки в пространстве, путем введения новых

dx _

OlJ.

_

v,

dz

_

"

переменных: dt -

и, dt

-

 

dt

-

w, приводится к нормальноn систе-

меДУ:

 

 

 

 

 

 

 

 

dx

'

 

 

 

dt

-

 

 

 

 

~

= v,

 

 

 

dt

 

 

 

 

dz

-w

'

 

 

 

dt

 

-

 

 

 

~~= F1(x;y;z;t;u;v;w),

~~= F2(x; у; z; t; и;v; w),

~~= Fэ(x;y;z;t;u;v;w).

Уравнение третьего порядка у111 = f(x; у; у'; у") путем замены у'= = р, у11 =р' = q сводится к нормальной системе ДУ

у' =р,

{рq'' =q,= f(x; у;р; q).

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

ных систем.

~Решением системы (52.1) называется совокупность из п функ­

ций У1, У2, ... , Уп, удовлетворяющих каждому из уравнений этой

системы.

Нача.л,ьные услови.я для системы (52.1) имеют вид

(52.2)

Задача Коши для системы (52.1) ставится следующим образом:

найти решение системы (52.1), удовлетворяющее начальным услови­

ям (52.2).

Условия существоlilания и единственности решения задачи Коши

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

368

Теорема 52.1 (Коши). Если в системе (52.1) все функции

fi(x; У1; · · ·, Yn)

непрерывны вместе со всеми своими частными производными по Yi

в некоторой области D ((п + 1)-мерного пространства), то в каждой

точке М00; у~; у~; ... ;у~) этой области существует, и притом един­

ственное, решение У1 = <р1(х), У2 = ЧJ2(х), ... , Уп = 1Рп(х) системы, удовлетворяющее начальным условиям (52 2)

Меняя в области D точку М0 (т. е. нача.пьныf' условия), получим

бесчисленное множество решений, которое можно записать в виде ре­

шения, -зависящего от п произвольных постоянных:

Э10 решение является общим, если по 3аданным начальным усло­

виям (52.2) можно однозначно опрсдРлить постоянные с1 , с2, ... , Сп иэ

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

{~_1_<~;-~1·;-~2.;."".. ;.~~~·~-~~:

<t'n(x; С1; с2; ... ; Сп) =У~·

Решение, получающееся из общего при конкретных значениях по­

стоянных с1 , с2, ... , Сп, называется 'Частным решением системы (52.1).

52.2. Интегрирование нормальных систем

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

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

(Обратная задача - переход от ДУ к системе - рассмотрена выше на примере.) Техника этого метода основана на следующих соображениях.

Пусть задана нормальная система (52.1). Продифференцируем по

х любое, например первое, уравнение:

d2 y1

8/i

8/i

dy1

8Ji

dy2

8/i

dуп

-

2

= -

+ -

. -

+ -

. -

+ ". + -

. - .

dx

 

ду1

dx

ду2

dx

дуn

dx

Подставив в это равенство значения производных ~, ~, ... , ~~ из

системы (52.1), получим

d2y1

дfi

8fi

8fi

8/i

·f п,

- 2

= -8

+ -8 ·/1

+ -8 ·/2

+ ···+ -

dx

Х

У1

У2

8Уп

 

369