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

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

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

2. Некоторые обобщения простейшей задачи вариационного исчисления.

 

51

3.

Задача со свободным концом. Фиксируем отрезок t1 6 t 6 t2, y(t),

 

 

 

 

F (y) =

t2

f

t; y(t); y0(t) dt;

 

 

 

 

Z

 

 

 

 

 

def

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t1

 

 

 

 

 

 

 

f

задана, определена при

t1 6 t 6 t2;

y;

y0

любые

.

 

 

" > 0

 

 

 

 

y(t) äàåò

 

 

 

F

 

 

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

 

 

 

 

минимум

 

, если существует

 

, такое, что

для любой допустимой функции y~(t), удовлетворяющей условию

 

 

 

maxt y(t) y~(t) + maxt

y0(t) y~0(t) < ";

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

имеет место

F (y) 6 F (~y):

Замечание: Данное определение дословно повторяет определение допустимой функции, данное в первом параграфе. Однако там допустимый подразумевало то, что функция удовлетворяет начальным условиям на обоих концах, а сейчас только на левом.

Определение: Пусть допустимая функция y(t) дает минимум F . Тогда y(t) удовлетворяет уравнению Эйлера

d

 

 

fy dtfy0

= 0:

(6)

и граничному условию

hfy0it=t2

= 0:

(7)

Доказательство: Берем y(t),

 

 

(t): (t1) = 0; число (малое): y~ = y(t) + (t); !( ) = F (y + )

!( ) имеет минимум при = 0, !0(0) = 0.

t2

 

 

 

 

 

 

 

!0(0) = tZ1

(fy + fy0 0)dt

 

 

 

 

t2

 

 

 

 

t1

= Z (fy dtfy0) dt + fy0

 

 

d

h

i

t2

t1

 

 

 

 

 

t2

(t) (t)dt + fy0 t=t2 ;

= fy dtfy0

= Z

 

 

 

h i

 

 

d

t1

 

 

 

 

 

 

Сначала возьмем (t2) = 0.

t2

Z

(t) (t)dt = 0;

t1

52

ГЛАВА 2. ЭЛЕМЕНТЫ ВАРИАЦИОННОГО ИСЧИСЛЕНИЯ

hi

откуда fy0 . Берем (t2) 6= 0. Получаем

t=t2

hi

fy0 = 0:

t=t2

Замечание: Условие (7) называется естественным граничным условием для F .

Некоторые итоги. Было рассмотрено несколько обобщений: когда функций много, но производные второго порядка, когда функция одна, но производных много, и когда граничное условие одно. Что делать, если встречается смешанная задача? Громоздко формулировать такую теорему, но метод решения таких задач ясен: он вытекает из приведенных выше рассуждений.

4. Обязательная задача.

Фиксируется отрезок t1 6 t 6 t2, y(t), z(t),

F (y; z) =

t2

f(t; y; y0; y00; z; z0; z00; z000)dt:

tZ1

def

 

 

Не задано никаких граничных условий. Написать систему уравнений Эйлера, вывести естественные условия. Рекомендуется это сделать, чтобы осознать, что решение первых трех пунктов дает метод для исследования более сложных случаев.

3. Условный экстремум

В этом параграфе все функции достаточно гладкие. 1. Изопериметрическая задача.

t1 6 t 6 t2, y(t),

t2

 

 

t2

 

 

F (y) = Z f t; y(t); y0

(t) dt; G(y) = Z g t; y(t); y0

(t) dt

t1

 

 

t1

 

 

F (y); G(y) заданные функции, f и g определены при t1 6 t 6 t2, y è y0 любые.

Определение: Функция y(t) называется допустимой, если y(t1) = y1; y(t2) = y2; G(y) = `;

ãäå y1, y2, ` заданы.

Определение: Допустимая функция y(t) дает минимум F , если существует " > 0 такое, что для всех допустимых функций y~(t), удовлетворяющих условию

maxt

y(t) y~(t)

+ maxt

y0

(t) y~0(t)

 

 

 

 

 

 

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

3. Условный экстремум

53

Теорема: Пусть допустимая функция y(t) дает минимум F , и y(t) не является экстремумом функционала G. Тогда существует число такое, что y(t) экстремаль функционала H = F + G, то есть удовлетворяет уравнению

d

hy dth = 0;

ãäå h = f + g.

Доказательство: Пусть y(t) решение задачи.

0 : 0(t1) = 0(t2) = 0;

t2

 

gy

t; y(t); y0(t) (t) + gy0

 

t; y(t); y0(t) 00

(t) dt 6= 0:

Z

 

 

t1

 

 

 

 

0(t) существует, так как y(t) не является экстремалью функционала G.

y(t) + 0(t) + (t);

(t1) = (t2), и числа (малые).

def

def

+ ):

!( ; ) = F (y + 0

+ ); '( ; ) = G(y + 0

= = 0 решение задачи: найти минимум ! при условии ' = `. Тогда:

9 :

@

= =0 = 0;

(

@

= =0 = 0:

 

@(! + ')

 

@

! + ')

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Из первого равенства следует:

 

 

 

 

 

t2

 

 

 

 

 

(

@

 

)

Z

fy 0 + fy0n00 + (gy 0 + gy0 00 ) = 0

= =0

@

! + ')

 

 

t1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

0 t2 fy 0

+ fy0 00

1 0 t2(gy 0 + gy0

00

) dt1

 

 

 

Z

 

 

 

,

Z

 

A

 

 

 

@t1

 

 

 

A

@t1

 

не зависит от выбора (t). Посмотрим, что дает второе равенство:

t2

(

@

= =0

)

Z

(fy + gy) + (fy0

+ gy0) 0 dt = 0:

@

! + ')

 

 

t1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

h = f + g; h не зависит от выбора .

t2

(hy + hy0 0) dt = 0 ) hy dthy0

= 0: (ñì. 1)

tZ1

 

 

d

 

54

ГЛАВА 2. ЭЛЕМЕНТЫ ВАРИАЦИОННОГО ИСЧИСЛЕНИЯ

Замечание: Мы получили уравнение Эйлера для функцонала H. В типичном случае мы получим решение в виде уравнения второго порядка, y будет функцией от t; c1; c2; ,

y(t1) = y1; y(t2) = y2; G(y) = `:

Задача: На плоскости заданы две точки A и B. Соединяем их отрезком, рассматриваем некоторую соединяющую их кривую . Надо выбрать ее так, чтобы длина кривой была постоянной, и площадь ограниченной фигуры была максимальна.

Выберем координаты так, чтобы отрезок AB лежал на оси Ox.

: y = y(x); y(x1) = y(x2) = 0; x1 < x2:

 

x2

 

 

 

x2

 

 

 

 

 

F (y) = Z

 

 

G(x) = Z

 

 

 

 

 

y(x)dx;

 

1 + (y0)2dx

 

 

x1

 

 

 

x1

p

 

G(y) = ` задано, F (y) max.

 

 

 

 

 

 

 

 

 

 

 

 

 

x2

 

 

 

 

 

 

h = y +

 

; H = Z

 

 

 

1 + (y0)2

y + 1 + (y0)2

dx

p

 

 

x1

 

 

p

 

Решая уравнение Эйлера, получаем:

(x + c1)2 + (y + c2)2 = 2

Если задача имеет решение, то это решение дуга окружности. Рассмотрим несколько слу- чаев:

1.` < x2 x1. Задача не имеет решения. (Множество допустимых функций пусто)

2.` = x2 x1. Допустимая функция только отрезок [x1; x2].

3. ` > x2 x1, ` < 2 (x2 x1). В этом случае получается дуга, которая является гладкой кривой.

3. Условный экстремум

55

4. ` = 2 (x2 x1). Получим полуокружность, jy0j = 1 на концах.

5. ` > 2 (x2 x1). Получим дугу окружности.

Кривая не представляется уравнением какой-нибудь функции.

Замечание: Под нашу теорию подходит только случай 3). Можно показать, что решение в случаях пунктов 3), 4), 5) действительно является искомым. Данная задача является изопери-

метрической в буквальном смысле. Эта конкретная задача дала название всему класса задач, который мы рассматриваем, несмотря на то, что функционал G может не давать длину, а

какое-нибудь другое значение. другое значение

56

 

 

 

 

ГЛАВА 2. ЭЛЕМЕНТЫ ВАРИАЦИОННОГО ИСЧИСЛЕНИЯ

2.

Задача Лагранжа

 

t1

6 t 6 t2,

z(t)

;

 

 

 

 

y(t)

 

t2

 

 

 

 

 

 

 

 

 

 

 

 

Z

F (y; z) = f t; y(t); z(t); y0(t); z0(t) dt:

t1

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

y(t)

Определение: Ïàðà z(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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q

p

2 2 0 ~0 2 0 ~0 2 max (y(t) y~(t)) + (z(t) z~(t)) + max (y (t) y (t)) + (z (t) z (t)) < ";

t t

имеет место F (y; z) 6 F (~y; z~). Геометрический смысл задачи Лагранжа:

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

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

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

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

3. Условный экстремум

z

 

 

 

 

 

 

 

 

 

 

57

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

да¼т минимум 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

 

 

 

 

 

 

 

 

>

 

 

 

 

 

 

 

 

 

 

:

 

 

 

 

 

 

 

 

 

 

 

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

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

 

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

 

 

 

 

 

 

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

58

ГЛАВА 2. ЭЛЕМЕНТЫ ВАРИАЦИОННОГО ИСЧИСЛЕНИЯ

f^ t; y; y0

 

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

 

 

 

 

t2

 

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

f^dt.

 

 

 

 

 

t1

 

 

 

 

^

d

^

 

 

 

 

fy

dt

fy0

= 0:

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

f^y = fy + fzay + fz0aty + fz0ayyy0

 

 

 

 

 

f^y0 = fy0 + fz0ay:

 

 

 

 

 

 

 

 

 

 

d

^

 

d

 

 

d

 

 

0

dt

fy0

=

dt

fy0

+

 

dt

fz0

ay +

fz0ayt + fz0ayyy

 

 

 

 

 

 

 

 

 

 

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

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

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

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

 

 

fy

 

d

 

 

fz

d

 

 

 

 

fy0

 

 

fz0

def

 

dt

 

dt

(t) =

 

 

 

 

 

=

 

 

 

 

 

 

; t1 6 t 6 t2:

 

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.

59

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

1.

Вспомогательные сведения

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1.

Рассмотрим матрицу A из функций apq,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A(z1; :::; zn) = 0 a: :11:

:: :: ::

a:

1:m:

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

@ an1

: : : anm

A

 

 

 

 

 

 

 

 

 

apq

= aqp(z1; :::; zr) где n; m фиксированные числа, z1; : : : ; zr

вещественные переменные,

apq вещественные или комплексные (нужно уточнить).

 

 

 

 

 

 

 

 

 

 

 

Определение: A(z1; :::; zn) непрерывна

def

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

, âñå apq непрерывны.

 

 

 

 

 

 

2.

 

 

 

 

 

 

 

0 y.1

1

 

 

~x = 0 x.1

 

1, òî

 

Определение: Если задана функция ~y(~x), ~y =

,

 

 

 

 

 

 

 

 

 

 

B yn

C

 

 

 

B xm

C

 

 

 

 

 

0

 

 

@

 

A

 

 

1

@

 

 

 

A

 

 

 

 

 

@y~1

 

:::

 

@y~1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

@~y def

@x~1

 

@x~m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

B

:::

 

:::

 

:::

 

 

C

 

 

 

 

 

 

 

 

 

 

 

 

@~x =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

@y~n

 

 

@y~n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

B

 

 

:::

 

 

 

 

C

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

@x~1

 

@x~m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

B

 

 

C

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

@

 

 

 

 

 

 

 

A

 

 

@~y

 

 

@~y @~x

 

 

 

 

 

 

 

 

~

 

 

 

 

 

 

 

 

 

Задача: Доказать, что если ~y = ~y(~x) , ~x = ~x(~z) (y) = (~x(~z)), то

 

 

 

=

 

 

 

 

 

@~z

@~x

@~z

 

0 a11 : : : a1m

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3.Пусть A = @ : : : : : : : : : : : : : A, aij вещественные или комплексные.

an1 : : : anm

Определение: Норма матрицы:

jAj = s

i;j

jaijj2

 

X

 

Свойства:

1.

Пусть A + B определено. Тогда

 

 

jA + Bj jAj + jBj

 

2.

Пусть AB определено. Тогда

 

 

jABj jAjjBj

 

 

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

 

 

C = AB; C = kcrsk; crs = Xp

arpbps

60

 

 

 

X

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

 

 

 

 

X

X

 

 

 

 

 

 

jcrsj2 = arpbps jarpj2

jbqsj2

 

 

jCj2 = r;s

 

 

 

p

p

! =

q

 

 

 

 

jcrsj2

 

r;s

jarpj2

jbqpj2

r;p

jarpj2

q;s

jbqsj2

= jAj2jBj2

X

 

 

X Xp

Xq

 

X

 

X

 

 

3. Если комплексное число, то j Aj = j jjAj

4. Теорема: Пусть f~ =

0 y.1

1

~x =

0 x.1

1, ~x вещественный,

 

B

C

 

B

C

 

@

A

 

@

A

yn xm

комплексная, и выполнены следующие условия:

1. Функция ~ f(~x) определена в замкнутой области G.

 

 

~

 

2.

~

@f

непрерывны в

 

 

f è

@~x

G.

3.

 

~

~

отрезок [~a; b] лежит в G

Тогда

j ( ) (

)j ~x

[~a;~b] @~x

j

 

j

 

 

 

 

 

~

 

 

~ ~

~ ~

 

 

@f

 

 

 

f b

f a

2

 

 

b

a

max

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~

f вещественная или

Доказательство: Рассмотрим

def ~ ~ ~ ~

 

0 t 1

!~(t) = f(a + t(b a));

 

1

1

@~x(~a + t(~b ~a))(~b ~a)!dt

f~(~b) f~(~a) = !~(1) !~(0) = Z !~0(t)dt = Z

 

 

~

 

 

@f

0

0

 

 

j~ ~ ~ ~ j f(b) f(a)

1

@~x

(~a + t(~b ~a))(~b ~a)!dt

Z

 

~

 

@f

 

1

~

 

Z

 

@f

(~a

@~x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+ (

) j

 

j

~x [~a;b~]

@~x

j

j

~

~

 

~

~

 

 

~

 

 

t b

a

 

 

b

a dt

2

 

 

 

b

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0 0

5. Теорема: (Лемма Гронуолла).

Пусть при 6 t 6 определена вещественная функция f(t). Пусть f(t) непрерывна и удовлетворяет на 6 t 6 соотношению

f(t) 6 A + B Z

t

 

f( )d ;

(1)

где A и B - вещественные числа, B > 0. Тогда при 6 t 6 имеет место f(t) AeB(t )