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

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

В этом параграфе рассмотрим Задачу Коши для линейной

системы с переменными коэффициентами:

 

y0 = A(t)y + f(t); t [¡T; T ]; 0 < T < 1;

 

(y(0) = y0 2 CN (или2RN ):

(1)

Здесь

A = (aij(t)); i; j = 1; N - матрица с непрерывными коэффициен-

тами aij0(t); t 21[¡T; T ]; f1(t)

f(t) = Bf.

N (t)C

- вектор-функция правых частей с непрерывны-

 

 

@

A

 

 

 

ми компонентами fi(t); i =

1; N

на отрезке [¡T; T ];

 

 

ByN (t)C

- вектор-функция искомых переменных,

y(t) = 0.y1(t) 1

 

 

@

A

 

 

 

y0

 

ByN0C

 

 

 

=

0.y10

1 - вектор начальных данных.

 

 

@

A

 

 

 

В связи с вышеизложенным существуют постоянные M; N > 0, такие, что

jjA(t)jj · M; jjf(t)jj · N на отрезке [¡T; T ]:

Предположим, что Задача Коши (1) имеет непрерывное и непрерывно-

1

Лекция №9, НГУ, ММФ, 2009

2

дифференцируемое решение y = y(t). Тогда

dtd [(y(t); y(t))] = dtd £jjy(t)jj2¤ = (y0(t); y(t)) + (y(t); y0(t)) =

=(A(t)y(t) + f(t); y(t)) + (y(t); A(t)y(t) + f(t)) =

=(A(t)y(t); y(t)) + (y(t); A(t)y(t)) + (f(t); y(t)) + (y(t); f(t)) =

= (B(t)y(t); y(t)) + 2Re(f(t); y(t)); где B(t) = A(t) + A¤(t) = B¤(t),

2Re(f(t); y(t)) = (f(t); y(t)) + (f(t); y(t)). p

Поскольку jjBjj = ¸max(B¤B) = maxfj¸min(B)j; j¸max(B)jg (см. Упражнение 1 к этому параграфу),

jjBjj = jjA + A¤jj · jjAjj + jjA¤jj = 2jjAjj

(см. Упражнение 2 к этому параграфу),

2Re(f; y) · 2j(f; y)j · 2jjfjj ¢ jjyjj · jjf(t)jj2 + jjy(t)jj;

то

dtd Z2 · jjBjjZ2 + L2 + Z2 · C1(M)Z2 + N2;

где Z = jjy(t)jj, L = jjf(t)jj, C1(M) = 2M + 1.

Умножив обе части полученного неравенства на eC1t и интегрируя полученное выражение (полагая, что 0 · t · T ), в итоге получим:

Z2

·

eC1tZ2

+

N

2

eC1t ¡ 1

; Z

 

= y

 

0

 

 

C1

0

jj 0jj

или более грубое неравенство:

Z2

·

eC1T Z2

+

N

2

eC1T ¡ 1

; 0

·

t

·

T:

(2)

 

0

 

 

C1

 

 

 

Используя прием из §2 с заменой ¿ = ¡t, получим аналогичную оценку и при ¡T · t · 0. Итак, при всех t 2 [¡T; T ]

Z2

·

eC1jtjZ2

+

N

2

eC1jtj ¡ 1

·

eC1T Z2

+

N

2

eC1T ¡ 1

;

t

j ·

T: (3)

C1

 

 

0

 

 

0

 

 

C1

j

 

Лекция №9, НГУ, ММФ, 2009

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

Легко видеть, что в случае Задачи Коши

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y0 = A(t)y + f(t); t · T;

 

 

 

 

 

 

(10)

 

 

(y(t0) = y0; t0

 

 

 

 

( T;j j

 

T )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 ¡

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

оценка (3) перепишется так:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Z2

·

eC1jt¡t0jZ2

+

N

2 eC1jt¡t0j ¡ 1

·

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

C1

 

 

 

 

 

 

 

 

·

e2C1T Z2

+

N

2

e2C1T ¡ 1

;

 

 

 

t

j ·

T:

 

 

 

(30)

 

0

 

 

 

 

C1

 

 

 

 

 

j

 

 

 

 

 

 

 

Поскольку p

 

· a + b при a; b > 0, то из (3); (30) следует:

a2 + b2

Z = jjy(t)jj · e21 C1jtjZ0 + Ns

 

 

 

 

 

 

 

 

 

 

 

 

 

e

C1

jCj1¡ 1 ·

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

· e21 C1T Z0 + Ns

e

C1T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

¡ 1

 

 

 

 

 

 

 

(4)

 

 

 

 

 

 

 

 

 

 

 

и

 

 

 

C1

 

 

 

 

 

 

 

 

 

 

 

+ Ns

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Z =· e2 C1jt¡t0jZ0

 

eC1

j

 

 

C1

¡

 

 

 

·

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t¡t0j

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

· eC1T Z0 + Ns

 

 

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

e2C1 ¡ 1

:

 

 

 

 

 

 

(40)

 

 

C1

 

 

 

 

 

 

Как и в §2, покажем, как с помощью полученных оценок доказать единственность решения Задачи Коши (1) (если оно существует). Предположим противное: одним и тем же векторам y0; f(t) отвечают два решения yI(t); yII(t). Их разность ¢(t) = yI(t) ¡ yII(t) есть решение следующей Задачи Коши:

(¢0 = A(t; jtj · T;

(5)

¢(0) = 0:

Лекция №9, НГУ, ММФ, 2009

4

Из оценки (4) для Задачи Коши (5) имеем:

 

jj¢(t)jj · e21 C1T ¢ 0 + 0 ¢ s

 

 

 

= 0;

 

e

C1¡ 1

 

 

 

C1T

 

т.е. yI(t) ´ yII(t) на отрезке [¡T; T ].

Рассмотрим теперь Задачу Коши для однородной системы (т.е.

f(t) ´ 0 при jtj · T ):

 

 

 

y0 = A(t)y;

t

T;

 

(y(t0) = y0:

j j ·

 

(10)

Предположим, что Задача Коши (10) однозначно разрешима при любом векторе y0. Тогда, как и в случае линейных систем с постоянными коэффициентами, можно утверждать, что:

а) Все решения линейной системы y0 = A(t)y образуют линейное пространство.

б) Размерность этого пространства равна N, т.е. размерности пространства начальных векторов y0.

в) N линейно-независимых решений системы y0 = A(t)y можно построить, решив N раз Задачу Коши (10) и взяв в качестве начальных данных N линейно-независимых векторов y0.

Пусть y[k](t); k = 1; N - какая-либо совокупность линейно-независимых

решений для системы y0 = A(t)y, при этом:

 

 

 

8

³[k]

 

 

´

0

 

 

 

 

 

 

 

y[k]

t

 

= A(t)y[k](t); t

T;

 

 

 

 

( )

 

 

[k]

;

j j ·

 

 

 

<y

 

(0) = y0

 

 

где

0

:

 

 

 

 

0

 

 

1.

 

y[k](t) =

.y1k(t) 1, y0[k] =

.y1k0

 

 

ByNk(t)C

 

 

 

ByNk0C

 

 

@

A

 

 

[k]

@

 

 

A

 

Составим из векторов y

 

(t) матрицу Y (t): Y (t) = (yij(t)); i; j =

1; N. Ясно, что как и в случае линейных систем с постоянными коэффициентами, если det Y (0) 6= 0, то и det Y (t) 6= 0, jtj · T . Точно

Y 0(t) = A(t)Y (t); t 2 [¡T; T ]; Y (0) = Y0; det Y0 =6 0:

Лекция №9, НГУ, ММФ, 2009

5

так же можно утверждать, что если det Y (0) = 0, то det Y (t) = 0 при всех t 2 [¡T; T ]. Очевидно, что матрицу Y (t); det Y (t) =6 0 можно назвать фундаментальной матрицей решений для системы y0 = A(t)y. Любое решение Задачи Коши (10) выражается так:

XN

y(t) = Y (t)C = Cky[k](t);

k=1

0C1 1

где C = B@. CA, при этом вектор C следующим образом выража-

CN

ется через вектор начальных данных y0: C = Y ¡1(0)y0, т.е. решение Задачи Коши (10) записывается так:

y(t) = Y (t)Y ¡1(0)y0:

(6)

Сама же фундаментальная матрица решений системы y0 = A(t)y

есть решение следующей Задачи Коши для матричной системы

(

(7)

Как и в случае линейной системы с постоянными коэффициентами справедлива формула (типа формулы (8) из §3), связывающая det Y (t) и det Y0. В самом деле, если мы снова вычислим производную от ¢(t) = det Y (t) (см. §3), то получим:

 

d¢(t)

= Ã

N

akk(t)!¢(t) = T rA(t)¢(t):

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

=1

 

 

 

 

 

 

 

 

Xk

 

 

 

 

 

 

 

 

 

 

t

T rA(s)ds

 

 

 

 

 

 

 

 

¡

Умножая это соотношение на e

0

, мы получим:

R

 

 

 

 

 

 

 

t

 

 

 

 

 

 

d

 

¡ T rA(s)ds

 

 

 

 

 

 

 

0

 

 

 

 

 

dt (e

R

 

¢(t)) = 0;

Лекция №9, НГУ, ММФ, 2009

6

т.е.

t

 

 

¢(t) = eR

T rA(s)ds

(8)

¢(0);

0

 

 

где ¢(0) = det Y0. Формула (8) носит название формулы Лиувилля. Она обобщает аналогичную формулу, полученную для систем с постоянными коэффициентами.

Вернемся к формулам (6). Поскольку (см. неравенство (4)):

Z = jjy(t)jj · e12 C1(M)jtjZ0;

то

и

jjY (t)Y ¡1(0)y0jj · e21 C1jtjjjy0jj

 

 

 

 

 

jjY (t)Y ¡1(0)y0jj

 

 

 

Y (t)Y ¡1

(0) =

sup

·

e21 C1jtj:

(9)

jj

jj

Z0=jjy0jj6=0

Z0

 

 

 

 

 

 

 

 

Используя фундаментальную матрицу Y (t), получим формулу для решения Задачи Коши (1). Умножим систему y0 = A(t)y + f(t) слева на Y ¡1(t). В итоге получим (см. Упражнение 3 к этому параграфу):

т.е.

©Y ¡1(t)y(t)ª0 = Y ¡1(t)f(t);

 

 

 

y(t) = Y (t)Y ¡1(0)y0 + Z0 t

 

Y (t)Y ¡1(s)f(s)ds:

(10)

Поскольку (см. формулу (9))

 

 

 

 

 

 

 

 

 

 

 

 

jjY (t)Y ¡1(s)jj · e21 C1jt¡sj;

 

 

 

 

 

то из (10) следует:

 

 

e21 C1jtjZ0 + ¯

 

 

e21 C1jt¡sj

 

 

f(s) ds¯

 

Z = y(t)

 

 

t

 

 

 

 

 

 

 

 

¯0

 

 

1

 

 

 

 

 

¯

 

jj

jj ·

 

 

¯

 

 

 

jj

 

jj

¯

 

 

 

¯Z

 

 

 

 

 

¯ ·

 

 

 

e2 C1jtjZ + 2¯

e

 

j j ¡ 1

:

 

¯

 

 

 

 

1

 

¯

 

 

2 C1

t

 

 

 

 

¯

 

 

·

 

 

0

N

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C1

 

 

 

 

 

 

Лекция №9, НГУ, ММФ, 2009

7

Это еще один вариант априорной оценки для решения Задачи Коши (1).

Кратко остановимся на случае одного линейного неоднородного уравнения высокого порядка с переменными коэф-

фициентами:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

8

Lx

 

x(N)

 

a

 

t x(N

¡

1)

 

:::

 

 

a

t x

0

x

 

F t ;

 

 

=

 

+

 

1

( )

 

+

 

 

+

 

 

1( )

 

+ aN (t)t

=

T;( )

(11)

>

 

 

 

 

 

 

 

 

 

 

 

(N

¡

1)

(0) = ®

 

 

j j ·

 

 

<x(0) = ® ; x0(0) = ® ; :::; x

 

 

 

 

 

;

 

 

 

:

i

 

1

 

 

 

2

 

 

 

 

 

 

 

 

 

N

 

 

 

 

>

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где a (t); i = 1; N; F (t) - непрерывные функции от t на отрезке [¡T; T ]. Как и в случае линейного уравнения высокого порядка с постоянными коэффициентами, Задача Коши (11) может быть сведена к Задаче Коши вида (1):

 

8

d

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y(t) = A(t)y(t) + f(t); jtj · T;

 

 

 

 

 

dt

 

 

 

(100)

 

>

 

 

 

 

 

 

®1

 

 

 

 

 

 

 

 

>

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

>

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

>

 

 

 

 

=

 

 

 

 

;

 

 

 

 

 

 

>y(0) = y0

0

 

1

 

 

 

 

 

 

<

 

 

 

 

 

.

 

 

 

 

 

 

 

>

 

 

 

 

 

B

®N C

 

 

 

 

 

 

 

>

 

 

 

 

 

@

 

A

 

 

 

 

 

 

 

>

 

 

 

 

 

 

 

 

 

 

 

 

где

>

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

>

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

:

 

1

 

 

 

0

 

 

 

 

 

 

1

 

 

x t

 

 

 

 

 

 

0

1 0 : : :

 

0

 

y(t) =

0x0((t))

 

; A(t) =

 

 

0

0 1 : : :

 

0

;

B.

 

 

C

B: : :0: : : : : :0: :0: : ::::::: : : : :1: : :C

 

B

 

 

C

 

 

 

B

 

 

 

 

 

 

C

 

 

Bx(1)

(t)C

 

 

 

B

 

aN (t) : : : : : : : :

 

a1

C

 

 

B

 

 

C

 

 

 

B

 

¡

(t)C

 

 

@

 

 

A

 

 

 

@¡

 

 

 

 

A

 

f(t) =

0

.0

1

:

 

BF (t)C

 

 

@

 

A

 

Понятно, что если Задача Коши (100) однозначно разрешима, то однозначно разрешима Задача Коши (11) (в силу эквивалентности этих задач). При F (t) ´ 0 (т.е. f(t) ´ 0) фундаментальная

0
где W (0) = det ©(0). Ясно, что решение Задачи Коши (11) при F (t) =6 0; jtj · T может быть найдено с помощью формулы (10).
8
< Zt
W (t) = W (0)exp :¡

Лекция №9, НГУ, ММФ, 2009

8

матрица решений для системы y0

= A(t)y такова:

 

 

 

©(t) = (y

 

(t); :::; y

 

(t)) =

B

x1

1)

: : :

xN

1)C

;

 

 

(N.

 

(N.

 

[1]

 

[N]

 

 

0

x10

 

: : :

xN0

1

 

 

 

 

 

 

 

Bx

¡

 

: : : x

¡

C

 

 

 

 

 

 

 

B

1

 

 

N

C

 

 

 

 

 

 

 

@

 

 

 

 

A

 

при этом любое решение уравнения Lx = 0 представимо так:

XN

x(t) = Ckxk(t):

k=1

Заметим, что вронскиан W (t) = det ©(t) =6 0, причем по формуле Лиувилля 9

=

a1(s)ds;; (12)

Лекция №9, НГУ, ММФ, 2009

9

Упражнения к §9

1. Доказать, что для эрмитовой матрицы A:

jjAjj = maxfj¸min(A)j; j¸max(A)jg:

2. Дана система y0 = A(t)y; t 2 [0; T ] с непрерывными коэффици-

ентами на отрезке [0; T ]. Пусть ¸max(B) · 0 на отрезке [0; T ], где B = A + A¤. Доказать, что для решения Задачи Коши (10) (если

оно существует): (y0 = A(t)y; t 2 (0; T ];

y(0) = y0

справедлива априорная оценка:

Z = jjy(t)jj · Z0 = jjy0jj; t 2 (0; T ]:

3. Рассмотрим следующие задачи:

(X0(t) = A(t)X(t); jtj · T; X(0) = IN

и

(Y 0(t) = ¡Y (t)A(t); jtj · T;

 

 

Y (0) = IN :

Показать, что Y (t) = X¡1(t).

4. Пусть любое решение Задачи Коши

(y0 = By; t 2 R1; y(0) = y0

где B - матрица с постоянными коэффициентами, таково, что jjy(t)jj = jjy0jj. Показать, что это справедливо тогда и только тогда, когда

B¤ = ¡B.

 

B, то

eB

¤

= e¡B =

eB

¡

1

, т.е. eB -

Замечание. Если B¤ =

 

 

унитарная матрица. ¤

¡

 

¡ ¢

 

 

¡ ¢

 

 

 

Соседние файлы в папке Диф.уравнения