Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
BORZENKOV_MATLAB ОДУ.pdf
Скачиваний:
184
Добавлен:
11.05.2015
Размер:
2 Мб
Скачать
☆

Предположим, что в задаче (1.11) коэффициент Q(t, ) 0, а P от не зависит: P(t, ) = P(t). Причем P(t) непрерывен при t 0; параметром считается y0 . Докажем, что если P(t) 0, то y(t, y0) непрерывна по y0 равномерно относительно t при t 0. Действительно, из формулы (1.12) имеем

t

~

 

 

 

 

 

 

 

 

~

 

 

 

 

 

 

 

 

P(t )dt

 

 

 

 

 

 

 

 

y(t,y0) y0 e0

.

Отсюда

сразу

 

получаем

следующее

представление:

 

 

 

t

~ ~

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P(t )dt

 

 

 

 

y y(t,y0 y0) y(t,y0) y0 e0

.

Поэтому, если

P(t) 0,

то

| y | | y0 | и, таким образом,

| y |

при | y0

| ( ) для всех t 0,

что

и требовалось доказать.

 

то утверждение о непрерывности y(t, y0 )

Докажем, что если P(t) 0,

по y0 , равномерной

относительно

t

при

t 0, неверно.

Рассуждаем

от

противного. Допустим, что имеет место утверждение, доказанное ранее, т.е. для

0

( ) такое, что

| y | для всех

t 0, если

| y0 | .

Возьмем

| y

0

 

1

 

 

 

t

~

~

 

 

 

 

 

 

|

 

 

( ). В силу условия P(t) 0 имеем exp P(t )dt exp( t)

и тем

 

2

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~ ~

при достаточно большом t становится больше

самым экспонента exp P(t )dt

 

 

 

 

 

0

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~ ~

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P(t )dt

 

 

 

 

 

 

 

4

. Но тогда из представления

y y0 e0

получаем

| y |

4

2 ,

 

 

 

 

 

 

 

 

 

 

 

y | .

 

 

2

 

 

 

что противоречит неравенству |

 

 

 

 

 

 

 

 

Замечание. Данный пример показывает, что теорема 2 не распространяется на бесконечный промежуток изменения t и для ее справедливости нужны еще некоторые дополнительные требования, например, условие на знак P(t).

§2. Дифференцирование по параметрам и начальным значениям

Кратко исследуем возможность дифференцирования решений по параметрам и начальным значениям. Не ограничивая общности, достаточно рассмотреть этот вопрос для какой-либо одной переменной . Запишем систему в виде

dyi

= fi(y(t, ),t, ) (i = 1, …, m).

(2.1)

 

dt

 

Будем считать, что выполнены условия теоремы 1. Тем самым решение начальной задачи для системы (2.1) существует и является непрерывной

функцией параметра при

| 0 |

c . Построим конечно-разностные

отношения – функции zi (t, ):

 

 

55

 

 

 

 

zi (t, ) =

yi (t, 0 ) yi (t, 0)

,

(2.2)

 

 

 

 

 

 

 

 

 

 

 

которые являются решениями следующей системы:

 

 

dzi

=

1

{ fi(y(t, 0

),t, 0 ) fi(y(t, 0),t, 0)}.

(2.3)

 

 

 

 

dt

 

 

 

 

Из (2.2) получаем следующее представление:

 

 

yi(t, 0 ) yi(t, ) zi(t, ) .

(2.4)

Предположим, что дополнительно к условиям теоремы 1 функции

fi (y,t, ) в

~

обладают непрерывными частными производными по

. Тогда

области D

пользуясь представлением (2.4) и тождеством Адамара, запишем (2.4) в виде

 

 

 

 

 

 

dzi

 

=

a

 

(t, )z

 

(t, ) (t, ),

(2.5)

 

 

 

 

 

 

 

 

 

 

m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ik

 

k

i

 

 

где

 

 

 

 

 

dt

 

 

k 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

aik (t, ) =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 f

i

(y1(t, 0), …,

 

 

 

 

 

 

 

 

 

 

ym (t, 0 ), 0 ) d ,

 

=

 

yk (t, 0) zk , …,

(2.6)

yk

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 f

i

 

 

 

 

 

 

 

 

 

 

 

i (t, )

 

 

(y1(t, 0),...,ym(t, 0), 0

)d .

(2.7)

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

fi

 

 

 

 

 

 

 

 

 

 

 

 

 

Здесь

 

равны нулю, если

 

 

является одним из

yi0 и отличны, вообще

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

говоря, от нуля, если является одним из k . Начальные условия для функций zi (t, ) также имеют различный вид в зависимости от того, является ли

параметр 0 каким-либо из начальных значений задачи или нет. Если 0

 

yj0

и поскольку

yi (t0, 0) yi0 ,

yi (t0, 0 ) yi0 , то в силу (2.4)

zi (t, ) 0.

Если же 0 =

yj0, то yi (t0, 0

) ( 0

) ij

и zi (t, ) ij ( ij

– символ

Кронекера).

 

 

 

 

 

 

 

 

 

 

 

 

 

Из (2.6), (2.7) видно, что ai,k (t, )

и i (t, )

– непрерывные функции t

и при | t t0

| H , | | c . Действительно, aik и i зависят от t и как

сложные функции: и непосредственно, и

через

y(t, 0 ).

Но в

силу

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

fi и

доказанной

выше

непрерывности

y(t, 0 ) по t и

(непрерывность

y(t, )

по

t и

0

эквивалентна непрерывности

y(t, 0

)

по

t

и

)

эти

сложные функции будут также непрерывными по t и .

 

 

 

 

 

 

Поэтому правые части (2.5) удовлетворяют условиям теоремы 1, из которой

следует, что

zi (t, ) – непрерывная функция

t

и

при

| t t0

| H ,

| | c . Это означает, в частности, что существуют предельные значения

56

limzi (t, )

zi (t,0) lim

 

yi (t, 0 ) yi (t, 0 )

 

yi

 

,

(2.8)

 

 

 

0

 

 

 

 

 

0

 

 

 

0

 

 

 

 

yi

 

 

 

 

 

 

 

 

 

 

 

т. е. производные

 

 

при 0. Эти производные удовлетворяют (2.5), где

 

 

aik aik (t,0),

 

 

 

 

 

 

 

 

 

 

 

 

 

i i (t,0). Из

(2.6), (2.7) видно,

что

 

эти

величины

представляют собой

 

 

fi

(y(t, 1,..., k ,y10,...,ym0),t, 1,..., s) и

fi

() при

 

 

 

 

 

 

 

 

 

 

 

 

yk

 

 

 

 

 

 

 

0. Поэтому, пользуясь той же теоремой 1, можно сделать заключение о

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

значений.

Теорема

 

3.

Если

функции

fi (y,t, )

непрерывны

вместе

с

частными

производными

по

y1,..., ym

1,..., s

в

~

то

существуют

 

производные от

D,

 

решения задачи (1.1), (1.2) по начальным значениям

y10,..., ym0

и параметрам

,...,

s

,

непрерывные при

|t t

c

| H ,

| y

i0

y0

|

i

|,

|

k

 

k0

| c

1

 

 

 

 

 

 

 

 

 

i0

 

 

 

 

 

k

(i 1,...,m;

k 1,...,s).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Замечание. В теореме сформулированы условия существования первых непрерывных производных по параметрам решения задачи Коши. Вопрос о существовании и непрерывности производных высших порядков по параметрам исследуется аналогично. Можно показать, что существование непрерывных частных производных до порядка k функций fi (y,t, ) является достаточным условием существования непрерывных частных производных

k -го порядка по параметрам y10,..., ym0 ,

1,..., s

решения задачи (1.1), (1.2).

 

 

 

 

 

 

 

 

 

k y

(i 1,...,m;

Существование непрерывных частных производных

 

i

 

k

 

 

 

 

 

 

 

 

 

 

k 1,...,r 1) позволяет при достаточно

малых 0

искать решение

задачи (1.1), (1.2) в виде асимптотического разложения

 

 

 

r

k

yi

 

 

( )

k

 

 

 

yi(t, ) yi(t, 0)

 

 

(t, 0)

 

O[( )r 1].

(2.9)

 

 

 

k!

 

k 1 k

 

 

 

 

 

 

 

Идея разложения решения по малому параметру , от которого правые части зависят регулярно, принадлежит выдающемуся французскому математику Анри Пуанкаре. На формуле (2.9) основаны многие методы асимптотического исследования и численного интегрирования дифференциальных уравнений, а также некоторые методы теории устойчивости.

57

§3. Уравнения в вариациях. Линеаризация

 

 

 

 

 

Предельные при 0

 

значения функций zi (t, ) – функции zi (t,0) =

 

yi

 

 

– являются решениями задачи Коши (для определенности рассмотрим

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

fi и равно yj0 ):

 

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

 

 

 

 

 

 

dz

i

 

m

f

i

 

 

 

 

 

 

 

 

 

 

 

 

(y,t)zk ,

zi (t0,0) ij ,

(3.1)

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

k 1

yk

перехода при 0

 

которую

получим в результате предельного

в (2.5).

Нетрудно видеть, эта же система может быть получена путем формального дифференцирования исходной задачи (1.1), (1.2) по параметру . Уравнения (3.1) часто называются системой уравнений в вариациях относительно рассматриваемого решения y(t, ).

Переход к уравнениям в вариациях связан с идеей линеаризации уравнений (1.1) в окрестности некоторого выбранного решения. Пусть, например, свойства

какого-то частного решения

yi

(i 1,...,m) нам известны и нужно исследовать

другое решение

 

 

 

 

 

этого же уравнения, близкое к yi

. Так, например,

yi

 

(i 1,...,m)

yi

и

 

 

 

 

могут различаться по начальным значениям и т.п.

 

Введем величины

 

yi

yi

 

 

 

yi (i 1,...,m). Предполагая, что правые части уравнений достаточно

yi

гладки, подставим

 

 

yi yi

в (2.1) и разложим fi

по степеням yi . Тогда

 

yi

получается уравнение вида

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d

 

 

 

m

 

 

 

 

fi

 

 

 

 

 

 

 

 

 

 

 

( yi) aik (t) yk i( y1 , …, ym,t),

aik (t)

(y,t).

(3.2)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

k 1

 

 

 

 

yk

 

 

Если параметр, значениями которого различаются

yi

и

 

,

явно не

 

yi

входит в

fi , то i

o( | y| ),

получаем линейную систему уравнений

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d

m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( yi) aik (t) yk ,

 

 

 

 

 

 

 

(3.3)

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k 1

 

 

 

 

 

 

 

 

которая представляет собой систему в вариациях, так как совпадает с точностью до обозначений с (3.1). Если рассмотрение ведется на конечном

промежутке [t0,T ], то в силу доказанных выше теорем yi и yi близки на всем [t0,T ] (так как их начальные значения или правые части уравнений отличаются

мало) и, зная yi , можно получить приближенное выражение для yi в виде yi yi , где yi определяется из линейных уравнений (3.4).

Если же, зная yi , мы хотим исследовать yi на неограниченном промежутке, то задача становится сложнее и относится к теории устойчивости. Линейные уравнения изучены значительно полнее, чем системы общего вида. Поэтому процесс линеаризации – замена (2.1) уравнениями вида (3.4) – приводит в тех случаях, когда эта операция законна, к серьезному упрощению задачи.

58

Замечание. В теореме 3 рассматривался вопрос о существовании производных решения по параметрам y10,...,ym0 , 1,..., s, входящем в начальные условия и правые части системы (1.1), (1.2). Аналогично может быть исследован вопрос

и о существовании производных

 

dyi

. Эти производные также являются

 

 

 

 

 

 

 

dt0

 

 

 

 

 

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

 

 

d

 

yi

 

m

fi

 

yk

 

 

 

 

 

 

(y,t)

,

(3.4)

 

 

 

 

 

 

 

 

t0

 

yk

t0

 

dt

 

k 1

 

 

 

получающейся дифференцированием системы (1.1) по t0 и является линейной

системой

относительно

yi

.

Чтобы получить

начальные

условия задачи,

 

 

 

t0

 

 

 

заменим

исходную систему

дифференциальных

уравнений

(1.1) системой

 

 

 

 

t

 

 

интегральных уравнений yi(t) yi0 fi(y, )d . Если подставить решение y(t)

 

 

 

 

t0

 

 

исходной системы, то мы получим тождества. Дифференцируя тождества по t0,

 

y

 

 

t

 

m f

i

 

y

k

 

 

 

 

 

получим

i

 

 

 

 

 

 

 

 

 

 

 

, откуда при

t t

 

в силу (1.2)

 

t0

f (y(t0),t0)

 

 

 

 

 

 

 

 

d

 

 

0

 

 

 

 

t

k 1 yk

 

t0

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

будем иметь

yi

 

 

fi (

y10, …, ym0 , t0). Эти условия и представляют собой

 

 

t0

 

 

 

 

t t0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

начальные данные для системы (3.4).

ИСТОРИЧЕСКАЯ СПРАВКА

Красовский Николай Николаевич. Выдающийся математик.

Ученик Е. А. Барбашина, И. Г. Малкина. Родился 7 сентября 1924 г. в Свердловске. В 1949 г. окончил Уральский политехнический институт (УПИ) (Свердловск). Специалист в области математики и механики. Избран членомкорреспондентом по Отделению механики и процессов управления (механика) в 1964 г., академиком по Отделению механики и процессов управления (теория устойчивости и регулирования) в 1968 г. Награжден золотой медалью им. А. М. Ляпунова за цикл статей «Работы по теории устойчивости и по теории оптимального управления» (1992 г.); золотой медалью им. М. В. Ломоносова за выдающиеся достижения в области математической теории управления и теории

дифференциальных игр (1996 г.). Основатель крупной научной школы. В настоящее время является главным научным сотрудником отдела динамических систем Уральского филиала института математики и механики РАН, членом президиума Российской Академии Наук. Среди его учеников академики: президент РАН Ю. С. Осипов, А. Б. Куржанский, А. И. Субботин, члены-корреспонденты РАН А. Г. Ченцов, В. Е. Третьяков и многие другие члены-корреспонденты РАН, доктора и кандидаты наук, инженеры и преподаватели.

59

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