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

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

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

1. Простейшая задача вариационного исчисления

41

Теорема: Пусть y(t) дважды непрерывно дифференцируема.

Тогда F (y) = 0 для любой допустимой вариации y тогда и только тогда, когда y(t) удовлетворяет уравнению (Эйлера)

d

 

 

fy dtfy0

= 0:

(2)

Здесь

 

d

- полная производная,

d

f

 

 

= f

 

 

+ f

 

y0 + f

y00.

 

 

 

 

 

 

y0

y0t

y0y

 

 

 

 

 

 

dt

 

 

 

 

dt

 

 

 

 

 

y0y0

 

 

 

 

Все производные существуют и непрерывны в силу наложенного заранее ограничения.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t2

 

 

 

 

 

 

 

Доказательство: Преобразуем !0(0) = tR1

(fy +fy0 0)dt. Для этого проинтегрируем по частям:

 

 

 

 

t2

 

fy dtfy0

dt +

 

 

fy0

t1

t2

(t) (t)dt;

(t) = fy dtfy0:

 

 

 

 

Z

 

 

 

= Z

 

 

 

 

 

 

 

 

d

 

h

 

 

 

i

t2

 

 

 

 

 

d

 

 

 

 

t1

 

 

 

 

 

=0

 

 

t1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

}

 

 

 

 

 

 

1)

(

: очевидно.

 

 

 

 

 

 

 

|

 

 

{z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2)

)

: от противного. Предполагаем, что

^

^

 

 

 

 

 

 

 

 

9t: t1

< t <

Берем

U

окрестность точки ^ такую, что функция

 

 

 

 

 

 

 

t

 

 

непрерывная функция.) Берем функцию (t):

 

(t) непрерывно дифференцируема на [t1; t2]

 

 

(t) > 0

íà

[t1; t2]

,

^

 

 

 

 

 

(t) 6= 0

 

 

 

(t) = 0

âíå U.

 

 

 

 

^ 6 t2; (t) = 0.

(t) сохраняет знак (напоминаем, что

На концах функция равняется 0. ( (t1) = (t2) = 0:)

t2

(t) (t)dt = Z

 

Z

(t) (t) 6= 0

t1

U

42

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

Мы пришли к противоречию, действительно выполняется уравнение Эйлера. Определение: Каждое решение уравнения (2) называется экстремалью функционала F .

Замечание: Экстремалью называется не решение задачи на экстремум, а решение уравнения Эйлера!

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

Получаем задачу

 

 

 

(2) ^ y(t1) = y1 ^ y(t2) = y2 :

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

Задача:

Äàíî: t1 6 t 6 t2, y(t), y(t1) = y1, y(t2) = y2.

t2

Z

p

F = 1 + (y0)2dt

t1

Нужно найти кратчайшую кривую, соединяющую эти точки. (Функционал, дающий минимум

F )

=0

=00

 

y=00

 

0

0

 

p

0

 

fy

fy t

f

yy0

 

fy

y

y00

= 0; f = 1 + (y

)2

|{z}

|{z}

| {z }

 

 

 

 

 

 

 

fy0y0y00 = 0:

 

fy0 =

1 2y0

 

 

=

 

 

 

y0

 

2

 

p

 

 

 

p

 

 

 

 

 

 

 

 

 

1 + (y0)2

1 + (y0)2

 

 

 

 

 

 

 

 

 

 

 

y0

 

 

 

 

 

 

 

 

1 p1 + (y0)2 y0

 

 

 

 

=

1 + (y0)2 (y0)2

 

 

 

 

 

 

 

fy0y0 =

 

1 + (y0)2

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( 1 + (y0)2)3

 

 

 

1 + (y0) p

 

 

Решением y00 = 0 является y = c1t + c2, c1, c2

 

 

 

 

 

 

 

p

константы. Таким образом, кратчайший путь

прямая (отрезок).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Замечание: Мы на самом деле ничего не доказали. Мы предположили, что этот путь задается однозначной функцией y = y(t), то есть отбросили варианты вида

1. Простейшая задача вариационного исчисления

43

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

3. Функции, непрерывно дифференцируемые, удовлетворяющие данным условиям, называются допустимыми.

 

 

 

F min:

 

 

(1)

y(t) решение (1) )

 

 

d

 

 

 

 

 

 

 

fy =

fy0 = 0

 

 

(2)

 

 

 

 

 

 

 

 

 

dt

 

 

 

d

fy0 = fy0

+ fy0yy0 + fy0y0y00

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

Случаи понижения порядка.

 

 

 

 

 

 

 

1. f не зависит от y, f(t; y; y0) = g(t; y0)

 

 

 

 

 

 

 

dg0

 

 

 

 

 

 

 

 

 

y

= 0

 

 

gy0(t; y) = c

 

const

(20)

dt

,

 

 

 

 

 

2.f не зависит от t, f(t; y; y0) = g(y; y0) Рассмотрим выражение g y0gy0. Найдем полную производную по t этого выражения.

 

d

 

(g y0gy0) = gy y0 +

gy0

y00

 

y00 gy0

y0gy0yy0 y0gy0y0y00 = y0 gy gy0yy0 gy0y0y00

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

gy

 

gy0yy0

 

gy0y0y00

= 0

 

(200)

 

 

 

 

d

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(g

 

y0gy0) = 0;

 

 

g

 

y0gy0

= c

 

const

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(200)

)

 

 

g

 

y0gy0 = const

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

g y0gy0 = const ^ y0 6= 0

) (200):

 

(a) y0 = 0 в отдельных точках. Тогда g

 

y0gy0

,

(200)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

44

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

(b) y0

= 0 на некотором отрезке. При y = y0 const выражение g y0gy0 равно нулю.

(200) выполняется, если gy(y0; 0) = 0. То есть нужно найти те y0, для которых gy(y0; 0) = 0.

4. Случай квадратичного функционала.

f многочлен второй степени относительно совокупности y; y0.

f(t; y; y0) = p(t)(y0)2 + q(t)y2 2m(t)y:

Замечание: Можно возразить, что выписаны не все члены. Но можно показать, что все функции можно привести к такому виду:

t2

a(t)yy0dt = 2

a(t)y2(t)

 

 

t1

2

t2

a0(t)y2(t)dt;

Z

i

Z

t1

 

1

h

 

t2

1

t1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t2

b(t)y0dt =

b(t)y(t)

 

 

 

 

t2

b0(t)y(t)dt;

 

Z

i

t1

Z

 

t1

 

 

 

h

t2

t1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t2

 

 

 

 

 

 

 

 

 

 

Z

c(t)dt не зависит от выбора y(t):

t1

Далее предполагаем, что p0(t), q(t), m(t) непрерывны на отрезке [t1; t2].

Вопрос: Как выглядит уравнение Эйлера в этом случае?

 

Ответ:

 

2qy 2m 2(py0)0 = 0

 

(p(t)y0)0 + q(t)y = m(t) , py00 p0y0 + qy = m

(3)

Теорема: Пусть p(t) > 0 и q(t) > 0 на [t1; t2]. Пусть y(t) решение уравнения (3), удовлетворяющее условиям y(t1) = y1, y(t2) = y2. Тогда y(t) дает минимум F .

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

Доказательство: Берем нашу функцию y(t), какую-то (t), такую, что (t1) = (t2) = 0, (t)непрерывно дифференцируемая, число (любое).

def

y~(t) = y(t) + y(t)

Приращение функционала представляем в виде суммы линейной части и членов более высо-

кого порядка.

t2

Z

F (~y) = F (y) + F + 2 p(t)( 0)2 + q(t) 2 dt > F (y)

|{z}

=0 t1

1. Простейшая задача вариационного исчисления

45

Мы доказали, что функция дает минимум. Важно: имеется ввиду глобальный минимум. Если(t) не тождественно нулевая функция, то F (~y) > F (y).

Замечание: Можно ли гарантировать, что существует решение уравнения (3), удовлетворяющее граничным условиям? Пока мы не можем ответить на этот вопрос. Во втором семестре мы докажем, что уравнение (3) имеет решение, и притом единственное.

5. Задача о брахистохроне.

Замечание: В 1696 году Иоганн Бернулли опубликовал работу задача о брахистохроне . Считается, что эта работа положила начало вариационному исчислению. Была рассмотрена

следующая задача: в вертикальной плоскости лежат точки A, B, соединенные некоторой кри-

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

Казалось бы, логичнее всего двигаться по отрезку прямой. Оказывается, это не так. Поместим точку A в начало координат, вертикально вниз направлена ось y, где-то задана точка B.

A(0; 0); B(x1; y1)

Ограничимся кривыми, задаваемыми уравнениями y = y(x); 0 6 x 6 x1.

v скорость, v = v(x; y).

mv2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= mgy

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x1

 

 

 

 

 

 

v = p2gy;

T = Z`

dt = Z`

 

v

 

 

s

y 0

 

 

dx

 

= p2g Z0

)2

 

 

 

 

 

 

 

x1

 

 

 

 

 

ds

1

 

 

 

1 + (y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F (x) = p2g

Z0

p

py 0

)2

 

dx;

F min

 

 

 

 

 

 

 

1

 

1 + (y

 

 

 

 

 

 

 

 

 

 

 

 

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

46

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

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

s

(y

)2

 

 

 

 

 

2y

 

 

 

 

1 + 0

 

y0

 

 

 

 

0

 

 

 

= c const:

y

 

2

 

p

y(1 +

(y

)2)

 

 

 

 

 

 

 

0

 

 

 

Сразу напишем ответ (легко проверить подстановкой). Рассмотрим функции, задаваемые уравнениями

 

 

x = R(t sin t) + Q

;

y = R(q cos t)

 

где R и Q const, t параметр на решение.

 

Эта кривая называется циклоидой. (Мы будем брать из нее кусочки) Возьмем Q = 0, тогда y(0) = 0. R =?

Нарисуем циклоиду, соответствующую R = 1, потом изменим нашу циклоиду (растяжением

AB

или сжатием) так, чтобы R = AC : R определена однозначно условием y(t1) = y1.

Оказывается, брахистохроной является часть циклоиды.

Замечание: Так как Бернулли жил до Эйлера, то он не мог воспользоваться тем приемом, которым пользовались мы, он делал это как-то по-другому. С другой стороны, мы тоже нестрого решили задачу, так как мы пользовались дополнительными предположениями. Путем кропотливой работы можно показать, что найденное нами решение действительно является мини-

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

47

мальным.

 

 

 

 

 

 

Rt3

Rt3

 

Замечание: При x близком к 0, 0, x

 

, y

 

.

 

6

2

 

x = o(y), касательная вертикальна в начальной точке.

[ На лекции было описано геометрическое построение R ] Дадим формулу для времени:

 

x = R(t sin t);

y = R(q cos t):

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

s

y

d = p2g Z0

0p2Rd = s g

0

t = Z0

 

 

(x)2

+ (y)2

 

1

 

 

 

 

R

Мы до сих пор рассматривали случай, когда B намного ниже точки A. Интересно рассмотреть

случай, когда x1 íå ìàëî, y1 мало. Точка падает, разгоняется, и доходит до точки B. Таким образом, под действием силы тяжести можно передвигаться в горизонтальном направлении.

AB = L; R =

L

; 0

= 2

2

ss

T =

L

2 =

2 L

2 q

g

Замечание: L = 600 êì, 600 ñåê 10 ìèí.

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

Замечание: В этом параграфе все используемые функции предполагаются достаточно гладкими. Если захотеть, можно четко посмотреть, где мы пользовались этими предположениями. Этого не делается, чтобы не отвлекаться от основной задачи.

48

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

1. Функционалы, зависящие от нескольких функций.

 

t1 6 t 6 t2;

~y(t) =

 

y1(t)

;

F (~y) =

t2 f

t; ~y(t); ~y 0(t) dt;

 

 

y (t)

 

 

Z

 

 

 

 

 

m

 

 

 

 

 

 

 

 

t1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

 

ft1 6 t 6 t2;

~y;

è ~y 0 любыеg:

 

 

~

~

~

 

~ ~ ~

 

 

 

 

y(t1) = a;

y(t2) = b; a; b заданы.

 

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

 

 

 

maxt

~y(t) ~y(t)

+ maxt

~y 0

(t) ~y 0

(t)

< ";

имеет место

 

 

 

 

 

 

 

 

 

 

s

 

 

 

 

F (~y) 6 F ( ~y):

 

 

 

y1(t)

 

 

 

 

 

 

m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P

Здесь j~uj = , ~u =

k=1

 

 

ym(t)

F min.

(1)

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

fyk

 

d

fy0 = 0; k = 1; : : : ; m:

(2)

 

 

 

dt k

 

dtd полная производная по t .

Доказательство: Доказательство аналогично тому, что происходит в математическом анализе. Фиксируется только одна переменная. Раз достигается минимум по совокупности, значит достигается минимум по данной конкретной переменной. Отсюда следует, что соответствующая производная равна нулю. Вместо слов переменная используется допустимая функция , и так далее.

Определение: Каждое решение системы (2) называется экстремалью функционала F .

Замечание: Не гарантируется, что экстремаль дает минимум.

2. Функционалы, зависящие от производных высших порядков. Фиксируем отрезок t1 6 t 6 t2, y(t), рассматриваем функционал

t2

 

 

F (y) = Z

f

t; y(t); y0(t); : : : ; y(n)(t) dt:

t1

 

 

f задана, определена при t1 6 t 6 t2;

y; y0; : : : y(m) любые.

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

49

Определение: Функция y(t) называется допустимой, если

 

y(t1) = a0;

y0(t1) = a1; : : : ;

y(m 1)(t1) = am 1;

 

y(t2) = b0;

y0(t2) = b1; : : : ;

y(m 1)(t2) = bm 1;

 

ãäå a0; a1; : : : ; am 1; b0; b1; : : : ; bm 1 заданы.

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

maxt

y(t) y~(t)

+ maxt y 0(t) y~0

(t)

< ";

имеет место

 

 

 

 

 

 

 

F (y) 6 F (~y):

 

 

Задача нахождения минимума:

 

 

 

 

 

 

F (y) min.

 

(3)

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

 

d

 

d2

m dm

 

fy

 

fy0

+

 

fy00 : : : + ( 1)

 

 

fy(m) = 0:

(4)

dt

dt2

dtm

dtd полная производная по t .

Доказательство: Мы взяли нашу функцию y(t). Берем какую либо гладкую функцию (t), удовлетворяющую условию

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

и составляем семейство функций

def

y~ = y(t) + (t); число (малое).

Составляем вспомогательную функцию

def 0

!( ) = F (y + ) = !(0) + ! (0) + o( ); !(0) = F (y)

50

 

 

 

 

 

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

!0(0) =

2dt Z f(t; y + ; y0 + 0; : : : ; y(m) + (m))dt3

 

 

 

t2

 

 

 

 

 

 

 

d

 

 

 

 

5 =0

 

4t2

t1

 

 

 

 

 

Z

 

 

 

 

 

 

 

 

= (fy + fy0 0 + fy00 00 + : : : + fy(m) (m))dt

 

t1

 

 

 

 

 

 

 

 

 

 

 

Используем обозначения fy = fy(t; y; y0; : : : ; y(m))

 

 

Затем интегрируем по частям:

 

 

 

t2

 

 

 

 

d2

dm

= tZ1

 

d

 

fy

 

fy0

+

 

fy00 : : : + ( 1)m

 

fy(m) +

dt

dt2

dtm

Когда интегрируем по частям, понижаем порядок , и,

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

hit2

+fy0

t1

t2

Z

t1

 

 

 

t2

d

 

 

t2

 

 

 

 

 

t2

d fy(m

 

1)

 

 

t2

+ hfy00 0it1 h

 

fy00 it1 + : : : + hfy(m) ((m 1) it1 h

 

 

 

(m 1)it1 + : : :

dt

dt

 

 

|

 

{z

 

} |

 

{z

 

}

d |

 

{z

 

 

} |

 

 

 

{z

 

 

}

 

 

 

 

 

 

 

dm

 

 

 

=0

 

 

 

 

=0

 

 

 

=0

 

 

 

 

 

 

=0

 

 

(t) (t)dt = 0;

 

 

(t) = fy

 

fy0 + : : : + ( 1)m

 

fy(m) :

 

 

 

 

dt

dtm

 

 

Покажем, что (t) = 0. Допустим, это не так,

^

^

:

(t) 6= 0; t1

< t < t2

 

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

 

(t)

сохраняет знак в окрестности

U

точки ^.

 

 

 

 

t

 

 

 

 

(t): (t) 6 0

íà

 

^

âíå

U:

 

 

[t1; t2]; (t) > 0; (t) = 0

 

(t) гладкая!

 

 

Z

 

 

 

t2

 

 

 

 

 

Z

 

 

 

 

 

 

(t) (t)dt =

(t) (t)dt 6= 0:

 

 

 

t1

 

 

U

 

 

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

Определение: Каждое решение уравнения (4) называется экстремалью функционала F .

Замечание: Эта задача разумная. Но не утверждается, что эта задача имеет решение, или даже единственное. Если удастся найти решение, отсюда не следует, что на решении достигается минимум. (Оно лишь является кандидатом). Проблема достаточности не затрагивается.