
Variatsionnoe_ischislenia
.pdf
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) |
|
Пусть экстремаль xˆ дает локальный минимум. По теореме 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 должна удовлетворяет начальным условиям
= 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 - окрестности графика экстремали xˆ, будет
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 , λ).
В итоге получаем |
|
|
dτ x(τ, λ) − D1x(τ, λ) |
|
||||
D1S(τ, ξ) = L(τ, ξ, u(τ, ξ)) + D3L(τ, ξ, u(τ, ξ)) · |
= |
|||||||
|
|
|
|
|
|
d |
|
|
|
|
|
|
|
|
|
= L(τ, ξ, u(τ, xi)) − D3L(τ, ξ, u(τ.ξ)) · u(τ, ξ), (V.10) |
|
здесь |
d |
x(τ, λ) = |
d |
ξ = 0. |
|
|
|
|
|
|
|
|
|
|
|||
|
dτ |
dτ |
|
|
|
|
||
Аналогично вычисляем |
|
|
|
|
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) = |
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). Сгладим эти разрывы, "скругляя"кривую |
|||||||||||||||
в малой ε-окрестности каждого разрыва подходящей дугой так, чтобы было xˆε C2, |
|
|
|
|
||||||||||||
|
|
xˆ |
xˆ |
6 ε, mes supp [ˆx(t) |
|
xˆ |
|
˙ |
|
˙ |
|
(t) 6 K |
|
|
|
|
|
k |
− |
(t)] 6 2N ε, sup xˆ(t) |
− |
xˆ |
|
|
|
|
|||||||
|
|
− ε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ˆ, например, если xˆ |
имеет разрывы производных). Значит, m сильный минимум |
и поэтому возможно xˆ |
квадратичного функционала является и его слабым минимумом. В частности, экстремаль, доставляющая слабый минимум, дает сильный минимум.

Глава 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