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

Pyrkova_--_Metody_reshenia_differentsialnykh_urav

.pdf
Скачиваний:
35
Добавлен:
27.03.2016
Размер:
2.03 Mб
Скачать

y = y

o

+ y

ч

= C e4x +C

2

sin x + C

3

cos x

 

 

1

 

 

 

 

 

 

 

 

 

x

 

1

 

4 x

 

 

x sin x 4x cos x +

 

 

e

 

.p

4

8

 

 

 

 

 

 

 

 

 

 

 

1.5. Задачи для самостоятельного решения

Найти все действительные решения уравнений:

1.(1-01) y IV + 4 y′′ = 8e2 x +8x2 .

2.(1-02) y IV y′′′=sin x + e x .

3.

(1-03)

y IV

y = 2e x + 5e x sin x .

4.

(1-04)

y IV

+ y′′′− y′′− y′= 2x + 2 sin x .

5.

(1-11)

y′′′−3y′′+ 3y′− y =6e x + x .

6.

(1-12)

y′′′− y′′− y′+ y = x 2 + 4e x .

7.(1-13) y′′′+5y′′+9 y′+5y = 2ex +5x +4 .

8.(1-14) y′′′+3y′′+ y′−5y =10e x 5x .

9.(1-21) y′′′− y′′+9 y′−9 y = 60 cos 3x +4xex .

10.(1-22) y′′′+ y′′+4 y′+4 y =8 sin 2x +20xe x .

11.(1-23) y′′′−9 y′′+ y′−9 y =164 sin x +`18xe3x + 1136 e9 x .

12.(1-24) y′′′+4 y′′+ y′+4 y = 34 sin x +` 34x +13 e4 x .

4

13.(1-31) y′′′−3y′−2 y =18(1x)ex +100 cos 2x .

14.(1-32) y′′′−3y′′+4 y′−12 y = (126x)e3x +30 cos 3x .

15.(1-33) y′′′+3y′′−4 y =18(1x)e2 x +40 cos 2x .

16.(1-34) y′′′+ y′′+16 y′+16 y = (34x 4)ex +30 sin x .

17.(1-41) y IV +8y′′+16 y =8 cos 2x .

18.(1-42) y′′ + y = sin x sin 2x .

14

19. (1-43) y′′−2 y′+ y = ex sin x +4e x . 20. (1-44) y′′+ 2 y′+ y = xex +cos x .

21.(1-51) y′′′− y′′+ y′− y = 4 cos x (4x 14)ex .

22.(1-52) y′′′+ y′′− y′− y = 4 sin x +8(x +1)e x .

23.(1-53) y′′′− y′′− y′+ y = 4 cos x +4(x 4)ex .

24.(1-54) y′′′+ y′′+ y′+ y = 4 sin x +2(x +1)e x .

25.(1-61) y′′+2 y′−3y = 2e x sin 2 x .

26.(1-62) y′′′−2 y′′+ y′ = 6xe x +4 sh x .

27.(1-63) y′′− y′−2 y =10e2 x cos2 2x .

28.(1-64) y′′′+2 y′′+ y′ =16 ch x 6xex .

1.6.Ответы:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

1

 

1

 

1.

y = C +C

2

x +C

3

cos 2x +C

4

sin 2x +

x 2

x

 

 

+

e2x

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6

 

 

 

 

2.

y = C +C

2

x +C

3

x2 +C

4

e x

 

+ xex

+ 1 (sin x cos x)

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3.

y = C ex +C

2

ex

+C

3

sin x +C

4

cos x + 1 xex

ex sin x

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

4.

y = C +C

2

ex +C

3

xex +C

4

e x

+ 2x x 2 +

(sin x +cos x)

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

y = C e x +C

 

xex

 

 

 

 

x2e x

+ x3e x x 3

 

 

 

 

 

 

5.

2

+C

3

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6.

y = C e x +C

2

xex

+C

3

ex +

 

 

x2 + 2x + 4 + x2e x

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

7.

y = C ex +C

2

cos xe2x +C

3

sin xe2x

+ xex

+ x 1

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

8.

y = C e x +C

2

cos xe2x +C

3

sin xe2x

+ xex

+ x + 1

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

7

 

 

 

 

 

 

 

 

9.

y = C cos3x +C sin 3x +C ex

 

 

 

+

 

ex x(3cos3x +sin 3x)

 

 

 

 

1

 

 

 

 

2

 

 

 

 

 

 

3

 

 

 

 

5

 

50

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

15

10.

y = C ex +C

2

sin 2x +C

3

cos 2x +

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

9

x

 

2

x(2 sin 2x + cos 2x)

 

 

 

 

 

 

 

2x

 

e

 

 

 

 

 

 

 

 

 

 

 

5

 

 

 

 

 

 

 

 

5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

11.

y = C sin x +C

2

cos x +C e9 x + 9x cos x xsin x +

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

+(0.3x +0.13)e3x

 

11

 

e9x

 

 

 

 

 

 

 

72

738

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

12.

 

 

 

4x

 

 

 

 

 

 

 

 

 

 

 

x

 

1

 

4 x

y = C1e

 

 

+C2 sin x +C3

cos x

x sin x 4x cos x +

 

 

e

 

 

 

4

8

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

13.y = (C1 x +C2 )ex +C3e2x + (x3 2x 2 )ex cos 2x 7 sin 2x

14.y = C1 cos 2x +C2 sin 2x +C3e3x + (x x 2 )e3x +cos 3x sin 3x

15.y = (C1x +C2 )e2x +C3ex + (x3 2x2 )e2x 2 cos 2x sin 2x

16.y = C1ex +C2 cos 4x +C3 sin 4x + x 2 e x cos x +sin x

17.y = (C1 +C2 x)cos 2x + (C3 +C4 x)sin 2x 14 x 2 cos 2x

18.y = C1 cos x +C2 sin x + 4x sin x + 161 cos 3x

19. y =

(C +C

2

x)ex +

2x 2 e x +

3 sin x +4 cos x

e x

 

 

1

 

25

 

 

 

 

 

 

20.y = (C1 +C2 x)ex + 12 sin x + 16 x3ex

21.y = C1e x +C2 cos x +C3 sin x x(cos x +sin x)+(x 2)ex

22.y = C1e x +(C2 +C3 x)ex + x2ex +cos x sin x

23.y = C1ex +(C2 +C3 x)e x + 12 (x 2 6x)ex +cos x sin x

24.y = C1ex +C2 cos x +C3 sin x x(cos x +sin x)+ 14 (2x 1)e x

25.

y = C1e

x

+C2e

3x

 

1

 

x

1

 

1

 

x

 

 

+

 

xe

 

+

 

cos 2x

 

 

sin 2x e

 

 

 

4

 

20

10

 

 

 

 

 

 

 

 

 

 

 

 

 

26. y = C + (C

2

+C

3

x)e x + (x3

2x2 )ex +

1

ex

 

1

 

 

2

 

 

 

 

 

 

 

 

 

 

 

16

 

 

 

27.

y = C1e

x

 

2x

 

5

 

2x

 

3

 

1

 

 

2x

 

+C2e

 

+

 

xe

 

+

 

sin x

 

cos x e

 

 

 

3

 

2

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

28.

y = C1 +(C2 +C3 x)ex

+ 2ex + (x3 x2 )ex

 

 

17

§ 2. ЛИНЕЙНЫЕ СИСТЕМЫ

2.1. Основные понятия

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

 

dxk

n

 

 

 

 

= akj x j , k =

 

,

 

 

(2.1)

 

1, n

 

 

 

dt

 

 

 

j=1

 

 

 

где akj

= const .

 

 

(t)

 

 

 

 

 

x

 

 

 

 

x

12

(t)

Вводя в рассмотрение вектор-функцию x

=

 

и мат-

 

 

 

 

 

 

...

 

 

 

 

 

xn (t)

рицу A = (akj ), уравнения (2.1) можно представить в векторной форме

 

 

d x

 

= A x .

(2.2)

dt

 

 

2.2. Общее решение однородной системы

Фундаментальной системой решений однородной системы дифференциальных уравнений (2.1) называется совокуп-

ность n линейно независимых решений

 

 

 

 

 

 

x

(t)

 

 

 

 

x

(t)

 

 

 

 

x

(t)

 

=

x1121 (t)

=

x1222 (t)

, ... ,

=

x12nn (t)

(2.3)

x

 

 

 

,

x

2

 

 

 

x

n

 

 

 

1

 

...

 

 

 

...

 

 

 

...

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xn1 (t)

 

 

 

 

xn2 (t)

 

 

 

 

xnn (t)

 

этой системы.

Общее решение векторного уравнения (2.2) представляется в виде

18

 

x

(t)=C1 x1

(t)+ C2 x2

(t)+K+ Cn xn (t),

(2.4)

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

2.3. Метод Эйлера

(Метод сведения решения системы к задаче отыскания собственных значений и собственных векторов матрицы сис-

темы).

 

Чтобы найти решения (2.2):

A , решив

1) Вычислим собственные значения матрицы

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

 

det(A λE )= 0 .

(2.5)

Обозначим λ1 ,λ2 ,...,λn корни (2.5), вообще говоря, ком-

плексные.

Для собственного значения λ , отвечающий ему собствен-

ный вектор h

определяется условием

 

Ah = λh ,

h 0 .

(2.6)

2.1) Корни характеристического уравнения (2.5) действительные, простые. Тогда существует базис из собственных

векторов матрицы A :

Ahm = λm hm ,

hm 0 ,

m =

1, n

.

 

→ →

 

 

 

 

 

 

Вектор-функции

xm = hm eλmt ,

m =

1, n

 

являются реше-

ниями (2.2).

Общее решение векторного уравнения (2.2) есть их произвольная линейная комбинация ( Cm – постоянные)

n

r

 

x

= Cm hm eλmt .

(2.8)

m=1

2.2) Корни характеристического уравнения (2.5) невещественные, простые.

Еще раз напомним, что для комплексного числа z = x + iy , где x, y R , его действительной и мнимой частью называют-

19

ся соответственно Re z = x , Im z = y . Кроме того, имеет место формула Эйлера e(α+iβ )t = eαt (cos βt +i sin βt).

Если среди корней характеристического уравнения (2.5) есть невещественный корень λ =α +iβ , то комплексно со-

пряженное ему число λ =α iβ также будет корнем этого

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

 

 

векторного уравнения (2.2) x

(t)= h eλt

и x

(t )= h eλt . По-

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(t)=

x

(t)+ x (t)

 

комбинации решений x

и x , а именно,

x

и

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

(t ), или x

(t)= Re x

(t) и x

 

(t)= Im x(t).

 

(t )= x

(t )x

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

1

 

 

2

 

 

 

 

 

 

 

 

2i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Замечание. В более общем случае, когда собственный век-

тор для числа λ берется не сопряженным с вектором h , действительная и мнимая части соответствующего комплекснозначного решения системы будут линейными комбина-

циями

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

 

и

x

(t)= Re h eλt

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

λ .

x2

(t)= Im h eλt , найденных для собственного значения

 

 

 

 

 

 

 

 

 

Т.о.,

если

α ±iβ - простые корни характеристического

уравнения (2.5), то компонента общего решения системы, соответствующая этой паре комплексных корней, записывается в виде

20

=C

(t)+ C

(t)

= C

 

 

+ C

 

 

 

x

x

x

Re h eλt

 

Im h eλt , (2.9)

 

1

1

 

2 1

 

1

 

 

 

 

2

 

 

 

где C1 , C2 - произвольные постоянные, h - собственный вектор, отвечающий собственному значению λ =α +iβ .

2.3) Корни характеристического уравнения действительные кратные. В этом случае матрица A может не иметь n линейно независимых собственных векторов. Тогда для построения общего решения (2.2) используется следующее понятие.

Жордановой цепочкой матрицы A , соответствующей соб-

ственному значению λ , называется система векторов r r r

h1 , h2 ,..., hp такая, что

 

Ah1

= λh1 ,

h1 0,

 

r

r

r

 

Ah2

= λh2

+ h1 ,

 

r

r

r

 

Ah3

= λh3

+ h2 ,

(2.10)

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

 

r

r

r

 

Ahp = λhp + hp1 .

 

Вектор h1 - собственный, а h2 , h3 ,..., hp - присоединенные векторы.

Равенства (2.10) можно записать также в виде

r

r

 

r

 

(A λE)h

= 0,

h 0,

 

r1

r

1

 

 

(2.11)

(A λE)hk = hk 1 , k =

2, p

.

 

Каждой цепочке

h1 , h2 ,..., hp

соответствует p линейно не-

зависимых решений x1 , x2 ,..., x p

векторного уравнения (2.2):

21

xr

 

= eλt h ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

λt t

 

 

 

r

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

= e

 

 

 

 

 

 

h

+ h

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1!

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

1

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

 

2

 

 

r

 

 

 

 

t

r

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

λt t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

= e

 

 

 

 

 

 

 

 

h

+

 

 

 

 

 

h

 

+ h

,

 

 

 

 

 

 

 

 

 

 

(2.12)

3

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

2!

 

 

1

 

 

1!

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

t

p1

 

 

 

 

r

 

 

 

t

p2 r

 

 

t

r

 

r

 

 

 

λt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

= e

 

 

 

 

 

 

 

 

 

 

 

 

h

 

+

 

 

 

 

h

 

+...+

 

 

 

h

 

+ h

.

p

 

 

(p 1)!

 

 

(p 2)!

2

 

 

 

p1

 

 

 

 

 

 

1

 

 

 

 

1!

 

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Замечание. Приведем правило запоминания формул

(2.12). Собственному вектору

h1 , соответствует решение

xr

= eλt h . Если везде отбросить

eλt , то каждая строка правой

1

1

 

части (2.12) получается интегрированием по t предыдущей строки, причем постоянную интегрирования надо взять равной следующему по порядку вектору серии.

Для кратного собственного значения λ может существовать несколько жордановых цепочек, содержащих линейно независимые собственные векторы матрицы A .

Компонента общего решения системы, соответствующая

действительному

собственному значению λ кратности p,

имеет вид

 

 

 

 

 

 

 

k

Cl(r )xrl(r )(t),

 

 

xr(t)= eλt r

 

 

 

r

l =1

 

 

kr

 

(r )

(r )

 

(r )

- произвольные постоянные,

= p .

где C1

, C2

,K, Ckr

r

Известно, что для любой квадратной матрицы A существует базис, составленный из ее жордановых цепочек, поэтому произвольная линейная комбинация решений вида (2.12) дает общее решение векторного уравнения (2.2).

22

2.4. Общее решение неоднородной системы

Решение неоднородной системы

 

 

dxr

r

r

 

 

 

= Ax

+ f ,

(2.13)

 

dt

можно найти методом вариации постоянных, если известно общее решение однородной системы (2.1) с той же матрицей A = (akj ). Для этого в формуле общего решения (2.4) одно-

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

Cm , m =

1, n

, на неизвестные функции C m (t ):

 

 

 

xr(t)= n

Cm (t )xrm (t ).

 

 

 

 

 

 

(2.14)

 

 

 

 

 

 

m=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Полученные

 

 

выражения

для

 

dxr(t)

 

n dCm

(t ) r

n

dxrm

(t)

 

 

 

 

 

 

 

=

 

 

 

 

 

 

xm

(t )+Cm (t )

 

 

 

 

 

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

 

dt

 

 

dt

 

dt

 

 

 

 

 

m=1

 

 

 

 

m=1

 

 

 

 

 

 

нородную

 

 

 

 

 

 

 

систему

 

 

 

 

(2.13).

Т.к.

n

Cm

(t )

dxrm

(t )

= n

Cm (t)Axrm (t ),

то

 

получаем

систему для

 

 

 

 

 

m=1

 

 

 

 

dt

 

 

 

 

m=1

 

 

 

 

 

 

 

 

 

 

определения

 

 

dCm (t)

 

, m =

 

:

 

 

 

 

 

 

 

 

 

 

1, n

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

n

dCm (t )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xrm (t )=

fv .

 

 

 

 

 

 

(2.15)

 

 

 

 

 

 

 

 

 

 

 

m=1

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Неизвестные функции C m (t ),

m =

 

, находим проинтег-

 

 

1, n

рировав, полученные при решении системы (2.15) функции

dCm (t), m =1, n . dt

Заметим, что если при нахождении функций Cm (t ) запи-

сывать всю совокупность первообразных, т.е. сохранять в записи выражений для Cm (t ) возникающие при интегриро-

вании произвольные постоянные, то (2.14) будет общим решением неоднородной системы. Частное решение неодно-

23

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