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

книги из ГПНТБ / Васильев Ф.П. Лекции по методам решения экстремальных задач

.pdf
Скачиваний:
15
Добавлен:
25.10.2023
Размер:
14.17 Mб
Скачать

200 ДИНАМ И ЧЕСКО Е ПРОГРАММИРОВАНИЕ. ПРОБЛЕМА СИНТЕЗА [Гл. А

принадлежит Д(;<c,t),

то

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(+Д*

 

 

 

 

 

 

 

 

 

В (х, t) <

J* (х,

и (т)) = .

J(х (х,

и*), ы* (т), т) dx +

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

+

B(x(t + At,u*),

t + At) =

t + A t

 

 

 

 

 

 

 

 

J /°(х(т,

u*),

w*(t),

+

 

 

 

 

 

 

 

 

*+Д<

 

 

 

 

(t), t) dx +

+

У*+А‘ (X (f +

Д£, и’), и**) <

j

/ °(*(T, и*). «*

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

+

Ji+At (x(t +

Ы, U * ) , U*) =

Jt ( X , «*) =

В (X ,

t).

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t+ A t

 

 

 

 

 

 

 

 

 

 

В (x, t) =

Jt (x,

 

u*) =

j

f° (x (т,

и*),

и* (t), t) dx +

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

~\- В {x(t -(- At, u*), t -j- At).

 

 

(8)

Далее,

 

возьмем

произвольное

w e V

и

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

w (-t)eA (x,£)

согласно предположению

(7).

Пусть этому управле­

нию соответствует траектория х(х,и)

из

(5).

При сделанных пред­

положениях

существует

управление

в * (т ) е А (х(£+Д£, и) Д +А О

такое,

что

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Д+д* (х (t +

At, и), и* (т)) — B(x(t -{- At, ti),

t - f

At).

Тогда управление

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

й ( т ) = { “ .

 

« ! «

 

+ * ,

 

 

 

 

 

 

 

 

 

 

1 о*(т), t + Д * < т < 7 \

 

 

 

принадлежит A(x,t),

поэтому

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t + A t

 

 

 

 

 

 

 

 

 

B(x,t) <CJ* (х, и(х))

=

J

 

f° (x(x, u),u,

x) d x B { x ( t

+ At,u), t + At).

Отсюда и из

(8) следует

 

 

 

 

 

 

 

 

 

 

 

 

t + A t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j

/° (х(т, и*), и*(х), x)dx-\-B{x{t +

At, u'),t +

At) ■

t + A t

— В (x, t) = 0 < j

f°(x(x,u),u, x)'dx]+

В (x(t + At, u),

t + At)— В (x, t).

§ 3]

 

 

 

Дифференциальное уравнение Р. Веллмана

 

 

201

Разделим

это

неравенство

на

>

0

и

перейдем к

пределу

при

At -*• +

0.

Получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(х, и* (t + 0), 0 +

d- -& £ +-0, u*h

 

= 0 <

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

< f° (X,

U,

t)

+

(*■(* + о, «0. О t и 6 V,

и < / <

Т.

 

(9)

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

производная

dB(x(t

 

0, и *), t)

с

учетом

уравнения (5) пред-

Полная

— v

 

ставляется

в виде

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d_B ( x ( t +

0 , u * ) , t )

=

 

^

0> /(Х)

 

+

0)> t ) ) + B (

{Xt

t)

 

 

 

at

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Аналогично

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d B ( x ( t + 0, u ),t)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

= (В х ( х , t ) , f { x , и , f ) ) + B t ( X , t ) .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Из (9) тогда имеем

 

 

 

 

 

 

 

 

 

 

 

 

 

f° (х,

и*(t +

0), t) +

(Bx (x,

t),

f (x,

и* {t +

0), t) + B t (x,

f) =

0 <

 

 

 

<

f°(x,

u, t) +

(Bx ix, t),

f ix, u,

t)) +

B( (x, t)

 

 

 

при всех

и 6 V.

Так

как

V замкнуто и и'

-f- 0) 6 V, то последнее

неравенство может быть переписано в виде

 

 

 

 

 

 

min [Bt ix,

t) +

(Bx ix,

t), f ix,

u, t)) +

(x, u, *)] = 0, t0 <

/ <

T.

(10)

u€V

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Полученное

соотношение

называется

уравнением

Веллмана.

К этому уравнению добавим начальное условие

 

 

 

 

 

 

 

 

 

В ix,

Т) =

Ф (х).

 

 

 

 

 

(11)

Таким образом, для определения функции Веллмана имеем задачу (10), (11), которую естественно назвать задачей К.оиш Веллмана [14— 18, 24, 27, 34, 195, 206, 233] и др.). Заметим, что присутствие знака min в левой части (10) весьма осложняет изу­ чениетаких уравнений, и на-сегодняшний день задача Коши — Веллмана исследована недостаточно. В тех случаях, когда удается найти решение задачи (10), (11), то нетрудно получить оптималь­ ные решения задачи (1) — (3) и (4) — (6), а также указать синте­ зирующую функцию для задачи (1) — (3). Об этом речь пойдет в следующем параграфе.

Упражнения. 1. Показать, что к задаче минимизации функции

 

N—1

J i u o> u i> ••• у u n ) —

J i i U i , Щ + 1) + J n iu ij)

(=0

202

ДИНАМ И ЧЕСКО Е ПРОГРАММИРОВАНИЕ. ПРОБЛЕМА СИНТЕЗА

[Гл. 4

векторных переменных [«„, ult ■■■, un при щ 6 У,- (i — 0, 1, . . . , N) можно применить метод динамического программирования.

У к а з а н и е .

Ввести функцию

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N—1

 

 

 

 

 

 

 

 

 

 

Bk (u) =

min [ J ]

/, (uh u;+i)]+ J N jv)],

 

 

 

 

 

 

 

 

 

 

 

i=k

 

 

 

 

 

 

 

 

где минимум берется по всем наборам

(uk — и, uk+u . . .

 

, uN), щ 6 V;

(i = k,

k + 1, . . .

,N),

и показать,

что

 

 

 

 

 

 

 

 

Bk (u )= min

[/*(«, о ) +

Дж

(и)],

(k =

.0, l , ..., N — 1),

 

 

 

 

vSVk n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

BN (и) = J N (и).

 

 

 

 

 

 

2. Написать и исследовать

уравнения

(1.13),

(2.5),

(3.10)

для

 

 

 

 

 

 

 

 

 

 

 

Г

 

 

 

 

 

задачи:

минимизировать

функционал

/ (и) = \x2(t)dt

при усло-

 

 

 

 

 

 

 

 

 

 

 

о

 

 

 

 

 

виях

х — А (t)x-{-K(t)u{t)+f(,t),

0s^.tz^T,

х(0) = х 0,

где

х =

{х\

..., хп), A(t), K(i),

f ( t ) — матрицы с

известными

непрерывными

элементами, х0, Т заданы;

управление u(t) — (u](t),

..., ur(t))

ку­

сочно-непрерывно

и

удовлетворяет

ограничениям

 

| ы '(/)[^1,

 

 

 

 

Г

 

 

 

 

 

 

 

 

 

 

 

 

£=1, 2,

г,-или

£

I

 

|2 < 1 .

 

 

 

 

 

 

 

 

3.

Решить

 

£=1

 

 

задачу,

заменив

функционал

на

предыдущую

J(u) =

 

(c,x(T)),

где

c = { c i ,

...,

сп) — заданный

вектор,

или

на

J ( u ) = x 2(T).

 

 

 

 

 

 

J(u)z=G>(x(T)),

'(f°(x, и, t) = 0 ) .

4.

Пусть в задаче

(1.1—4)

Показать, что схема Веллмана позволяет решить проблему синтеза для этой задачи. Провести аналогичные исследования для схемы Моисеева.

5. Найти функцию Веллмана В (х, t) для задачи: минимизиро­ вать функционал

т

J{u) = $ u 2(t)dt + x2(T)

о

при условиях x ( t ) = u , д:(0)=л:о, ( п = г = 1 ) .

6. Вывести дифференциальное уравнение Веллмана для зада­ чи быстродействия: минимизировать J ( u ) = T 10 (t0 задано) при

условиях x = f ( x , и d), t o ^ t ^ /Г, x(t0) = x о, x ( T ) = x h u(t)^.V, t^ .t0.

7. Показать, что функция Веллмана B(x;t) для задачи из при­ мера 3.2.2. не является непрерывно дифференцируемой.

§ 4]

Проблема синтеза. Оценка погрешности

203

§ 4. ПРОБЛЕМА СИНТЕЗА ДЛЯ СИСТЕМ С НЕПРЕРЫВНЫМ ВРЕМЕНЕМ. ОЦЕНКА ПОГРЕШНОСТИ

Проблема синтеза для задачи (1.1—4) с непрерывным време­ нем заключается в построении функций u— u(x,t), представляю­ щей собой значение оптимального управления при условии, что в момент t объект (1.2) находится в точке х фазового пространства. Как отмечалось в § 1, такая функция u(x,t) называется синтези­ рующей, и с ее помощью дальнейшее оптимальное движение объ­

екта определяется

условиями

х(х) = f ( x ( x ) , и(х, т), т),

 

x ( t ) = x .

 

 

 

 

 

 

Решение проблемы синтеза для

задачи

(1.1—4)

экви­

валентно решению

следующих

задач:

определить управление

и(х) = и * { х , t, т), доставляющее функционалу

 

 

 

т

 

 

 

 

 

J* (х, и (х)) = | /° (т), и (т), х) dx + Ф (х (Т))

(1)

 

t

 

 

 

 

 

минимальное значение при условиях

 

 

 

х (х) — f (х (т), и(х),

т),

t < х <

Т\ х (0

= х,

(2)

 

x(x)EG(x),

/ < т < Г ,

 

(3)

и = и (т) 6 V (т), t < т < Г ; управление и (х) кусочно-непрерывно, (4)

для каждого фиксированного x^G(t) и для каждого фиксирован­ ного момента времени /, t0^ t c T . Зная функцию и*(х, t,x), синте­ зирующую функцию можно получить сразу: u ( x ,t ) = u * ( x ,t ,t ) для

всех x^G(t) и всех

Т],

для которых задача (1)— (4) имеет

решение.

 

 

При решении проблемы

синтеза для дискретных систем важ ­

ную роль сыграло уравнение Веллмана (1.13). По аналогии попро­ буем использовать дифференциальное уравнение Веллмана для ре­ шения проблемы синтеза для задачи (1.1—4) с непрерывным вре­ менем. Для этого отвлечемся от тех ограничительных условий, при

которых выводилось

уравнение

Веллмана

(3.10), и рассмотрим

следующую задачу Коши — Веллмана:

 

 

 

«екщ

^

^

+

В{

*)] =

°»

(5)

 

x£G(t),

t0 < t < T ,

 

 

 

 

В ( х , Т ) =

Ф(х), xEG (T)

 

 

(6)

для определения функции В (х, t) . Если ввести функцию

 

 

R (х, и, t) г (Вх (х, 0 ,

f (х, и, t)) + Bt {х,

t ) + f ° (х, и,

t),

(7)

204

ДИНАМ И ЧЕСКО Е

ПРОГРАММИРОВАНИЕ. ПРОБЛЕМА СИНТЕЗА

\Гл. 4

то задача (5), (6) перепишется в виде

 

 

 

inf

R(x,u,t) =

0,

xeG {t),

В (х ,Т )

= Ф (х), x £ G {T ) . (8)

uev(0

 

 

 

 

 

 

 

Оказывается, решение задачи

Коши — Веллмана

(5),

(6), (или

(8))

равносильно решению проблемы синтеза для задачи

(1.1— 4).

Точнее, справедлива

 

 

 

 

 

 

Т е о р е м а

1.

Пусть функция B(x,t)

кусочно-непрерывна,

кусочно-гладка

и удовлетворяет условиям (5), (6),

(или

(8)) при

x ^ G ( t ) ,

 

Пусть для

каждой пары (х (т ),ы (т )),

удовлет­

воряющей условиям

(2)— (4),

функция Л (х (т ),т )

переменной т

непрерывна и кусочно-гладка на отрезке

 

при всех x e G ( t )

и t, t ( ^ t < . T . Пусть u = u ( x ,t )

кусочно-непрерывная функция, на

которой реализуется нижняя грань в (б)

или (8).

Тогда

функция.

u — u(x,t)

представляет собой решение проблемы

синтеза

для

за ­

дачи

(1.1— 4).

 

 

 

Возьмем некоторый момент t, t<ys^.t<T,

Д о к а з а т е л ь с т в о .

точку

x e G (f)

и через

D£t, Т\ обозначим

множество всех

пар

(х(т), ы (т)),

 

 

удовлетворяющих условиям

(2) — (4). Р ас­

смотрим задачу Коши:

 

 

 

 

 

 

 

 

 

 

X(т) =

/ (т), и (х (т), т), т),

t < т

<

Т; х (t) = х.

 

(9)

Пусть эта

задача

имеет

решение

х *(т ),

 

определенное

на

от­

резке

 

 

 

и

пусть

л* ( т) ё

6 ( т),

 

 

Положим

« * ( т )= ы ( х * ( т ),г ),

 

 

 

 

Очевидно,

(х *(т ),

и * (т )) eZXJf, Т].

Покажем,

что пара (х* (т),

и* (т )) является

оптимальным

решени­

ем задачи

(1) — (4)-. Для этого прежде всего выведем формулу

 

 

 

 

 

 

 

т

 

 

 

 

 

 

 

 

 

 

J 1(х, и(т)) =

J R (x( т), и (т), r)dx-\- В{х,

0 ,

 

(Ю)

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

справедливую при всех

(х (т), ы (т))еД Д ^, Т] и при всех

фиксиро­

ванных x ^ G ( t ) ,

t o ^ t c T - , здесь R(x,u,t)

взята из

(7).

 

 

Всюду

на отрезке

 

 

за исключением, быть может, ко­

нечного числа точек,

имеем

 

 

 

 

 

 

 

 

 

 

dB - . (T),T)

= R ( x (т), и (т), т) — /°( х (т), и (т), т).

 

 

 

 

dr

 

 

 

 

 

 

 

 

 

 

 

Так как по условию В ( х (г ), г)

непрерывна на \t, Т], то, интегрируя

последнее тождество с учетом условия (6), получим

 

 

 

 

 

 

 

т

 

 

 

 

 

т

 

 

 

Ф ( х ( Т )) — В(х,

t) = f

R(x(r), и (г), г)dr—f /°(х(т), и (г), г) dr,

j

<

§ 4] Проблема синтеза. Оценка погрешности 205

что равносильно (10). Далее, из условия (8) имеем

 

R (х,

и, т) >

inf

R (х , и,х) =

R (х, и (х, т),

т) = 0,

 

 

 

u£V{т)

 

 

 

 

где и е У (т ), x g G(t),

Отсюда следует, что для любой

пары (х (х) , и (г )) eZ)JY, 7]

справедливо неравенство

 

Я (х (т), и(т), т) >

0 е= # ( х*(т), и* (г), х), * < т < 7 \

(11)

А тогда с помощью

формулы (10)

и условия

х* (t) = x ( t )

= х

сразу получим

 

 

 

 

 

 

 

 

т

 

 

Jt{x ,u (т)) — Д (х, и* (т)) = J

[R (х (т), и (т), т) —

 

 

 

 

t

 

 

 

 

 

— R (х* (т), и* (т), т)] dx > 0

 

 

при всех х (т ),

и(х) ) g O JI, Т\. ±

 

 

 

Кстати, из

(10),

(11)

следует, что

 

 

 

Д (х, и* (т)) = inf Д (х, и (т)) == В (х, t).

DXV.T1

Как видим, в рассматриваемом случае В (х, t) есть функция Вел­ лмана задачи (1.1— 4). Теперь нетрудно .получить решение задачи

(1 .1 - 4 ) .

Т е о р е м а 2. Пусть

функции B (x ,t), u ( x ,t )

при x£ G (t),

t0 Kt<^.T, удовлетворяют

всем условиям теоремы 1, и

точка х* оп­

ределена из условий

в

( * ; , g = r i n f в ( * . *о). К £ х о>

х£Хо

где Х 0 = { х : х £ G (t0), Dx [t0, Т ] ф 0 } .

 

Тогда траектория х* (т) € G (х),

/0 < т

определяемая

из

(9) при

t — t0,

х =

х*

и

управление

ц* (т) =

и (х* (т), т), t0 < т; <

Т,

представляют

собой

оптимальное [ре­

шение задачи (1 .1 — 4).

 

 

 

 

 

 

 

 

Д о к а з а т е л ь с т в о .

 

Возьмем

произвольные

х е ! 0 и

(х(т),

u{x).)<=Dyt[tQ, Г ]. Неравенство

(11) при t = t 0, очевидно, так­

же справедливо. .Поэтому с учетом

определения

хо

из формулы

(10) имеем

 

 

 

 

 

 

 

 

. J (и) J («’) ЕЕ*До (х, и (т)) -

До (х;,

и* (т)) -

| [R (X (х), и (т), т) -

 

 

 

 

 

 

io

 

 

 

- R (х* (т), и* (т), т)] dx + в (X, g - В (х‘ , g > 0

при всех допустимых ы(т) из (1.2— 4). А

206

ДИНАМ И ЧЕСКО Е ПРОГРАММИРОВАНИЕ. ПРОБЛЕМА СИНТЕЗА

[ Г л , ./

2.Заметим, что теоремы 1—2 не работают в тех случ

когда нижняя грань в (5) не достигается и приходится иметь дело

с функциями,

 

приближенно реализующими эту нижнюю

грань.

Аналогично

при решении

задачи с подвижным

левым

концом

inf В (х, t0)

 

также может не достигаться. В таких случаях могут

*ех0

 

 

 

 

 

 

 

 

 

 

 

 

оказаться полезными следующие две теоремы.

 

 

 

 

Т е о р е м а

3. Пусть функция В (х, t) удовлетворяет

всем ус­

ловиям теоремы 1, и пусть

имеется

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

функций

uk (x,x)

( k = l ,

2, ...), таких,

что:

 

 

 

 

 

 

 

1)

Uk(x,т) определена,

кусочно-непрерывна и ии(х,т ) е К ( т )

при всех x<=G(т),

 

 

 

 

 

 

 

 

 

2)

решение хк(х) задачи Коши:

 

 

 

 

 

 

 

х(х)

 

= f ( x ( x ) , u k (x{т ),т ),т ),

t <

т < Г; л:(t) =

х

 

(12)

существует и хк (х) в G (х) при ^ <

т <

Г ;

 

 

 

 

 

 

г

 

 

 

 

 

 

 

 

 

 

 

3)

П т Г R (xk (т), uk (т), т) dx =

0,

 

 

 

 

 

 

 

Л**»оо J

 

 

 

 

 

 

 

 

 

 

 

где uk ( x ) = u ( x k (x),x), t ^ x ^ T , (й = 1 , 2,

...,).Тогда последователь­

ность

(хь(т),

uh(x))

{ k = l ,

2, ...) является минимизирующей для

задачи

(1) — (4) и,

таким образом,

функция uh(x,

t),

k = \,

2, ...,

осуществляет приближенный синтез.

 

 

 

 

 

 

Д о к а з а т е л ь с т в о .

Возьмем

произвольную

 

пару

(x(x),u(x))^D^t,T\

при некоторых

фиксированных

t,

to^t<^T,

x ^ G ( t ) . Согласно формуле

(10)

и условиям Xk{t)— x { t ) = x

тогда

имеем

 

 

 

 

 

 

 

 

 

 

 

 

т

j ‘ (х, и (х)) — Д (х, uk (х)) = j* [R (х ), и (x),x)— R (xk (т), 'ик (т), х)] dx. t

По условию теоремы правая часть этого равенства имеет предел

при

k-+oo, поэтому и

левая

часть имеет,

предел. Поскольку в

силу

(8)

 

 

 

 

R (х (х), и (х), х) >

inf

R (х (х), и, х) =

0, t < т < Т,

 

'

иек(т) .

 

то

при любых (jt (т), и (х)) 6 Дс А Г ]. Следовательно,

= inf J* (х, и(х)) > lim J*{x, uk (т)).

D x [ t,T ]

Л->оо

§ 4 ]

 

 

 

Проблема синтеза. Оценка погрешности

 

 

 

207

С другой стороны,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

J<(x,uk ( x ) ) > J * ‘,

k =

1, 2, . . . ,

 

 

 

 

поэтому

lim J* (х, uk (х)) =

J u .

 

 

 

 

 

 

 

 

 

k-+00

 

 

 

 

 

 

 

 

 

 

 

 

 

Т е о р е м а

4.

Пусть

функция В (х,

t)

' довлеткоряет

всем

ус­

ловиям теоремы

1, и пусть последовательность {xft}

такова,

что

 

 

 

xk e X 0 и lim

В (xk, t0) =

inf

В (х, t0).

 

 

 

 

 

 

 

 

 

 

 

k->oo

 

х £ Х о

 

 

 

 

 

Кроме

того,

пусть

имеется

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

функций

uh=

uh

(х, т), обладающая

свойствами

1)— 3) из теоремы 3

при

t— t0, х = хи .

Тогда

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

пар

(хй(т),

ик(т)),

t o ^ x ^ T ,

где uh (x) =

uk (xh(x),x),

а хк(х) — решение задачи

Коши

(12)

при

t=<to,

x=Xh,

является

минимизирующей для

 

задачи

(1 .1 - 4 ) .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Доказательство этой теоремы .опирается на формулу

(10)

при

t= tQи полностью аналогично доказательству теоремы 3.

 

 

 

 

3.

 

Умение решать проблему синтеза крайне важно в

разли

ных прикладных задачах оптимального управления. В самом деле, пусть управляемый процесс описывается условиями (1.1—4), и пусть известна синтезирующая функция u(x,t). Тогда техническое осуществление оптимального хода процесса обычно производится по схеме: с измерительного прибора, замеряющего в каждый мо­ мент t фазовое состояние x(t), на ЭВМ или какое-либо другое вы­ числительное средство подается величина x(t), вычисляется зна­

чение управления u ( t ) = u ( x ( t ) , i ) , после чего

оптимальное в дан­

ном положении управление u(t) передается

на исполнительный

механизм, непосредственно регулирующий требуемое течение уп­

равляемого процесса.

 

 

 

Как вытекает

из вышеизложенных теорем 1—4, для решения

проблемы синтеза

для

задачи (1.1—4) достаточно решить задачу

Коши — Веллмана

(5),

(6). Однако явное аналитическое выраже­

ние решения B(x,t) этой задачи

и функции u(x,t), на

которой

может достигаться

нижняя грань

в (5), удается лишь в

редких

случаях, поэтому задачу (5) , (6) приходится решать приближенно. Наиболее удобными и эффективными при решении задачи (5), (6), по-видимому, являются методы, изложенные выше в § 1, 2, ибо рекуррентные соотношения Веллмана (1.13) и Моисеева (2.5) по существу представляют собой некоторую дискретную аппроксима­

цию задачи (5),

(6) (В (х, th) &Вц{х) ttC k(x)).

Однако здесь воз­

можны и иные

подходы. В тех случаях, когда

удается получить

явное выражение для и=и (х, t, B x) ^ V ( t ) , на котором достигается нижняя грань в левой части (5), то, подставив такое и в (5), при­ дем к обычной задаче Коши для квазилинейного уравнения с част­

208

Д ИНАМ И ЧЕ СКО Е ПРОГРАММИРОВАНИЕ. ПРОБЛЕМА СИНТЕЗА

[,Гл. 4

ными производными первого порядка, и здесь тогда можно вос­ пользоваться известным арсеналом численных методов (разност­ ные методы, метод характеристик, ме*год прямых и т. п.) [20].

В. Ф. Кротов [142] предлагает искать решение B (x,t) задачи

(5)— (6) в виде многочлена по переменным х\ х2....... хп с неоп­

ределенными коэффициентами, зависящими от времени:

т,

т,

тп

 

 

Е

■ • •i

( 0 № ( * * ) '• •••№ •

* 1= 0

*1=0

£„=0

 

Если подставить это выражение для B(x,t) в (5) — (6), то для оп­

ределения

 

,{л (t)

получится система дифференциальных урав­

нений следующего вида:

 

 

 

 

 

 

 

т,

тп

 

 

 

 

B t ( х , 0 =

J . . .

£

(0 (я1/1 . . .

=

 

 

 

 

£.= 0

£„=0

 

 

 

 

=

inf

F (“Фо....О(0,

■••, фт,...тп (0.

U, t),

(13)

 

 

u6V(i)

 

 

 

 

 

 

 

 

x e G ( t) , t0 < t < c T ,

 

 

с начальным условием

 

 

 

 

т,

 

т п

 

 

 

 

 

Е

•••

Е

 

 

=

х е о ( Т ) .

(14)

£,=0

 

£„=0

 

 

 

 

 

Если Ф (х),

inf

F являются

многочленами относительно хх.......

хп,

иек(<)

то, приравняв коэффициенты при одинаковых степенях в (13)— (14), получим задачу Коши для системы дифференциальных уравнений относительно %,...£„ (0> записанной в нормальной форме Коши (см. ниже пример 1). Для ее решения можно использовать различ­ ные численные методы решения задачи Коши (методы Адамса,

Рунге— Кутта и т. д.)

'[20] .

 

Если же Ф (х) и л и

inf Т7

не являются многочленами ОТНОСИ-

 

ueV(£)

тельно хх, ..., хп, то условия

(13), (14) не могут быть, вообще гово­

ря, удовлетворены во всей области G(t), t0^.t^iT, ни при каком

выборе

N = (m i+ 1) (m2+ l ) - ( ^ n + l )

коэффициентов

%,...£„(0-

В этом случае в [142] предлагает задать в области G(t)

N кривых

gi(/), ...,

init)

и рекомендуется определять %,...£„ (0

из условия

удовлетворения

(13), (14) не всюду в

G(t), а лишь на этих кри­

вых. Этот подход перекликается с известными методами коллокации и интегральных соотношений и приводит к задаче Коши для

§ 4] Проблема синтеза. Оценка погрешности 209

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

шенных

относительно

производной

(0 (эти

производные в

уравнения будут входить линейно). Кривые £ i ( 0,

.... £jv(0> во-пер­

вых, должны

образовать

достаточно

густую

сетку в

области

G(t),

 

и, во-вторых, иметь достаточно простое аналити­

ческое выражение (например,

семейство

прямых, параллельных

осям координат, семейство парабол и т. п.).

 

 

 

 

П р и м е р 1.

Пусть требуется минимизировать функционал

 

 

 

 

т

 

 

 

 

 

 

 

 

J (и) — I* и1 (t) dt + kx2 (Т)

 

 

 

 

 

 

 

о

 

 

 

 

 

 

при условиях х — и, х (0) =

х0; и = u(t) — кусочно-непрерывная ска­

лярная функция;

числа Т,

х0 заданы, К =

const > 0.

 

 

Задача (5), (6) здесь имеет вид

 

 

 

 

 

inf [Вх (х,

t) и +

В ( (х, t) +

u2] = ------В\ (х, t) - f B t (х, 0 = 0 ,

х £ Е 1г

«e£i

 

 

 

 

4

 

 

 

 

 

 

 

0 < t < Т ; В (х, Т) = Хх2, х £ Е х.

 

 

 

Функцию B (x ,t)

будем искать в

виде многочлена

В-{х,

t) =

фо(0 +

+ Фг (0 х + Ф2 (0 * 2- Подставим в

предыдущие уравнения

_

 

---- —(Фг + 2ф2х)2 + Фо + Фг* + Фг*2 =

 

х ^Е1у 0 < ^ < Т ;

4

 

 

 

 

 

 

 

 

 

 

 

Фо (Т) + Ф1 (Т) х + ф2 (Т) х2 = Хх2, х Ех.

 

 

Приравнивая коэффициенты

при

одинаковых

степенях х, придем к

следующей задаче Коши:

 

 

 

 

 

 

 

 

ф 0----1-ф| = 0, Ф1 —

ф 1ф 2 = 0 , ф 2—

Ф| = 0,

0 < ^ < Т ,

 

 

Фо {Т) = 0, ф! (Т) = 0, фа ( Г )= Я .

 

 

 

Отсюда находим

 

 

 

 

 

 

 

 

 

 

Фо (0 = Ф ж (0

=

° ,

Фа(*) =

 

Я,

 

 

 

 

%(t_ T) + l

 

 

Таким образом, функция Веллмана здесь имеет вид Яг3

В (х, 0 =

— X(t — T) + 1

синтезирующей является функция

Хх

и(х, 0 — ----- 2 ~ВХ=

X(t Т) — 1

Соседние файлы в папке книги из ГПНТБ