Файл: Васильев Ф.П. Лекции по методам решения экстремальных задач.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 11.04.2024

Просмотров: 209

Скачиваний: 1

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

296 МЕТОДЫ

М И Н И М И З А Ц И И В Ф У НКЦ И О НАЛ ЬНЫ Х ПРОСТРАНСТВАХ [Гл. 6

 

Т

Т

 

 

 

Т

 

 

 

х{1, t)u(t)dt<Cv

 

t)di +

^ и2 (t)dt.

 

 

О

0

 

 

 

0

 

 

Подставляя это неравенство в

(8), сразу

придем

к

требуемой

оценке (6).

 

является частным

случаем задачи

(5.2). Что­

Задача (1) — (5)

бы убедиться в этом, достаточно в

(5.1)

положить

 

 

B =

L 2[0, Т], Х = Ь 2[0,

l],x{u ) = x (s, Т, и), х (0) =

 

=

a-(s, Т, 0) =

0,

Ku = x(s, Т, и),

 

 

причём оператор К,

действующий

из L2[0,

Т] в L2[0,

/], очевидно,

•линеен, а его ограниченность равносильна оценке (6) с ||/(||<

Из результатов § 5 тогда вытекает, что функционал (1) при усло­ виях (2) — (5) является выпуклым и дважды непрерывно диффе­ ренцируемым на U, причем градиент его удовлетворяет условию Липшица с константой L=2||/(||2^ v . Градиент этого . функцио­ нала в точке u = u ( t ) ^ U может быть представлен в виде

 

=

 

 

 

0 < г < 7 \ *

 

 

(9)

где ф(5, ^ ) = tJ5(s ,

t, и) является

решением

следующей

краевой

задачи:

 

 

 

 

 

 

 

 

 

 

 

дт|>

 

d2ip

 

(S, t) е Q,

 

 

 

( 10)

 

dt

 

ds2

 

 

 

dV(o,t)

=0| J W

J L

=

— Vip(itt)t о < t < T ,

( П )

ds

ds

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

il)(s, T) = c (s), 0 < s < / ,

 

 

 

( 12)

где c ( s ) = 2 [ x (5,

T, u) — y (s )],

O ^ s ^ l

Можно доказать

сущест­

вование и единственность

решения

ф(5, t)==ty(s,

t, и)

задачи

(10) — (12), которое удовлетворяет условиям

(10),

(12) в

класси­

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

(11)

— в

слабом

смысле,

причем

имеет смысл говорить о значениях ф(/,

t ) ^ L 2[0,

Г ], являющихся

пределом в Ь2[0, Т] значений ф(«, t):

 

 

 

 

 

 

 

s-lim*l—0 I|tJ)(s ,

 

t) |2 dt

 

 

 

 

Выведем формулу (9)

для градиента. Пусть u(t),

v(t) = u {t) +

+ h ( t ) ^ U , и пусть x{s, t,

u )= x (s , t), x ( s ,t ,v ) = x ( s ,t )

+Ax(s, t) —


§ 6]

Оптимальное управление процессом нагрева

стержня

297

соответствующие этим управлениям решения

краевой

задачи

(2) — (4). Приращение функционала (1) тогда запишется в виде.

i i

J (и + h) J (и) = 2 ^ [х (s, Т) у (s)]Ах (s, T ) d s + $ \ b x (s, Т) |2 ds.

 

 

 

 

о

 

 

 

 

 

о

(13)

 

 

 

 

 

 

 

 

 

 

 

 

Заметим, что функция Ax(s,

t) является решением краевой задачи

(2) — (4)

при u ( t ) = h ( t ) , и,

следовательно,

справедлива

оценка

 

 

 

 

i

 

 

 

т

 

 

 

 

 

 

 

| | Д * (8 ,Г )| « * < - ^ | А '(*)| * Я ,

 

(14)

 

 

 

 

о

 

 

 

о

 

 

 

 

вытекающая из

(6)

 

 

 

 

 

 

 

 

Учитывая это обстоятельство и условия

(10) — (12),

имеем

 

 

 

 

 

 

 

 

 

i

т

 

 

j с (s) Ax (s, Т) ds = J

ф(s, Т) Ax (s, Т) ds = j

( j

(фДх) dt^j ds =

= Я ( ‘f ^ + 16 ^

 

 

= Я h ■ S - + ^

“ =

 

=

\

( г 1 7 А*

+

1т г ) С

'

=

[l' t ) M '(/’ t] +

 

 

0

4

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

T

 

 

 

 

 

•+ ф (/, t)v[h(t) — &x(l, t)]}d t= v ^ф(/, t)h(t)dt.

(15)

 

 

 

 

 

 

 

 

 

0

 

 

 

Так

как

c (s )= 2 [x (s ,

T, u) — y (s)],

то

формула

(13)

с . учетом

(14),

(15)

запишется в виде

 

 

 

 

 

 

 

 

J ( u + ti) j

(«) = J

Vij?(/, t)h{t)dt +

0 ( |

| Л .

 

 

 

 

 

 

 

о

 

 

 

 

 

 

Сравнивая

полученное

выражение для приращения функционал.а.

с (1.1),

находим, что

градиент J'(u)

выражается

формулой (9),

а дифференциал

(/'(u),

К) имеет вид

 

 

 

 

 

 

 

 

т

 

 

 

 

I

 

 

 

( J 1(«), h)Lt[0,T] =

v

 

и) h(t)dt =

2 j

[х (s, Т, и) — у (s)] Ах(s, Т) ds.

 

 

 

 

о

 

 

 

о

 

 

 

 

( 1 6 )


I

298 МЕТОДЫ М И Н И М И З А Ц И И В ФУНКЦ И ОНАЛ ЬНЫ Х ПРОСТРАНСТВАХ [Гл. б

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

условиях (2) — (5)

требуется

последовательно

решить

краевые

задачи

(2) — (4) и

(10) — (12).

Для решения этих краевых задач

можно

использовать различные приближенные

методы,

в част­

ности расчеты можно вести по неявной разностной схеме [20, 207,

227];

интегралы (1),

(16) вычисляют с помощью квадратурных

формул [19].

 

 

 

 

Заметим, что приведенное выше доказательство равенств

(7),

(15)

носит формальный характер, поскольку законность всех

про­

деланных при выводе

(7), (15) операций осталась необоснованной;

строгий вывод (7), (15)

см. в работе [35].

 

 

В

силу теоремы

1.4

оптимальное решение u = u *(t)

задачи

(1) — (5) существует,

а согласно теореме 2.1.3 и формуле

(16)

для

оптимальности u*(t)

необходимо и достаточно, чтобы

 

 

 

т

 

 

 

 

J ф (I, t, и") [и (t) — и* (t)] dt > 0,

j (x(s, Т, u*) — y{s)) (.x(s, T, и) — x(s, Т, и*)) ds > 0

о

при всех u = u (t)^ U .

Опираясь на результаты § 5, опишем метод условного гради­

ента для приближенного

решения задачи

(1 )— (5). Пусть

un(t),

О ^ ^ Т ( п ^ О ) известно.

Определим

un(t)

из условия

 

 

f Ф(I,

“«)

(0 dt = min

Г ф(I, t, ип) и (t) dt,

 

 

где ф(в, t, ип) — решение задачи

(10) — (12) при и = и п.

В

рас­

сматриваемом случае

un{t) легко выписывается в явном виде

 

 

 

1,

если ф(/, t,

ип) < 0 ,

 

 

“л Н

- >

;

если ф(£, t,

ил) > 0 , ' 0 ;< < < Т .

 

 

Далее, г полагаем

ип+х (t)

=

ип (t) + а„ [ип(t) ип(if)]' 0 < t <

Т,

где

«л = m in {l; а * } > 0 ,

 

 

 

 

 

 

 

 

 

 

г

ф)(/, t, a)[un {t) — un {t))dt

 

 

 

 

v J

 

 

 

 

2 J

|дс (s, Г , un) — x ( s ,T ,

un) |a ds

 

 

о


§ 6)

Оптимальное

управление процессом нагрева

стержня

299

 

f [X (s, Т ,

и,) у (s)] [л: (s, Т, ип) — x (s, Т ,

ип)] ds

 

------------------------------------------------------- --------------------------

\i

[ |х (s, Т , ип) — x (s . Г , ип) |2ds

о

(если а„ = 0 или x(s, Т, un) = x(s, Т, ип),

то un(t) = u*(t), /0 < ^ < 7 \

— оптимальное решение задачи (1) — (5),

и итерации на этом

закан­

чиваются). Сходимость этого метода следует из теоремы 2.3.

 

Кратко остановимся также на одном из вариантов

метода

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

станты

Липшица L

для

градиента

1'{и)

здесь

можно взять

L —v. Тогда

метод

проекции

градиента для

задачи

(1) — (5)

будет заключаться в построении последовательности

{un.(t)}czU

по правилу

 

 

 

 

 

 

^

 

 

 

 

 

 

ип+ 1(t) =

Рц («„ (0 — a„vi|>(l,t,

ип)) =

 

 

 

' ип(t) — a„vij) (l,t,

ип),

если |ип (t) — «„vife (/, t,

ип) |< 1,

=

1,

 

 

 

если

ип (t) — a„vaj)](/, t, ип) >

1,

 

 

— 1,

 

 

 

если

ип(t) — a„vtj3 (/,

 

 

 

где a„

выбирается из условия

ех -< а„ <

v +

•;

elt

е — параметры

 

 

 

 

 

 

 

 

2e

 

 

8> 0.

алгоритма, выбираемые

вычислителем,

0 < е х < ;-

 

 

 

 

 

 

 

 

 

 

 

 

v + 2 e

В частности,

возможно е =

-,

е1 =

а п =

 

 

 

 

Более общие задачи оптимального управления для параболи­ ческих уравнений и формулы для градиентов см. в работе [35].

Упражнения. 1. Пусть дан функционал

гI

J(u ) = $ J|u(f/|*tf + 8 j > ( s , 7 )-i/ (s )| 2ds

оо

при условиях (2) — (4), где

р,

6 = const^sO,

y(s)

заданная не­

прерывная функция. Требуется:

 

 

 

 

 

 

а)

доказать,

что при р > 0

функционал J

(и)

сильно выпуклый

в L2[0,

Г ];

 

 

 

 

 

 

 

 

ет

б)

убедиться,

что градиент этого функционала в L2[0,

Г] име­

вид

 

 

 

 

 

 

 

 

 

 

 

 

 

J' («) = 2ры(t) +

П и ) ,

£0 <

* <

7\

 

 

где

ф(5,

t, и)

решение

задачи

(10) — (12)

при

c(s) =

= 2 6 [x (s,

Т, u)—y (s)],