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

Лекции по дифференциальным уравнениям (Абрамов А.А

.).pdf
Скачиваний:
142
Добавлен:
12.03.2016
Размер:
3.56 Mб
Скачать

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

~

~

åñëè j 1

2j < . Таким образом, обе функции являются "-приближ¼нными. Одна из них

является точным решением. Вс¼ происходит в замкнутом выпуклом множестве, функции удовлетворяют одним и тем же начальным условиям, поэтому

 

 

 

 

 

~

~

 

6 2 "e

K

 

 

 

 

 

 

 

 

 

 

 

 

 

j~y(t; 1) ~y(t; 2)j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~

~

6

Замечание: Непрерывность следует из неравенства треугольника: j~y(t1; 1)

~y(t2; 2)j

~

~

 

~

~

2 "e

L

+ Mjt1 t2j Теорема доказана.

 

 

j~y(t1; 1)

~y(t1; 2)j + j~y(t1

; 2)

~y(t2; 2)j 6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~

 

 

 

2. Теорема: Пусть в дополнение к условиям предыдущей теоремы

@f

непрерывна в

 

~

 

 

 

G.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

@

 

 

 

 

~

 

 

 

 

 

 

 

 

 

 

 

 

@~y

 

 

 

 

 

Тогда функция ~y(t; ) имеет непрерывные производные V =

 

и функция V удовлетворяет

~

уравнению

 

 

 

 

 

 

 

 

 

 

 

 

 

@

 

 

 

 

 

 

dt =

 

@~y t; ~y(t; ~ ); ~ !V

+ @~ t; ~y(t; ~ ); ~

 

 

 

 

(4)

 

 

 

 

 

 

 

 

 

dV

 

~

 

 

~

 

 

 

 

 

 

 

 

 

 

 

@f

 

 

 

 

@f

 

 

 

 

 

 

 

 

 

и начальному условию

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

V (t0) = 0

 

 

 

 

 

 

 

(5)

Замечание: Уравнение (4) называется уравнением вариации для уравнения (1). Оно получа-

~

ется из (1) формальным дифференцированием по . (Это всего лишь способ запомнить это

уравнение. Вообще, никто не давал разрешения своевольно переставлять порядок дифференцирования.)

Замечание: V прямоугольная матрица частных производных, которая имеет m столбцов и n строк. (n m).

Доказательство: Рассмотрим уравнение (4), с начальными условиями (5). Это уравнение,

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

~

^ (5)

V (t; ) решение, единственное задачи (4)

То, что это решение распространяется на весь отрезок [t0 ; t0 + ], будет доказано позднее, при изучении линейных уравнений с переменными коэффициентами.

~ ~

 

 

 

 

 

 

 

 

 

 

 

 

 

 

@~y

jt t0j 6 ; j 0j 6 . Хотим доказать, что эта функция является частной производной

~

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

@

Áåð¼ì ~ ~

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

; , и составляем функцию

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~ ~

def

~

~

~

 

~

 

~

 

 

 

 

 

 

 

 

~z(t; ; ) = ~y(t; + ) ~y(t; ) V (t; )

 

 

 

 

 

~

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Мы докажем, что jj~zjj = o(jj jj). Здесь и далее o( ); O( ) равномерные в Q оценки.

 

 

~

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Если доказать, что ~z = o(jj jj), то вс¼ будет следовать.

 

 

 

 

 

 

 

 

 

 

Напишем дифференциальное уравнения для ~z.

 

 

 

 

 

 

 

 

 

 

 

 

 

d~z

 

 

 

 

 

~

 

 

 

~

 

 

 

 

 

 

 

 

 

@f

 

 

 

@f

 

 

 

 

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

~ + ~

f~ t; ~y(t; ~); ~

 

 

 

V ~

 

 

~

 

 

 

dt

@~y

@~

 

 

 

 

 

 

 

 

|

 

 

{z

 

 

}

 

 

 

 

 

 

 

 

 

 

~

f

72

ГЛАВА 3. ИССЛЕДОВАНИЕ ЗАДАЧИ КОШИ

Оценим эту величину. Отметим, что из соотношения (3), использованного при доказательстве предыдущей теоремы, следует, что " = O( ). Дело в том, что мы предположили, что функ-

ция дифференцируема. Тогда разность в неравенстве (3) не превосходит разницы частных производных, умноженной на ?

Поэтому (см. доказательство предыдущей теоремы)

~ ~ ~ j ~ j ~y(t; ) ~y(t; ) = O( )

Преобразуем:

 

 

d~z

~

 

 

~

 

~

~

~

 

 

 

 

 

@f

 

@f

@f

 

 

 

 

 

 

= f

 

 

 

~y

 

 

+

 

 

~z

 

 

 

 

dt

@~y

~

@~y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

@

 

 

 

 

 

Поскольку ~

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f имеет непрерывные частные производные,

 

 

 

j ~ j

f~ @~y

~t @~ ~ = o

j ~yj + j ~ j = o

 

~

 

~

 

 

 

 

 

 

 

 

 

 

 

@f

 

 

@f

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

def

 

@f~

~ ~

 

 

 

 

 

 

 

 

A(t)

=

 

 

(t; ~y(t; ); )

 

 

 

 

 

 

 

 

 

@~y

 

 

 

 

 

d~z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~

 

 

dt

= A(t)~z + ~(t);

ãäå j~(t)j = o(j j)

 

~

~z(t0) = 0.

j ~ j

Надо доказать, что ~z = o( ). Введ¼м вспомогательную функцию

def ~

'~(t) = 0

Эта функция, подставленная в наше дифференциальное уравнение

да¼т погрешность

Теорема доказана.

Осталось доказать, что

d dt

A(t)'~(t) ~(t) = j~(t)j = o(j ~ j);

 

'~(t)

 

 

 

 

 

 

 

 

 

~

K

~

j~z(t) '~(t)j 6 2o( ) e

 

= o(j j)

j j j ~ j

~z(t) = o( ):

dV

dt (t0) = 0. Так как мы не меняем начальных данных, то так и есть.

3. Зависимость решения от начальных данных.

Замечание: Это очень просто исследовать

Сделаем замену ~y = y~0 + , t = t0 + . Тогда уравнение (1) становится эквивалентным

d~

def

 

 

 

= f~(t0 + ; y~0 + ~; ~ ) = ~g( ; ~; ~ )(t0

; y~0)

(10)

d

5. Уравнения, не разреш¼нные относительно производной. Особые решения.

73

,~

(2)~(0) = 0

Начальные данные фиксировали, а в преобразованные данные стали параметрами.

4.Несколько дополнительных замечаний

Обязательная задача. Для уравнения (стрелочек нет)

y(n) = f(t; y; y0; : : : ; y(n 1))

пользуясь тем, что оно сводится к нормальной системе y0 = z1, z10 = z2; : : : ; zn0 2 = zn 1; zn0 1 = f(t; t; y; z1; : : : ; zn 1) сформулировать (и доказать) результаты, соответству-

ющие теоремам пп.1-3.

Мы привыкли к тому, что система n уравнений имеет общее решение, зависящее от n произвольных постоянных.

[ Тут нет рисунка, нарисуйте его сами ]

t0 фиксированное, y0 = C параметр.

~

~y = ~y(t; C)

Наложив более ж¼сткие ограничения на функцию непрерывность высших производных ~y по ~.

~

f, можно доказать существование и

Замечание: Представьте, что дифференциальное уравнение, описывающее, например, пол¼т самол¼та, разрывно зависит от некоторых параметров (например, веса). Малые изменения параметра влекут большие изменения решения. Поэтому теорема пункта 1 имеет фундаментальное значение для использования в практических целях.

5. Уравнения, не разреш¼нные относительно производной. Особые решения.

1. В этом параграфе мы рассмотрим уравнение (1), где функция f(t; y; p) определена в области Gt;y;p и непрерывна там вместе с первыми частными производными.

f(t; y; y0) = 0

(1)

Напомню геометрическую интерпретацию этого уравнения.

В тр¼хмерном пространстве

(t; y; p) рассматриваем множетво , для которого f(t; y; p) = 0, и рассматриваем кривые `, заданные уравнениями y = y(t); p = p(t).

74

ГЛАВА 3. ИССЛЕДОВАНИЕ ЗАДАЧИ КОШИ

` график решения (1) , ` 2 Q^ dy = p dt íà `.

Задача Коши для (1): (различные формулировки)

1.Íà Q задана точка. Найти то решение уравнения (1), график которого проходит через эту точку.

2.Заданы числа t0; y0; p0, удовлетворяющие условию f(t0; y0; p0) = 0. Найти y(t) то решение уравнения (1), для которого выполняется (2):

y(t0) = y0; y0(t0) = p0

(2)

Замечание: Мы зада¼м три числа, чтобы исключить возможность протыкания прямой поверхности в нескольких точках.

Все величины, которые рассматриваются в этом параграфе, вещественные.

d

Теорема: Пусть f(t0; y0; p0) = 0, dpf(t0; y0; p0) 6= 0. Тогда существует решение задачи (1)^(2). Если y1(t); y2(t) решения задачи (1) ^ (2), то y1(t) = y2(t) в некоторой окрестности точки t0.

[ Вопрос на понимание: окрестность равномерная? ]

Доказательство: В силу наших ограничений по теореме о неявных функция в некоторой окрестности точки t0 соотношение f(t; y; p) = 0 эквивалентно

p = g(t; y);

где g непрерывна вместе с первыми частными производными. Тем самым, в некоторой окрестности нашей точки уравнение f(t; y; y0) = 0 эквивалентно

y0 = g(t; y)

(10)

Начальные данные приобретают следующий вид:

y(t0) = y0;

, y(t0) = y0

(20)

y0(t0) = p0

Тем самым, в некоторой фиксированной окрестности точки (t0; y0; p0) используем теорему 3.

5. Уравнения, не разреш¼нные относительно производной. Особые решения.

75

Замечание: Из доказательства видно, что окрестность фиксированная, то есть не зависит от выбора функции.

2.

Определение: Решение уравнения (1) называется особым, если на каждом интервале графи-

ка этого решения существуют точки, в которых нарушается единственность решения задачи Коши.

Теорема: Пусть y = y(t) особое решение (1). Тогда в каждой точке графика этого решения

имеет место f = 0 и @f@p = 0. (от противного, если бы производная не равнялась нулю, то можно было бы воспользоваться предыдущей теоремой существования и единственности решения задачи Коши).

Определение: Множество точек, в которых f = 0 ^ @f@p , называется p-дискриминантным множеством уравнения p = 0. Стандартное обозначение: . Тогда предыдущую теорему можно

переформулировать так:

Теорема: График особого решения лежит на p-дискриминантном множестве этого уравнения.

Замечание: В типичном случае множество двумерное множество в тр¼хмерном пространстве. p-дискриминантное множество зада¼тся двумя уравнениями, поэтому имеет вид кривой. Кандидатами в особое решение являются какие-либо из этих кусков кривых.

Если кусок кривой является решением, то отсюда не следует, что это решение особое. Вообще говоря, задача формулировки достаточного условия особого решения выходит за рамки программы МФТИ. Поэтому когда мы применяем эти при¼мы к исследованию уравнения, необходимо всегда помнить, что эти требования являются обязательными. Для доказательства надо провести дополнительное исследование.

3. Примеры

1. y0 = p3

y2

, t; y

любые. Этот пример мы решили на самой первой лекции. Решениями

являются y =

 

3 + C

3

; y 0.

 

 

 

 

t

 

76

ГЛАВА 3. ИССЛЕДОВАНИЕ ЗАДАЧИ КОШИ

Этот пример не подходит под нашу теорию. Для того, чтобы он удовлетворял условию, возвед¼м в куб и преобразуем:

(y0)3 y2 = 0; f(t; y; p) = p3 y2:

Если взять y0 6= 0, то решение задачи Коши единственно (локально).

Ïðè y0 = 0 решение задачи Коши не единственно.

 

3p2 = 0

 

,

y = 0

:

p3 y2 = 0

 

 

p = 0

y = 0 особое решение.

5. Уравнения, не разреш¼нные относительно производной. Особые решения.

77

2. (y0)2

= t2, t; y любые. Это уравнение легко решается: y = t, y =

t2

 

+ C, C const

2

Ïðè p0 = t0 решение единственное, хотя через точку проходит два графика решения. (Касательные различны)

Ïðè t0 = 0 через каждую точку проходит 4 решения

Тут имеет место нарушения единственности решения задачи Коши.

p-дискриминантное множество:

 

2p = 0

 

,

t = 0

p2

t2 = 0

 

 

p = 0

Особых решений нет, так как графиком дискриминантного множества является вертикальная прямая, которая не соответствует никакому решению.

3.Пример, который показывает, что все приведенные условия не являются достаточными для выделения особых решений. (y0)2 = 0, y0 = 0, y = C const.

78

ГЛАВА 3. ИССЛЕДОВАНИЕ ЗАДАЧИ КОШИ

Совокупность графиков решений множество горизонтальных прямых. Дискриминантное множество:

:

p2 = 0

2p = 0

Каждая точка плоскости принадлежит дискриминантному множеству. С другой стороны, особых решений нет.

4. Уравнение Клеро.

y = ty0 + f(y0);

где f определена на интервале ( ; ), и f00 непрерывна на ( ; ), f00 > 0.

Замечание: Клод Алексис Клеро. Работал в 1 половине 18 века. Первую математиче- скую работу (исследование геометрических свойств алгебраических кривых 4 порядка) написал в 12 лет. В 16 лет изучил кривые двоякой кривизны. В 18 лет стал адъюнктом Парижской Академии Наук. Вв¼л криволинейные интегралы, понятие полного дифференциала функции нескольких переменных. Кроме того, занимался теорией формы Земли, движения Луны, и т.д.

y = tp + f(p); dy = p dt. Решаем уравнение в переменных t; p.

p dt + t dp + f0(p) = p dt

(p + f0(p)) dp = 0

Надо рассмотреть две возможности:

(a) dp = 0 , p = C const. y = Ct + f(C). Мы получили семейство прямых `(C), зависящих от параметра C. Подстановкой в исходное уравнение убеждаемся, что действительно решение.

(b) На решении t + f0(p) = 0.

t = f0(p); y = f0(p)p + f(p) (параметрическая запись решения.)

Òàê êàê ddpt = f00 < 0, то оно не обращается в 0. По теореме о неявной функции, соотношение t = t(p) , p = p(t).

5. Уравнения, не разреш¼нные относительно производной. Особые решения.

79

Из нашего выражения для t следует

dp

=

1

< 0, y00 =

dy0

=

dp

< 0.

dt

dt= dp

dt

dt

 

 

 

 

 

Обозначим эту кривую . Она выпукла вверх. Мы нашли семейство решений, и ещ¼ одно решение. Как они связаны друг с другом? Докажем, что совокупность этих `(C) совокупность касательных к нашей кривой .

Доказательство: Çàäà¼ì t0 = f0(p0), получаем y0 = f0(p0)p0 + f(p0).

dt p0

= p p0

= p0

 

dy

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

`(C); C = p0; y = p0t + f(p0); t0 = f(p0); y0 = f0(p0)p0 + f(p0); y00 = C = p0:

Действительно, это та же самая точка, касательная в этой точке совпадает с нашей кривой.

:

y tp f(p) = 0

=

 

t f0(p) = 0

 

Общая картина такая. Есть выпуклая вверх кривая , и берутся всевозможнные каса-

тельные к этой кривой. Множество касательных в объединении с кривой есть решение нашей задачи.

Замечание: существенно было то, что мы предположили, что f00 > 0. Хорошая задача: исследовать случаи

80 ГЛАВА 4. АВТОНОМНЫЕ СИСТЕМЫ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ

4.Мы рассмотрели одно уравнение, где есть только одна искомая функция. Рассмотрим

уравнение

~ 0 ~ f(t; ~y; ~y ) = 0;

Все функции достаточно гладкие. Предположим, задана точка t0; y~0; p~0 :

nf~(t0; y~0; p~0) = ~0o

^ ndet

~

(t; y~0; p~0) 6= 0o

 

@f

 

 

@~p

Тогда (см 3)

0

~

 

 

0

 

 

 

 

è ò.ä.

~

,

~y

= ~g(t; ~y)

f(t; ~y; ~y

) = 0

 

 

Глава 4. Автономные системы дифференциальных уравнений

1. Основные определения и простейшие свойства

1.

Определение: Система дифференциальных уравнений

åñëè ~

 

 

 

 

 

 

 

 

 

f(t; ~y) не зависит от t.

 

 

 

 

 

 

 

 

 

Далее будем рассматривать систему

 

 

 

 

 

 

 

 

 

 

 

 

d~y

~

 

 

 

 

 

 

dt

 

= f(~y);

f.1

ãäå

~y =

y.1

 

; f~ =

 

 

 

 

 

 

 

 

 

 

 

 

yn

 

 

fn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d~y

dt

;

~

= f(t; ~y) называется автономной,

(1)

 

 

 

~

 

~

~

~

@f

 

 

t; ~y; f вещественные. f определена в области G, f и

@~y

непрерывна в G.