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

Линейные дифференциальные уравнения и системы линейных уравнений. Учебное пособие

.pdf
Скачиваний:
0
Добавлен:
12.08.2026
Размер:
1 Мб
Скачать

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

x 4x y 5t 1,

6.y x 2y t 1.

Найти все действительные решения систем уравнений, построив частные решения по виду вектора F :

x 2 y 3t,

7.y 2x 4.

8.x y t2 6t 1,y x 3t2 3t 1.

9.x 3x y e2t ,y x 5y et .

10.x 2x 4 y cost,

y x 2 y sin t.

11.x x y 2e t ,

y 6x 4 y 4e t .

x x y t2 ,

12.y y z 2t,z z t.

Ответы

1.

x

C

cost

e2t C

2

 

 

sin t

e2t

 

 

 

1

 

 

 

 

 

 

 

 

 

y

 

cost

sin t

 

 

 

cost sin t

 

 

e2t t

cost sin t ln

 

cost

 

 

 

cost sin t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

2t

t sin cost ln

 

 

 

 

 

 

 

 

2e

 

cost

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

111

Нежельская Л. А. Линейные дифференциальныеуравнения

 

 

 

 

 

 

 

 

 

2

t2

t 7

2t e3t

 

 

 

 

 

 

 

 

 

 

 

2.

x

C

1 e3t C

 

t

1

 

 

35

 

 

 

 

 

 

 

2

0 e3t

 

 

 

 

 

 

.

 

 

1

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

y

 

1

 

 

1

 

1

 

 

t

3

te

3t

 

 

 

 

 

 

 

 

 

 

 

35

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x cost sin t sin t ln tgt

3.y C1 sin t C2 cost cost ln tgt .

4.

x

1

2t

1

3t

 

t2

 

 

C1 e

 

C2

e

 

 

 

.

 

y

1

 

 

2

 

t

2

x cost sin t t cost sin t ln sin t

5.y C1 sin t C2 cos t t sin t cost ln sin t 1 .

6.

x

C

1 e3t

C

t

1

 

 

0 e3t

 

t .

 

 

1

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

1

 

 

 

 

1

 

1

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5

 

 

 

7.

x

C

cos 2t

C

 

sin 2t

 

 

 

4

 

 

 

2

 

 

.

 

 

1

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

y

 

sin 2t

 

 

 

 

cos 2t

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

x

 

cost

 

 

 

 

sin t

 

 

 

 

3t

2

 

 

 

 

 

8.

 

C2

 

 

 

 

 

t 1

 

C1

 

 

 

 

 

 

 

 

 

 

t2

2

.

 

y

 

sin t

 

 

 

 

cost

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 e 4t C

 

t

 

 

 

 

 

 

1

et

 

7

e2t

 

 

 

 

 

 

 

 

 

 

 

 

 

9.

x

C

 

2

1

 

 

1 e 4t 25

 

 

 

36

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

4

 

t

 

 

1

 

2t

 

y

 

 

1

 

 

 

 

 

1

 

0

e

 

e

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

25

 

36

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

10.

x

C1

 

4

 

0

 

 

C2

 

2

3sin t 2cost

 

 

 

 

t

 

 

 

 

 

 

 

 

2sin t

 

 

 

.

 

 

 

y

 

 

2

 

1

 

 

 

 

1

 

 

 

 

 

 

 

.

112

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

11.

x

C

1

et C

1

 

e t

 

 

 

 

 

 

 

 

 

 

2

 

e2t

 

.

 

 

 

 

 

 

 

1

 

 

 

 

 

2e t

 

 

 

 

 

 

y

 

 

 

2

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

1

 

 

 

0

0

 

 

1

 

0

 

 

12.

 

C1

 

 

2

 

 

 

 

 

 

t

 

 

 

 

t

 

y

t

 

0

t 2

0 e

 

C2 t 0

 

1 e

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

0

 

 

 

0

2

 

 

0

 

0

 

 

 

 

1

 

 

 

t2

3t 3

 

 

 

 

 

 

 

 

C

0 e t

 

 

t

 

 

.

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

1

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3.6.Общие методы интегрирования систем

3.6.1.Метод исключения

Метод интегрирования нормальной системы n дифференциальных уравнений

dxi

a

x

 

f

 

t , i 1, n ,

(3.1)

 

n

 

 

 

 

 

 

 

 

dt

j 1 ij

 

j

 

i

 

 

 

 

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

В качестве примера применения метода рассмотрим решение системы дифференциальных уравнений в научной задаче, связанной с изучением MAP-потока (Markovian Arrival Process) – дважды стоха-

стического потока событий, являющегося математической моделью реальных информационных потоков запросов (сообщений) в телекоммуникационных сетях и системах связи [4, стр. 29–33].

Пример 3.10. Рассматривается поток событий, сопровождающий случайный процесс которого t является кусочно-

постоянным с двумя состояниями – S1 и S2 : если t i , то

113

Нежельская Л. А. Линейные дифференциальныеуравнения

имеет место i -е состояние ( Si ) процесса t , i 1,2 ,

1

2 0 .

Длительность пребывания процесса

t в состоянии

Si

опреде-

ляется случайной

величиной с

функцией распределения

F t 1 e it , t 0 ,

i 1,2 . В момент окончания состояния S

i

i

 

 

 

 

процесса t возможны следующие ситуации, каждая из которых протекает мгновенно: 1) наступает событие потока, и процесс t

переходит из состояния Si

в состояние S j ; совместная вероят-

ность описанной ситуации

P1 j

 

i , i, j 1,2 ; 2) не наступает

 

событие потока, и процесс t переходит из состояния Si

в со-

стояние S j ; совместная вероятность этой ситуации P0 j

 

i ,

 

i, j 1,2 , i j . При этом

 

 

 

 

 

P0 j i P1 j i P1 i i 1, i, j 1,2 , i j .

Вариант возможного поведения случайного процесса t и MAP-потока событий представлен на рис. 1, где S1 и S2 – состояния процесса t ; t1,t2 ,...,tk ,... – моменты наступления событий в потоке.

S1

S2 ...

...

...

)

P1( 1| 1)

 

)

 

 

 

 

 

 

2

 

 

 

 

 

2

 

 

 

 

 

 

|

 

 

 

 

 

|

 

 

 

 

 

1

 

 

1

 

 

1

 

1

 

 

 

 

(

 

 

 

 

(

 

 

 

 

 

 

 

 

)

 

P1( 2| 2)

 

 

)

 

 

 

 

 

 

 

2

 

 

 

2

 

 

 

 

0

 

 

|

 

 

1

 

|

 

 

 

P

 

 

 

P

 

...

 

 

 

(

 

 

(

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

1

 

 

 

 

 

 

 

P

 

 

 

P

 

 

 

 

 

 

 

 

 

 

 

 

Процесс (t)

t

 

 

 

 

 

 

 

 

 

 

...

t

 

t1

t2

t3 t4

t5

Рис. 1. Формирование MAP-потока событий

114

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

Одной из характеристик случайного процесса t является априорная вероятность состояния S1 или S2 , в котором может находиться процесс t в произвольный момент времени t .

Необходимо получить систему дифференциальных уравнений для априорных вероятностей состояний процесса t для MAP-

потока событий, проинтегрировать её и выписать явный вид априорных вероятностей.

Решение. Обозначим

i t

 

t0 – априорная вероятность того,

 

что в момент времени t

значение процесса t i , i 1,2 , при

условии, что функционирование MAP-потока

(процесса t )

началось в момент времени t0 . Таким образом,

в предлагаемой

задаче время t

независимая переменная, i t

 

t0 – искомая

 

функция.

 

i t t

 

t0 того,

 

 

 

 

Найдём вероятность

 

что в момент времени

 

t t процесс t

принимает значение

i , т.е. имеет место со-

стояние Si , i 1,2 ,

случайного процесса

t . Опишем возмож-

ные ситуации на полуинтервале t,t t , связанные с поведением процесса t :

1)

в момент времени t

значение процесса t i с вероят-

ностью

i

t

 

t0 , и за промежуток времени длительности

t не

 

закончится

 

 

состояние

Si

процесса t с вероятностью

e i t 1 i t o t ,

и с вероятностью единица процесс

t

останется в состоянии Si ; вероятность этой ситуации есть

 

i t

 

t0 1 i t o t i t

 

t0 i i t

 

t0 t o t ;

 

 

 

 

 

115

Нежельская Л. А. Линейные дифференциальныеуравнения

2) в момент времени t значение процесса t i с вероят-

ностью i t

 

t0 , и за

промежуток времени длительности t за-

 

кончится

состояние

Si

процесса t

с

вероятностью

1 e i t

i t o t ,

наступит событие

потока, и процесс

t

 

перейдёт из состояния Si в состояние

Si

с вероятностью

P1 i

 

i ; вероятность описанной ситуации равна

 

 

 

i t t0 i t o t P1 i i i P1 i i i t t0 t o t ;

3)в момент времени t значение процесса t j с вероят-

ностью j t t0 , j 1,2 , j i , и за промежуток времени длитель-

ности t закончится состояние S j процесса t с вероятностью

1 e j t j t o t , наступит событие

потока, и процесс

t перейдёт из состояния S j в состояние

Si с вероятностью

P1 i j ; вероятность рассмотренной ситуации равна

j t t0 j t o t P1 i j j P1 i j j t t0 t o t ;

4)в момент времени t значение процесса t j с вероят-

ностью j t

 

t0 ,

j 1,2 ,

j i , и за промежуток времени длитель-

 

ности t закончится состояние S j процесса t с вероятностью1 e j t j t o t , не наступит событие потока, и процесс

t

перейдёт из состояния S j в состояние Si с вероятностью

P0 i

 

j ; вероятность ситуации есть

 

116

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

j t t0 j t o t P0 i j j P0 i j j t t0 t o t .

Других возможностей не имеется.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Относительно вероятностей

1 t

 

t0

и

 

2 t

 

 

 

t0 в момент вре-

 

 

 

мени t t имеем систему равенств:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 t t

 

t0 1 t

 

 

t0 1 1 t

 

 

t0 t 1P1 1

 

1 1 t

 

t0 t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

0

t o t ,

 

 

 

 

 

 

 

 

 

 

2 P1 1

 

 

2 P0 1

 

 

 

2

2

 

t

 

 

 

 

 

 

 

 

2 t t

 

t0 2 t

 

t0 2 2 t

 

 

t0 t 2 P1 2

 

2 2 t

 

t0 t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 t

 

 

 

 

0

 

t o t .

 

 

 

 

 

 

 

 

 

 

 

 

 

1 P1 2

 

1 P0 2

 

 

1

 

t

 

 

 

 

 

 

 

 

 

В результате несложных преобразований находим

 

 

 

 

 

 

 

 

 

 

 

 

 

1 t t

 

t0 1 t

 

 

 

t0

 

 

 

 

 

 

 

 

 

 

t o t ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 t

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 1

P1 1

 

1

 

t

 

 

2 1 P1 2

 

2 2

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 t t

 

t0 2 t

 

t0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P

 

 

 

 

 

 

 

 

 

 

 

t

 

t

0

 

 

 

P

 

 

 

 

 

 

 

t

 

 

t

0

t o t .

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

2

 

1

 

2

 

 

2

 

 

 

2

 

 

 

 

1

 

1

1

 

1

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Разделив оба соотношения на t и выполнив в них предельный переход при t 0 , получим однородную систему линейных дифференциальных уравнений относительно априорных вероятно-

стей i t t0 состояний Si , i 1,2 , процесса t :

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 t

 

t

 

1 1 P1 1

 

1

 

1 t

 

t

 

 

2 1 P1 2

 

2

 

2

t

 

t

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 t

t

 

 

1 1

P1 1

1

1 t

t

 

2 1

P1 2

2

 

2 t

t

 

,

 

1 t

 

t0

 

2 t

 

t0 1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Обозначим

 

 

 

 

 

, b

 

1 P

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

1 P

|

 

2

2

|

2

.

 

 

 

 

 

(3.27)

 

 

 

 

 

 

 

 

 

1

 

1

1

1

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

117

Нежельская Л. А. Линейные дифференциальныеуравнения

Тогда система запишется в виде

 

 

 

 

 

 

 

 

 

 

 

 

1

t

 

 

 

t0 a 1 t

 

t0 b 2 t

 

 

t0 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

t0 a 1 t

 

t0 b 2 t

 

t0 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

t0 2 t

 

t0 1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 t | t0

Из

первого

 

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

исключим

 

 

1

1 t

 

t0

a

1 t

 

 

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

 

 

 

b

 

 

 

 

 

 

 

 

b

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

1

t

 

t0

1 t

 

t0 a 1 t

 

t0 1 t

 

t0 a 1 t

 

t0 ,

 

 

 

 

 

 

 

 

 

 

 

b

 

 

 

 

 

 

 

 

b

 

 

 

 

 

 

 

 

 

 

 

 

0 или 1 t

 

 

 

a b 1 t

 

 

 

0 .

1

1 t

 

t

0

 

 

 

a

1 t

 

t

0

 

t

0

 

t

0

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

b

 

 

 

 

 

 

 

 

 

b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Решим полученное линейное однородное уравнение с постоянными коэффициентами заменой 1 t t0 ekt . Подставляя

1 t

 

t0 kekt ,

1 t

 

t0 k2ekt в последнее уравнение и сокращая

 

 

на множитель ekt 0 , приходим

к

 

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

уравнению

 

 

 

a b .

k2 a b k 0 , откуда k

0 ,

k

2

1

 

 

 

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

1 t t0 1 , 1 t t0 e a b t .

Общее решение уравнения есть 1 t t0 C1 C2e a b t . Тогда

2 t t0 1b dtd C1 C2e a b t ab C1 C2e a b t ba C1 C2e a b t .

Значения произвольных постоянных C1 и C2 в общем решении

118

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

1 t

 

t0 C1 C2e a b t ,

2

t

 

t0 a C1 C2e a b t

(3.28)

 

 

 

 

 

 

 

 

b

 

системы уравнений определим из начальных условий: в момент

времени t0 априорная вероятность

1 t0

 

t0 , априорная

 

вероятность 2 t 0

 

t0 1 .

 

С

 

учётом

 

сказанного (3.28)

 

 

 

 

принимает вид

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C C

2

e a b t0 ,

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

a C

C

 

e a

b t

0

.

 

 

 

 

1

 

 

 

 

 

 

 

b

1

 

2

 

 

 

 

 

 

Решая последнюю систему, находим

Cb , C b e a b t0 .

1a b 2 a b

Сучётом вычисленных C1 , C2 , а также введённых

обозначений, решение (3.28) запишется в виде

 

1

t

 

 

 

t0

 

2 1 P1( 2

| 2 )

 

 

 

 

 

 

 

 

 

 

 

 

1 1 P1( 1 | 1) 2

1 P1( 2 | 2 )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

1 P1( 2 | 2 )

 

 

 

 

0

 

 

 

 

 

 

 

 

e a b t t

 

,

 

 

 

 

 

 

 

 

 

 

 

1 1 P1( 1

| 1) 2 1 P1( 2 | 2 )

 

 

 

 

 

 

 

 

 

 

 

 

2

t

 

t0

 

1 1 P1( 1

| 1)

 

 

 

 

 

 

 

 

 

 

 

 

1 1 P1( 1 | 1) 2

1 P1( 2 | 2 )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

1 P1( 2 | 2 )

 

 

 

 

0

 

 

 

 

 

 

 

 

e a b t t

 

,

 

 

 

 

 

 

 

 

 

 

 

1 1 P1( 1

| 1) 2 1 P1( 2 | 2 )

 

 

 

 

 

 

 

 

 

 

 

где a и b определны в (3.27).

Методом исключения можно проинтегрировать и неоднородную систему порядка n 2 .

Пример 3.11. Найти все решения системы

119

Нежельская Л. А. Линейные дифференциальныеуравнения

x 5x 2 y 40et ,y x 6 y 9e t ,z x 10et .

Решение. Из первого уравнения системы выразим

y12 x 5x 40et

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

12 x 5x 40et x 3 x 5x 40et 9e t , x 11x 28x 280et 18e t .

Из третьего уравнения системы выразим

(3.29)

(3.30)

x z 10et

(3.31)

и подставим в (3.30).

Запишем уравнение третьего порядка относительно неизвест-

ной функции z t :

 

z 11z 28z 120et 18e t .

(3.32)

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

неизвестной функции z t .

1) Решим соответствующее однородное уравнение с постоянными коэффициентами. Характеристическое уравнение

k3 11k 2 28k 0

имеет простые действительные корни k1 0 , k2 4 , k3 7 .

Корням характеристического уравнения соответствуют линейно независимые решения

z1 1 , z2 e 4t , z3 e 7t .

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

z0 C1 C2e 4t C3e 7t .

120

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