Скачиваний:
0
Добавлен:
12.02.2026
Размер:
963.13 Кб
Скачать

Глава III

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

§ 8. Линейные ОДУ n-ого порядка с переменными коэффициентами.

8.1. Эквивалентность линейного ОДУ n-ого порядка системе n уравнений первого порядка.

Рассмотрим линейное ОДУ n-ого порядка в общем виде:

 

 

x(n) (t) a

(t)x(n 1) (t) ... a

n 1

(t)x (t) a

n

(t)x(t) b(t) ,

 

 

 

 

 

1

 

 

 

 

где a j , b C a, b

 

 

 

 

 

 

для любого j 1, n .

 

 

 

Задача Коши: Найти решение x x(t ) ,

x Cn a, b уравнения (51),

удовлетворяющее следующим условиям:

 

 

 

x(t

0

) x1

,

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

2

,

 

 

 

 

 

 

x (t0 )

x0

t0 a, b ,

 

 

 

 

 

 

 

 

 

 

,

 

 

 

 

 

. . . . . . .

 

 

 

 

 

 

x(n 1) (t

0

) xn ,

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

(51)

(52)

здесь x1

, x2

, ..., xn – заданные постоянные числа из .

0

0

0

 

 

 

Произведем следующие замены:

y 1 (t ) x(t ) ,

y 2 (t)

тогда получим следующую систему:

 

 

 

 

 

 

 

 

dy 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

]в силу (53)[ x (t) y 2 (t ) ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dy 2

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x (t ) y

 

(t ) ,

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

dyn 1

x(n 1) (t ) y n (t ) ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dy n

 

x(n) (t ) a (t )x(n 1) (t ) ... a

 

 

(t )x (t) a (t )

 

 

dt

 

 

 

 

 

1

 

 

 

n 1

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

(t ) y n (t ) ... a

 

 

(t ) y 2 (t ) a

 

(t ) y

 

 

 

 

 

 

n 1

n

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

x (t) , …, y n (t ) x(n 1) (t ) , (53)

(21)

x(t ) b(t )

1 (t ) b(t )

 

dy (t )

 

 

 

или в векторной форме

A(t ) y

(t ) f (t ) ,

(21 )

 

 

dt

 

 

 

 

 

 

 

y 1 (t)

 

 

x(t)

 

 

 

 

 

 

 

2

 

 

 

x (t )

 

где y

 

(t )

 

y

 

(t)

]в силу (53)[

 

,

 

 

. . . . .

 

. . . . .

 

 

 

 

 

y

n

 

 

 

(n 1)

 

 

 

 

 

 

 

 

(t)

x

 

(t)

 

 

0

 

 

 

1

0

 

0

 

 

 

 

 

0

 

 

 

0

1

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A(t )

 

 

 

 

 

 

 

 

 

,

 

 

0

 

 

 

0

0

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

an (t )

 

an 1

(t ) an 2 (t)

a1 (t)

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

f

(t )

0

.

 

 

 

 

 

 

b(t )

Вывод 8.1. Любое решение уравнения (51) вместе с его производными до ( n 1 )-ого

порядка в совокупности задают решение системы ( 21 ), и, наоборот, первая координата

любого решения системы ( 21 ) является решением уравнения (51). Тогда в этом смысле

 

 

 

 

 

 

 

 

 

 

 

 

 

будем понимать, что уравнение (51) является эквивалентным системе ( 21 ) .

 

Аналогично перепишем начальные условия:

 

 

 

 

 

 

 

 

x(t0 )

 

 

x 10

 

 

 

 

y

(t

 

 

x (t0 )

 

 

x 2

 

y

.

(22 )

 

)

 

 

 

 

0

 

 

 

0

 

. . . . .

 

. . . . .

 

0

 

 

 

 

 

(n 1)

(t0 )

 

 

n

 

 

 

 

 

 

 

x

 

 

 

x 0

 

 

 

 

Теорема 8.1. (ТСЕ) Пусть функции a j , b C a, b

 

 

 

для любого j 1, n .

Тогда для любого t0 a, b

и для любых x 10, x 02, ..., x 0n

существует единственное

решение задачи Коши (51),(52)

класса Cn a, b , определенное на a, b .

 

 

 

 

# В силу вывода 8.1 задача Коши (51),(52) эквивалентна задаче Коши ( 21 ),( 22 ) ,

причем из условия a j , b C a, b для любого

 

 

 

 

 

 

 

j 1, n следует, что вектор-функция

 

 

 

 

 

 

 

 

 

 

f (t ) и матрица A(t) являются непрерывными на a, b .

 

 

 

 

 

 

 

 

1 a, b задачи

Тогда по теореме 4.1 существует единственное решение y

C

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Коши ( 21 ),( 22 ) . Следовательно, y j C1 a, b

для любого

j 1, n , в том числе и

для

y n x(n 1) C1 a, b , т.е. x Cn a, b .

 

 

 

 

 

 

 

 

Итак, существует единственное решение x x(t ) класса Cn a, b

задачи Коши

(51),(52) на отрезке a, b .

#

Замечание 8.1. Решение задачи Коши (51),(52) является глобальным, т.к. оно определено на всем отрезке a, b , где решается данная задача Коши .

8.2. Однородное линейное ОДУ n-ого порядка, его ФСР n общее решение.

Рассмотрим однородное линейное ОДУ n-ого порядка вида:

 

 

 

x(n) (t) a (t)x(n 1) (t) ... a

 

 

(t)x (t) a (t)x(t) 0 ,

 

 

(51 )

 

 

 

 

1

 

 

 

 

n 1

 

n

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где a j C a,b для любого j

1, n

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

При этом, как и в пункте 8.1, уравнение ( 510 ) сводится эквивалентным образом к

 

 

 

 

 

 

dy (t )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

системе:

 

 

A(t ) y (t) .

 

 

 

 

 

 

 

 

 

 

(21 0 )

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

Лемма 8.1. Если решение уравнения ( 510 ) обращается в ноль одновременно со своими

производными до ( n 1 ) порядка хотя бы в одной точке из a, b

, то оно тривиально,

 

 

т.е. x(t ) 0 на a, b .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

)

# Утверждение леммы 8.1 следует из эквивалентности уравнения ( 510 ) и системы ( 210

и леммы 5.1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

#

 

 

Пусть H – множество решений уравнения ( 510 ), которое, очевидно, не является

 

 

пустым множеством, т.к. x(t ) 0 есть решение уравнения ( 510 ) .

 

 

 

 

 

 

 

 

 

Для элементов множества H введем операции сложения элементов этого множества

 

и умножения элементов на число.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Лемма 8.2.

Для множества H справедливы следующие свойства :

 

 

 

 

 

 

 

1о

H – линейное пространство (относительно введенных операций).

 

 

 

 

 

 

 

2о

dim H n .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

# Свойство

о

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

, а также

 

1 следует из эквивалентности уравнения ( 510 ) и системы ( 210 )

 

из леммы 5.2.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

о

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Свойство 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

следует из эквивалентности уравнения ( 510 ) и системы ( 210 ) , а также из

 

теоремы 5.1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

#

 

 

 

 

 

 

 

 

 

 

 

 

j

 

j 1

 

 

 

 

 

 

 

 

 

 

Определение 8.1. Система функций

 

x

 

(t) n

называется:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

n

 

 

1) линейно зависимой на a, b , если существуют числа j j 1

такие, что

 

j

0 ,

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а

j x j (t) 0 (*) для любых t a, b

;

 

 

 

 

 

 

 

 

 

 

 

 

 

j 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2) линейно независимой на a, b , если равенство (*) выполняется в том и только в том случае, когда j 0 для любых j 1, n .

 

 

 

 

j

j 1

 

 

 

 

 

 

 

 

Лемма 8.3. Система функций

 

x

 

(t)

n

линейно зависима (линейно независима) на

a, b

тогда и только тогда, когда линейно зависима (линейно независима) на a, b

 

 

 

 

 

 

 

 

 

 

 

 

x j (t )

 

 

 

 

 

n

 

 

 

 

 

 

 

x j (t )

 

 

система вектор-функций y j (t) j 1 ,

 

 

y j (t )

 

 

 

где

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

x (n 1) (t )

 

 

 

 

 

 

j 1

 

 

 

 

 

j

 

 

#

 

 

j

 

линейно зависимая система на a, b ,

Необходимость. а) Пусть

 

x

 

(t)

n

 

n

 

 

 

 

 

 

n

 

 

 

 

n

 

т.е. существуют числа j j 1

 

такие, что

 

j

0 , но

j x j (t) 0 (*)

на a, b .

 

 

 

 

 

 

 

 

 

 

 

 

 

j 1

 

 

j 1

 

Продифференцируем (*) по t

k 0, n 1 . Следовательно,

 

 

 

 

n

Итак, система y j (t) j 1

 

 

 

j

j 1

б) Пусть система

 

x

 

(t) n

n

( n 1 ) раз, тогда j x(jk ) (t) 0 на a, b

j 1

n

j y j (t) 0 для любого t a, b .

j 1

является линейно зависимой на a, b .

является линейно независимой на a, b , т.е.

для любого

n

j x j (t) 0 j 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

на a, b тогда и только тогда,

когда j

0

для любых j 1, n , тогда j y j (t) 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

для любого t a, b в том и только в том случае, когда j 0 для любых j 1, n .

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Следовательно, система y j (t) j 1

является линейно независимой на a, b .

 

 

 

 

 

 

 

 

 

 

 

Достаточность, Пункт а) очевиден, т.к.

x j (t ) для любых j 1, n являются первыми

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

координатами вектор-функций y j (t) ,

 

j 1, n .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

Пункт б) докажем от противного. Предположим, что система y j (t) j 1

является

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j

j 1

 

 

 

 

 

 

 

 

 

 

 

 

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

 

a, b

 

, но система

 

x

 

(t) n

будет линейно зависимой на

 

a, b

 

,

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

тогда в силу пункта а) система y j (t) j 1

должна быть линейно зависимой на a, b ,

 

 

что противоречит предположению.

 

#

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Определение 8.2. Фундаментальной системой решений (ФСР) уравнения ( 510 )

называется любая линейно независимая на a, b система из n решений уравнения ( 510 ) .

Теорема 8.2. Фундаментальная система решений (ФСР) уравнения ( 510 ) обладает следующими свойствами:

1о ФСР уравнения ( 510 ) всегда существует и может быть найдена решением n задач Коши для уравнения ( 510 ) и начальных условий:

x(t

0

) 1

,

 

 

 

 

 

 

 

 

 

 

 

 

 

j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

x (t0 ) j

,

,

 

 

 

 

для любого t0 a, b ,

 

j 1, n ,

(52j )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x(n 1) (t

0

) n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j

 

 

 

 

det 1

 

 

0 .

где векторы-столбцы

j

 

 

таковы, что

, ..., n

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j

 

 

 

 

 

 

 

 

2о Если x j (t) nj 1 ФСР уравнения ( 510 ), то общее решение уравнения ( 510 ) имеет вид:

n

xоо (t) С j x j (t) , (54)

j 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

С j – произвольные постоянные для любого

j 1, n .

 

 

 

 

 

 

 

 

 

 

#

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1о ): В силу теоремы 8.1 для любого j 1, n

задача Коши ( 51

),( 52

j

) имеет

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

единственное решение x j C n a, b .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(t) n

 

задач Коши ( 51 ),( 52

 

 

 

 

 

 

 

Рассмотрим систему полученных решений

x

 

 

 

),

j 1, n .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j

 

 

j 1

0

 

 

 

j

 

 

 

 

Эта система линейно независимая на a, b также, как система вектор-функций

 

 

 

 

 

 

 

 

 

 

x j (t )

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x j (t )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(t )

 

 

 

 

,

ибо в противном случае существуют числа j

 

n

 

 

такие,

y j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j 1

 

 

 

 

 

 

 

 

 

 

x (n 1) (t )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j

 

 

 

j 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

n

 

 

 

 

для любого t a, b ; в том числе и для t

 

 

 

 

 

что

j

0 , но

j y j

(t) 0

t0 , т.е.

 

 

 

j 1

 

 

 

 

j 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

j y j (t0 ) 0

или j j 0 . Это означает, что система j j 1

 

будет линейно

j 1

 

 

 

 

 

 

 

 

 

 

j 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

зависимой, и det

1 , ..., n 0 , что противоречит условию задачи Коши.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j

 

 

 

j 1

 

 

 

 

 

 

 

 

 

 

 

 

 

Итак, в силу леммы 8.3 система решений

 

x

 

 

(t)

n

является линейно независимой

на a, b , а значит она является ФСР уравнения ( 510 ) .

 

 

 

 

 

 

 

 

 

 

 

( 2о ): Поскольку

 

x

 

(t)

n

– ФСР уравнения (

51

) и dim H n (см. лемму 8.2), то система

 

 

 

 

j 1

 

 

 

j

 

j 1

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

(t) n

является базисом пространства H .

 

 

Тогда для любого решения x(t ) H

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

что x(t) С j x j (t) .

 

 

 

 

 

существуют произвольные постоянные С j j 1

 

 

такие,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j 1

 

 

 

 

 

 

 

 

 

n

Итак, x(t) С j x j (t) является общим решением уравнения ( 510 ) . #

j 1

8.3. Определитель Вронского системы функций и его свойства.

 

Определение 8.3. Определителем Вронского (вронскианом) системы функций

 

j

j 1

x1 (t )

xn (t )

 

x

(t )

x

(t )

 

 

 

n

 

x (t)

называется определитель W ( x(t )) det

1

 

 

n

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(n 1)

 

(n 1)

 

 

 

 

 

x1

 

(t )

xn

 

(t )

 

 

 

j

j 1

 

 

 

 

0

 

 

 

 

Теорема 8.3. Пусть

 

x

 

(t)

n

система решений уравнения ( 51

) на

 

a, b

 

,

W ( x(t )) вронскиан этой системы решений.

 

 

 

 

 

 

 

Тогда для любых t, t0 a, b

справедлива формула Лиувилля:

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

W ( x(t )) W ( x(t0 ))exp

 

 

a1

( )d .

 

 

 

 

(55)

 

 

 

 

 

 

 

 

t0

 

 

 

 

 

 

 

#Формула (55) следует из теоремы 5.4, поскольку уравнение ( 510 ) эквивалентно

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

#

 

 

 

 

 

 

),

и Sp A(t) a1 (t)

 

(см. пункт 8.1) .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

системе ( 210

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Теорема 8.4.

 

Определитель Вронского системы функций обладает следующими

 

 

 

свойствами :

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j

 

j 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1о

Если

 

z

 

(t)

n

является линейно зависимой системой функций на

 

a, b

 

, то

 

 

 

 

 

 

 

W ( x(t )) 0 на a, b .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j

j 1

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2о

Система

 

x

 

(t) n

решений уравнения ( 51

) является линейно зависимой на

 

a, b

 

 

 

тогда и только тогда, когда существует t0 a, b такое, что W ( x(t0 )) 0 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j

 

j 1

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3о

Система решений

 

x

 

(t) n

 

уравнения ( 51

) является линейно независимой на

 

a, b

 

 

тогда и только тогда, когда W ( x(t )) 0 для любого t a, b .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4о

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j

 

 

n

 

 

 

0

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

Если для системы решений

 

x

 

(t)

 

уравнения ( 51

) существует t

 

 

 

a, b

 

,

для

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

которого W ( x(t1 )) 0 ,

то W ( x(t )) 0 для любого t a, b .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

#

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j

j 1

является линейно зависимой на a, b

 

 

 

(1о ): Пусть система функций

 

z

 

(t) n

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z j (t )

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z j (t )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(t )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

тогда в силу леммы 8.3 эквивалентная система z j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z (n 1) (t )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j

 

j 1

 

 

 

 

 

 

 

 

 

 

 

 

также будет линейно зависимой на a, b , а тогда в силу свойства 1о теоремы 5.5

Wz (t) 0 на a, b , но Wz (t ) W (z(t )) на a, b .

Свойства 2о , 3о , 4о следуют из свойств 2о , 3о , 4о теоремы 5.5 в силу эквивалентности

 

 

 

#

 

 

) , а также на основание леммы 8.3 .

уравнения ( 510 ) и системы ( 210

8.4. Неоднородные линейные ОДУ n-ого порядка.

Рассмотрим неоднородное уравнение вида:

x(n) (t) a (t)x(n 1) (t) ... a

n 1

(t)x (t) a

n

(t)x(t) b(t) ,

(51)

1

 

 

 

где a j , b C a, b для любого j

1, n

.

 

 

 

Теорема 8.5. Для любого t a, b

общее решение неоднородного ОДУ (51) имеет

следующую структуру :

 

 

 

 

 

xон (t ) xоо (t) xчн (t ) ,

 

 

(56)

где xоо (t ) общее решение однородного уравнения ( 510 ), а xчн (t ) частное решение неоднородного уравнения (51).

# Пусть xоо (t ) – общее решение однородного уравнения ( 510 ), а xчн (t ) – частное решение неоднородного уравнения (51).

1) Рассмотрим функцию x(t) xоо (t) xчн (t) и подставим ее и ее производные до (n 1) -ого порядка включительно в ОДУ (51).

x(n) (t ) a1 (t)x(n 1) (t) ... an 1 (t)x (t) an (t)x(t) xоо (t) xчн (t) (n) a1 (t) xоо (t)xчн (t) (n 1) ... an 1 (t) xоо (t) xчн (t) an (t) xоо (t) xчн (t) 0 b(t) b(t) .

Следовательно, функция x(t ) есть решение неоднородного уравнения (51) при любых постоянных, входящих в xоо (t ) .

2) Пусть теперь x(t ) есть произвольное частное решение уравнения (51), тогда рассмотрим разность x(t ) xчн (t ) . Подставив эту разность в уравнение (51), получим

x(t) xчн (t) (n) a1 (t) x(t) xчн (t) (n 1) ... an 1 (t) x(t) xчн (t) an (t) x(t) xчн (t)

b(t ) b(t) 0 , т.е. эта разность является частным решением уравнения ( 510 ). И его

можно найти из общего решения ОДУ ( 510 ) подбором соответствующих констант

 

 

 

 

из xоо (t ) .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Итак, из 1) и 2) следует, что x(t ) xоо (t) xчн (t ) является общим решением

 

 

 

 

 

неоднородного уравнения (51).

 

 

 

 

 

 

 

 

 

 

 

#

 

 

 

 

 

 

 

 

j

j 1

0

 

 

 

 

 

 

 

 

 

 

 

Следствие 8.1. Если

 

x

 

(t) n

ФСР уравнения ( 51 ) на

 

a, b

 

,

то общее решение

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

уравнения (51) может быть записано в виде: xон (t) С j x j (t) xчн (t) .

 

 

 

(57)

 

 

 

 

 

 

 

 

 

 

j 1

 

 

 

 

 

 

 

 

 

 

 

# Формула (57) следует из (56) в силу равенства (54).

 

 

 

 

 

 

 

 

 

#

Метод вариации постоянных. Для поиска частного решения xчн (t ) уравнения (51)

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

 

 

 

 

 

 

 

(t) n

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

Пусть

x

 

– ФСР уравнения ( 51 ), тогда в силу формулы (54)

x (t)

 

С

 

x

 

(t) .

 

 

j

j 1

 

 

 

 

 

0

 

 

 

 

оо

 

j

 

j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j 1

 

 

 

 

 

На основании этой формулы частное решение xчн (t ) будем искать в виде

n

 

 

 

x(t) q j (t)x j (t) ,

где q j (t ) – искомые функции для любого j 1, n .

j 1

 

 

Поскольку уравнение (51) эквивалентно системе ( 21 ) , то как и в пункте 6.2 рассмотрим

 

 

 

 

равенство X (t)q

(t)

f

 

(t) (*) , где

 

 

 

 

x1 (t )

xn (t )

 

q (t )

 

 

 

x (t )

x (t )

 

 

 

 

1

(t )

 

 

X (t )

1

n

,

q

(t )

q2

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x(n 1) (t )

x(n 1) (t )

 

q

n

(t )

 

 

 

1

n

 

 

 

 

 

 

 

Распишем равенство (*) поэлементно

 

 

 

 

 

(t

x1 (t )q1

(t ) x2 (t )q2

x (t )q

(t ) x

(t )q

(t

 

1

1

2

2

 

 

(n 2)

 

 

(n 2)

x1

(t )q1 (t ) x 2

 

 

x(n 1) (t )q (t ) x(n 1)

 

1

 

1

2

 

 

 

 

)

x

n

(t )q

(t )

 

 

n

 

)

x

(t )q

(t )

 

 

n

n

 

(t )q

(t )

x(n 2)

2

 

 

 

n

(t )q

(t )

x( n 1)

2

 

 

 

n

 

 

 

 

0

 

 

 

 

 

 

 

f

 

(t )

0

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b(t )

 

0 ,

 

 

 

 

 

0 ,

 

 

 

 

 

 

 

 

 

,

 

(**)

 

 

 

 

 

 

(t )q

(t) 0 ,

 

 

 

 

n

 

 

 

 

(t )q

(t ) b(t )

 

 

 

 

n

 

 

 

 

Поскольку

 

x

 

(t)

n

– ФСР уравнения ( 51 ), то det X (t) W ( x(t)) 0 для любого

 

 

 

 

 

 

j

 

j 1

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Wj ( x(t))

 

 

t a, b . Тогда по теореме Крамера для любого j 1, n q j (t)

, где W j ( x(t ))

W ( x(t))

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

получен из W ( x(t )) заменой j -ого столбца матрицы X (t ) столбцом f (t) .

 

 

 

 

 

 

 

 

 

 

 

t

W

( x( ))

 

 

 

 

 

 

 

 

 

 

 

 

Итак, q j (t )

 

j

 

d

для любого j 1, n .

 

 

(58)

W ( x( ))

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Вывод 8.2.

 

Общее решение уравнения (51) представимо в виде :

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

n

t

W ( x( ))

 

 

 

 

 

 

 

 

 

 

 

 

xон (t) С j x j (t) x j (t)

j

d для любых t, t0 a, b

,

(59)

 

 

 

 

W ( x( ))

 

 

 

 

 

 

 

 

 

 

j 1

 

j 1

t0

 

 

 

 

 

 

 

 

 

 

где

 

x

 

(t)

n

 

 

ФСР уравнения ( 51

) , а С

 

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

для

 

 

j

j 1

 

 

 

 

 

 

0

 

j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

любого

 

j 1, n .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

§ 9. Линейные ОДУ n-ого порядка с постоянными коэффициентами.

Рассмотрим линейное ОДУ n-ого порядка с постоянными коэффициентами в виде

оператора:

A x(t) a0 x(n) (t) a1 x(n 1) (t) ... an 1 x (t) an x(t) b(t) ,

(61)

где a0 0 ,

 

 

 

 

a j – в общем случае комплексные числа для любого j 0, n .

 

9.1. Линейные ОДУ n-ого порядка с постоянными коэффициентами.

Неоднородному уравнению (61) соответствует однородное вида:

 

A x(t) a0 x(n) (t) a1 x(n 1) (t) ... an 1 x (t) an x(t) 0 ,

(610 )

или кратко A x 0 (*) .

 

Замечание 9.1. Решение уравнения ( 610 ) всегда существует и определено на .

 

Будем искать решение уравнения ( 61 ) в виде x(t) e t , при этом требуется

 

0

 

определить. Из равенства (*) найдем A

e t (a0 n a1 n 1 ... an 1 an )e t 0 ,

т.е.

A e t P ( )e t

,

 

 

 

 

(62)

где P ( ) a

n a n 1

... a

n 1

a

n

называется характеристическим многочленом.

0

1

 

 

 

 

 

 

 

 

Определение 9.1. Уравнение P ( ) 0 или a n a n 1 ... a

n 1

a

n

0 (63)

 

 

 

 

 

 

0

1

 

 

называется характеристическим уравнением ОДУ ( 610 ).

 

 

 

 

Вывод 9.1. Функция x(t) e t

является решением уравнения ( 61 ) тогда и только

 

 

 

 

 

 

 

 

0

 

 

 

тогда, когда число есть корень характеристического уравнения (63) .

Рассмотрим все возможные варианты корней характеристического уравнения (63).

Первый случай. (Корни характеристического уравнения простые).

Теорема 9.1. Пусть j n

 

простые корни характеристического уравнения (63),

 

 

 

 

 

 

 

 

j 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

т.е. j k для любых

j k ,

 

 

 

 

 

 

 

 

 

 

 

 

j, k 1, n .

 

 

 

 

 

Тогда e j t n

 

ФСР уравнения ( 610 ).

 

 

 

 

 

 

 

j 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

# Пусть j n

 

– простые корни характеристического уравнения (63), тогда в силу

 

 

j 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(t ) e j t

является решением уравнения ( 61

 

вывода 9.1 для любого

j 1, n функция

x

j

).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

Рассмотрим систему x j (t) e j t n

 

 

и построим для нее вронскиан W ( x(t )) . Вычислим

 

 

 

 

 

 

 

 

 

 

 

 

j 1

 

 

 

 

 

 

 

 

этот вронскиан при t 0 :

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e 1t

 

 

e nt

 

 

 

 

1

1

 

 

 

 

1

e 1t

 

 

n

e nt

 

 

 

 

 

1

n

 

 

W ( x(0)) det

 

 

 

 

 

 

 

 

 

 

 

 

det

 

 

]это определитель

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n 1

 

1t

 

n 1

 

nt

 

 

 

 

n 1

n 1

 

 

 

 

1

 

 

e

 

 

n

 

e

 

 

 

t 0

 

 

 

 

1

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(*)

Вандермонда[ ( j k ) 0 .

k j

Тогда в силу свойств 4о и 3о теоремы 8.4 система x j (t) e j t n

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

j 1

 

а значит является ФСР уравнения ( 610 ) .

#

Упражнение 9.1. Доказать равенство (*) самостоятельно.

Второй случай. (У характеристического уравнения есть кратные корни).

Лемма 9.1. Пусть Pm (t ) многочлен степени m 0 .

 

 

 

 

 

 

Тогда производная

d k

Pm (t)e t Qm (t)e t , k

,

где 0 ,

а Qm (t )

 

 

k

 

 

 

d t

 

 

 

 

 

 

 

 

 

 

 

 

 

некоторый многочлен степени m .

 

 

 

 

 

 

 

 

 

 

 

# Рассмотрим производную

d

 

 

t

 

 

 

t

 

Pm (t) Pm 1 (t) e

t

 

 

 

Pm (t)e

 

Pm (t) Pm (t) e

 

 

,

d t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где Pm 1 (t) – многочлен степени m 1 , если m 0 и

Pm 1 (t) 0 ,

если m 0 .

 

 

Далее по индукции получаем, что deg Qm (t) m .

 

 

 

 

#

 

 

Теорема 9.2. Пусть 1 , 2 ,..., m

корни характеристического уравнения (63)

 

 

кратности k1 , k2 , ..., km соответственно, при этом k1 k2

... km n .

 

 

x

1

(t ) e 1t , x

2

(t ) te 1t , ..., x

(t ) t k1 1e 1t

Тогда система

 

 

 

 

k1

 

 

 

 

 

 

 

 

x

n km 1

(t ) e m t , ........... , x

 

(t ) t km 1e m t

 

 

 

 

 

n

 

(64)

будет ФСР уравнения ( 610 ) .

# а) Возьмем произвольный номер p 1, m и рассмотрим корень p характеристического уравнения (63) кратности k p , т.е.

P (

 

) 0, P (

 

 

) 0, ...,

 

 

P

k p 1

(

 

)

0

(65) ,

но P

k p

(

 

 

) 0 .

p

p

 

 

 

 

 

 

 

 

p

 

 

p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

подействуем на t r e t оператором A :

Обозначим

p

 

, тогда для любого

r 0, k

p

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d r

 

 

 

 

 

 

 

 

d r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d r

 

 

 

 

 

A t r e t A

 

 

 

 

 

e t

 

 

 

 

 

A e t ]в силу формулы (62)[

 

 

 

P ( )e t

 

 

 

 

r

d

r

d

r

 

 

d

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

]по формуле Лейбница[ Сrj (e t )(r j )P ( j ) ( ) ]по формуле (65)[ 0 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Следовательно, A e pt 0 ,

 

A te pt

0 , …, A t k p 1e pt

0

 

 

 

 

для любого p 1, m .

Итак, функции из формулы (64) являются решениями уравнения ( 610 ).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j

 

j 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

б)

Докажем, что система

 

 

x

 

 

(t)

n

 

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

 

 

. От противного.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j

 

 

j 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Предположим, что система

 

 

x

 

(t)

n

будет линейно зависимой, т.е. существуют числа

 

 

n

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j j 1 такие, что

 

j

 

 

0 , но j x j (t) 0 (*)

на .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j 1

 

 

 

 

 

 

 

 

 

 

 

j 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Соседние файлы в папке Ещё лекции