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

08 ЛДУ 2

.pdf
Скачиваний:
12
Добавлен:
12.02.2015
Размер:
439.98 Кб
Скачать

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

Уравнение вида

y n a1 y n 1 a2 y n 2 an 1 y an y 0 ,

(1)

в котором a1 , a2 , , an – действительные или комплексные числа, называется линейным

однородным дифференциальным уравнением порядка n с постоянными коэффициентами.

В дальнейшем будем предполагать, что a1 , a2 , , an – действительные числа, и

находить действительные решения уравнения (1).

В соответствии с методом Эйлера решение уравнения (1) будем искать в виде

 

 

 

 

y e x ,

 

 

 

(2)

где – некоторое число. Последовательно дифференцируя функцию (2), находим

 

 

y e x , y 2e x , , y n ne x .

(3)

Подставляя функцию

y e x

и ее производные (3) в уравнение (1), получим равенство

ne x a n 1e x a

e x

a

e x 0 ,

из которого после сокращения на множитель

1

n 1

n

 

 

 

 

 

 

 

e x 0 вытекает соотношение

 

 

 

 

 

 

 

 

 

n a n 1 a

n 1

a

n

0 .

(4)

 

 

 

1

 

 

 

 

Уравнение (4) называется характеристическим уравнением для линейного однородного дифференциального уравнения (1).

Из курса линейной алгебры известно, что число различных корней j характеристического уравнения (4) не больше n , но сумма кратностей всех корней равна n . В

случае

действительных коэффициентов a1 , a2 , , an каждому комплексному корню

j j

i j соответствует сопряженный корень

 

j j

i j той же кратности, что и

 

корень j .

 

 

Каждому действительному корню кратности r

соответствуют r

линейно не-

зависимых решений уравнения (1):

 

 

 

e x , xe x , x2e x , , xr 1e x ,

 

 

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

соответствуют s

пар линейно

независимых решений:

 

 

e x cos x, xe x cos x, , xs 1e x cos x , e x sin x, xe x sin x, , xs 1e x sin x .

1

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

корней характеристического уравнения (4) составляют фундаментальную систему решений y1 x , y2 x , , yn x уравнения (1). Тогда общее решение уравнения (1) запишется в виде

y x C1 y1 x C2 y2 x Cn yn x ,

(5)

где C1 , C2 , , Cn – произвольные постоянные.

Пример 1. Найти общее решение уравнения y y 6y 0 .

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

 

3 2 6 0

 

или 2 3 0 .

 

Его корни 1 0, 2

2, 3

3 – действительные числа кратности r1 r2

r3 1. Им

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

y e0 x ,

y

2

e 2 x , y

3

e3x , поэтому в соответствии с формулой (5)

 

1

 

 

 

 

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

y C1 C2e 2 x C3e3x . ▲

Пример 2. Найти общее решение уравнения

y 5 5y 4 8y 8y 7 y 13y 0 .

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

 

5 5 4 8 3 8 2

7 13 0

или

 

 

1 2 1 2 6 13 0

 

имеет

корни 1 1,

2,3

i ,

4,5 3 2i .

Этим

корням соответствуют

решения

y e

x , y

2

cos x ,

y

3

sin x ,

y

4

e3x cos 2x ,

y

5

e3x

sin 2x . Тогда общее решение ис-

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ходного уравнения имеет вид

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y C e x C

2

cos x C sin x e3x

C

4

cos 2x C sin 2x . ▲

 

 

 

 

 

1

 

 

 

 

 

3

 

 

 

 

 

5

 

 

Пример 3.

Пусть

1

2 3

2 3i ,

 

4

5 6 2 3i , 7

8 4 ,

9 1 – корни характеристического уравнения. Записать формулу общего решения

соответствующего дифференциального уравнения.

 

 

 

 

Корням

2 3i

кратности

s 3

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

решения

e 2 x cos 3x ,

xe 2 x cos 3x ,

x2e 2 x cos 3x , e 2 x

sin 3x ,

xe 2 x sin 3x ,

x2e 2 x sin 3x . Корню 4 кратности

r 2 соответствуют решения

e4 x ,

 

xe 4 x

и корню 1 соответствует решение e x .

Тогда общее решение исходного уравнения имеет вид

 

 

 

 

y C C

2

x C

x2 e 2 x cos 3x C

4

C

x C

6

x2

e 2 x sin 3x C

7

C x e4 x C

e x . ▲

1

3

 

 

5

 

 

 

8

9

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

В частном случае, линейное однородное дифференциальное уравнение 2-го по-

рядка

 

 

 

y a1 y a2 y 0

(6)

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

 

 

 

2 a a

2

0 .

(7)

1

 

 

В зависимости от корней 1 и 2 характеристического уравнения (7) получаем общее решение уравнения (6) в виде:

1.y C1e 1x C2 e 2 x , если корни действительны и 1 2 ;

2.y (C1 C2 x)e 1x , если корни действительны и 1 2 , другими словами имеется один корень кратности 2;

3.y (C1 cos x C2 sin x)e x , если корни 1,2 i – комплексные числа;

где C1 , C2 – произвольные постоянные.

Пример 4. Найти решение задачи Коши

y

 

5y

 

4y 0

,

 

1.

 

 

 

y(0) y (0)

 

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

2 5 4 0 имеет

корни 1 1,

2 4 .

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

y C1ex C2e4 x .

Для определения произвольных постоянных C1 и C2 найдем производную

y C1ex 4C2e4 x

и используем начальные условия:

y(0)

C

C

,

 

 

C

C

 

1,

 

C 1,

 

1

2

 

 

 

1

 

2

1,

 

1

y (0)

C1

4C2 ,

 

C1 4C2

 

C2 0.

Подставляя найденные постоянные C1 и C2 в общее решение, находим частное реше-

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

y ex . ▲

 

 

 

 

 

 

 

Пример 5. Найти решение задачи Коши y

 

 

 

y 0 ,

y(0) 2 ,

 

1.

 

2y

y (0)

Характеристическое уравнение 2 2 1 0

имеет корень 1 кратности

r 2 . Таким образом, общее решение дифференциального уравнения имеет вид

3

y (C1 C2 x)e x .

Для определения произвольных постоянных C1 и C2 найдем производную

y (C2 C1 C2 x)e x

и используем начальные условия:

y(0)

C ,

 

 

C 2,

 

C 2,

 

1

 

 

1

1,

 

1

y (0)

C2 C1 ,

 

C2

C1

 

C2 1.

Подставляя найденные постоянные C1 и C2 в общее решение, находим частное реше-

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

y (2 x)e x . ▲

 

 

 

 

 

Пример 6. Найти решение задачи Коши y

 

 

5y 0 ,

y(0) 1,

 

5 .

 

2y

y (0)

Характеристическое уравнение 2 2 5 0

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

женных корня 1,2 1 2i . Следовательно, общее решение дифференциального урав-

нения имеет вид

y (C1 cos 2x C2 sin 2x)ex .

Для определения произвольных постоянных C1 и C2 найдем производную

y (C1 2C2 )cos 2x (C2 2C1 )sin 2x ex

и используем начальные условия:

y(0)

C ,

 

 

C 1,

 

C 1,

 

1

2C2 ,

 

 

1

 

1

y (0)

C1

 

C1

2C2 5,

 

C2 2.

Подставляя найденные постоянные C1 и C2 в общее решение, находим частное реше-

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

y (cos 2x 2sin 2x)ex . ▲

4

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

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

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

(8)

Как следует из теоремы о структуре общего решения неоднородного уравнения, общее решение уравнения (8) имеет вид

y y * C1 y1 C2 y2 Cn yn ,

где y * – частное решение уравнения (8), а y1, y2 , , yn – фундаментальная система решений соответствующего (8) однородного дифференциального уравнения

y n a1 y n 1 a2 y n 2 an 1 y an y 0 ,

(9)

способ нахождения которой уже рассмотрен. Поэтому основной задачей при нахождении общего решения является задача отыскания частного решения y * уравнения (8).

Рассмотрим частные случаи для правых частей f x специального вида, когда отыскать частное решение y * уравнения (8) наиболее просто, а именно, рассмотрим

метод подбора частного решения (или метод неопределенных коэффициентов).

Пусть

 

 

x A ,

 

f x e x A xn A xn 1

A

(10)

0

1

n 1

n

 

где A0 0 .

Если не является корнем характеристического уравнения дифференциального уравнения (9), то частное решение y * дифференциального уравнения (8) с правой частью (10) ищется в виде

y* e x B xn B xn 1

B

x B ,

(11)

0

1

n 1

n

 

где B0 , B1, , Bn – неизвестные постоянные коэффициенты.

 

Если – корень кратности r

характеристического уравнения дифференциаль-

ного уравнения (9) (резонансный случай), то частное решение

y * дифференциального

уравнения (8) с правой частью (10) ищется в виде

y* xr e x B0 xn B1 xn 1 Bn 1 x Bn .

5

В частном случае формулы (10), при 0 , то есть в случае f x A0 xn A1 xn 1 An 1 x An

частное решение y * следует искать в виде

y* B0 xn B1xn 1 Bn 1x Bn ,

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

y* xr B0 xn B1 xn 1 Bn 1 x Bn ,

если число 0 является корнем кратности r характеристического уравнения.

Подставляя соответствующее представление частного решения y * в уравнение

(8), сокращая полученное на e x и приравнивая коэффициенты при одинаковых степенях x в левой и правой частях сформированного равенства, находим систему линейных

алгебраических

уравнений

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

неопределенных

коэффициентов

B0

, B1, , Bn 1

, Bn .

Решение

этой

системы

определяет

вид

коэффициентов

B0

, B1, , Bn 1

, Bn

и частное решение уравнения (8).

 

 

 

 

 

Пример 7. Найти общее решение уравнения y 2y 8x 2 .

 

 

 

Характеристическое уравнение 2 2 0 имеет корни 1 0 , 2

2 . Тогда

общее

решение

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

однородного

уравнения

y 2y 0

имеет вид

y0

C

C

e 2 x . Поскольку правая часть

f x 8x 2 есть многочлен первой степени, а

 

1

2

 

 

 

 

 

 

 

 

 

 

число 0 является корнем кратности

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

ное решение y * ищем в виде

y* x B0 x B1 B0 x2 B1 x .

Отсюда

 

 

 

 

2B0 x B1 ,

 

2B0 .

y*

y*

Подставляя

 

,

 

 

 

 

y *

y * в исходное уравнение, получаем равенство

2B0 4B0 x 2B1

8x 2 или

 

4B0 x (2B0 2B1 ) 8x 2

Отсюда

 

 

 

 

 

 

 

 

 

 

 

4B 8,

 

B 2,

 

 

 

x1

 

 

 

 

x0

 

0

 

0

 

 

 

2B0 2B1 2,

B1 3.

Таким образом, искомое частное

решение исходного уравнения имеет вид

y* 2x2 3x , а общее решение определяется выражением

y y0 y* C C

e 2 x 2x2 3x . ▲

1

2

 

 

6

 

Пример 8. Найти общее решение уравнения y 2y y (4x 12)e x .

Так как характеристическое уравнение 2 2 1 0

имеет один корень 1

кратности r 2 , то общим решением однородного уравнения

y 2y y 0 является

функция y0 (C C

x)ex . Правая часть исходного

уравнения есть функция

1

2

 

 

 

f x (4x 12)e x

вида (10), где 1. Поскольку число

1 не является корнем

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

 

 

 

 

 

 

 

 

y* e x B x B .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

1

 

 

 

 

 

 

 

Отсюда

 

 

B0 x B1

e

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

x

B0

e

x

(B0 B1 B0 x) ,

 

e

x

(B0 x

B1 2B0 ) .

y* e

 

 

 

y*

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

0

, получаем равенст-

 

и y *

в исходное уравнение и сокращая на e

 

 

во

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(B0 x B1 2B0 ) 2(B0

B1

B0 x) (B0 x B1 ) 4x 12

или

4B0 x (4B1 4B0 ) 4x 12 ,

откуда получаем систему линейных алгебраических уравнений

 

 

 

 

 

 

 

 

x1

 

 

4B 4,

 

B 1,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x0

 

 

 

0

 

0

2.

 

 

 

 

 

 

 

 

 

 

4B1 4B0 12,

B1

 

 

 

 

 

Итак, y* e x x 2 , и, таким образом, общее решение исходного уравнения имеет вид

y y0 y* (C1 C2 x)ex x 2 e x . ▲

Рассмотрим более общий случай, пусть правая часть дифференциального уравнения (8) имеет вид

f x e x ( A xn A )cos x (B xm B

)sin x ,

(12)

0

n

0

m

 

 

где A0 0 , B0 0 .

Если i не являются корнями характеристического уравнения дифференциального уравнения (9), то частное решение y * дифференциального уравнения (8) с правой частью (12) ищется в виде

y* e x (C0 xk Ck )cos x (D0 xk Dk )sin x ,

где k max( n, m) , C0 , C1 , , Ck , D0 , D1 , , Dk – неизвестные постоянные коэффициенты.

7

Если i – корни кратности r характеристического уравнения дифференциального уравнения (9) (резонансный случай), то частное решение y * дифференциального уравнения (8) с правой частью (12) ищется в виде

y* xr e x (C0 xk Ck )cos x (D0 xk Dk )sin x .

Подставляя соответствующее представление частного решения y * в уравнение

(8), сокращая полученное на e x и приравнивая коэффициенты при одинаковых функциях в левой и правой частях сформированного равенства, находим систему линейных алгебраических уравнений относительно неопределенных коэффициентов C0 , C1 , , Ck , D0 , D1 , , Dk . Решение этой системы определяет вид коэффициентов

C0 , C1 , , Ck , D0 , D1 , , Dk и частное решение уравнения (8).

Пример 9. Найти общее решение уравнения y 3y 2y 8cos 2x 4sin 2x .

Так как корнями характеристического уравнения 2 3 2 0 являются числа 1 1 и 2 2 , то общим решением однородного уравнения y 3y 2y 0 яв-

ляется y0 C e x C

e 2 x . По виду правой части находим числа

i 0 2i , которые

1

2

 

 

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

y* C cos 2x Dsin 2x .

Отсюда

 

2C sin 2x 2Dcos 2x ,

 

4C cos 2x 4Dsin 2x .

y*

y*

Подставляя

 

и

 

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

y *

y *

( 4C cos 2x 4Dsin 2x) 3( 2C sin 2x 2Dcos 2x) 2(C cos 2x Dsin 2x) 8cos 2x 4sin 2x

или

( 2C 6D)cos 2x ( 6C 2D)sin 2x 8cos 2x 4sin 2x .

Отсюда, приравнивая коэффициенты при cos 2x и sin 2x в обеих частях равенства соответственно, получаем систему уравнений

2C 6D 8,

C 1,

 

 

6C 2D 4,

D 1.

Таким образом, частное решение исходного дифференциального уравнения есть y* cos 2x sin 2x , а его общее решение имеет вид

y y0 y* C1e x C2e-2 x cos 2x sin 2x . ▲

8

Пример 10. Для уравнения

y y 2x 3 e2 x x2 1 x3 cos x

указать вид ча-

стного решения.

 

 

 

 

 

 

 

 

 

 

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

уравнение

3 2

1 0 имеем

корни

0 ,

 

 

 

 

 

 

 

 

 

1

 

2,3 i . Правая часть дифференциального уравнения состоит из трех слагаемых:

 

f x 2x 3,

f

2

x e2 x x2

1 ,

f

3

x x3 cos x .

 

 

1

 

 

 

 

 

 

 

 

По принципу суперпозиции частное решение ищется в виде

 

 

 

 

y* y y

y

,

 

 

 

 

 

 

 

1

2

3

 

 

 

 

 

Где

 

 

 

y e2x Cx2 Dx E ,

 

 

y x Ax B ,

 

 

1

 

 

 

2

 

 

 

 

 

 

y x Fx3 Gx2 Kx L cos x Mx3 Nx2

Px Q sin x . ▲

 

3

 

 

 

 

 

 

 

 

 

 

9

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

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

y

n

a1 x y

n 1

a2 x y

n 2

 

an x y f x ,

(13)

 

 

 

an 1 x y

если известна фундаментальная система решений соответствующего однородного уравнения

y

n

a1

x y

n 1

a2 x y

n 2

 

an x y 0 .

(14)

 

 

 

an 1 x y

Пусть y1, y2 , , yn

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

(14). Тогда по теореме о структуре общего решения однородного уравнения его общее решение имеет вид

y x C1 y1 x C2 y2 x Cn yn x ,

(15)

где C1, C2 , , Cn – произвольные постоянные.

В соответствии с методом вариации произвольных постоянных будем искать общее решение неоднородного уравнения (13) в виде решения (15), но, считая, что C1, C2 , , Cn есть некоторые дифференцируемые функции переменной x , то есть

n

x yi

 

 

y C1 x y1 x C2 x y2 x Cn x yn x Ci

x ,

(16)

i 1

Выберем неизвестные функции Ci x , i 1, n таким образом, чтобы функция (16) удовлетворяла неоднородному уравнению (13). Для этого потребуем, чтобы производные Ci x , i 1, n этих функций являлись решением системы линейных алгебраических уравнений

 

x y1 x

 

x y2 x

 

 

x yn x

0,

 

 

C1

C2

... Cn

 

 

 

 

 

x

 

 

 

x

0,

 

 

 

 

 

C1

x y1 x

C2

x y2

... Cn

x yn

 

 

 

(17)

 

 

n 2

 

n 2

 

 

 

n 2

 

 

 

 

x y1

x y2

 

 

x yn

 

x 0,

 

 

C1

x C2

 

x ... Cn

 

 

 

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

 

 

 

 

 

 

x y n 1 x f x .

 

1

1

2

2

 

 

n

 

n

 

 

 

 

 

10