Файл: Васильев Ф.П. Лекции по методам решения экстремальных задач.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)], |
|
|
|
|
|
|
|