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

Математический_анализ_2_семестр

.pdf
Скачиваний:
23
Добавлен:
24.05.2015
Размер:
2.05 Mб
Скачать

Заметим, всегда, применяя метод вариации, надо делить на коэффициент при старшей производной, т.е. приводить уравнение.

Пусть найдено решение однородного уравнения yoo x C1 y1 x ... Cn yn x .

Варьируем произвольные постоянные, ищем решение неоднородного уравнения в виде

yx C1 x y1 x ... Cn x yn x .

Дифференцируем это соотношение

yx C1 x y1 x ... Cn x yn x C1 x y1' x ... Cn x yn' x .

Потребуем, чтобы

C1 x y1 x ... Cn x yn x 0 .,

тогда yx C1 x y1' x ... Cn x yn' x .

Дифференцируем еще раз

yx C1 x y1' x ... Cn x yn' x C1 y1'' x ...Cn yn'' .

Потребуем, чтобы

C1 x y1 x ... Cn x yn x 0 .,

тогда yx C1 x y1'' x ... Cn x yn'' x .

Вновь дифференцируем и т.д., в результате, после n-2 дифференцирования получим

C

x y n 2 x ... C

x y

n 2 x 0 .

 

 

 

 

 

 

 

 

 

 

 

1

 

1

 

 

 

 

n

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

y

n 1 x C x y n 1 x ... C x y n 1

x .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

1

 

 

 

 

n

n

 

 

 

 

 

 

 

 

 

 

Дифференцируем и подставляем

 

 

 

 

 

 

 

 

 

 

 

y

n

x C y n

... C

n

y n + C x y n 1 x ... C

x y

n 1 x .

 

 

 

 

 

 

 

 

1

1

 

 

n

 

1

1

 

n

 

n

 

 

 

 

 

 

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

 

yn Ly f x .

 

 

 

 

 

 

 

 

 

C y n ... C

n

y n

+ C x y

n 1 x ... C

x y

n 1 x = L C y

...C

n

y

n

f x

1 1

 

 

 

n

1

1

 

 

 

 

n

n

 

1 1

 

 

 

 

Так

как

y1 ,... , yn -

 

 

решения

однородного

 

уравнения,

то

ykn Lyk

0, k 1,...n .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Получим C ' y n 1 ... C

' y

n 1 f x .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

1

 

 

 

n

n

 

 

 

 

 

 

 

 

 

 

 

Это – последнее уравнение системы для определения варьированных

констант. Соберем все уравнения в систему для определения констант.

C1 x y1 x ... Cn x yn x 0 ., C1 x y1 x ... Cn x yn x 0 .,

........................................................

C1' y1n 1 ... Cn' ynn 1 f x .

Так как определитель системы – определитель Вронского, не равный нулю в силу линейной независимости решений, то функции C1 x ,...Cn x определяются из этой системы однозначно.

61

Теперь общее решение неоднородного уравнения определяется по формуле yx C1 x y1 x ... Cn x yn x .

Лекции 17-18. Линейные дифференциальные уравнения с постоянными коэффициентами.

Начнем изучение линейных дифференциальных уравнений с постоянными коэффициентами с однородных уравнений второго порядка. Дело в том, что в приближенных инженерных расчетах, в инженерной практике, в исследовании процессов и систем все часто строится на анализе систем, моделями которых служат линейные дифференциальные уравнения с постоянными коэффициентами первого и второго порядка. Вспомним, например, что вся механика строится на втором законе Ньютона, который можно записать в виде дифференциального уравнения второго порядка. Основные элементарные функции являются решениями уравнений первого и второго порядков. Экспонента является решением уравнения x ax , sin x, cosx, shx, chx - решения

уравнения x x 0 .

Рассмотрим линейное однородное дифференциальное уравнение с

постоянными коэффициентами второго порядка y py qy 0 .

Будем

 

искать

его решение в виде y e kx .

Подставляя y в

дифференциальное уравнение, получим

 

ekx k 2

pk q 0. Так как ekx 0, то имеем

 

k 2 pk q 0

-

характеристическое уравнение.

Решая его, получим

корни

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p

 

 

p

2

 

 

k1,2

 

 

 

 

 

 

q .

 

2

 

 

 

 

 

 

 

2

 

 

 

Возможно три случая:

 

1) k1 , k2

 

действительны и различны,

 

2)k1 i , k2 i - комплексно сопряженные корни,

3)k1 k2 - действительный кратный корень.

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

y ek1x ,

y

2

ek2 x .

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Для того, чтобы доказать, что решения составляют фундаментальную

систему решений и общее решение записывается в виде

 

 

y

oo

C ek1x C

2

ek2 x ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

надо проверить линейную независимость

 

 

y1 , y2 . Составим определитель

Вронского

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y2

 

 

ek1x

ek2 x

 

e k1 k2 x

 

1

1

 

k k

 

e k1 k2 x 0 , так как

W

y1

 

 

 

 

 

2

 

 

 

y'

y'

 

 

 

k ek1x

k

 

ek2 x

 

 

 

k k

 

 

1

 

 

 

 

 

 

 

2

 

 

 

2

 

 

 

 

 

 

 

1

2

 

 

 

1

 

 

 

 

 

1

 

 

 

 

 

k1 k2 .

62

Заметим, что для уравнения второго порядка проверять линейную

независимость можно проще. Надо показать, что y1 m const . Тогда столбцы y2

определителя Вронского линейно независимы и W 0 . В нашем случае

k1 i ,

k2 i

при k1

k2 .

 

 

 

 

В случае комплексно

сопряженных корней

k1 i ,

k2 i ,

применяя

формулу

Эйлера eiz cos z i sin z,

 

получим

комплексно

сопряженные решения

 

 

 

e x cos x i sin x . Так

y e x cos x i sin x , y

2

 

 

1

 

 

 

 

как линейная комбинация решений линейного однородного уравнения тоже

 

 

1

 

 

 

 

 

 

 

 

1

 

 

 

 

является решением, то

y

 

y

y

2

e x cos x,

y

2

 

 

 

y

y

2

e x sin x

 

 

 

 

1

2

1

 

 

 

 

 

2i

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

являются решениями. Они линейно независимы, так как

 

y1

 

ctg x 0 .

 

y2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

yoo e x C1 cos x C2 sin x .

В случае кратного действительного корня

k1 k2

k одно из решений

можно выбрать в форме

 

y

ekx .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

u x ekx . Подставим в

Второе

решение

 

будем выбирать

 

в

виде

y

2

дифференциальное уравнение, чтобы определить u x .

 

 

 

 

 

'

 

 

 

kx

 

 

kx

e

kx

u

 

ku ,

 

 

 

 

 

 

 

 

 

 

 

y2

u e

 

kue

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y2''

ekx

u ku ku k 2u ekx u 2ku k 2u ,

 

 

 

 

 

e

u

 

2ku

 

k u p u

 

ku qu e

 

u

 

u

p 2k u k pk q 0

 

kx

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

kx

 

 

 

 

 

 

 

 

2

Так как k - корень характеристического уравнения, то k 2 pk q 0 . Так

как k

еще и кратный корень, то по теореме Виета

 

k1 k2 k k 2k p .

Поэтому

 

p 2k 0 .

 

Для определения u x

имеем уравнение

u 0 , отсюда

u x ax b . Выберем a 1, b 0 , получим u x x .

 

 

 

 

 

Следовательно,

 

 

y

2

u x ekx xekx .

Решения

y , y

2

линейно независимы,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

так как y2 x m . y1

Поэтому общее решение линейного дифференциального уравнения с постоянными коэффициентами в случае кратного корня можно записать по формуле

yoo ekx C1 C2 x .

Примеры.

 

 

 

1) y y 0, k 2 1 0, k1,2 1, y1 ex , y2 e x ,

yoo

C1ex C2e x

2) y 5 y 6 y 0,

k 2 5k 6 0, k1 2, k2 3,

yoo C1e2 x C2 e3x

3) y y 0,

k 2 1 0, k1,2 i, y1 cosx, y2 sin x,

yoo C1 cosx C2 sin x

63

4) y 2y 5y 0,

 

k 2 2k 5 0, k1,2 1 2i, y1 e x cos2x, y2

e x sin 2x,

y

oo

e x C cos 2x C

2

sin 2x .

 

 

 

 

1

 

 

 

 

 

5) y 2y y 0,

k 2 2k 1 0, k1,2

1, y1 e x , y2 xe x ,

 

 

y

oo

e x C C

2

x .

 

 

 

 

 

 

1

 

 

 

 

 

Рассмотрим теперь линейное однородное дифференциальное уравнение

n - го порядка с постоянными коэффициентами. y n a1 y n 1 ... an 1 y an y 0 .

Будем искать его решение в виде y e kx . Дифференцируя и подставляя в

дифференциальное уравнение, получим характеристическое уравнение k n a1k n 1 ... an 1k a .

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

Если среди корней характеристического уравнения есть простой

действительный корень

k

1

,

то ему соответствует частное решение

y

ek1x в

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

фундаментальной системе решений и слагаемое

C e k1x в y

oo

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

Если все корни характеристического уравнения k1 ,...kn действительны и

различны,

 

то

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

 

 

им

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

будут

равны

y ek1x ,... y

n

ekn x .

Покажем, что эти решения линейно независимы. Составим

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

определитель Вронского

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ek1x

 

ek2 x

 

...

ekn x

 

 

 

 

1

...

1

 

 

 

 

 

 

 

 

 

1

 

 

W x

 

k ek1x

 

k

2

ek2 x ...

k

n

ekn x

e k1

... kn x

k

 

k

 

..

k

 

 

 

 

1

 

 

 

 

 

 

 

 

1

 

 

2

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

...

 

 

...

 

 

...

 

 

...

 

 

...

 

... ... ...

 

 

 

 

k n 1ek1x

 

k n 1ek2 x ...

k n 1ekn x

 

 

k n 1

 

k n 1

...

k n 1

 

 

 

 

1

 

 

2

 

 

 

 

 

n

 

 

 

 

1

 

2

 

 

n

 

 

e k1 ...kn x

 

ki

k j

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 j i n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Полученный определитель известен в алгебре как определитель Вандермонда, он равен нулю только, когда какие-либо из корней совпадают.

Так как корни различны, то определитель Вронского не равен нулю, следовательно, решения y1 ek1x ,... yn ekn x линейно независимы и составляют

фундаментальную систему решений. Поэтому

yoo C1ek1x ... Cn ekn x .

Если среди корней имеется действительный корень k кратности r, то ему соответствуют частные решения

y

 

ekx , y

2

xekx , y

3

x 2 ekx , ... y

r

x r 1ekx и группа слагаемых в общем

1

 

 

 

 

 

 

 

 

 

 

 

 

решении

...ekx C

 

 

 

 

 

 

 

 

x r 1 ...

y

oo

C

2

x C

3

x 2

... C

r

 

 

 

1

 

 

 

 

 

 

 

 

Если среди корней имеется простая пара комплексно сопряженных

корней

 

k1 i , k2

 

i ,

то им

 

соответствуют частные решения в

64

фундаментальной системе решений

 

y e x cos x, y

2

 

e x sin x

 

и

группа

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

слагаемых в общем решении

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

00

 

...e x (C cos x C

2

sin x) ...

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Если

 

среди

 

 

корней

 

имеется

пара

комплексно

 

 

сопряженных

корней

k1

i ,

 

k2 i ,

кратности r,

то им соответствуют частные решения

в

 

 

фундаментальной

 

 

 

системе

 

решений

 

 

y e x

cos x,

y

2

e x sin x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

3

 

e x x cos x,

 

 

y

4

e x x sin x ... y

2r 1

e x x r 1 cos x,

 

 

 

y

2r

e x x r 1 sin x

 

 

и

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

группа слагаемых в общем решении

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

00

...e x [(C

1

cos x C

2

sin x) .x C

3

cos x C

4

sin x

... x r 1 C

2r 1

cos x C

2r

sin x ] ...

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Примеры.

 

 

 

 

k 4 1 k 1 k 1 k i k i 0, k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y 4

 

y 0,

 

 

 

1

1, k

2

 

1, k

3

i, k

4

i,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

1

 

e x , y

2

e x , y

3

cos x, y

4

sin x,

 

 

 

y

oo

C e x C

 

e

x

 

C

3

cos x C

4

sin x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y 6

 

y 4

0,

 

k 6 k 4 k 4 k 1 k 1 0, k

 

 

0, k

5

1, k

6

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1,2,3,4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y 1, y

2

 

 

x, y

3

x

2 , y

4

 

x3 , y

5

e x , y

6

e x

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

oo

C C

2

x C

3

x

2 C

4

x3

C

5

e x C

e x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Линейное неоднородное дифференциальное уравнение n –го порядка с

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y n

 

a1 y n 1 ... an 1 y an y f x .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Теорема о наложении частных решений.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пусть y1 x - решение неоднородного уравнения с правой частью

 

f1 x ,

 

 

y2

x

-

 

 

решение

 

 

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

с

 

правой

 

 

частью

 

f 2 x . Тогда

y1

x y2 x - решение неоднородного уравнения с правой частью

 

f1

x f 2 x .

Доказательство. Подставим y1 x y2 x в неоднородное уравнение:

y1 x y2 x n a1 x y1 x y2 x n 1 ... an x y1 x y2 x

y1 x n a1 x y1 x n 1 ... an x y1 x

y2 x n a1 x y2 x n 1 ... an x y2 x f1 x f2 x .

По

теореме

о

структуре решения неоднородного

уравнения

yон x yoo

x yчн

x .

Общее решение однородного уравнения

мы строить

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

Метод подбора формы частного решения.

65

 

 

Рассмотрим сначала уравнение второго порядка

 

 

 

y py qy f

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1)

Пусть

 

правая

часть

 

представляет собой

 

квазиполином

 

 

 

f x P x e x .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ищем частное решение в виде y

ч

x Q x e x . Здесь P

x

- полином n-ой

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

степени, Q x - полином, степень которого надо определить.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

yч

x Q x e

 

Q x e

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

yч

x Q x e

 

Q x e

Q

x e

 

Q x

e

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y py qy e

x

Q x (2 p)Q (

2

p q Q x ) e

x

Pn x

 

 

 

 

 

 

 

а) Если - не корень характеристического уравнения, то 2

p q 0 ,

и многочлен Q x

надо выбирать той же степени, что и Pn x , т.е. степени n.

 

 

б)

Если

 

 

 

-

 

простой

корень

 

характеристического

уравнения, то

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

p q 0,

 

2 p 0 . В этом случае многочлен Q x

надо выбирать той

 

 

же

степени, что

и

Pn x ,

т.е.

степени

n. Тогда степень

многочлена надо

выбирать равной

n+1. Однако

при

 

дифференцировании

Q x

производная

свободного члена (постоянной) равна нулю, поэтому Q x можно выбирать в

виде Q x = xQn x .

 

 

 

 

 

 

 

 

 

 

 

 

 

в) Если

 

 

- кратный

корень

характеристического уравнения, то

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

p q 0,

 

2 p 0 . В этом случае многочлен Q x надо выбирать той

 

 

же степени, что и

 

Pn x , т.е. степени n.

Тогда степень многочлена Q x надо

выбирать равной n+2. Однако при двукратном дифференцировании Q x

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

только

свободного

члена

 

равна нулю, но и производная

линейного

 

члена

 

равна

нулю.

Поэтому

Q x

можно выбирать в виде

Q x = x 2Q

n

x .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример. y y x e x

 

 

 

 

 

 

 

 

 

 

k 2 1 0 ,

k

1

1, k

2

1,

y

oo

C e x C

e x

 

 

 

 

 

 

 

 

 

 

 

 

1

 

2

 

 

 

 

f1 x x,

Pn

x x, 0 ,

 

0

-

не

корень характеристического

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

часть,

y Ax B,

y A,

y 0 . Подставляем в неоднородное уравнение с

 

 

1

x

1

1

 

правой частью f1

x .

 

 

Ax B x B 0, A B 1

 

y1

x x .

 

 

 

 

f

2

x e x ,

P x 1,

1 .

Корень 1 содержится один раз среди

 

 

n

 

 

 

корней характеристического уравнения, поэтому частное решение ищется в

виде y

2

Dxex ,

y

 

Dex 1 x ,

y

 

Dex 2 x .

 

 

 

 

 

 

2

 

 

2

 

 

 

x e x .

Подставляем в неоднородное уравнение с правой частью

f

2

 

 

 

 

 

 

 

 

 

 

 

66

Dex 2 x x e x D 12 y2 x 12 xe x .

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

yчн x 12 xe x .

Общее решение неоднородного уравнения будет

y

 

C e x C

e x x

1

xex .

он

 

 

1

2

2

 

 

 

 

 

 

 

f x e x M x cos x N x sin x

2) Правая часть имеет вид

a)

Если

i не корни

характеристического уравнения, то частное

решение ищется в том виде, в котором задана правая часть: yч e x U m x cos x Vm x sin x ,

где U m x ,Vm x - полиномы степени m – максимальной из степеней полиномов M x , N x .

б) Если i - пара корней характеристического уравнения, то частное

решение ищется в виде

yч xe x U m x cos x Vm x sin x ,

Пример. y y sin x

k 2 1 0, k

1,2

i,

y

oo

С cos x C

2

sin x

 

 

 

1

 

 

f x sin x,

0, 1,

M x , N x

степени 0

Пара корней i = i - пара корней характеристического уравнения. yч x Acos x B sin x

yч Acosx B sin x Axsin x Bxcosx,

yч 2Asin x 2B cosx Axcosx Bxsin x

Подставляем в неоднородное уравнение, получаем

2Asin x 2Bcosx sin x , откуда B 0, A 12 yч 12 x cosx , yон C1 cosx C2 sin x 12 x cosx

Рассмотрим неоднородное уравнение n-го порядка, покажем, как в нем применять метод подбора формы частного решения.

Здесь ситуация сложнее, так как в характеристическом уравнении n корней, действительные корни и комплексно сопряженные, простые и кратные корни.

1) Пусть правая часть неоднородного уравнения имеет вид f x e x Pn x

67

a)Если не является корнем характеристического уравнения, то частное решение неоднородного уравнения ищется в том же виде, что

и правая часть

y

ч

e x Q

n

x .

 

 

 

 

b) Если - корень характеристического уравнения r-ой кратности, то частное решение неоднородного уравнения ищется в виде yч x r e x Qn x .

2)Пусть правая часть неоднородного уравнения имеет вид f x e x M x cos x N x sin x

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

y

ч

e x U

m

cos x V

m

sin x , где степень m многочленов –

максимальная из степеней многочленов M x , N x .

b)Если пара комплексно сопряженных корней является корнями характеристического уравнения r-ой кратности, то частное решение неоднородного уравнения ищется в виде

yч x r e x U m cos x Vm sin x .

Пример.

y 5 y x sin x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k 5 k 2 k 2 k 1 k 2 k 1 0, k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0, k

 

1, k

 

 

1

i

 

3

,

 

1,2

 

3

4,5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

y

 

C

C

 

x C

e x e 2

C

 

cos

 

 

x C

 

sin

 

 

x .

 

 

 

 

 

 

 

 

 

oo

2

4

 

 

 

5

 

 

 

 

 

 

 

 

 

 

 

 

1

 

3

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f1 x x, 0, Pn x x, n 1 .

 

 

 

0

 

 

 

 

содержится

 

в

корнях

характеристического уравнения 2 раза,

поэтому

 

yч1

x 2 Ax B . Подставляя

это частное решение в неоднородное уравнение с правой частью

f1 x x ,

получим 6Ax 2B x A

1

, B 0,

y

 

 

 

 

 

1

x3 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ч1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6

 

 

 

 

 

 

 

 

6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f2 x sin x,

0, 1,

i i . Корни i не содержатся в корнях

характеристического уравнения,

поэтому

 

 

 

yч,2 D cos x E sin x . Подставляя это

частное

решение

в неоднородное уравнение с

 

правой

частью

 

f2 x sin x ,

получим E D cosx E D sin x sin x, E D 0, E D 1 E D

1

.

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

D cos x E sin x . y

 

 

 

1

(cosx sin x) .

 

 

 

 

ч,2

ч,2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

1

x3 +

1

 

(cosx sin x) .

 

 

 

 

 

 

 

 

ч

 

 

 

 

 

 

 

 

 

 

 

 

6

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример. y 5

y x sin x

 

 

 

 

 

 

 

 

k 5 k 3 k 3 k i k i 0, k

0, k

4,5

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1,2,3

 

 

 

 

 

 

 

y

oo

 

C

C

2

x C

x 2 C

4

cos x C

5

sin x .

 

 

 

 

 

 

 

1

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

f1 x x, 0, Pn

x x, n 1

 

0

содержится

в

корнях

характеристического уравнения 3 раза, поэтому yч1 x3 Ax B .

68

f2 x sin x,

0, 1, i i .

Корни

i (пара

корней)

содержатся в корнях характеристического

уравнения

один раз,

поэтому

yч,2 x(D cos x E sin x) . Неопределенные коэффициенты определяются, как и

выше, подстановкой в уравнение и сравнением коэффициентов при одинаковых степенях x, при sinx, cosx, xsinx, xcosx.

Лекции 19-20. Нормальные системы дифференциальных уравнений.

Система дифференциальных уравнений – это система уравнений относительно независимой переменной x, функций этой переменной и их

производных

y , y

,

y

...y(m1 ) ...y

n

, y

n

...y mn

. Система может быть записана в

 

 

1

1

 

1

1

 

 

 

 

 

n

 

общем виде

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F ( x,

y , y , y

 

...y(m1 ) ...y

n

, y

 

...y mn )=0

 

1

1

1

1

 

1

 

 

 

n

 

n

 

....................................................................

 

F ( x,

y , y

, y

...y(m1 ) ...y

n

, y

...y mn )=0

 

n

1

1

1

 

1

 

 

 

 

n

 

n

 

Порядок этой системы равен m1 ... mn .

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

 

 

y

m1

 

1

( x, y

, y

 

, y

...y

(m1 1)

...y

n

, y

 

...y

mn 1 )

 

 

 

 

 

1

 

 

 

 

1

 

1

 

1

1

 

 

 

 

 

n

 

n

 

 

 

 

 

 

 

..................................................................................

 

 

 

 

 

y

mn

 

 

 

 

( x, y , y

 

, y

 

...y

(m1 1)

...y

 

, y

 

...y

mn 1

)

 

 

 

 

 

n

 

 

1

 

 

 

 

 

n

n

n

 

 

 

 

 

 

 

 

 

 

 

1

 

1

 

1

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Теорема.

Любое

 

 

дифференциальное

 

уравнение,

разрешенное

относительно

 

 

старшей

 

 

производной,

 

 

можно

свести

к

системе

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

 

 

 

 

 

 

 

Доказательство. Рассмотрим дифференциальное уравнение n-ого порядка

 

 

y

n

f x, y,y

 

y

n 1

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Обозначим

 

 

 

 

...

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

0

x y x , y

0

x y ,

y

x y

2

,..y

 

 

y

n 1

.

Дифференциальное

уравнение

 

 

 

 

 

 

 

 

 

1

 

1

 

 

 

 

 

n 2

 

 

 

 

 

 

 

 

n-ого порядка

удалось

свести к

 

системе n

 

дифференциальных

уравнений

первого порядка

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

x y ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

x y

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

....................

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

yn 2

 

yn 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f

 

x, y0 , y1

 

yn 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

yn 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

69

y10 x y1 x

 

 

 

 

y

y

 

 

 

 

10

11

 

 

 

 

y

y

 

 

 

 

11

12

 

 

 

 

................

x, y10 ,... y1 m1 1 ,.... yn0 ,... yn mn 1

 

 

 

 

y1,m1 1 f1

.........................................................................................

yn 0

x yn

x

 

 

 

 

 

yn 1

 

 

 

 

yn 0

 

 

 

 

.................

 

x,

y10 ,... y1 т1 1 ,... yn 0 ,... yn mn 1

 

 

 

 

 

yn mn 1 fn

Получена система

 

из

m1 ...mn дифференциальных уравнений первого

порядка.

Удобнее нормальную систему дифференциальных уравнений (систему в нормальной форме) записывать в виде:

y

 

f

1

x,

 

y ,... y

n

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

..................................

 

(покоординатная форма)

y

 

 

f

n

x,

y ,... y

n

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

или в виде

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

f

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

1

 

 

y

f

x, y

, где y

 

 

... ,

f

f 2

 

(векторная форма).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

yn

 

 

...

 

 

 

 

 

 

 

 

 

 

y

 

sin y

cos y

 

 

 

 

 

 

Пример.

 

1

 

 

 

 

 

1

 

 

 

2

 

Эти

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

xy

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

1

 

 

 

2

 

 

 

 

 

 

 

системе

( y10 y1 )

y10 y11

y11 sin y11 cos y21

( y20 y2 )

y20 y21

y21 xy10 y21

Оказывается, не только дифференциальное уравнение n- ого порядка сводится к системе n дифференциальных уравнений первого порядка – нормальной системе, но и нормальная система может быть сведена к одному дифференциальному уравнению.

70