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

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

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

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

Добавлен: 11.04.2024

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

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

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

210

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

[ Г л . 4

4.Предположим, что, пользуясь тем или иным методом, уд

лось

получить

некоторое приближенное решение B (x,t)

задачи

(5),

(6). Если

эта функция получена разностным методом

(напри­

мер,

методами

§ 1, 2) на

какой-то дискретной сетке точек, то до­

определим ее

(например,

интерполяцией) во всех точках

области

x^.G(t), t o ^ t ^ T , до некоторой непрерывной, кусочно-гладкой функции B(x,t). Тогда функции u — u(x,t), на которой точно или

приближенно реализуется

inf R ( x ,u ,t )

можем

принять

в

каче-

 

 

H£V«)

 

 

 

для

задачи

стве приближенного решения проблемы синтеза

(1.1—4). Это значит, чтТГ приближенное решение

(х(т),

ц (т))

зада­

чи (1) — (4)

будем определять из условий

 

 

 

 

 

 

*

ф = / (X(т), и (х (т), г), т), t <

т < Т;

х(() =

х,

 

 

 

-

-

 

Т.

 

 

 

(15)

х (т) 6 G (т),. и (т) =

и (х (т), т), t < т <

 

 

 

 

Приближенное решение исходной задачи (1.1—4) находится аналогично: сначала определим точку ЯоеЛ'о, на которой точно или приближенно реализуется inf В {х, t0), затем, решая задачу Коши

 

 

 

х £ Х о

 

 

 

_

 

 

(15) при_ t = t 0,

х = х 0,

найдем траекторию х(х)

и управление

и (т) = и(х,(т), т ) ,

 

 

Найденную пару (х (т),

и(т )) примем

в качестве приближенного решения задачи

(1.1— 4).

Спрашивает­

ся, какая при этом будет допущена погрешность?

 

 

Т е о р е м а 5.

Пусть

функция В (х, t) кусочно-непрерывна, ку­

сочно-гладка при x ^ G ( t ) ,

t o ^ t ^ T ,

и

В (х , Т )= ф ( х )

x ^ G (T )

(подчеркнем, что B(x,t)

не обязана

быть

решением

уравнения

(5 )). Пусть для

каждой

 

пары (х(т),

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

функция

Б (х (т ),т )

переменной т

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

при всех x ^ G ( t ) , to^.t<.T. Тогда для любой фиксирован­

ной пары

{х(х), и(х) ) e D JY , 7] справедлива оценка

 

 

где

 

o < j ' ( * , « ( 0 ) - • / * '<

М «),

 

(16)

 

 

 

 

 

 

 

 

 

 

 

 

J*‘ =

inf

J ‘ (x ,u ( т))

 

 

 

 

 

 

 

DXU.T]

 

 

 

 

 

 

_

Г

_

 

_

 

 

 

 

 

 

ei (и) = J [Я (* (*),

и (т),

т) — Rmin (X)) dx,

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

Япцп(т )=

inf

inf

R(x,u, т);

 

 

 

 

 

 

y 6 G ( t ) u G K ( t )

 

 

 

 

функция R(x, и, т) определена равенством (7). Кроме того, для любой пары (х(х), u ( x )) ^ D Xo[t0, Т] справедлива оценка

0 J*° (хи (т)) J

е„ (х0, и),

(17)


§ 4]

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

211

где

 

 

 

 

J* =

inf inf До (лг, и (()),

е0 (х0, и) =

f [R (х (т),

и (т), г) —

 

х & Х а D x[ t „ T ]

 

£

 

 

— ^min W ] dx +

в (xQ, g —

inf В (x, g

,

X 0 = { x : x 6 G ( g , £>* [ g T] Ф 0 }.

Д о к а з а т е л ь с т в о . При выполнении условий настоящей теоремы формула (10), очевидно, остается справедливой. Возьмем произвольную пару {х{%),u(x))<^DJ[t,T]. Пользуясь формулой (10),

имеем

 

 

 

 

т

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

J ‘ (х,

и (/)) — J t (х,

и(х)) =

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

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

R {х (т), и (т), т)] dx < г( (и).

 

 

 

 

Так как вг(й) не зависит от

(x(x),u (x))^ D jt,T ],

то отсюда

сразу

получаем оценку (16).

 

 

 

 

 

 

 

 

Аналогично,

взяв

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

допустимую

 

пару

(jc(t ) , u (t ) ) g D J 10,7 ]

при некотором х е 1 0,

с помощью

формулы

(10) получаем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(х0, и (т)) — До (х, и (т)) =

 

 

 

 

т.

 

 

 

 

 

 

 

 

 

 

 

= § [R (x (т),

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

и (т), т)] dx +

В (х0, t0) —

 

to

 

 

— В(х, g < е 0(х0, и),

 

 

 

 

 

 

 

 

 

 

 

 

 

откуда немедленно следует оценка

(17).

Д,

 

 

 

 

 

В тех случаях, когда конструктивное описание множества Х0

затруднительно, в выражении для

во(х,й)

вместо Х0 часто

берут

G(t0). Очевидно, при такой замене величина ео(х0,й)

может лишь

увеличиться.

 

 

 

 

 

 

 

 

 

 

 

Из полученных

оценок (16),

(17)

видно,

что

чем

ближе

Я (х (т),гГ (т),т)

к

# т 1п(т), В (х0, g

— к inf В (х,

t0),

тем

меньше

 

 

 

 

 

 

x£Xq

 

 

 

 

погрешность. Отсюда ясно, что при практических вычислениях же­ лательно выбирать функцию В(х, t), которая удовлетворяет урав­

нению (5) по возможности

точнее, а функция u(x,i)

должна

по­

точнее реализовать

in fi?(x,

u,t). Тогда на парах (х (т ),« (т )),

оп-

ределяемых из (15

u6V(l)

 

Т\ (в частно­

)>при различных x ^ G (t) и f e [ g

сти, х = х 0, t = t 0),

погрешность в определении минимального зна­

чения функционала будет небольшой.

 

 


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

Если полученная точность неудовлетворительная, то надо ис­ кать новые функции B (x,t), u(x,t), используя, например, упомя­ нутые выше методы решения задачи Коши — Веллмана с дальней­ шим дроблением разностной сетки или увеличением степени мно­

гочлена и числа кривых

{£*(£) } e G ( i ) и т. п. Постепенно уточняя

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

вообще говоря, можно построить последо­

вательности B m(x,t), Wm(x,t), для которых погрешность монотонно убывает. Если окажется, что при т-*-оо погрешность стремится к нулю, то функцию Um{х, t) при достаточно больших т можно при­ нять в качестве приближенной синтезирующей функции.

Различные аспекты проблемы синтеза рассмотрены в работах

[13, 14— 18, 24, 27, 160, 171, 195, 229] и др.

Упражнения. 1. Решить проблему синтеза для задачи миними­ зации функционала

т

J (и) — ^ (х* + ы2) dt

о

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

2. Найти решение задачи Коши — Веллмана для задачи м нимизации функционала

 

 

 

т

 

 

 

 

J {и) =

т

I [{а {t)>

х {t)) + ъ (“ W* 0] dt + (с, -V[Г])

 

 

 

о

 

 

 

при условиях

x = A ( t ) x + C ( u ( t ) ,t ) , О ^ г ^ Г ;

x(Q )=x0,u(t)<=V(t);

u(t)

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

где

известными

предполагаются мо­

мент

Т, матрица

A{t) порядка

пХп, л-мерные вектор-функцин

C(u,t), a(t), скалярная функция b(u,t), л-мерные векторы х0, с и

множества V ( t ) ^ E r, 0^.tsS^T.

У к а з а н и е . Функцию B(x,t)

искать в

виде

многочлена

первой степени переменных х = (х\

..., хп) :В ( х ,

0 =

(ф(^),х).


Г л а в а 5

Достаточные условия оптимальности

При решении задач оптимального управления часто возникает следующий вопрос: будут ли на самом деле оптимальными те уп­ равления и соответствующие им траектории, которые мы нашли, используя какие-либо точные или приближенные методы решения таких задач? Такой вопрос, например, естественно возникает, когда управление и траектория найдены из краевой задачи прин­ ципа максимума, поскольку принцип максимума выражает собой необходимое условие оптимальности, не являясь в общем случае достаточным для оптимальности. Известные к настоящему времени достаточные условия оптимальности либо опираются нд^свойство выпуклости данных задачи, либо тесно связаны с динамическим программированием [23, 24, 59, 141, 142, 161, 206]. .Здесь.ограни­ чимся изложением подхода В. Ф. Кротова [142].

§ 1. ДОСТАТОЧНЫЕ УСЛОВИЯ ОПТИМАЛЬНОСТИ ДЛЯ ЗАДАЧ

СЗАКРЕПЛЕННЫМ ВРЕМЕНЕМ

1.Рассмотрим следующую задачу оптимального управления: минимизировать функционал

J (и) = J ( x ( t ) , u(t)) —

= f /° (X (t),

и (/), t) dt +

Фх (x (T)) +

ф0 (X (*„))

■ (1)

^0

 

 

 

 

при условиях

 

 

 

 

X (t) =

f ( x (t), и (t), t), t0< t < T ,

(2)

 

x(t)eG (t),

t0 < t < T ,

 

(3)

u = u (t)eV {t),

f0 < f < 7 \

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

(4)

где моменты tQ, T будем считать

известными,

Ф0(х) и Ф](д:) — за ­

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

в §

3.1; граничные условия при t— t0 и t = T включены в ограниче­

ния

G(t).

 

 

 

 

 

Пару (x(t), u{t)) назовем допустимой на отрезке

 

если управление u { t)^ V {t), а функция x(t)

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

гладка на

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

(2), (3). Множество

всех допустимых

пар

(x(t), u(t)),

обозначим

через

D[tQ, Г]. В пространстве

Е пу,Етвведем множество Dt точек

(х, и)