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

УМК

.PDF
Скачиваний:
91
Добавлен:
17.03.2015
Размер:
3.93 Mб
Скачать

 

y2

= α = const; т.е. решения y1 и y2 линейно зависимы, что противоречит

 

 

 

y1

 

 

 

 

 

 

 

 

предположению об их линейной независимости.

 

 

 

 

 

Допустим далее, что y1=0 в точках x1,x2,…,x k, принадлежащих отрезку [a; b].

Рассмотрим интервал (а,x1). На этом интервале

y1 0 . Следовательно, на

основании только что доказанного следует, что на интервале (а,x1)

 

 

 

 

 

y2

= λ = const, или y

2

= λ y .

 

 

 

 

 

 

 

 

 

 

 

y1

1

 

 

 

 

 

 

 

 

 

 

Рассмотрим функцию y = y2 − λy1. Так как y2 и y1 есть решения уравнения

(1.23),

то y = y2 − λy1 - решение уравнения (1.23)

и y 0 на интервале (a;x1).

Следовательно, на основании замечания в начале

доказательства следует, что

 

y = y

2

− λy 0 на отрезке [a; b], или

y2

= λ на

[a; b], т.е. y2 и y1 - линейно

 

 

 

 

 

1

 

 

 

 

 

 

y1

зависимы.

Но это противоречит предположению о линейной независимости решений y2 и y1. Мы доказали, что определитель Вронского не обращается в ноль ни в одной из точек отрезка [a; b].

Теорема 1.5 (о структуре общего решения линейного однородного уравнения). Если y1,y2,… yn - линейно независимые частные решения однородного

линейного

уравнения

a0 y(n ) + a1y(n 1)

+ ... + a n y = 0 ,

то

y = C1y1 + C2 y2 + ... + Cn yn - общее

решение

этого уравнения (С1,C2,…,C n

-

произвольные постоянные).

 

 

 

 

Доказательство. Докажем теорему на

примере

линейного однородного

уравнения 2-го порядка (n=2). Сначала покажем, что y=C1y1+C2y2 - решение уравнения (1.23), т.е. y′′ + a1y′ + a 2 y = 0 . Подставим функцию y=C1y1+C2y2 и её

производные

 

 

y′ = C y

+ C

2

y

,

 

 

 

 

 

y′′ = C y′′ + C

2

y′′

в

уравнение

(1.23):

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

1

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

1

1

 

 

 

 

 

 

2

 

 

 

 

 

(C y + c

2

y

2

)+ a

1

(C y + c

2

y

2

)

+ a

2

(C y + c

2

y

 

) = 0,

 

 

 

 

1

1

 

 

 

 

 

 

 

 

 

 

1

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

1

 

 

 

 

 

2

 

 

 

 

 

 

 

C y′′ + C

2

y′′

+ a C y′ + a C

2

y

 

 

+ a

2

C y + a

2

C

2

y

2

= 0,

 

 

 

 

 

1

1

 

 

2

 

1

 

1

1

1

 

 

 

2

 

 

 

 

 

 

1

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C

(y

+ a

y

+ a

 

y

)+ C

2

(y

 

 

+ a

1

y

 

+ a

2

y

2

)= 0.

 

 

 

 

 

 

 

 

1

1

 

 

 

1

 

1

 

 

 

2

 

1

 

 

 

 

 

2

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Т.к.

 

 

 

y1

 

 

 

 

и

 

 

 

y2

+ a

 

 

-

+ a

частные

 

 

 

 

 

 

 

решения

 

(1.23),

тогда

y′′

+ a y

 

+ a

2

y = 0 и

y′′

 

y

2

y

2

= 0 .

Имеем

C 0 + C

2

0 = 0

- верное

1

 

1

1

 

 

 

 

 

1

 

 

 

 

2

 

 

 

 

1

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

равенство, таким образом y

 

 

= C1y1 + C2 y2

 

 

 

-

 

 

решение

уравнения (1.23) для

любых С1 и С2.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Теперь

 

 

 

докажем,

 

что

 

 

каковы

 

 

бы

 

 

 

ни

 

были

начальные

условия

 

= y0 , y

 

x =x 0

= y0 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

y

x =x 0

 

 

 

 

можно

 

 

 

так

 

 

подобрать

значение

произвольных

постоянных

C1

и

 

C2,

 

чтобы

 

 

 

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

 

 

частное

решение

C1y1 + C2 y2

удовлетворяло заданным начальным условиям.

 

 

 

 

 

 

 

 

 

 

Подставляя начальные условия в равенствоy = C1y1 + C2 y2 , имеем

 

 

 

 

 

 

 

 

 

 

 

 

y

0

 

= C (y ) + c

2

(y

2

)

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

1

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1.24)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= C (y)

+ c

 

(y)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

0

2

0

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

1

 

0

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

где

 

обозначено

 

y

 

x = x 0

=

(y

 

)

0

,

 

 

 

y

2

 

 

= (y

2

)

0

,

y

 

x = x 0

= (y) ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

1

 

 

 

 

 

 

 

 

x = x 0

 

 

1

 

1 0

 

y2

 

x = x 0 = (y2 )0 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Из системы (1.24) можно определить c1 и c2, т.к. определитель этой системы

 

(y1 )0

(y2 )0

 

 

= (y

)

(y)

 

(y

 

)

 

(y)

 

есть определитель Вронского при x=x0 и,

 

 

 

 

 

 

(y) (y

)

 

 

1

0

2

0

 

 

2

0

 

 

1

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

0

2

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

следовательно, не равен нулю (в силу линейной независимости решений y1 и y2).

Частное решение, которое получится из семейства

y = C1y1 + C2 y2 при

найденных значениях C1 иC2 , удовлетворяет начальным

условиям. Теорема

доказана.

 

Теорема 1.6 Если известно одно частное решение y1 линейного однородного уравнения (1.23), то второе его решение y2, линейно независимое с первым, можно найти интегрированием первого по формуле

y

2 = y1

e

a1dx

dx

(1.25)

 

y2

 

 

 

1

 

 

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

Формула (1.25) дает возможность интегрировать линейные однородные уравнения 2-го порядка сразу, не прибегая к понижению порядка.

ПРИМЕР 1.19 Записать общее решение уравнения y ′′+ 2 y′ + y = 0, если x

известно его частное решение y1 = sin x . x

Решение. Найдем второе частное решение y2, линейно независимое с первым

 

 

 

 

 

2

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ln

1

 

 

 

 

 

 

 

 

 

 

 

 

e2 ln

 

x

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y2

=

sin x

e

x

dx =

sin x

 

 

 

 

dx =

sin x

 

e

x

 

 

 

dx =

 

 

 

 

 

 

x

sin x 2

x

sin x 2

x

sin x 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

x

 

 

 

x

 

 

 

 

= sin x x

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx =

sin x

 

dx

=

sin

(ctgx) = −

cos x

.

 

 

 

x 2

sin2 x

x

sin 2 x

x

 

 

 

 

x

 

 

 

x 2

 

 

 

 

 

 

 

 

 

Запишем общее решение:

Y = C y

+ C

y

 

= C

sin x

C

 

 

cos x

 

2

 

2

 

1

1

2

 

1

x

 

x .

 

 

 

 

 

 

 

 

1.11 ЛИНЕЙНЫЕ ОДНОРОДНЫЕ ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ n-ГО ПОРЯДКА С ПОСТОЯННЫМИ КОЭФФИЦИЕНТАМИ

ОПРЕДЕЛЕНИЕ 1.29 Линейным однородным дифференциальным уравнением (ЛОДУ) n-го порядка с постоянными коэффициентами

называется уравнение вида

y(n ) + a1y(n 1) + a2y(n 2) + ... + an 1 y+ an y = 0,

(1.26)

где a1,a2,…,a n - постоянные.

 

Общее решение уравнения (1.26) имеет структуру

 

 

 

= C1 y1 + C2y 2 + ... + Cn yn ,

 

 

Y

(1.27)

где y1,y2,…,y n – линейно независимые частные решения однородного уравнения

(1.26) или функции, вронскиан которых не равен нулю.

ОПРЕДЕЛЕНИЕ 1.30 Совокупность n решений ЛОДУ n - го порядка,

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

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

ПРИМЕР 1.20 Дано уравнениеy′′′ - y= 0 . Составляют ли функции ex , e-x , ch x фундаментальную систему решений?

Решение. Для проверки линейной независимости этих решений вычислим

вронскиан:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

W(y , y

 

 

 

) =

 

y1

y

2

y3

 

 

ex

ex

chx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

, y

3

 

y

y

y

 

=

ex

ex

shx

= 0 .

(1.28)

1

 

 

 

1

 

2

3

 

 

 

 

 

 

 

 

 

 

 

 

 

y1′′

y

2′′

y3′′

 

 

ex

ex

chx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Определитель равен нулю, так как соответствующие элементы первой и третьей строк равны. Так как W=0, то данные функции ex , e-x , ch x - линейно зависимы и не составляют фундаментальную систему решений данного ДУ.

Решим задачу о нахождении общего решения ЛОДУ (1.26). Будем искать частные решения в виде функции

y=ekx, где k=const,

(1.29)

тогда y′ = k ekx , y ′′= k2 ekx , ...,

y(n ) = k n ekx .

 

 

 

Подставляя полученные выражения в уравнение (1.26), получим

 

ekx (k n + a1 k n1 + a 2 k n2 + ... + a n1 k + a n )= 0 .

 

Так как ekx ≠ 0 , то

 

 

 

 

 

k n + a1 k n 1 + a 2 k n 2 + ... + a n 1 k + a n = 0.

 

(1.30)

ОПРЕДЕЛЕНИЕ

1.31

Уравнение

вида

(1.30)

называется

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

Если число k удовлетворяет характеристическому уравнению (1.30), то

функция ekx

будет решением однородного уравнения (1.26).

 

 

 

 

 

Для простоты сначала рассмотрим частный

случай

(n=2):

линейное

однородное

 

 

дифференциальное уравнение 2-го порядка

с постоянными

коэффициентами

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y′′ + a1y+ a2y = 0.

 

 

 

 

 

(1.31)

Соответствующее характеристическое уравнение имеет вид

 

 

 

 

 

 

 

 

 

 

 

 

k 2 + a1k + a2 = 0.

 

 

 

 

 

(1.32)

Возможны следующие случаи.

 

 

 

 

 

 

1.

D > 0 . Характеристическое уравнение имеет действительные различные

корни

k

1

≠ k

2

. Тогда частные решения уравнения (1.31)

y

= ek1x ,

y

2

= ek 2 x .

 

 

 

 

 

 

= C1y1 + C2 y2 уравнения (1.31) имеет вид

1

 

 

 

Общее решение

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= C1ek1x + C2ek 2 x .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Y

 

 

 

 

 

(1.33)

Убедимся, что функции y1 и y2 образуют фундаментальную систему решений (линейно независимы). Для этого покажем, что их вроскиан не равен нулю.

 

W(y , y ) =

y1

 

y2

=

ek1x

ek 2 x

 

= k ek1x ek 2 x k ek1x ek 2 x

= (k k )e(k1 + k 2 )x

1

 

 

2

 

 

 

y

 

y

 

k1e

k1x

k2e

k 2 x

 

2

 

 

 

 

 

 

 

 

1

 

 

 

 

2

1

 

 

 

 

 

1

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

W(y1, y2 ) 0 , т.к. k1

k 2 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2.

 

D = 0 .

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

корни

k , кратность m=2. Частные решения y

= ekx , y

2

= x ekx

. Общее решение

 

 

= C1y1 + C2 y2

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

уравнения (1.31) имеет вид

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= C1ekx + C2 x ekx .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Y

 

 

 

 

 

 

 

 

 

 

 

(1.34)

 

 

 

 

Покажем,

что

y2 = x ekx является решением уравнения (1.31). Для этого

подставим y2 в уравнение (1.31)

 

(x ekx )+ a1 (x ekx )+ a 2 (x ekx )= 0 . После

несложных

 

 

 

 

 

вычислений

 

получим

 

 

 

(k 2 + a1k + a 2 )x + (a1 + 2k) = 0 .

k 2 + a1k + a 2 = 0 ,

т.к.

 

k корень

 

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

уравнения

(1.32);

a1 + 2k = 0 ,

т.к.

по условию

D = 0 , следовательно,

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

b ±

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

D

 

 

уравнения вычисляется по формуле k = −

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

k1,2

 

 

2a

 

 

. Убедимся, что

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

и

 

 

 

y

2

 

 

линейно

 

 

независимы,

 

 

 

т.е.

 

y1

≠ const .

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

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

ekx

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

=

 

 

 

 

 

=

 

 

 

const .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xekx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y2

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3.D<0. Характеристическое уравнение имеет пару комплексно -

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

1

= α + βi, k

2

= α − βi.

Частные решения y

= eαx cos βx ,

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

y

 

= eαx sin βx . Общее решение

 

= C y + C

 

y

 

уравнения (1.31) имеет вид

2

y

2

2

 

 

 

 

 

 

 

 

 

1

1

 

 

 

 

 

 

 

= C1eαx cos βx + C2eαx sin βx.

(1.35)

 

 

Y

Убедимся, что y1 и y2 - решения уравнения (1.31).

Так как функции ek1x и ek 2 x - частные решения (1.31) (см. формулу (1.29)), то по теореме (1.5) их линейная комбинация так же является решением уравнения

(1.31). Покажем, что y

и y

2

есть линейная комбинация функций ek1x и ek 2 x .

1

 

 

 

(ek1x + ek 2 x )=

 

(e(α+βi)x + e(α−βi)x )=

Действительно, y = eαx cos βx =

1

1

 

 

1

 

 

2

2

 

 

 

 

 

= 1 (eαx (cos βx + i sin βx) + eαx (cos βx i sin βx))= 1 eαx (2 cos βx). Аналогично,

2

 

1

(ek1x ek 2 x ).

 

2

 

 

y2

= eαx sin βx =

 

 

 

 

 

 

 

 

 

 

 

 

 

2i

 

 

 

 

 

y1

 

 

Теперь убедимся, что y

и y

2

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

const .

 

 

 

 

 

 

 

 

1

 

 

y2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

=

eαx cos βx

= ctg βx const .

 

 

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

 

1

 

 

 

 

y2

eαx sin βx

 

 

 

 

 

 

 

 

 

 

Вернемся к линейным однородным дифференциальным уравнениям n-го порядка (1.26).

Сформулируем необходимые утверждения (уже без доказательств), а затем рассмотрим примеры.

Для корней характеристического уравнения (1.30) возможны следующие случаи.

1.Уравнение (1.30) имеет n простых (кратности 1) действительных корней.

Каждому корню ki (i =

 

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

yi = ek i x . Общее

1, n

решение ЛОДУ (1.26) согласно (1.27) имеет вид

 

 

 

= C1ek1x + C2ek 2 x + ... + Cn ek n x .

 

 

Y

(1.36)

2.Уравнение (1.30) имеет m действительных кратных корней k. Каждому

корню k кратности m соответствует m линейно независимых частных решений ekx , x ekx , x 2 ekx , ..., x m 1 ekx .

3.Уравнение (1.30) имеет комплексно сопряженные корни. Каждой паре

простых комплексно - сопряженных корней

k1 = α + βi

и k 2 = α − βi

соответствуют функции eαx cos x

и eαx sin βx .

 

 

4. Уравнение (1.30)

имеет

комплексно

- сопряженные

кратные корни

k1 = α + βi и k 2 = α − βi

кратности m > 1. Каждой такой паре соответствует 2m

частных решений

 

 

 

 

 

 

eαx cos βx ,

x eαx cos βx ,

x2 eαx cos βx ,

… ,

xm 1 eαx cos βx ;

eαx sin βx ,

x eαx sin βx ,

x2 eαx sin βx ,

 

 

(1.37)

… ,

xm 1 eαx sin βx .

ПРИМЕР 1.21 Найти общее решение уравнения

 

 

 

 

y′′′ − y′′ − 4y′ + 4y = 0.

 

 

Решение. Запишем характеристическое уравнение. Для этого заменим

формально функцию y и ее производные y, y′′ и y′′′

соответствующими степенями

k . y заменяем на k0 = 1, yна k1 , y′′ на k 2 , y′′′ на k3 . Тогда получим k3 k2 4k + 4 = 0, откуда, раскладывая левую часть уравнения на множители,

имеем

k 2 (k 1) 4(k 1) = 0, или (k 1)(k + 2)(k 2) = 0.

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

k1 = 1, k2 = −2, k3 = 2.

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

действительные и различные, то общее решение данного дифференциального уравнения имеет вид Y = C1ex + C2e2x + C3e2x .

ПРИМЕР 1.22 Найти общее решение уравнения y′′′ − 5y′′ + 3y′ + 9y = 0.

Решение. Составим характеристическое уравнение k3 5k2 + 3k + 9 = 0 или k3 + k 2 6k2 6k + 9k + 9 = 0, (k +1)(k2 6k + 9)= 0, (k +1)(k 3)2 = 0 . Таким

образом,

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

уравнение

имеет один

простой k1 = -1 и

двукратный корень k 2,3

= 3 . Им соответствует фундаментальная система решений

y = ex ,

y

2

= e3x , y

3

= x e3x .

Общее решение дифференциального уравнения

1

 

 

 

 

 

 

 

 

 

 

 

 

по формуле (1.27) запишется как их линейная комбинация:

 

 

 

 

 

 

 

 

= C ex + (C

 

+ C

 

x)e3x .

 

 

 

 

 

 

 

Y

2

3

 

 

 

 

 

 

 

 

 

1

 

 

 

 

ПРИМЕР 1.23 Найти общее решение уравнения

 

 

 

 

 

 

 

yIV + y ′′′+ 4y ′′+ 4y′ = 0

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Решение. Этому уравнению соответствует характеристическое уравнение

k 4 + k3 + 4k2 + 4k = 0

или

k 2 (k 2 + 4)+ k(k2 + 4)= 0 ,

k(k +1)(k 2 + 4)= 0 .

Характеристическое уравнение имеет два действительных k1 = 0 , k 2 = −1 и два комплексных корня k3,4 = ±2i . (все они являются простыми). Таким образом,

общее решение данного дифференциального уравнения запишется в виде

Y = C1 + C2ex + C3 cos 2x + C4 sin 2x.

ПРИМЕР 1.24 Найти общее решение уравнения yIV +18y ′′+ 81y = 0.

 

Решение. Данному уравнению соответствует характеристическое уравнение

k 4

+18k 2 + 81 = 0 или (k 2 + 9)2 = 0,

имеющее двукратные комплексные корни

k1,2,3,4 = ±3i .

Им соответствуют

четыре частных решения y1 = cos 3x ,

y2

= sin 3x ,

y3 = x cos 3x , y4 = x sin 3x . Тогда общее решение однородного

уравнения запишется по формуле (1.27) как их линейная комбинация

Y = (C1 + C2 x)cos 3x + (C3 + C4 x)sin 3x.

1.12ЛИНЕЙНЫЕ НЕОДНОРОДНЫЕ ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ n - ГО ПОРЯДКА С ПОСТОЯННЫМИ

КОЭФФИЦИЕНТАМИ

Рассмотрим линейное неоднородное дифференциальное уравнение (ЛНДУ) n-го порядка

y(n ) + a1 y(n 1) + a 2 y(n 2) + ... + a n 1 y′ + a n y = f (x).

(1.38)

Теорема 1.7 (о структуре общего решения ЛНДУ). Общее решение Y ЛНДУ

(1.38) представляет собой сумму какого-нибудь частного решения Y* этого

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

Y соответствующего однородного уравнения

y(n ) + a1 y(n 1) + a

2 y(n 2) + ... + a n 1y′ + a n y = 0,

(1.39)

т.е.

 

 

 

 

 

Y =

 

+ Y*.

(1.40)

 

Y

Доказательство. Для простоты, доказательство теоремы проведем для ЛНДУ 2 - го порядка

y′′ + a1y+ a1y = f(x).

(1.41)

Докажем сначала, что функция (1.40) есть решение уравнения (1.41).

Подставляя сумму Y + Y* в уравнение (1.41) вместо y, имеем

(Y + Y* )+ a1 (Y + Y* )+ a 2 (Y + Y* )= f (x) или

 

 

 

 

+ a

 

 

 

 

*

+ a1Y

*

+ a

*

= f (x).

 

 

 

Y

+ a1 Y

2 Y

+ Y

 

 

2Y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Так как Y есть решение однородного уравнения

y′′ + a1y+ a 2 y = 0 ,

(1.42)

(1.43)

то выражение, стоящее в первых скобках, тождественно равно нулю. Так как Y* есть решение уравнения (1.41), то выражение, стоящее во вторых скобках, равно

f(x). Следовательно, равенство (1.42) является тождеством, функция Y + Y* является решением уравнения (1.41). Первая часть теоремы доказана.

Докажем теперь, что (1.40) есть общее решение уравнения (1.41), т.е.

докажем, что входящие в решение Y = Y + Y* = C1y1 + C2 y2 + Y* произвольные постоянные можно подобрать так, чтобы удовлетворялись начальные условия

y (x0 ) = y0

( ) = ′ (1.44)y x0 y0 .

Продифференцировав функцию

Y = C1y1 + C2 y2 + Y*

(1.45)

и подставив начальные условия (1.44) в функцию (1.45) и её производную, получаем систему уравнений

C y (x

0

) + C

2

y

2

(x

0

) = y

0

Y(x

0

),

 

 

1

1

 

 

 

 

 

 

 

 

(1.46)

 

 

y

(x

 

) + C

 

y

(x

 

) = y

(Y(x

 

 

))

C

0

2

0

0

,

 

1

1

 

 

 

2

 

 

0

 

 

 

 

где y0 = y(x0 ), y0 = y(x0 ), с неизвестными С1 и С2.

Определителем этой системы является определитель Вронского W(x0) для функций y1(x) и y2(x) в точке x=x0. Так как эти функции по условию линейно независимы (образуют фундаментальную систему решений), то W(x0) 0.

Следовательно система (1.46) имеет единственное решение: C1 = C10 и C2 = C02 .

Решение Y = Y* + C10 y1 (x) + C02 y2 (x) является частным решением уравнения (1.41). Теорема доказана.

1.13 МЕТОД ВАРИАЦИИ ПРОИЗВОЛЬНЫХ ПОСТОЯННЫХ

Рассмотрим универсальный метод нахождения частного решения Y*

неоднородного уравнения (1.38) на примере уравнения 2-го порядка (1.41).

Метод вариации произвольных постоянных (метод Лагранжа) состоит в

следующем: пусть

 

= C1 y1 (x) + C2 y2 (x) - общее

 

Y

решение однородного

уравнения (1.43).

 

Будем искать частное решение Y* неоднородного уравнения (1.41) в виде

Y* = C1 (x) y1 (x) + C2 (x) y2 (x).

(1.47)

Найдем производную

(Y )= C1(x) y1 (x) + C1 (x) y1(x) + C2 (x) y2 (x) + C2 (x) y2 (x).

Подберем функции C1(x) и C2(x) так, чтобы выполнялось равенство

 

 

 

 

 

 

 

 

C(x) y (x) + C

(x) y

2

(x) = 0 .

 

 

 

 

 

 

 

 

 

 

(1.48)

 

 

 

 

 

 

 

 

 

 

1

 

1

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Тогда (Y* )

= C

(x) y(x) + C

2

(x) y

(x),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

1

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(Y* )= C

(x) y (x)

+ C (x) y

(x) + C

 

(x) y

(x) + C

2

(x) y

 

(x).

 

 

 

 

 

1

 

 

 

1

 

 

 

 

1

 

 

 

1

 

 

2

 

 

2

 

 

 

2

 

 

 

Подставляя выражения для Y* , (Y* )

, (Y* )в уравнение (1.41), получим:

 

 

 

C

(x) y (x) + C (x)y (x) + C

(x)

y

(x)

+ c

2

(x) y

 

(x) +

 

 

 

 

 

 

 

1

 

 

 

1

 

1

 

 

 

1

 

 

 

 

2

 

 

 

 

2

 

 

 

 

2

 

 

 

 

 

 

 

 

 

+ a

1

[C

(x) y (x) + C

2

(x) y

(x)]+ a

2

[C (x)y (x) + C

2

(x)y

2

(x)]

= f (x)

 

 

 

 

 

1

1

 

 

 

 

 

2

 

 

 

y

 

 

1

 

1

 

 

 

y(x) + a

 

y

 

(x)]+

или C (x) [y′′(x) + a

1

y(x) + a

2

(x)]+ C

2

(x) [y′′ (x) + a

1

2

2

 

 

 

1

 

1

 

 

 

1

 

 

 

 

1

 

 

 

 

 

 

2

 

 

 

2

 

 

 

 

+ C

 

(x) y (x) + C

(x) y

(x) = f (x).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

1

 

 

2

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Так как y1(x) и y2(x) - решения уравнения (1.43), то выражения в квадратных скобках равны нулю, а потому

C1(x) y1(x) + C2(x) y2(x) = f (x).

(1.49)

Таким образом, функция Y* (1.47) будет частным решением уравнения (1.41), если функции C1(x) и C2(x) удовлетворяют системе уравнений, составленной из уравнений (1.48) и (1.49):

C1(x) y1

(x) + C2

(x) y2

(x) = 0,

 

C

(x) y

(x) + C

(x) y

(x) = f (x).

(1.50)

1

1

2

2

 

 

 

y1

(x) y2

(x)

 

0 , так как это определитель

 

 

Определитель системы

y

(x) y

(x)

 

 

1

2

 

 

 

Вронского для фундаментальной системы частных решений y1(x) и y2(x)

уравнения

(1.43).

Поэтому

система (1.50) имеет единственное решение:

C

(x) = α

(x) иC

(x) = α

2

(x).

Проинтегрировав функции α (x) иα

2

(x).,

1

1

2

 

 

1

 

находим С1(x) и С2(x), и по формуле (1.47) записываем частное решение неоднородного уравнения (1.41).

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]