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

Rozova_Maximova

.pdf
Скачиваний:
153
Добавлен:
15.04.2015
Размер:
714.37 Кб
Скачать

3. Найти допустимые экстремали, т.е. допустимые решения уравнения Эйлера для лагранжиана L(x) на векторе множите-

лей Лагранжа λ = (λ0 ,λ1 ,..., λm ) , не равном нулю.

При этом бывает полезно отдельно рассмотреть

случаи λ0=0 и λ00. Во втором случае можно положить λ0 равным единице или любой другой, отличной от нуля константе.

Пример. Найти экстремаль в следующей задаче:

1

 

1

 

 

2

dt extr, xdt = 0, x(0)

= 0, x(1) = 1.

x

 

0

 

0

 

Решение. Составляем лагранжиан

1

2

+ λ1 x)dt .

 

L(x) = (λ0 x

 

0

 

 

Выписываем уравнение Эйлера для лагранжиана

λ1 dtd λ0 2x = 0,

λ1 2λ0 x = 0.

Если λ0=0, то λ1=0, то все множители Лагранжа равны нулю, что противоречит правилу множителей Лагранжа.

Рассмотрим λ0 0. Пусть λ0=0,5, тогда

x = λ1 , x = λ1t + C1 ,

x = λ1 t 2 + C1t + C2 . 2

Подставляем краевые условия и получаем

C

2

= 0,C = 1

λ1

.

 

 

 

 

1

2

 

 

 

 

 

 

 

 

 

 

t 2

 

λ

x = λ

 

+ (1

 

1 )t.

2

 

 

 

1

 

2

 

41

Подставляем полученное решение в изопериметрическое условие

1

t 2

 

λ

(λ1

 

+ (1

1

)t)dt = 1.

2

2

0

 

 

Отсюда находим λ1=6, С1=2.

Тогда допустимая экстремаль имеет вид xˆ = 3t 2 2t.

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

Пусть задан функционал от функции нескольких переменных

F(x1

b

,..., xn , x1

,..., xn )dt.

 

,..., xn ) = f (t, x1

 

 

 

 

 

 

 

a

 

 

 

Найти экстремум функционала

 

 

 

b

 

 

(6.3)

F(x1,..., xn ) = f (t, x1,..., xn , x1,..., xn )dt extr

a

с краевыми условиямиx(a) = xa , x(b) = xb –заданные n-мерные векторы при наличии связей

gk(t,x1,…,xn)=0, k<n.

(6.4)

Для нахождения экстремума в данной задаче используется следующая теорема.

Теорема. Пусть функции x1(t),…,xn(t) реализуют экстремум функционала (6.3) при наличии связей (6.4), а функции f, gk дважды непрерывно дифференцируемы по совокупности переменных.

Тогда существует ненулевой вектор λ (t) = (λ1 (t),...,λm (t)) такой, что x1,…,xn – стационарные кри-

вые лагранжиана:

42

b

,..., xn , x1

m

L(x, λ ) = ( f (t, x1

,..., xn ) + λk (t)gk (t, x1 ,..., xn ))dt .

 

 

 

a

 

k =1

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

Найти экстремум функционала

b

F( y, z) = f (t, y(t), z(t), y(t), z(t))dt extr

a

с краевыми условиями

y(a) = ya , y(b) = yb , z(a) = za , z(b) = zb

при наличии связи

g(t,y(t),z(t))=0, y,z C1[a;b].

Теорема. Пусть xˆ = ( yˆ, zˆ) – решение задачи и функции f (t, y, z, y, z), g(t, y, z) – дважды непрерывно дифференцируемы по совокупности переменных в точке ( yˆ, zˆ) . Функционал F(y,z) дифференцируем по Фреше в точке xˆ = ( yˆ, zˆ) ,

т.е. существует F (xˆ)[h] = F ( yˆ, zˆ)[ξ ,η ].

 

 

Пусть

также grad g(t, y, z) 0,

t [a, b], т.е.

g

и

 

 

 

 

 

 

 

 

 

y

 

g

одновременно в ноль не обращаются, т.е.

на поверхно-

 

z

 

 

 

 

 

 

 

сти g(t,y(t),z(t))=0 нет особых точек.

λ (t) 0

 

 

 

 

 

Тогда

существует функция

такая,

что

xˆ = ( yˆ, zˆ) – стационарная кривая лагранжиана

b

L(x) = ( f + λ (t)g)dt , где x(t)=(y(t),z(t)).

a

Док-во. Выпишем приращение функционала и воспользуемся формулой Тейлора:

43

 

 

 

 

 

b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F( y, x) = [ f (t, y + ξ , z +η,

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

y, z, y, z)]dt =

 

 

 

 

+ ξ , z +η ) f (t,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

(t,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f (t,

 

 

 

 

 

 

 

 

 

= [

f

y, z, y, z)

ξ +

 

f (t, y, z, y, z)

ξ +

y, z, y, z)

η +

 

 

 

a

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f (t, y, z, y,

z)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

η ]dt + R( y, z,ξ ,η ).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Заметим, что

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

 

f

(t, y,

 

 

 

 

 

 

 

 

 

 

 

 

 

f (t, y, z,

y,

z)

ξdt =

z, y,

z)

dξ

=

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f (t, y, z,

 

 

 

 

 

b

 

d

f (t, y,

 

 

 

 

 

 

 

 

 

b

 

d

 

f (t, y, z,

 

 

 

 

 

 

 

 

 

 

 

 

 

=

y,

z)

ξ

ba

(

 

z, y, z)

)ξdt = −

(

 

y, z)

)ξdt.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

a

dt

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

a

 

dt

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f (t, y, z, y, z)

 

 

 

b

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Слагаемое

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ξ

 

a

 

равно нулю в силу закреплен-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ных концов x(a) = xa , x(b) = xb траектории.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

 

 

d

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Аналогично

 

f

(t, y, z, y, z)

ηdt = −

(

 

 

f (t, y, z, y, z)

)ηdt .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

dt

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таким

 

 

образом,

 

 

 

 

 

 

 

 

 

 

по

 

 

 

 

 

принципу

 

Ферма

 

 

 

 

 

 

 

 

 

 

b

 

f

 

 

 

 

 

d

 

 

f

 

 

 

f

 

 

 

d

 

f

 

 

 

 

 

 

 

 

 

 

 

F ( yˆ, zˆ)[ξ ,η ] = [(

 

 

 

 

 

 

)ξ +(

 

 

)η ]dt = 0 .

 

 

 

y

 

 

dt y

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

dt z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Однако допустимыми являются не всякие функции ξ (t) и η (t) , а только удовлетворяющие ограничению

g(t,y+ξ,z+η)=0, т.е. ξ (t) и η (t) связаны между собой.

Выявим эту связь:

 

 

 

 

 

g(t,y+ξ,z+η)=0 dg = g(t, y, z) ξ +

g(t, y, z)

η = 0 .

 

y

 

z

 

Пусть t* [a,b] такая точка, что

g

 

(t*, y, z) 0 .

 

z

 

 

 

 

 

44

Тогда, учитывая непрерывность всех производных и функций, мы можем сказать, что gz (t, y, z) 0 и в некоторой ок-

рестности Sδ(t*) точки t*.

Выразим явно η через ξ в этой окрестности:

g(t, y, z)

η (t) = ( ,y, ) ξ (t) , t Sδ(t*). g t y z

z

Вне Sδ(t*) положим η и ξ

равными нулю. Подставим полу-

ченное выражение в производную, получим:

g

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F ( yˆ, zˆ)[ξ ,η ] =

 

f

 

 

d f

f

 

d f

y

[(

 

 

 

 

 

) +(

 

 

 

 

)(

 

)]ξdt = 0 .

y

 

 

 

z

 

 

g

Sδ (t*)

 

 

dt y

 

dt z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Теперь ξ не имеет ограничений, и воспользовавшись леммой Лагранжа, получим:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

g

f

 

d f

(

f

 

d f

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

)

 

 

= 0.

y

 

 

 

 

 

 

z

 

 

 

 

 

g

 

 

 

dt y

 

 

 

 

dt z

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

g

 

f

d f

= (

f

 

d f

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

)

 

.

 

y

 

 

 

 

 

 

z

dt

 

 

 

g

 

 

 

dt y

 

 

 

 

z

 

 

 

z

Откуда по правилу пропорции получаем:

 

f

d f

 

 

f

d

 

f

 

 

 

 

 

 

 

 

 

 

 

 

y

dt y

 

 

 

not

 

 

 

 

 

 

 

 

 

 

 

 

 

=

z

 

dt z

 

=λ (t).

 

 

g

 

 

g

 

 

 

 

 

 

 

 

 

 

y

 

 

 

z

 

 

Таким образом, имеем систему

45

 

f (t, y, z, y, z)

d

 

f (t, y, z, y, z)

λ (t) g(t, y, z) = 0,

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

dt

 

y

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(6.5)

 

 

 

 

 

 

 

 

g(t, y, z)

 

 

f (t, y, z, y, z)

d f (t, y, z, y, z)

λ (t)

= 0,

 

 

z

 

dt

 

z

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

при t Sδ(t*).

Так как grad g(t, y, z) 0 на всем отрезке [a,b] , аналогичные рассуждения верны и в некоторой окрестности любой точки t [a,b] (эти окрестности в совокупности образуют покрытие отрезка [a,b] , из которого, в силу компактности,

можно выбрать конечное подпокрытие). Таким образом, для функционала

b

L(x) = ( f (t, y, z, y, z) + λ (t)g(t, y, z))dt

a

равенство L(xˆ)[ξ ,η ] = 0 h имеет место на всем промежутке[a,b] . Полученная система (6.5) является системой уравнений Эйлера для лагранжиана L(y,z). Теорема доказана.

6.5.Правило множителей Лагранжа

вобщем случае

Перед тем как рассмотреть задачу, напомним некоторые факты из функционального анализа.

F(x): XY – оператор, X, Y – линейные нормированные пространства.

Опр. 1. Пусть U – окрестность точки xˆ X. Говорят, что

оператор F: XY дифференцируем по Фреше в точке xˆ , если существует такой линейный непрерывный оператор L:

46

XY, что для всех h X, для которых xˆ +h U, справедливо равенство

F( xˆ +h) F( xˆ )=Lh+r(h)

с оценкой остатка

ε>0 δ>0: h X <δ r(h) Y < ε h X .

Оператор L называется производной по Фреше оператора F в точкеxˆ и обозначаетсяFФ(xˆ)[h] . Это можно записать следующим образом:

F(xˆ + h) F(xˆ) L(xˆ)[h] Y < ε h при h X < δ (ε , xˆ) .

Опр. 2. Оператор F(x) строго дифференцируем по Фреше в точке xˆ , если существует такой линейный непрерывный оператор L( xˆ ), что ε>0 δ(ε): x1 , x2 Sδ (xˆ) вы-

полняется условие

F(x1 ) F(x2 ) L(xˆ)[x1 x2 ] Y < ε x1 x2 X .

Опр. 3. Пусть M – некоторое линейное подпространство в линейном нормированном пространстве X. Аннулятором подпространства M в пространстве X называется множество линейных непрерывных функционалов, обладающих

следующим свойством: Λ M : ΛM = 0.

Замечание. M { }, так как нулевой оператор всегда принадлежит аннулятору.

Замечание. Если X,Y – линейные нормированные

пространства, а Z=X×Y также линейное нормированное пространство, то норма Z – это сумма норм

z Z = x X + y Y , z = (x, y) Z .

47

Лемма об аннуляторе. Если M – замкнутое собствен-

ное (MX) подпространство в X, то его аннуляторM всегда содержит хотя бы один ненулевой элемент.

Лемма. Пусть X,Y,Z – линейные нормированные про-

странства. Если A(x): XY, B(x): XZ, C(x)=(A(x),B(x)):

XY×Z, ImA(Ker B) – замкнутое подпространство в Y, а ImB

– замкнутое подпространство в Z, то ImC – замкнутое подпространство в Y×Z.

Опр. 4. Норма линейного функционала F: XR вводится как F = sup F(x) .

x≤1

Теперь вернемся к задаче Лагранжа. Пусть f(x): XR,

F(x): XY. Задача состоит в том, чтобы найти экстремум функционала f(x) при наличии ограничений F(x)=0.

Сформулируем необходимое условие экстремума данной задачи.

Пусть xˆ – решение поставленной задачи.

Теорема. Пусть X,Y – линейные нормированные Банаховы пространства, U( xˆ ) – открытое множество в X. Отображения f(x), F(x) строго дифференцируемы по Фреше в

точке xˆ . f( xˆ ): U( xˆ ) R , F( xˆ ): U( xˆ ) Y. Пусть

Im F (xˆ) –

замкнутое подпространство в Y.

 

Тогда существуют λ0 R и y* Y*

( λ0 , y * –

множители Лагранжа) одновременно не равные нулю, такие, что для функционала

L(x) = λ0 f (x) + y*, F(x)

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

Lx (xˆ, λ0 , y*, h) = λ0 f (xˆ)[h] + y*, F (xˆ)[h] = 0 , h X,

48

т.е. экстремальxˆ является стационарной точкой лагранжиана.

Замечание. Рассмотрим отображение L(x): XR. Найдем производную L(x).

L= L(x + h) L(x) =

=λ0 f (x + h) + y*, F(x + h)λ0 f (x) y*, F(x) =

=λ0 [ f (x + h) f (x)] + y*, F(x + h)y*, F(x) =

=λ0 [ f (x + h) f (x)] + y*, F(x + h) F(x) =

 

 

= λ0 f + y*, F = λ0 f (xˆ)[h] + y*, F (x)[h] + R(x, h) =

 

 

 

 

 

= L(x)[h] + R(x, h).

 

 

 

 

 

 

x

 

 

~

Док-во теоремы.

Рассмотрим следующее отображе-

ние

(x) = ( f

(x) f (xˆ), F

(x)),

~

F

F(x) : X R×Y=Z.

~

 

 

 

 

 

 

 

F(xˆ) = (0,0) .

 

 

~

 

 

 

 

Очевидно,

что

(x)

строго дифференцируемо по

 

 

F

Фреше в точке xˆ

~

(xˆ) = ( f (xˆ), F (xˆ)) – линейное непре-

и F

рывное отображение.

 

 

 

 

 

Возможны~два случая:

 

 

 

I.

Im F (xˆ) Z.

 

 

 

 

II.

~

(xˆ) =Z.

 

 

 

Im F

 

Рассмотрим каждый отдельно.

Случай I.

По условию, Im F (xˆ) – замкнутое линейное подпространство в Y, F (xˆ) – линейное непрерывное отображение, а следовательно, Ker F (xˆ) ={h X: F (xˆ) [h]=0} линейное подпространство в X. Тогда f (xˆ) [Ker F (xˆ) ] – линейное под-

пространство в R (в силу линейности и непрерывности производной f (xˆ) ).

49

~

(xˆ)

– замк-

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

В R всего два линейных подпространства

{0}

и R, и

оба они замкнуты. Таким образом, f (xˆ) [Ker F (xˆ) ] замкну-

то. Im F (xˆ) – замкнутое линейное подпространство в Y, f (xˆ) [Ker F (xˆ) ]– замкнутое линейное подпространство в R,

нутое подпространство в Z. По лемме об аннуляторе

~

~

(xˆ) )=0.

Λ [ImF(x)]ˆ

:ZR, Λ≠0: Λ(Im F

Далее, поскольку линейный непрерывный функционал на прямом произведении линейных нормированных пространств представим в виде суммы линейных непрерывных функционалов на соответствующих пространствах, имеем

Λ(x, y) = x*, x + y*, y .

~

(xˆ) )=0,

следует, что

Тогда из того, что Λ(Im F

Λ( f (xˆ)[h], F (xˆ)[h]) = λ0 f (xˆ)[h] + y*, F (xˆ)[h] = 0 , h X.

Так как Λ≠0, то λ0 и y* одновременно в ноль не обращаются.

L

(xˆ)[h] = λ

0

f (xˆ)[h] + y*, F (xˆ)[h] = 0, h X:

 

 

x

 

 

~

 

L

(xˆ)[h] : XR, следовательно, в случае, когда

(xˆ) Z,

Im F

x

утверждение теоремы доказано.

Случай II. ~( ˆ)

Im F x =Z.

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

Теорема о неявной функции. Пусть X,Z – Банаховы

пространства,

U( xˆ )

окрестность

в

X

и

~

 

. Пусть

~

 

 

 

F(xˆ) :U (xˆ) R × Y = Z

F(x) строго дифференци-

руемо в точке

~

(xˆ) =Z. Тогда k>0 и Sρ(z0) Z, где

xˆ , Im F

~

 

 

 

~

(ϕ (z)) = z

z0= F(xˆ) : z Sρ(z0), найдется ϕ(z)=x, такая, что

F

и выполнена оценка:

50