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

Diff

.pdf
Скачиваний:
20
Добавлен:
03.06.2015
Размер:
529.83 Кб
Скачать

Дана функция g(t; y; z). Заданные функции f и g определены при t1 6 t 6 t2, остальные

аргументы любые.

z(t)

допустимая, если

Определение: Пара

 

 

y(t)

 

 

 

 

 

 

 

 

 

1.g t; y(t); z(t) = 0 ïðè t1 6 t 6 t2.

2.y(t1) = y1, z(t1) = z1, y(t2) = y2, z(t2) = z2, ãäå y1, z1, y2, z2 заданы (разумеется g(t; y1; z1) = 0, g(t2; y2; z2) = 0.)

Определение: Допустимая пара

z

 

дает минимум функционала F , если существует " > 0

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ïàðû

 

y

 

, удовлетворяющей условию

 

 

такое, что для любой допустимой

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~

 

 

~

 

 

maxt

p(y(t) y~(t))

2

 

 

 

 

2

 

 

 

2

2

 

 

+ (z(t) z~(t))

+ maxt

q(y0(t) y0

(t))

 

+ (z0(t) z0(t))

 

< ";

имеет место F (y; z) 6 F (~y; z~).

Геометрический смысл задачи Лагранжа:

[ Тут будет рисунок ]

Задана поверхность , рассматриваются кривые, соединяющие данные две точки, и лежащие

на данной поверхности. Среди этих кривых надо найти ту, которая дает минимум функционала F .

Замечание: Если мы можем выразить, скажем, y через t и z, то надо это сделать и полу-

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

Теорема: Пусть допустимая пара z

да¼т минимум F . Пусть в каждой точке кривой y =

y(t), z = z(t), t [t1; t2

] имеет место

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

_

 

 

 

 

 

 

 

 

@y

6= 0

@g

6= 0 :

 

 

 

 

@y

@z

 

 

 

 

 

 

 

 

 

 

 

 

 

z

является экстремалью функ-

Тогда существует функция (t), t 2 [t1; t2] такая, что пара

ционала H(y; z), где

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

H(y; z) =

 

t2

 

 

 

 

 

 

 

 

 

 

 

 

tZ1

h(t; y; z; y0; z0)dt; h = f + g:

 

def

 

 

 

 

 

 

 

 

 

 

 

 

 

То есть, удовлетворяется система уравнений

 

 

 

 

 

 

 

 

 

 

 

 

8 hy

dthy0 = 0;

 

 

 

 

 

 

> hz

 

 

d

 

 

 

 

 

 

 

 

 

 

 

d

hz

 

= 0:

 

 

 

 

 

 

<

 

 

 

 

0

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

>

 

 

 

 

 

 

 

 

 

:

 

 

 

 

 

 

 

 

 

 

31

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

 

 

 

 

 

 

 

 

[ Тут будет рисунок ]

 

 

 

Поверхность задана уравнением

 

. Берем любую точку ^ на кривой, удовлетворяющую

 

 

 

 

 

 

 

g = 0

 

 

 

 

 

t

 

условию ^

 

. Затем проведем построения в некоторой окрестности этой точки.

 

 

t 2 (t1; t2)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

^

 

 

^

 

 

 

 

 

 

 

 

 

 

 

 

y(t) = y;^ z(t) = z:^

 

 

 

Пусть в этой точке частная производная

@g

^

 

не равна нулю. Существует окрестность

 

 

@z (t; y;^ z^)

 

 

 

. (Это следует из теоремы о

точки

U

точки ^

 

, для которой

g(t; y; z) = 0

, z = a(t; y)

 

 

t; y;^ z^

 

 

 

неявной функции.)

 

 

 

 

 

 

 

 

 

 

 

Возьмем ^

^ ^

 

 

 

 

 

 

 

 

 

 

 

 

 

t1

< t < t2.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t2

 

 

 

t^1

 

t^2

 

 

t2

 

 

 

Z

f(t; y; z; y0; z0)dt = Z ( ) + Z

f^ t; y; y0

 

dt + Z ( ); ãäå

 

 

 

t1

 

 

 

t1

 

t^1

 

t^2

 

 

 

 

 

 

f^ t; y; y0

 

= f t; y; a(t; y); y0; at(t; y) + ay(t; y)y0

 

 

 

 

 

 

 

 

 

t^2

 

 

 

 

 

На отрезке [t^1; t2^1] функция y(t) экстремаль Z^

f^dt.

 

 

 

 

 

 

 

 

 

 

 

 

 

t1

 

 

 

 

 

Выразим производные:

f^y

dtd f^y0 =

f^y dtd f^y0 = 0:

= fy + fzay + fz0aty + fz0ayyy0

f^y0 = fy0 + fz0ay:

d

fy0

+

 

d

fz0

ay + fz0ayt + fz0ayyy0

 

 

dt

 

dt

 

 

 

 

 

По теореме о равенстве смешанных производных, подчеркнутые слагаемые взаимно уничтожаются. Следовательно,

0 = fy

d

 

fz

d

 

 

 

ay

 

 

fy0

+

 

 

fz0

dt

dt

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

 

 

gy

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ay =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

gz

 

 

 

 

Значит, последнее равенство можно записать в виде

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d

 

 

 

 

d

 

gy

 

fy

dt

fy0 fz

 

dt

fz0

 

gz

= 0:

 

fy

d

 

fz

d

 

 

 

 

 

 

fy0

 

 

fz0

 

 

dt

 

dt

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

gy

 

 

 

gz

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

32

Вопрос: Если gy = 0?

Ответ: Можно трактовать равенство следующим образом: хотя бы один из знаменателей отличен от нуля; если оба знаменателя отличны от нуля, то это обычное равенство. Если один из знаменателей равен нулю, то соответствующий числитель тоже равен нулю.

Введем функцию (t):

 

 

fy

d

 

 

fz

d

 

 

 

 

fy0

 

 

fz0

 

def

 

dt

=

dt

t1 6 t 6 t2:

(t) =

 

 

 

 

 

 

 

 

 

;

 

gy

 

 

 

gz

 

def

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

h = f + g.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

hy = fy + g + y; hy0 = fy0

 

 

hy

d

= fy

d

 

 

 

 

 

 

hy0

 

fy0

+ gy

= 0

 

dt

dt

Первое уравнение удовлетворяется. Аналогично показываем, что удовлетворяется и второе. Теорема доказана.

Замечание: Здесь y = y(t; c1; c2; c3; c4; (t)), z = z(t; c1; c2; c3; c4; (t)). Находим параметры из граничных условий и из равенства g(t; y; z) 0.

Глава 3. Исследование задачи Коши

Лекция 07.02.2012 пропущена!

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

Замечание: В этом параграфе все величины вещественные. Далее будем иметь дело с набо-

z1

рами чисел вида k . k

def

~

 

 

= Z, n фиксированное

zn

 

 

 

 

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

 

 

d~y

 

 

 

 

~

 

 

 

dt

= f(t;

f определена в замкнутой области G.

для всего параграфа.

~y);

(1)

Определение: Функция ~y(t) решение (1), если:

1.~y(t) определена на некотором промежутке [ ; ], < ,

2.~y0(t) непрерывна на [ ; ],

3.

d~y(t)

~

dt

 

= f(t; ~y(t)) íà [ ; ]

33

0

4. Для каждого < t < точка (t; ~y (t)) внутренняя точка G.

Замечание: Что означает условие 4.

[ Тут будет рисунок ]

. Мы отбрасываем случаи, когда внутренние точки графика решения лежат на границе нашей области. Допускаются только случаи, когда они лежат внутри этой замкнутой области.

 

~

1. Теорема: Пусть ~

@f

непрерывны в

f,

@~y

G. Пусть

1.существует > 0, что на отрезке [t0 ; t0 удовлетворяющее условию

~y(t0) =

(t0; y~0) внутренняя точка G. Тогда:

+ ] определено решение уравнения

(1),

y~0

(2)

2.Åñëè y~1(t); y~2(t) какие-либо решения задачи (1) ^ (2), то y~1(t) = y~2(t) на пересечении отрезков их определения.

Замечание: во первых фиксируется, что при заданных предположениях задача Коши имеет решение. Пункт второй утверждает единственность этого решение. Обратите внимание: утверждение 1 носит локальный характер. Решение существует в окрестности точки. Пункт 2 же носит глобальный характер.

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

 

 

 

 

 

 

1. Возьм¼м U какая-либо окрестность (t0; y~0), такая, что U

G.

 

 

 

 

 

 

Возьм¼м Qfjt t0j 6 ; j~y y~0j 6 g, 0 < ; 0 < ; M 6 ; Q 2 U, где M максимум

~

 

 

 

 

 

нормы M = maxU jfj.

[ Тут будет рисунок ]

 

 

 

 

 

Обозначим ! = [t0 ; t0 + ], K = maxQ @~y

: Мы будем строить решение на отрезке !

 

 

~

 

 

 

 

@f

 

 

 

следующим образом. Бер¼м какое-то положительное

h > 0, и строим на нашем отрезке

 

 

 

 

 

 

сетку:

 

 

 

 

 

 

[ Тут будет рисунок ]

 

 

 

tk = t0 + kh;

k = 0; 1; 2; : : :

 

Для зафиксированной сетки строим функцию '~(t):

(далее вс¼ при t > t0). Åñëè íàäî

провести рассуждения при меньших t, надо сделать замену переменной, поменяв знак.

~

'~(t) = y~0 + (t t0)f(t0; y~0) ïðè t0 6 t 6 t1,

~

'~(t) = '~(t1) + (t t1)f (t1; '~(t1)) ïðè t1 6 t 6 t2

~

'~(t) = '~(t2) + (t t2)f (t2; '~(t2)) ïðè t2 6 t 6 t3

Замечание: Мы строим такой отрезок от t0 äî t1, который касается графика решения в точке (t0; y0). Мы доходим до точки t1. Бер¼м график решения, если он существует,

34

в этой точке, и доходим до t2. Задавать вопрос ¾почему решение существует¿, не нужно. Мы просто строим последовательность значений по этим формулам. Ни на какие соображения мы не ссылаемся.

t

 

't~ '~(t ) = tZ

'~0(t) dt

 

 

[ Тут будет рисунок ] '~0 кусочно-непрерывная функция, поэтому можно воспользоваться формулой Ньютона-Лейбница. Поэтому

'~(t ) '~(t )

6 jt t jM:

 

 

 

В частности,

'~(t) y~0 6 jt t0jM 6 M 6 :

Мы получили кривую, которая через боковые стенки нашего цилиндра не может выйти за пределы нашего цилиндра.

[ Тут будет рисунок ] Следовательно, '~(t) определена при jt t0j 6 .

Теперь докажем, что '~(t) есть "-приближ¼нное решение (1), " ! 0 при h ! 0. Давайте оценим результат подстановки

 

 

 

(t) =

 

'~(t)

f~(t; '~(t))

d dt

 

 

 

 

 

 

 

 

 

 

в точках непрерывности '~0(t).

tr < t < tr+1 : d'~(t) dt

~

= f(t2; '~(t2)).

~

(t) = f(t2; '~(t2))

~

f(t; '~(t))

j j , j j ~ . ) t t2 6 h '~(t) '~(t2) 6 Mh, f равномерно непрерывна в Q

(t) 6 ", ãäå " ! 0 ïðè h ! 0.

h1; h2; : : : ; hs > 0; hs ! 0 ïðè s ! 1.

'~1(t) "s- приближ¼нное решение "s ! 0 ïðè s ! 1. j'~s(t) '~p(t)j 6 2" ek ;

" = max("s; "p). Ссылаемся на теорему, доказанную во втором параграфе.

'~s(t0) = '~p(t0) = y~0:

Используя критерий Коши, получаем

'~1(t) íà ! ~(t) ïðè s ! 1:

35

~(t) непрерывна, ~(t) = y~0:

d'~s(t)

dt

в точках непрерывности '~s0(t).

~

= f(t; '~s(t)) + s(t)

t

t

ZZ

 

~

'~s(t) = y~0 + f( ; '~s( )) d + s( ) d

t0

t0

~

~

~

(t))

f(t; '~s(t)) íà ! f(t;

 

~(t) = y~0 + tZ0

t

 

 

f~( ; ~( )) d

(3)

(4)

d~(t)

~

~

 

~

 

dt

 

= f(t;

 

(t));

 

(t0) = y~0

Осталось только проверить, что всякая внутренняя точка графика является внутренней точкой области.

jt t0j < ) j'~(t)j < M

2. Доказательство единственности решения. Докажем от противного. Пусть существуют

два решения задачи

(1)

^ (2) y~1(t)

è

y~2

(t)

,

y~1

^

^ .

 

 

 

(t) 6= y~2

(t)

Возьм¼м

 

 

 

 

^.

 

 

 

 

 

^

 

S = ft: y~1(t) = y~2(t)g \ [t0; t]

 

 

 

S

 

самая правая точка замкнутого ограниченного множества

þ

t0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

^

^

(t) 6= ~y2

(t)

ïðè ^

^

 

 

 

 

y^1(t0) = y^2(t0); ~y1

t0 < t 6 t:

^

^

 

(îíà æå

^

^

) внутренняя точка

 

. Для этой точки строим

(t0); y^1(t0)

 

(t0); y^2

(t0)

 

 

G

 

 

 

[ Тут будет рисунок ]

^ ^ ^

Q; M; K.

 ^

(t)

è

~y2

(t)

0-приближ¼нные решения (точные решения),

~y1

^

^

.

Q ~y1

 

 

(t)

6 2^ 0 ek^^

 

(t0) = ~y2

(t0)

 

Поэтому

~y1(t) ~y2

= 0. Противоречие.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2. Теорема: (о продолжении решения) Пусть выполняются условия теоремы п.1. Пусть

G

ограничено. Тогда решение задачи (1) ^ (2) может быть продолжено так, чтобы концы его

графика лежали на границе

G.

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

j~j

M = max f ; M > 0 (для M = 0 очевидно)

G

(t0; y~0); Q = fjt t0j 6 ; j~y y~0j 6 M g.

p

: 1 + M2 = ;

где расстояние от (t0; y~0) от границы G.

[ Тут будет рисунок ]

36

Далее t > t0. Мы построили решение на отрезке [t0; t0 + ] (см конструкцию, фигурирующую в доказательстве предыдущей теоремы).

Если правый конец достиг границы, то продолжать конструкцию незачем. Но пуст правый конец, то есть точка (t + ; ~y(t + )) не является граничной. Возьм¼м эту точку на рисунке.

Она является внутренней, поэтому мы можем продолжить решение направо. Бер¼м точку (t + ; ~y(t0 + )) в качестве начальной, [t0 + ; t0 + + 1], ãäå 1 выбирается из следующих

соображений: p

1 : 1 1 + M2 = 1;

ãäå 1 расстояние от точки (t + ; ~y(t0 + )) до границы G, и так далее.

В итоге получаем последовательность отрезков, на которую распространяем наше решение.

[t0; t0 + ]; [t0; t0 + + 1]; [t0; t0 + + 1 + 2]; : : :

k расстояние от точки t0 + + : : : + k 1; ~y(t0 + + : : : + k 1) .

Докажем, что эти точки стремятся к границе.

p X X X

1 + M2 k = k; k сходится:

k k

Следовательно, k ! 0 ïðè k ! 1.

В утверждении теоремы сказано, что точка должна лежать на границе, а не просто стремиться. Мы получили решение на [t0; t0 + ^), ãäå ^ = + 1 + 2 + : : :. На этом полуотрезке имеет место формула

 

t

 

 

 

 

~y(t) = y~0 + tZ0

f~( ; ~y( )) d

Обратим внимание на обстоятельство:

 

 

t t

M

Z t

f~( ; ~y( )) d

6

 

 

 

 

 

 

 

 

 

 

 

 

t0

 

 

 

 

 

По критерию Коши существует

t

Z

~

lim f( ; ~y( )) d

t!t0+^

t0

Естественным образом определяем решение по формуле:

0 + ^) = 0

+ t!t0+^ tZ0

t

~

~y(t y~

lim

f( ; ~y( )) d

~y(t) определена при t0 6 t 6 t0 + ^;

0 ~

непрерывна, и на отрезке [t0; t0 + ^] имеет место ~y (t) = f(t; ~y(t)).

Замечание:

37

1. Если имеем уравнение вида

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

(5)

то заменой z1 = y0, z2 = y00, : : :, zn 1 = y(n 1) уравнение сводится к системе y0

= z1; z10 =

z2; : : : ; zn0 1 = f(t; y; z1; : : : ; zn 1). Обязательная задача. Используя результаты пунктов 1 и 2, сформулировать соответствующие теоремы для уравнения (5)

2. Эта теорема может быть доказана очень многими способами. Мы выбрали метод ломаных Эйлера.

[ Тут будет рисунок ]

Эта построенная ломаная называется Ломаной Эйлера. Таким образом был предложен приближ¼нный метод решения дифференциального уравнения. На каждом отрезке вместо решения строим отрезок, который в начальной точке имеет направление, заданное нашим дифференциальным уравнением.

Затем мы доказывали, что наше приближ¼нное решение при h ! 0 становится точным.

Задача: составить программу приближ¼нного решения задачи Коши, используя метод Эйлера приближ¼нного решения. Можно попробовать решать задачи из задания таким способом.

3.В одной из теорем мы объединили два формально разных вопроса: существования и единственности решения задачи Коши. В более тонких исследованиях эти два вопроса разделяются. Оказывается, что если правая часть уравнения только непрерывна (раньше требовалось дифференцируемость), то хотя бы одно решение задачи Коши существует.

4.Конечно, условие в этой теореме это достаточное условие, и вовсе не является необхо-

димым. Рассмотрим 2 примера.

p

y0 = 3 y2

p

y0 = 3 y2 + 1

Контрольный вопрос: Почему? Какие условия теоремы не выполнены?

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

В этом параграфе все величины вещественные. Мы рассматриваем уравнение

 

 

d~y

 

= f~ t; ~y; ~y0

 

 

dt

y1

 

 

 

 

f1

~y = k .

~

 

. k; ~ = k

k; f = k

yn

 

 

 

 

fn

~ 0

f(t; ~y; ~y ) определена в замкнутой области Gt;~y;~ .

(1)

1

. k;n

38

Тогда мы получим, что решение зависит от параметра и начальных данных.

~y(t0) = y~0

(2)

В этой параграфе мы будем изучать функцию y, как функцию всей совокупности переменных.

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

 

 

 

 

~

 

 

 

 

 

 

 

 

 

 

 

6

 

 

 

 

 

 

 

 

1. Теорема: Пусть ~

 

 

df

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f è

 

 

 

; , такие, что

 

 

 

t

t0

 

 

 

~

~0

 

6

 

 

 

 

 

dt непрерывны в G. Пусть

t0; y~0; ~0

внутренняя точка G. Тогда

существуют положительные числа

 

 

 

 

ïðè

j

 

j

 

è

j

 

j

 

определено

решение (единственное) задачи (1) ^ (2) и функция ~y(t; ~) непрерывна.

 

 

 

 

 

 

M = maxU

 

f . (Можно

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t0; y~0; ~0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

G.

Доказательство: Мы бер¼м нашу точку

, и такую окрестность U, что U

Обозначим

 

 

j

~

 

 

показать, что такое обозначение корректно, максимум

существует).

 

 

j

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Возьм¼м такой цилиндр Qfjt t0j 6 ; j~y y~0j 6 ; j~ ~0j 6 g (бицилиндр), что

 

 

 

 

 

 

 

 

 

 

 

6

 

 

 

 

 

 

=

 

Q

@~y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~

 

 

 

 

 

 

 

0 < ; 0 < ; 0 < ; M

 

; Q

 

U:K

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

max

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Мы бер¼м два значения ~1; ~2; ~y(t; ~1); ~y)(t; ~2) решение (1) ^ (2).

 

 

 

 

 

 

 

 

Как следует из конструкции построения этих решений, приведенных в предыдущем параграфе, эти решения определены на вс¼м отрезке.

Посмотрим, будут ли эти решения близки, если ~1 близко к ~2.

 

 

 

t; ~

1)

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

 

 

 

 

d~y(

 

 

 

 

dt

 

 

 

~

 

 

 

~

 

(3)

 

jf(t; ~y(t; ~

1); ~1) f(t; ~y(t; ~2); ~1)j < ";

åñëè

j~1 ~2j < ; ãäå " ! 0 ïðè ! 0;

 

 

 

òàê êàê ~

 

 

 

 

 

 

 

 

f равномерно непрерывна в Q.

 

 

 

 

Тогда

 

 

 

 

 

 

 

 

t; ~

 

 

< ";

 

 

d~y(dt

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

åñëè j~1 ~2j < . Таким образом, обе функции являются "-приближ¼нными. Одна из них является точным решением. Вс¼ происходит в замкнутом выпуклом множестве, функции удовлетворяют одним и тем же начальным условиям, поэтому

j~y(t; ~1) ~y(t; ~2)j 6 2 "eK

Замечание: Непрерывность следует из неравенства треугольника: j~y(t1; ~1) ~y(t2; ~2)j 6 j~y(t1; ~1) ~y(t1; ~2)j + j~y(t1; ~2) ~y(t2; ~2)j 6 2 "eL + Mjt1 t2j Теорема доказана.

39

~

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

@~

G.

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

@~y

 

и функция V

удовлетворяет

@~

уравнению

 

 

 

 

 

 

 

 

 

 

 

@~y

 

t; ~y(t; ~); ~

!V

+ @~

t; ~y(t; ~); ~

(4)

 

dt =

 

 

dV

 

~

 

 

 

~

 

 

 

 

 

@f

 

 

 

 

 

@f

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

V (t0) = 0

 

 

 

 

 

 

 

(5)

Замечание: Уравнение (4) называется уравнением вариации для уравнения (1). Оно получа- ется из (1) формальным дифференцированием по ~. (Это всего лишь способ запомнить это

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

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

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

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

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

То, что это решение распространяется на весь отрезок [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

 

 

 

dt

=

 

 

(

t; ~

+

~

)

f~

 

 

 

 

 

@~y

 

~

 

@~

 

~

 

 

f t; ~y

 

 

; ~ + ~

 

 

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

 

 

 

V

 

 

 

 

 

 

 

|

 

 

 

 

 

 

 

{z

 

 

 

 

}

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

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

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

~y(t; ~ ~)

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

40

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