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

Variatsionnoe_ischislenia

.pdf
Скачиваний:
22
Добавлен:
18.03.2015
Размер:
899.48 Кб
Скачать

V.3 Необходимые условия минимума в простейшей вариационной задаче

V.3.1 Условие Лежандра

Общую теорему о необходимом условии минимума мы перенесем на простейшую вариационную задачу

b

Z

J : x → L(t, x(t), x˙ (t))dt, x : [a, b] → R, x(a) = xa, x(b) = xb

a

Отметим, что необходимые условия слабого локального минимума будет одновременно необходимым условием слабого глобального минимума и сильного локального и глобального минимумов.

Мы предполагаем достаточную гладкость функции L(t, x, v). Пусть x = xˆ(t)- экстремаль:

 

 

 

 

 

˙

d

˙

(V.4)

 

 

DJ (ˆx) = 0, D2L(t, xˆ(t),

xˆ(t)) −

dt

D3L(t, xˆ(t), xˆ(t)) = 0

Вычислим вторую производную функционала.

 

 

 

 

 

d2

 

 

 

 

 

D2J (ˆx)(h, h) =

 

J (ˆx + sh)|s=0 =

 

 

 

 

ds2

 

 

 

 

 

 

b

 

 

 

 

 

 

 

= Z

D22L(t, xˆ(t), xˆ(˙ t)) · h2

(h) + 2D2D3L(t, xˆ(t), xˆ(˙ t)h(t)h˙ (t))+

 

 

 

a

h

 

 

+D32L(t, xˆ(t)xˆ(˙ t)) · h˙ 2(t)i dt

 

 

 

 

 

 

 

(V.5)

 

 

 

 

˙

 

ˆ

 

Далее для краткости будем обозначать:L(t, x(t), xˆ(t)) через

L(t), а если x = xˆ(t), то L(t)

 

Пусть экстремаль дает локальный минимум. По теореме 24 необходима неотрицательность квадратичной

формы

2

 

(ˆx)(h, h)

h C(a, b).

D

J

2

 

0

 

 

 

 

 

0 6 D

J (ˆx)(h, h) =

 

 

 

 

 

 

b

hD32Lhˆ ˙ 2(t) + D2Lˆ(t)h2(t) + D2D3Lˆ(t)

d

(h2(t))i dt =

= a

dt

Rb

nD32Lˆ(t)h˙ 2(t) + hD22Lˆ(t) −

d

D2D3Lˆ(t)i h2(t)o dt.

= a

dt

Возьмем приращение h

специального вида

R

 

 

 

hε(t) = ϕ( t − t0 )

ε

с фиксированной функцией ϕ C0, t0 (a, b) и достаточно малым ε, таким чтобы supphε (a, b). На таких

функциях имеем

D2J (ˆx)(hε, hε) =

Za

b

3

 

ε

2

 

ε2

 

 

=

1

D2Lˆ(t)ϕ2

(

t − t0

) + D2Lˆ

(t)

 

 

 

 

d

dt D2D3

 

 

ε

 

Lˆ(t) ϕ2

(

t − t0

) dt =

 

 

 

=

1

 

ϕ 2

(τ )dτ + o(1)

при ε → +0

 

 

ε D32Lˆ(t0) Z

 

 

 

 

 

 

 

Таким образом, из необходимого условия минимума функционала теоремы 24 следует

 

2 ˆ

0

для всех t [a, b]

 

 

 

 

 

Теорема 27. D3 L(t) 1

 

 

 

 

 

Это неравенство называется неравенством Лежандра 2) (усиленным условием Лежандра называется строгое

2 ˆ

при t [a, b]).

неравенство: D3 L(t) > 0

2) Legendre A. M.(Лежандр Адриен Мари, 1752 - 1833) французский математик, исследовал "преобразования Лежандра"в 1789

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

V.3.2 Условие Якоби

Получим еще одно следствие необходимого условия теоремы 24 неотрицательности квадратичной формы второй

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

˙

производной функционала на экстремали. Сначала введем сокращенные обозначения P (t) = D3 L(t, xˆ(t), xˆ(t)), Q(t) =

2

˙

d

 

 

˙

 

 

 

 

 

 

 

 

 

 

D2 L(t, xˆ(t), xˆ(t)) −

dt

D2D3L(t, xˆ(t), xˆ(t)).

 

 

 

 

 

 

 

 

 

 

 

 

 

2

b

h

 

 

 

Q(t)h

2

(t)i

2

 

 

 

 

 

 

a

 

 

 

 

 

В этих обозначениях D

 

J (ˆx)(h, h) =

P (t)h

(t) +

 

 

 

dt.

 

 

 

 

 

 

R

 

 

 

 

 

 

 

 

 

 

 

Уравнение Эйлера для квадратичного функционала h → D J (ˆx)(h, h) есть

 

 

 

 

 

 

 

 

d

(P h) + Qh = 0

(V.6)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

Рассмотрим его на функциях с начальными данными ˙ . Определим точку как сопря- h h(a) = 0, h(a) = 1 t > a

женную с точкой a, если введенная выше функция h обращается в нуль в точке t .

˙

Теорема 28. Пусть выполнено усиленное условие Лежандра. Для того, чтобы экстремаль x = xˆ(t) давала слабый локальный минимум функционалу J (x) необходимо, чтобы интервал (a, b) не содержал точек, сопряженных с a.

Отсутствие сопряженных точек на интервале (a, b) называется условием Якоби3). Усиленным условием Якоби называется отсутствие сопряженных точек на полуинтервале (a, b]. Само уравнение Эйлера для квадратичного функционала, образованного второй производной функционала J на экстремали(V.6), называется уравнением Якоби функционала J .

Доказательство. Ранее было показано, что неотрицательность квадратичной формы второй производной функционала на экстремали является необходимым условием минимума. Значит, D32J (ˆx)(h, h) > 0. Рассмотрим про-

b ˙ 2

¨

имеет общее решение

стейший квадратичный функционал K : h → a h

(t)dt. Его уравнение Эйлера h(t) = 0

R

 

 

h(t) = αt + β. Очевидно, у функционала K вообще нет сопряженных точек. Образуем линейную комбинацию функционалов D2J (ˆx) и K

Ms : h → sD2J (ˆx) + (1 − s)K =

b

Z

n

˙2

=[sP (t) + (1 − s)]h (t)

o

+ sQ(t)h2(t) dt (V.7)

a

Для всех s [0, 1) квадратичный функционал Ms положительно определен на гильбертовом пространстве

W21,0[a, b] = {h|h : [a, b] → R, h(a) = 0,

s

 

 

 

0 dt[h˙ 2

(t) + h2(t)] ≡ khk < ∞}.

 

 

 

 

 

1

 

 

Покажем это.

 

 

 

 

R

 

 

b

b

t

dτ h˙ (τ ))2 6 Z

t

 

 

 

 

Z

dth2(t) = Z

dt(Z

dτ h˙ 2(τ ) =

 

 

a

a

a

a

b

Z

˙ 2

= dth (τ )

a

Поэтому

b

 

 

 

(b − a)2

 

Za

dt(t

a) =

(V.8)

2

 

 

 

1

 

Ms(h, h) = (1 − s) Z0

dth˙ 2(t) >

> (1 − s)

 

2

b

dth˙

2(t)|2

·

(b

a)2

 

b

dth2(t)

>

 

 

 

Z

Z

 

 

 

 

 

1

 

 

1

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

a

 

 

 

 

b

 

 

 

 

 

 

 

 

 

 

 

1

 

 

1

} · Z

dt h˙ 2

(t) + h2(t)

≡ α · khk2

 

 

 

 

 

> (1 − s)min{

 

,

 

 

 

 

 

 

2

(b

a)2

 

 

 

 

 

 

 

 

 

 

 

 

 

a

h

 

i

3)Jacobi С.П.(Якоби Карл Густав, 1804 - 1851) немецкий математик, известны: Уравнение Якоби, Условие Якоби, определитель

-Якобиан

Рис. V.2:

При s = 1, Ms(h, h) совпадает с D2J (ˆx)(h, h) и является вещественной неотрицательной квадратичной фор-

мой. Уравнение Эйлера функционала Ms есть

 

d

˙

dt

[(sP + (1 − s)) h]J + sQh = 0

Пусть h(t, s) решение этого уравнения с начальными данными Коши h(a, s) = 0, D1h(a, s) = 1. Эта функция непрерывна по (t, s) в прямоугольнике Π = [a, b] × [1, 0], равна t − a при s = 0 и переходит в h(t) при s = 1.

Заметим, что если в некоторой точке(t0, s0) Π h(t0, s0) = 0, то необходимо D1h(t0, s0) 6= 0. Иначе, по теореме единственности решения задачи Коши для дифференциальных уравнений было бы h(t, s0) ≡ 0, и в частности, 0 = D1h(a, s0) = 1

Теперь будем рассуждать от противного. Пусть на интервале (a, b) есть сопряженная с a точка t , то есть

h(t , 1) = 0. Множество l = {(t, s)|h(t, s) = 0} есть некоторая замкнутая кривая в

. На ней всегда D1h(t, s) 6= 0

и по теореме о неявных функциях она представима в окрестности любой

своей точки в виде, разрешенном

 

Q

относительно t : t = T (s), c непрерывно диференцируемой функцией T. Эта кривая l начинается в точке (t , 1).

Но где кончается? См.рис.(V.2)

l не может окончится внутри Q - по теореме о неявных функциях.

lне может окончится на стороне s = 1, так как по пути прошла бы через точку (t0, s0) с D1h(t0, s0) = 0.

lне может окончится на стороне t = a, так как всей этой стороне h(a, s) ≡ 0 и в точке встречи l с этой стороной

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

lне может окончится при s = 0, так как h(t, 0) = t − a не обращается в нуль для t > a.

Покажем что l не может окончиться и на стороне t = b, например в точке (b, s), 0 6 s < 1. Действительно в этом случае уравнение Эйлера квадратичного функционала Ms имело бы ненулевое решение с нулевыми граничными условиями h(a, s) = h(b, s) = 0. И мы получаем

b

dt [(sP + (1 − s))h˙

] + sQh hdt =

 

 

 

 

 

0 = Za

 

 

 

 

 

 

 

d

 

 

 

 

 

 

 

 

 

 

a

 

(sP + (1 − s))h˙

 

b

 

 

 

 

 

= Z

n

2 + Qh2

dt > (1 − s) Z

h˙ 2(t)dt,

 

 

 

 

b

 

o

a

 

Что влечет h(t) ≡ const = 0-противоречие.

V.4 Достаточные условия минимума простейшей вариационной задачи

V.4.1 Достаточное условие слабого локального минимума

Мы конкретизируем сформулированные ранее общие достаточные условия слабого локального минимума на

b

примере простейшей вариационной задачи для функционала J (x) = R Ldt. Накладываемые ограничения, есте-

a

ственно, должны включать все необходимые условия, и, как увидим, будут достаточно близки к ним.

Теорема 29. Пусть x = xˆ(t) является экстремалью, то-есть решением уравнения Эйлера

d

D2L − dt D3L = 0

и для этой экстремали выполняются усиленные условия Лежандра и Якоби. Тогда функция x = xˆ(t) реализует слабый локальный минимум функционала J (x).

Доказательство. Так как полуотрезок (a, b] не содержит сопряженных с a точек, а P (t) > 0 на [a, b], то можно указать больший отрезок [a, b + ε], на котором нет сопряженных точек и P (t) > 0. Рассмотрим квадратичный

b

˙ 2

2

2

b ˙ 2

dt и его уравнение Эйлера.

функционал Nα : h →

(P h

+ Qh

)dt − α

h

R

 

 

 

R

 

a

 

 

 

a

 

d2 ˙

dt (P − α )h + Qh = 0

По свойству непрерывной зависимости решения дифференциального уравнения от параметров, при достаточно малых α и ε будет P (t)−α2 > 0 на [a, b] и решение h(t, α) написанного уравнения Эйлера с начальными данными

h(a − ε) = 0,

 

˙

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

h(a − ε) = 1 не будет иметь сопряженных с (a − ε) точек на отрезке [a, b]. Это позволяет нам

привести функционал Nα к виду

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

 

 

 

w

 

 

 

 

 

 

 

 

 

 

 

Nα : h → Za

(P − α2)(h˙ +

 

 

 

 

 

 

 

 

 

 

 

 

h)2dt.

 

(V.9)

 

 

 

 

 

 

 

 

(P − α2)

Действительно, достаточно, чтобы неопределенная пока функция w удовлетворяла условию

˙

2

 

1

 

2

 

2

 

d

2

 

 

 

 

 

 

 

 

 

2whh + w

 

·

 

· h

 

= Qh

 

+

 

(w · h

) =

 

 

 

 

 

 

 

 

 

P − α2

 

 

dt

 

2

 

˙

 

 

 

 

 

 

 

 

 

 

 

 

 

= (Q + w˙ ) · h

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+ 2whh

 

 

То есть w должно быть каким-то решением дифференциального уравнения w˙ =

1

 

w2 − Q, Это дифферен-

P −

α

циальное уравнение Риккати4).Его можно свести к линейному уравнению заменой переменных по формуле w = −uu· (P − α2). При этом для u получаем уравнение dtd (P − α2)u+ Qu = 0 Это как раз уравнение Эйлера для функционала Nα и за u мы можем взять h(t, α). Оно не обращается в нуль на [a, b] и замена w на u законна.

Установленная для Nα формула (V.9) означает неотрицательность Nα

 

b

b

 

 

Za

(P h˙ 2 + Qh2)dt − α2 Za

h˙ 2dt > 0,

Отсюда и по оценке (V.8) получаем

 

 

D2J (ˆx)(h, h) > α2 · Z

b

 

 

dth˙ 2(t) >

 

 

a

α2

> 2

b

Z

˙2

dth (t)

a

+ (b

 

b

dth2(t) >

2

 

min 1, (b

2a)2

b

a)2

 

dt[h˙ 2(t) + h2(t)]

2

 

Z

 

α2

 

 

 

Z

 

a

 

 

 

 

 

a

Теперь по формуле (II.12)для приращения функционала имеем

 

 

 

b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Z

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

h

 

 

 

 

 

 

 

 

 

 

 

 

2

2

2

˙ 2

 

 

 

2

(t) >

 

 

 

 

 

 

J (ˆx + εh) − J (ˆx) = ε

D

J (ˆx)(h, h) + o(ǫ

) · dt h

(t) + h

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

α

2

2

b

 

 

 

 

· min 1;

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

>

ε

 

Z dt

h˙ 2

(t) + h2

(t)

 

· (1 + o(1) > 0

 

 

 

 

2

 

(b

a)2

 

 

 

 

 

 

 

a

 

h

 

i

 

 

 

 

V.4.2 Достаточное условие сильного минимума.

Локальный сильный минимум – это минимум среди кусочно непрерывно дифференцируемых функций, взятых из малой окрестности, определяемой метрикой непрерывных функций. Эта окрестность может быть очень большой в метрике кусочно непрерывно дифференцируемых функций. Например, x(t) = ε sin nt при малых ε близко

к нулю в норме kxkc =

max ε sin nt

 

 

Однако

 

1

-норма таких функций может быть сколь угодно большой

t

 

|

 

 

| = ε

 

 

C

 

при достаточно большом n

k

x

1

=

max

{|

εsinnt , εn cos nt

|}

= εn.

 

 

kC

 

t

| |

 

 

 

4) Rikatti J.(Рикатти Я.) Впервые исследовал уравнение Рикатти в 1723 г.

Каждый сильный минимум является одновременно и слабым. Поэтому естественно искать достаточные усло-

вия сильного минимума, "взяв за основу"полученные ранее условия слабого минимума.

a

Мы рассматриваем простейшую вариационную задачу J (x) = R L(t, x(t), xˆ(t))dt → min, x(a) = xa, x(b) = xb

b

Теорема 30. Пусть L C4(a, b]×R2), функция x = xˆ(t) -класса C3[a, b] является экстремалью функционала J , Выполнены усиленные условия Лежандра и Якоби . Пусть дополнительно функция L(t, x, v) является строго выпуклой по v R1 в некоторой окрестности v R2 графика экстремали = {(t, x)|x = xˆ(t), t [a, b]}. Тогда эта экстремаль дает локальный сильный минимум в поставленной задаче.

Пример 19. Пример сильного минимума функционала.

 

 

 

 

T0

 

 

 

 

 

 

 

 

 

 

 

Z0

(x˙ 2 − x2)dt → inf,

x(0) = 0, x(T0) = ξ, T0 > 0.

 

 

 

 

Решение:

2

 

 

есть x¨ + x

= 0. Его общее решение x(t) = c

 

sin t + c

 

cos t. Условие

L(t, x, p)

p2

, уравнение Эйлера

 

 

 

 

− x

 

2

 

1

 

2

 

Лежандра, даже усиленное, выполнено:

D3 L ≡ 2 > 0. Проверим условиe Якоби. Так как наш функционал

квадратичный, то уравнение Якоби совпадает с уравнением Эйлера. Его ненулевое решение h(t) с начальным условием h(0) = 0 есть h(t) = A sin t с A =6 0. Сопряженные точки удовлетворяют условию sin t = 0, ближайшая к нулю сопряженная точкa t0 = π. Поэтому, если T0 < π, то xˆ(t) = ξ sin t/ sin T0 дает сильный минимум.

Перед доказательством теоремы введем понятие поля экстремалей. Запишем уравнение Эйлер функционала

J

D32L(t, x, x˙ )¨x + D2D3L(t, x, x˙ )x˙ + D1D3L(t, x, x˙ ) − D2L(t, x, x˙ ) = 0.

2

˙

3

экстремали

Так как D3 L(t, xˆ(t), xˆ(t)) > 0, то по непрерывности в некоторой окрестности U R

 

˙

{(t, x, y) | t [a, b], x = xˆ(t), y = xˆ(t)}

продолжает соблюдаться условие D32L(t, x, y) > 0. Поэтому в этой окрестности уравнение Эйлера эквивалентно

системе, разрешенной относительно производных

x˙ = y, y˙ = Φ(t, x, y) =

1

[D2L(t, x, y)−

 

D32L(t, x, y)

 

 

− D1D3L(t, x, y) − D2D3L(t, x, y) · y].

Благодаря предположенной гладкости функции L записанная система удовлетворяет условиям известных тео-

рем локальной однозначной разрешимости задачи Коши и продолжимости решения. То есть найдутся такие числа ε > 0 и δ > 0, что решение x = xˆ(t) продолжимо на отрезок [a − ε, b + ε] и для любого λ R1, |λ| < δ, на отрезке [a−ε, b+ε] однозначно определено решение x = x(t, λ) задачи Коши с начальными данными x(t ) = xˆ(t ),

˙

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x˙ (t ) = xˆ(t ) + λ, где t некоторая точка интервала (a − ε, a). По теореме о дифференцируемости решения

по параметрам функция (t, λ) → x(t, λ) непрерывно дифференцируема.

2Будем говорить, что экстремаль xˆ(t)

окружена полем экстремалей x(t, λ), если существует окрестность G R

графика экстремали

 

 

 

 

 

 

 

 

ˆ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= {(t, x) | t [a, b], x = xˆ(t)}

 

 

 

 

 

 

такая, что для любой точки (τ, ξ)

G

имеется единственная

экстремаль семейства

{

x(t, λ)

, проходящая через

 

 

C

1

 

 

}λ

ξ)

 

λ = λ(τ, ξ));

эту точку. Точнее,

существует функция λ = λ(τ, ξ) класса

 

(G) такая, что (x(τ, λ) =

 

n

 

d

 

 

 

 

 

 

 

 

 

 

 

 

 

функция u: G → R

 

, u(τ, ξ) =

 

x(t, λ(τ, ξ))

t=τ , называется функцией наклона поля. Если существует такая

 

dt

точка (t , x ), что x(t , λ) = x

для всех λ,

то говорят, что xˆ(t) окружена центральным полем экстремалей.

 

 

 

 

 

 

 

семейство x(t, λ) центральным полем экстремалей.

 

Точка (t , x ) называется центром поля, а

 

 

 

 

 

 

 

 

 

 

 

 

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

Лемма 1. Экстремаль x = xˆ(t) окружена центральным полем экстремалей.

Доказательство. Обозначим Dλx(t, λ)|λ=0 ≡ H(t, t ).

Дифференциальное уравнение для H получается дифференцированием по параметру λ уравнения Эйлера:

Dλ[−dt D3L(t, x(t, λ), x˙ (t, λ)) + D2L(t, x(t, λ), x˙ (t, λ))] λ=0 =

 

d

 

 

 

 

 

 

 

d

2

˙

˙

= −dt [D3 L(t, xˆ(t), xˆ(t))H (t, t ) + + D3D2

˙

D2D3L(t, xˆ(t), xˆ(t))H(t, t )]+

˙ ˙ 2 ˙

L(t, xˆ(t), xˆ(t))H (t, t ) + D2 L(t, xˆ(t), xˆ(t))H(t, t ) = 0.

H(t , t ) = Dλx(t , λ)|λ=0

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

= Dλxˆ(t ) = 0,

D1H(t , t ) = Dλx˙ (t , λ)|λ=0 = Dλ(x˙ (t ) + λ) = 1.

Пусть решение уравнения Якоби с условиями , ˙ . Выполнение усиленного

H(t, t0) H(a, a) = 0 H(a, a) = 1

условия Якоби означает, что H(t, a) 6= 0 при любом t (a, b]. По непрерывности H(t, t ) 6= 0 для любого t [a, b] при достаточной близости t к a.

Теперь возьмем отображение

(t, λ) 7→ψ(t, λ) ≡ (t, x(t, λ))

¯

¯

 

 

 

 

в некоторой точке (t, 0), t [a, b]

x˙ (t,¯

 

 

 

 

det Dψ(t,¯ 0) = det

0)

Dλx(t,¯

0)

 

 

1

 

0

 

= H(t,¯ t ) =6 0.

¯ ¯

¯

¯

 

Значит, определено обратное отображение в некоторой окрестности точки (t, xˆ(t))

vδ {(t, x) | |t−t| < δ, |x−xˆ(t)| <

δ}, то есть существует единственное λ = λ(τ, ξ) такое, что

 

 

 

(ψ(τ, λ(τ, ξ)) = (τ, ξ)) x(τ, λ(τ, ξ)) = ξ.

 

 

 

ˆ

¯

J

,

В силу компактности можно выделить конечное покрытие соответствующее конечному числу точек {tj }j=1

и построить единую для всей эктремали x = xˆ(t) функцию λ = λ(τ, ξ) C2. Построение центрального поля экстремалей завершено. Остается положить функцию наклона поля равной u(τ, ξ) = x˙ (τ, λ(τ, ξ)).

Доказательство теоремы 30. Нам надо показать, что для всех кусочно гладких кривых x = x(t), соединяющих точки (a, xa) и (b, xb) и лежащих в малой R2 - окрестности графика экстремали , будет

a

a

 

0 6 J (x) − J (ˆx) ≡ Zb

L(t, x(t), x˙ (t))dt − Zb

L(t, xˆ(t), xˆ(˙ t))dt

Выписанные интегралы можно записать в виде криволинейных интегралов второго рода(см. учебник Фихтенгольц Г.М. Основы математического анализа, т. II), взятых по кривым (t, x(t)) и (t, xˆ(t)) :

b

 

(b,xb )

Za

L(t, x(t), x˙ (t)dt =

Z

L(t, x(t), x˙ (t)|x=x(t)dt;

 

 

 

(a,xa)

 

 

b

b,xb

 

 

Za

L(t, xˆ(t), xˆ(˙ t) =a,xZa

L(t, x(t), x˙ (t)|x=ˆx(t)dt.

Идея доказательства состоит в том, чтобы записать второй интеграл в виде интеграла второго рода по той же кривой, что и первый интеграл. Тогда неотрицательность разности интегралов будет следствием неотрицательности разности подынтегральных выражений.

Для обоснования этих преобразований рассмотрим сначала функцию

τ

Z

S(τ, ξ) = L(t, x(t, λ), x˙ (t, λ)|λ=λ(τ,ξ)dt

t

Очевидно, она непрерывно дифференцируема, а интеграл второго рода от ее дифференциала не зависит от

 

 

 

 

 

(b,xb )

выбора пути интегрирования, а только от начальной и конечной точек этого пути. Интеграл

 

 

dS(t, x)

 

 

 

 

 

 

a,x

)

называется инариантным интегралом Гильберта. Вычислим дифференциал функции S,

(

R a

 

dS = D1S · dτ + D2S · dξ

 

 

 

 

 

 

τ

 

 

 

 

 

 

 

tZ

 

 

 

 

 

 

 

D1S(τ, ξ) = L(τ, x(τ, λ), u(τ, ξ)) + [D2L(t, x(t, λ), x˙ (t, λ))D2x(t, λ) · D1λ+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+D3L(t, x(t, λ), x˙ (t, λ))D2 x˙ (t, λ)D1λ]|λ=λ(τ,ξ) dt = L(τ, ξ, u(τ, ξ))+

 

 

 

 

t=τ

 

 

 

 

 

 

+ D3L(t, x(t, λ), x˙ (t, λ))D2 x˙ (t, λ)D1λ|t=t

+

 

 

 

 

t

D2L − dt D3L x=x(t,λ) · D2x(t, λ) · D1λdt.

+

Z

 

d

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Слагаемое с интегралом пропадает, так как x = x(t, λ) есть экстремаль, D2x(t , λ) = 0 из-за того, что x(t , λ) = xˆ(t , λ).

В итоге получаем

 

 

x(τ, λ) − D1x(τ, λ)

 

D1S(τ, ξ) = L(τ, ξ, u(τ, ξ)) + D3L(τ, ξ, u(τ, ξ)) ·

=

 

 

 

 

 

 

d

 

 

 

 

 

 

 

 

= L(τ, ξ, u(τ, xi)) − D3L(τ, ξ, u(τ.ξ)) · u(τ, ξ), (V.10)

здесь

d

x(τ, λ) =

d

ξ = 0.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Аналогично вычисляем

 

 

 

 

D2S(τ, ξ) =

τ

Z

=[D2L(t, x(t, λ), x˙ (t, λ)) · D2x(t, λ) · D2λ+

t

+ D3L(t, x(t, λ), x˙ (t, λ)) · D2x˙ (t, λ) · D2λ]dt =

τ

Zd

=(D2L − dt D3L)|x=x(t,λ) · D2x(t, λ) · D2λdt+

t

 

 

 

 

 

t=τ

=

 

 

+ D3L(t, x(t, λ), x˙ (t, λ))D2x(t, λ)D2λ|t=t

 

= D3L(τ, ξ, u(τ, ξ)).

 

=2Dξ x(t, λ(τ,|ξ)) = Dξ ξ = 1

 

D x(t, x)D2λ t=τ =

 

 

 

 

 

 

 

 

 

 

 

 

Так как криволинейный интеграл второго рода, взятый от дифференциала гладкой функции, не зависит от выбранного пути интегрирования, имеем

(b,xb )

 

 

 

(b,xb )

 

 

Z

dS(t, x)|x=x(t) = S(b, xb) − S(a, xa) =

Z

dS(t, x)|x=ˆx(t) =

 

(a,xa )

 

 

 

(a,xa)

 

 

 

 

 

(b,xb )

 

 

 

 

 

 

Z

 

 

 

 

 

=

{[L(τ, ξ, u(τ, ξ)) − D3L(τ, ξ, u(τ, ξ))u(τ, ξ)]dτ + +

 

 

 

 

(a,xa)

 

 

 

 

b

 

+D3L(τ, ξ, u(τ, ξ))dξ} |ξ=ˆx(τ ) =

 

 

 

 

 

 

 

 

= Za

[L(τ, xˆ(τ ), xˆ(˙ τ )) − D3L(τ, xˆ(τ ), xˆ(˙ τ ))xˆ(˙ τ ) + D3L(τ, xˆ(τ ), xˆ(˙ τ )) · xˆ(˙ τ )]dτ =

 

 

 

 

 

b

 

 

 

 

 

 

= Za

L(t, xˆ(t), xˆ(˙ t))dt ≡ J (ˆx).

Итак, мы получили

 

(b,xb )

 

(b,xb )

 

 

 

 

 

 

 

 

J (x) − J (ˆx) =

Z

L(t, x, x˙ (t)|x=x(t)dt −

Z

dS(t, x)|x=x(t) =

 

 

 

 

 

 

(a,xa)

 

(a,xa )

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

Z

 

 

 

 

 

 

 

 

 

 

 

=

{L(t, x(t), x˙ (t)) − [L(t, x(t), u(t, x(t))) − D3L(t, x(t), u(t, x, (t))) · u˙ (t, x(t))]−

 

b

 

 

 

 

 

 

 

 

 

 

 

 

− D3L(t, x(t), u(t, x(t))) · x˙ (t)}dt =

 

 

u)]

 

 

 

 

=

b [L(t, x, y)

L(t, x, u)

D3L(t, x, u)(y

dt

 

 

 

 

Z

 

 

 

y = x(t)

 

 

 

 

a

 

 

 

 

 

 

u = u(t, x(t))

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Теперь заметим, что мы предполагали выпуклость функции y → L(t, x, y), определенной для всех y при всех

из малой окрестности графика экстремали ˆ . Поэтому неотрицательна и стоящая под

(t, x) = {(t, x)|x = xˆ(t)}

интегралом функция это взятая в точке y разность между графиком выпуклой функции и касательной к ней в точке u.

Замечание 5. Стоящая под интегралом функция E(t, x, y, u) ≡ L(t, x, y) −[L(t, x, u) + D3L(t, x, u)(y −u)] называется функцией Вейерштрасса функционала J (x). Соотношение E(t, x, y, u) > 0 называется основной формулой

Вейерштрасса.

Если на экстремали x = xˆ(t) для всех u R1 функция Вейерштрасса5) строго больше нуля, то она будет

неотрицательной в некоторой окрестности графика экстремали ˆ что достаточно для выполнения неравенства

,

J (x) − J (ˆx) > 0

Задача 22. Неотрицательность функции Вейерштрасса на экстремали необходима для локального сильного

минимума. Обосновать самостоятельно.

Замечание 6. Теорема о достаточном условии локального сильного минимума обобщается на многомерный случай: x Rn, L : R2n+1 → R

 

1

 

 

Пример 20. Рассмотрим функционал J (x) =

3

(t)dt, x(0) = 0, x(1) = 1. L(t, x, v) = v3

, уравнение Эйлера:

2

R

 

 

 

0

 

 

Dt(x˙ (t)) = 0. имеет единственное решение xˆ(t) = t. Достаточное условие слабого локального экстремума заклю-

 

2 ˆ

¨

чается в выполнении усиленных условий Лежандра и Якоби.D3L(t) = 6ˆx(t) = 6 > 0, уравнение Якоби h(t) = 0.

˙

дают h(t) = t и, значит вообще нет сопряженных с нулем точек. Таким

Начальные условия h(0) = 0, h(0) = 1

образом xˆ(t) = t доставляет слабый локальный минимум. Установленные нами достаточные условия сильного локального минимума не выполнены, так как функция L = v3 не является выпуклой. Проверим необходимые

условия сильного локального минимума -условие Вейерштрасса(см. задачу 22 приведенную выше).

E(t, x, x,˙ u) = L(t, x, x˙ ) − [L(t, x, u) + D3L(t, x, u) · (x˙ − u)] =

= 1 − u3 − 3u2(1 − u) = 1 + u2(2u − 3) > 0? u R

Это условие не выполнено, например, для u = −1. Значит функция xˆ(t) = t не доставляет сильного локального

минимума.

Замечание 7. Общее понятие поля и связь его с известным методом прогонки численного решения дифферен-

циальных уравнений изложены в книге Гельфанд И.М., Фомин С.В. Вариационное исчисление, М., Физматгиз 1961г.

V.4.3 Совпадение слабого и сильного минимума вариационной задачи в случае простейшего квадратичного функционала

b

Мы рассматриваем вектор-функции x: R → Rn и функционал J : x 7→R L(t, x(t), x˙ (t))dt, x(a) = xa, x(b) = xb.

a

Функция L(t, x, v) является однородным многочленом второго порядка от переменных (x, v). Тогда он должен

иметь вид

b

Z

J (x) = [hA(t)x˙ (t), x˙ (t)i + 2hC(t)x˙ (t), x(t)i + hB(t)x(t), x(t)i]dt,

a

xa = xb = 0,

где A(t), B(t), C(t) квадратные матрицы размера n ×n. Коэффициенты A, B, C предполагаем непрерывными функциями от t.

Пусть xˆ(t) экстремаль, отвечающая его локальному сильному минимуму, то есть минимуму на множестве KC1 функций непрерывных, с кусочно непрерывными производными, могущими иметь конечное число разры-

вов первого рода, а малая окрестность экстремали выделяется ограничением в C-норме: t [a,b]

|

x(t)

|

 

 

 

 

 

 

 

 

 

 

 

 

max

xˆ(t)

 

< δ,

N

 

 

 

 

 

 

˙

 

 

 

 

 

 

 

 

 

{tj }j=1

множество точек разрыва первого рода производной xˆ(t). Сгладим эти разрывы, "скругляя"кривую

в малой ε-окрестности каждого разрыва подходящей дугой так, чтобы было ε C2,

 

 

 

 

 

 

6 ε, mes supp [ˆx(t)

 

 

˙

 

˙

 

(t) 6 K

 

 

 

 

 

k

(t)] 6 2N ε, sup xˆ(t)

 

 

 

 

 

 

εkC

 

ε

t |

 

 

ε

|

 

 

 

 

с некоторой постоянной K. Тогда

N

Z

J (ˆxε) = J (ˆx) + X

{[hAxˆ˙ ε, xˆ˙ ε − xˆ˙ i + 2hCxˆ˙ ε, xˆε − xˆi+

j=1|t−tj |<ε

˙˙ ˙ ˙

+hBxˆε , xˆε − xˆi] + [hA(xˆ − xˆε), xˆi + 2hC(xˆ − xˆε), xˆi + hB(ˆx − xˆε), xˆi]}dt.

5)Weirshtrass K.G.(Вейерштрасс Карл Гюстав, 1815-1897) - немецкий математик, один из создателей современной теории функций

Легко видеть, что 0 6 J (ˆxε ) − J (ˆx) 6 K· ε с некоторой постоянной K, или J (ˆx) 6 J (ˆxε) 6 J (ˆx) + Kε,

J (ˆx) = m сильный минимум функционала J (x) является и

inf J (x). Наш квадратичный функционал, очевидно,

 

 

 

x C2

 

непрерывен на C

1

 

ˆ

1

 

и его инфинум является минимумом (но достигается он на некоторой экстремали xˆ(t) KC

 

ˆ

6= xˆ, например, если

имеет разрывы производных). Значит, m сильный минимум

и поэтому возможно

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

Глава VI

Преобразование Лежандра

VI.1 Общее определение

Рассмотрим J : B → R функционал, определенный на множестве M B, непрерывно дифференцируемый (по Фреше). Oтображение u 7→DJ (u) называется градиентным, пусть оно является взаимно однозначным отображением M на N B , где B сопряженное к B пространство (то есть B это пространство линейных непрерывных функционалов на B).

Определение 9. Преобразование Лежандра

это отображение

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L :

 

{u, J } 7→ u{ , J },

 

 

 

 

 

 

 

 

 

 

 

 

 

 

заданное правилом u = DJ (u) B ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

J : u 7→(u , ((DJ )−1)(u )) − J ((DJ )−1(u ))

 

 

 

 

 

 

 

 

Здесь имеется в виду DJ : u 7→DJ (u) и

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(DJ )−1 :

v 7→ (DJ )−1(v)

v = DJ (u).

 

 

 

 

 

Значит,

J (u ) = (DJ (u), (DJ )−1(DJ (u))) − J ((DJ )−1(DJ (u)))|DJ(u)=u =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= hu , ui − J (u)|DJ(u)=u .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример 21.

Пусть f :

Rn

R1

. Для

f (x) =

|x|p

c

p (1, ∞) будет f

 

(y)

=

 

|y|p

, где

1 +

1

= 1. Действительно,

 

 

 

 

 

 

 

 

 

 

 

 

p

 

 

 

 

 

p

 

p

p

1

 

Df (x) = |x|p−2 · x ≡ y. Из этого соотношения выразим x через y: |x|p−1 = |y| |x| = |y|

p−1

, sgn x = sgn y.

Поэтому

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

1

 

 

p−2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x = |x| · sgn x = |y|

 

· sgn y = |y|

 

 

−1 · y = |y|p−1

· y ≡ ((Df )−1)(y)

 

 

 

 

p−1

p−1

 

 

 

 

 

 

 

 

 

p−2

 

1

 

 

p

p−2

 

1

 

 

 

p

 

1

 

 

 

 

 

 

 

f (y) = hy, |y|p−1 · yi −

 

 

· |y|

p−1

= |y|2− p−1

 

 

|y|

p−1

=

 

|y|p .

 

 

 

 

 

 

p

p

p

 

 

 

 

Теорема 31. Если J C2, то L2 = I (такое свойство называется инволютивностью оператора L). Доказательство. L : {x, J (x)} 7→p,{ J (p)} означает p = DJ (x), x = (DJ )−1(p) и

J (p) = (p, x) − J (x).

Так как J является дважды непрерывно дифференцируемым функционалом, то мы имеем для вычисления DJ (p) формулу

q ≡ DJ (p) = (DJ )−1(p) + Dp((DJ −1)(p))T · p − −Dp((DJ )−1(p))T · DJ ((DJ )−1(p)) = x.

Поэтому

J (p) = (DJ (p), p) − J (p) = DJ (p) = q = x = (DJ )−1(x) = p

= (x, p) − (p, x) + J (x) = [(x, p) ≡ (p, x)] = J (x).

50

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