Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
ODU_4.pdf
Скачиваний:
76
Добавлен:
15.04.2015
Размер:
725.58 Кб
Скачать

Фундаментальная система решений (ФСР) для линейной однородной системы ОДУ. Существование ФСР и их взаимосвязь. Общее решение линейной однородной неоднородной системы.

def: Линейно независимая система решений y1,, yn системы (2) называется фундаментальной системой решений (ФСР) системы (2), а матрица

y11(x)

Y (x)=

y1n (x)

y1n (x)

ynn (x)

называется фундаментальной матрицей решений (ФМР) системы (2).

Теорема:

ФСР системы (2) существуют. Их бесконечно много. Все они могут быть получены из одной по формуле Z (x)=Y (x) S , где Y (x) – фиксированная ФМР, а S – матрица

n ×n , с условием det S 0 .

 

 

Доказательство:

 

 

◄ Фиксируем произвольную точку x0 = [a,b]. Пусть B – матрица n ×n

с условием

det B 0 , тогда по теореме о ! решения задачи Коши j =

 

 

! решение

y j (x) систе-

1,n

мы (2) ( y′ = A(x)y (2)), удовлетворяющее начальным условиям

ykj (x0 )= bkj . В этом слу-

чае y1,, yn – ФСР, так как W (x)= detY (x0 )= det B 0 , поскольку матриц B с условием det B 0 бесконечно много, то и ФСР бесконечно много. Пусть теперь Y (x) – фиксиро-

ванная ФМР, а Z (x) – произвольная ФМР системы (2). Положим S =Y

1(x

) Z(x ), тогда

d(Y (x)S ) =

dY (x) S = (A(x) Y (x)) S = A(x)(Y (x) S ),

0

0

 

 

dx

dx

 

 

то есть Y (x)S – матрица из решений системы (2), причем

 

 

Y (x0 )S =Y (x0 )(Y 1(x0 ) Z(x0 ))= (Y (x0 )Y 1(x0 )) Z(x0 )= Z(x0 ).

 

По теореме о единственности решения задачи Коши (для каждого столбца) получим

Y (x)S Z (x) на [a,b]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Теорема:

 

 

 

 

 

 

 

 

y(x)= C1 y1(x)+…+Cn yn (x) (3), где

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

системы

(2)

имеет вид

y1,, yn – ФСР системы (2), а C1,,Cn

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

Доказательство:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

◄ ( ) y(x) из формулы (3) является решением системы (2) в силу ее линейности и од-

нородности.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~

 

 

 

 

 

 

 

 

 

 

1

,, y

n

– ФСР системы (2).

( ) Пусть y(x) – произвольное решение системы (2), а

 

y

 

Фиксируем точку x0 [a,b] тогда система относительно C1,,Cn :

 

 

 

 

 

 

 

 

 

1

 

n

 

~

 

 

 

 

 

 

 

 

 

 

C1 y1 (x0 )+…+Cn y1

(x0 )= y1(x0 )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

n

 

~

 

 

 

 

 

 

 

 

 

 

C1 yn

(x0 )+…+Cn yn

(x0 )= yn (x0 )

 

 

 

 

 

имеет единственное решение

~

~

, то есть ее определитель равен W (x0 )0 . Поло-

C1,,Cn

~ 1

 

~

n

(x), тогда

y(x) – решение системы (2) в силу ее линейности

жим y(x)= C1 y

(x)+…+ Cn y

 

и однородности, причем

 

 

 

~

(x0 ) по построению. По теореме о единственности ре-

y(x0 )= y

 

 

~

 

 

 

 

~

~

1

(x)+…

 

~

n

(x). ►

шения задачи Коши y(x)y(x) на [a,b], то есть

y

(x)= C1 y

 

+ Cn y

 

15

Теорема:

Общее решение неоднородной линейной системы (1) y′ = A(x)y + f (x) имеет вид: y(x)= yчастное (x)+C1 y1(x)+…+Cn yn (x),

где yчастное (x) – частное решение системы (1), y1,, yn – ФСР, соответствующая однородной системе (2), а C1,,Cn – произвольные постоянные. Это следует из теории линейных пространств и предыдущей теоремы.

Резольвента линейной системы ОДУ и ее свойства.

def: Резольвентой системы (2) называют отображение R : [a,b]×[a,b]GL(n) – множество невырожденных матриц порядка n, удовлетворяющих следующим условиям:

1) x0 [a,b]: R(x, x0 ) – ФМР системы (2),

2) R(x0 , x0 )= E =

1

0

.

 

 

 

0

1

 

Свойства:

dR(x, x0 ) = A(x)R(x, x0 ) – следует из определения. dx

2° Пусть

y(x)

 

решение системы

(2) с начальными условием

y(x0 )= y0 ,

тогда

y(x)= R(x, x0 )y0 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Доказательство:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

◄ Имеем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d(R(x, x )y

 

)

=

dR(x, x

)

y

 

â ñèëó 1o

)y

 

– решение системы (2),

 

 

 

 

 

0

 

0

 

 

0

 

0

=

A(x)R(x, x

0

 

 

 

 

dx

 

 

 

 

 

 

dx

 

 

 

 

0

 

 

 

 

 

 

 

 

причем

 

R(x0 , x0 )y0

Ey0 = y0 .

 

По

теореме о

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

 

задачи

Коши

R(x, x0 )y0 y(x) на [a,b]. ►

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

′ ′′

 

 

 

 

 

′′

) (полугрупповое свойство)

 

 

 

 

 

 

 

 

 

R(x, x

)R(x , x

)= R(x, x

 

 

 

 

 

 

 

 

 

Доказательство:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

◄ Пусть

 

y0 – произвольный вектор из Cn , тогда R(x, x′′)y0 – значение в точке

x

того

решения системы (2), которое в точке

x′′ равно y0

(свойство 2°). R(x, x)R((x, x′′)y0 )

значения точки x того решения, которое в xравно R(x, x′′)y0 (свойство 2°), то есть того

решения

 

(в силу

 

существования

единственности),

которое

в точке

′′

равно

y0 .

 

 

x

R((x, x′′)y0 ) – значение

в точке

x

того решения, которое в

точке

x′′

равно

y0

 

R(x, x)R(x, x′′)y0 =

R(x, x′′)y0 , в силу произвольного y0 получаем требуемое. ►

 

 

[R(x, x0 )]1 = R(x0 , x). Это следует из 3°, так как R(x, x0 ) R(x0 , x)= R(x, x)= E (свойство 3°). Слева аналогично.

dR(x1, x)= −R(x1, x)A(x) dx

Доказательство:

16

 

R(x1, x + x)R(x1, x)

= R(x , x + x)

R(x, x1 )R(x + x, x1 )

R(x

, x)

=

[R(x + x, x )]1 ×

 

 

 

 

 

 

x

 

 

 

 

 

1

x

 

 

 

 

1

 

 

 

1

 

R(x + x, x )R(x, x )

 

 

 

1

 

 

dR(x, x )

 

 

 

 

×

1

1

R(x

, x)

[R(x, x )]

 

1

R(x

, x)

=

R(x , x)×

 

 

 

 

x

 

 

1

 

x0

1

 

 

dx

 

1

 

 

1

 

 

 

 

 

 

 

 

 

свойство 3î

×(A(x)R(x, x1 )) R(x1, x)= −R(x1, x)A(x)E

 

 

 

 

 

 

 

 

 

 

свойство 1î

свойство 3î

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Построение линейной, однородной системы по известной ФСР. Формула Лиувилля.

Теорема:

Пусть y1(x),, yn (x) – произвольные n-мерные векторы-функции класса C1([a,b]), причем W (x)0 x [a,b], тогда ! система вида (2), для которой y1,, yn являются

ФСР. Доказательство:

◄ Рассмотрим систему

 

 

 

dyk

dy1k

dykn

 

 

 

1

 

dx

dx

dx

 

 

 

 

 

1

n

 

 

 

 

 

 

 

 

 

 

 

y1

y1

y1

= 0 , k =1, n .

W (x)

 

 

 

yn

y1n

ynn

 

 

 

Это система вида (2) (надо разложить по первому столбцу)

dyk = n akj (x)y j ,

dx j =1

где akj (x) – известные непрерывные функции на [a,b]. Заметим, что y1,, yn – решение

этой системы (так как при подстановке y j вместо y получаем два одинаковых столбца).►

Единственность.

 

 

~

Пусть Y (x) – ФМР системы

y′ =

 

A(x)y и y′ = A(x)y . Тогда

 

Y

(x)

~

 

A(x)Y (x)A(x)Y (x).

 

 

 

Умножаем справа на Y 1(x)

(Y 1(x) – существует, так как detY (x)=W (x)0

~

 

 

 

Получим, что A(x)A(x) на [a,b].

 

Формула Лиувиля.

Пусть y1,, yn – ФСР системы (2), тогда W (x)= J1 +…+ Jn , где

 

 

y11

y1n

 

y11

y1n

 

 

y11

 

 

 

 

 

 

 

n

 

n

n

линейность

n

 

 

 

1

 

определите ля

 

Jk

=

dyk

dyk

=

aks (x)y1s

aks (x)ysn

=

aks (x)

y1s

 

 

dx k я строка

dx

 

s=1

s=1

по строкам

s=1

 

 

 

 

 

 

 

 

 

y1n

 

 

y1n

ynn

 

y1n

ynn

 

 

 

 

 

 

 

 

 

 

 

 

k ya

 

 

 

 

= akk (x) W (x), stroka

так как при s k получаем две одинаковые строки

x [a,b]).

y1n

ysn =

ynn

17

Рассмотрим систему:

 

n

 

 

 

ÎÄÓ ñ

ìèñÿ

 

 

(x)=

разделяющи

W (x)= akk (x)W

=

Tr A(x)W (x)

 

k=1

 

 

 

переменным

è ñëåä A(x )

 

 

 

x

 

 

 

 

 

 

 

 

 

– Формула Лиувиля.

 

W (x)=W (x )exp

 

Tr A(t)dt

 

0

x0

 

 

 

 

 

 

 

 

 

Нахождение частного решения линейной неоднородной системы ОДУ методом вариации постоянных. Формула Коши.

y′ = A(x)y + f (x) (1)

y′ = A(x)y (2) – собственная однородная система.

Будем искать частное решение системы (1) в виде:

 

y = C1

(x)y′+…+Cn (x)yn =Y C(x) (3)

 

 

C

(x)

 

 

 

1

 

 

 

C(x)=

 

,

 

 

Cn

(x)

 

 

 

 

где y1,, yn – ФСР собственной однородной системы (2), а C1(x),Cn (x) – функции класса C1([a,b]). Подставим y из (3) в (1):

(Y C(x))=YC(x)+YC(x) = A(x)(Y C(x))+ f (x) YC(x)= f (x)

C(x)=Y 1(x)f (x)

(Y 1(x) – существует, так как detY (x)=W (x)0 , x [a,b])

C(x)= x Y 1(t)f (t)dt = [Y (x)= R(x, x0 )S, Y 1(t)= (R(t, x0 )S )1 = S 1 R(x0 ,t)]=

x0

= x R(x, x0 )S S 1R(x0 , x)f (t)dt = x R(x,t)f (t)dt .

x0

 

 

 

 

x0

 

y′ = A(x)y +

f (x)

имеет вид:

Итак, решение задачи Коши y(x

)= y

 

 

 

 

0

 

0

 

 

 

y(x)= R(x, x0 )y0 + x R(x,t)f (t)dt

 

 

 

 

 

 

x0

Линейные системы с постоянными коэффициентами и методы их решения Случай нормализуемой системы.

где

qi =

l11(p)y1 +…+ l1n (p)yn = f1

 

(1),

ln1

(p)y1 +…+ lnn (p)yn = fn

 

 

p = dxd , lij (p) – многочлен от p с постоянными коэффициентами из C. Положим

max deglij , N = q1 +…+ qn – порядок системы (формальный).

1in

18

 

l11(p)

l1n (p)

 

L(p)=

,

 

ln1(p)

lnn (p)

Lij (p)

– алгебраическое дополнение элемента lij (p) в матрице L(p). (p)= det L(p)

(дифференциальный оператор, а не число). “Умножим” i-ое уравнение из (1) на Lis (p)

слева и просуммируем по i, от 1 до n. (Если умножать элементы столбца матрицы на соответствующие алгебраические дополнения элементов того же столбца, то получим определитель матрицы, а если умножать на алгебраические дополнения другого столбца, то получим 0, на этом принципе построена обратная матрица.)

n

В итоге получим: (2) (p)ys = fs (x)= Lis (p)fi (x), s =1,n . (В дальнейшем будет

i =1

обосновано, что такое “умножение” возможно при достаточно гладких fi (x)).

Система (2) распадается на n независимых ОДУ, отличающихся только правыми частями. Заметим, что каждое (достаточно гладкое) решение системы (1) (или (1одн.)) является решением системы (2) (или (2одн.)), но не наоборот. Из прошлого семестра заключа-

ем, что общее решение системы (2одн.) имеет вид: λ1,λs – корни характеристического уравнения

y(x)= p1(x)eλ1 x +…+ ps (x)eλs x (3), где

(λ)= 0 , кратности r1,rs , а p j (x)

векторный многочлен, каждая компонента которого является многочленом степени rj 1. Для нахождения общего решения системы (1одн.) можно было бы подставить

y(x) из (3) в (1одн.), найти соотношение между коэффициентами многочленов p1(x),, ps (x), выразить базисные неизвестные через свободные и обозначить эти свободные через C1,C2 ,. Однако это слишком сложно. Для упрощения докажем лемму.

Лемма:

Пусть y(x) – решение системы (1одн.) (как показано выше, y(x) охватывается формулой (3)). Тогда каждое слагаемое в записи (3) также является решением системы (1одн.) (то есть можно привести предыдущие рассмотрения для каждой λj в отдельно-

сти, что значительно проще). Доказательство:

 

 

 

 

 

 

 

 

 

по теореме

◄ Подставим y(x) в (1одн.): L(p)y(x)= L(p)p1(x)eλ1 x +…+ L(p)ps (x)eλs x

=

= eλ1 x L(p + λ1 )

 

(x)+…+ eλs x L(p + λs )

 

(x)= 0 . Поскольку все функции xk eλi x

о сдвиге

 

 

 

p1

ps

линейно

независимы при различных k и различных λj (“ФСР в случае кратных корней” – из про-

шлого семестра), то все коэффициенты всех многочленов в этом соотношении равны нулю. То есть

 

 

 

 

 

по теореме

 

 

 

 

 

L(p + λj )p j (x)0 L(p)(p j (x)eλj x )

=

eλj x L(p + λj )p j (x)0 p j (x)eλj x – реше-

 

 

 

 

 

о сдвиге

 

 

 

 

 

ние системы (1одн.). ►

Таким образом, общая схема решения системы (1одн.) изложена. Осталось ее обосно-

вать. Для этого положим:

l

(p)= a pq j +

младшие степени p, некоторые a

могут рав-

 

 

 

 

 

 

ij

 

 

ij

 

 

 

ij

 

няться нулю.

 

 

 

 

 

 

l11(p)

l1n (p)

 

 

 

ω =

 

a11

a1n

 

,

(p)=

 

 

=ω pN + младшие степени

p

 

 

 

 

 

 

an1

ann

 

 

 

 

 

ln1(p)

lnn (p)

 

(N =q1 +…+qn )

 

 

 

 

 

(p)=

 

 

 

 

 

 

 

 

 

 

 

ω pN + + младшие

 

степени p .

 

19

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