Ф.П. Васильев - Методы решения экстремальных задач (1981) (1158203), страница 65
Текст из файла (страница 65)
Наряду с задачей (1) — (5) рассмотрим расширенную задачу быстродействия. Для этого при каждом е)0 введем множество ()е(Т) всех управлений ие и(!), 1)1„ которые удовлетворяют условиям (2) и для которых существует момент ((и), (О<((и) =Т, такой, что сама последовательность [и») слабо в ~„'[г„ Т'! сходится к некоторому управлению о, ~ !Г (Т).
Согласно (1.12) тогда зцр !х(С и») — х(С и„)/=)»»- О при л-»-оо. ь<~<т Отсюда и из (21) следует, что х(С 0 ) еБ 6» ' "» (г), х((„о,) а=У'»" "», й=1, 2, ... (22) Кроме того, согласно оценке (1.11) lх(г, о») — х(Т», о»)(~С»(Т» — 1») =!)» 1»((=Т». поэтому х(С о») е= 6»» "» а»(1»), 1»(1(Т» (23) х(Тм о„) ен У'»+ "»+а», юг=1, 2, ... (24) Далее, учитывая, что множество 6 (г) непрерывно по Х а усдорф у и й (6 (1) * 6 (Т»)) ( ыо (Т» — 1) ( е»о (Т» — !») = у» при всех С (» 1(Т„имеем 6(Г») ен 6 "»(1) при 1, 1» (1=- Т,.
Отсюда и из (23) получаем включение х((, о )ен6'»+ "»+а»+~»(() при 1» =.1-=Т,. С учетом первого из включений (22) тогда будем иметь, что х(1 о ), 6»»+е»+а,+т»(() при всех 1, (»(1 =То, 1=1, 2, " (25) Поскольку е»+р»+й»+у»-» О при л — ~со, а множества 6((), У замкнуты, то из (24), (25) при й-» со получим, что х(С о,) ен 6((), 1,==(==Т„х(Т„о„) ен У.
Таким образом, о енУ(Т) и время первой встречи 1(о,) траектории х(1, о,) с множеством г" таково, что („(((о,) ( ( Т,. Отсюда имеем неравенство 1,„~ Т,. Выше было установлено, что Т, ( 1,. Следовательно, („ = Т, = =!!ш („(е). Лемма 3 доказана. е 0 Таким образом, если принять У(и)=((и), иенУ= =(/(Т), (л([и)м)=1л([и7м), [и)» вне=(/л(Т), то из лемм 1 — 3 вытекает выполнение условий 1) — 3) теоремьг 2.4. Отсюда следует справедливость утверждения,теоремы 2. 4. Таким образом, показано, что при выполнении условий теоремы 2 для приближенного решения задачи 357 быстродействия (1) — (5) может быть использована последовательность разностных аппроксимирующих задач (7)— (11).
В свою очередь, для решения разностной задачи (7) — (11) при каждом фиксированном гУ можно рассмотреть следующее семейство задач: минимизировать функцию 7го ([и]л, у) = ( хг ([и]л) — у,' (26) при условиях (7), (8) и х,.[[и] ) ен гг~", г= О, 7', у е=- У, (27) где 1 — фиксированный номер, последовательно пробегающий значения г =О, 1, ..., тл., здесь номер тм определяется условием 1,,л, ( Т ( 1 .„, го, момент Т взят из теоремы 2.
Для решения задачи (26), (27), (7), (8) при каждом фиксированном / могут быть использованы известные методы минимизации функций конечного числа переменных или дискретные аналоги методов из главы 1. Обозначим ргл, =1п1 7л ([и]го, у), где нижняя грань берется по всем ([и]л, у), удовлетворяющим условиям (7), - (8), (27). Может слУчитьсЯ, что Ргл ) О пРи всех 1', О~ -=)ь(гг, а Р„и=Π— это значит, что (м„=1ол. Если же Рул) О пРи всех /=О, пгго, то Ясно, что 1л,„)Т)1 гол.
Отсюда следует, что, взяв номер У достаточно большим, согласно теореме 2 в принципе можно получить достаточно точное значение оптимального времени задачи (1)— (5). Однако нужно заметить, что такой подход к решенизо задачи быстродействия на практике может оказаться не очень удобным, поскольку с ростом У растет число задач вида (26), (27), (7), (8) и, следовательно, вообще говоря, растет и объем вычислений. Поэтому желательно иметь другие более удобные методы решения задачи (1) — (5), не требующие перебора всех задач вида (26), (27), (7), (8).
5. Остановимся на одном из таких методов. Для простоты ограничимся рассмотрением задачи быстродействия (1) — (5) при дополнительных предположениях, когда фазовые ограничения отсутствуют, множество У состоит из одной точки, а множество )г(г) не зависит от времени, т. е. 6(()=Е", Ъ'(1)=1' при (о(Г~Т, У=[у). (28) Как и выше„будем предполагать, что матрицы А (г), В(1), 7(1) кусочно непрерывны на любом конечном отрезке 358 [(,, а), множество У выпукло, замкнуто и ограни- чено.
Возьмем некоторое достаточно большое число Т)(„ зафиксируем 1, 1,< Г ~ Т, и рассмотрим задачу У (и, () = ( х (1, и) — у !'-»- (п1, х (т) = А (т) х (т) + В (т) и (т) + [(т), (о~т~( х((а) =ха, (30) и = и (т) ен Ж' = ЪГ ( Т) = [ и (т) е= й~ [(и Т): и (т) е= Ъ' почти всюду на [1„ТД. (31) (29) Заметим, что значения управлений и(т) при (~т(Т на задачу (29) — (31) не влияют. Но тем не менее мы здесь рассматриваем множество (31), так как в дальнейшем нам будет удобно считать, что управления доопределены на всем отрезке [(„Т'1. Обозначим р(()= 1п1 у(и, 1), 1,<( =Т; (32) ием при 1= 1, положим р (1,) = ~ х, — у,".
Будем считать, что р(1,))0, так как при р(1,) =О= ~х,— у!' задача (1) — (5), (28) становится тривиальной: 1„= 1,. Так как множество К слабо компактно в Л,'[1,, Т'1 и функция У(и, 1) слабо непрерывна на )Р', то в (32) нижняя грань достигается, т. е. существует управление и = = и,е=(Р' такое, что р(1)=У(а„()=,'х(1, и,) — у)'. Отсюда ясно, что для того, чтобы момент („был оптимальным временем задачи (1) — (5), (28), необходимо и достаточно, чтобы р(1) О, р(1)>0 при 1, =.1<г„, Р(1)=~х(1, и,) — У1', Р(т)=1х(т, и,) — У~', иь и,.ен)1'. 359 т. е. 1,— минимальный корень уравнения р(1)=0.
Это значит, что для поиска 1, могут быть использованы известные методы решения уравнений. В частности, здесь может быть использован метод, описанный в п. 4 ~ 16 гл. 5 книги [4). Этот метод был описан в предположении, что функция р(1) удовлетворяет условию Липшица. Покажем, что в рассматриваемой задаче зто условие выполняется. Пусть Тогда из определения (32) функции р(1) с учетом оценок (1.8), (1.11) имеем р(г) — р(т)((х((, и,) — у~' — (х(т, и,) — у('» =.2(Со+!у()!х((, и,) — х(т, ит)(-=2(Со+!у!)С~,'1 — т(, р(г) — р(т) ) ~х(г, ид — у/' — (х(т, и,) — у/') ) — 2(С,+ /у/) /х(1, и) — х(т, и)() ) — 2 (Со -(- / у () С~ (1 — т (.
Следовательно, )р(1) — р(т)( й~г — т( (э~1, т=аТ, Е.=2(С,+/у/)С,. (ЗЗ) Для поиска минимального корня 1 уравнения р(г) =0 на отрезке ((„Т1 может быть использован итерационный процесс (~~,=(~+р(1~)(Т., Й=О, 1, ... (34) По условию р(Г,))0. Сделаем индуктивное предположение: пусть р(())0 при (,(1((м Тогда при всех 1, 4а(Г<(эем имеем Р (() ~ Р ((а) — 1. (1 — (а) ) Р (1~) — 1. ((„„., — („) = О. Следовательно, р (Г) ) 0 при (, ( Г < Гь„. Может случиться, что р(1„„,)=0. Тогда Гх„=Ä— искомый корень уравнения р(1)=0. Если р((ьм))0, то процесс (34) продолжаем дальше. Здесь имеются три возможности: 1) Процесс (34) закончится определением момента Гэ такого, что (,<(а<.Т, р(1;))О, 1=0, й — 1, р(гх) =О.
Тогда (а=(, — задача решена. 2) р(Г,))0, 1„<Т, прн всех й=О, 1, ... Так как последовательность (г~) монотонна, то существует 1пп 1„= = („( Т. Из (34) при А -э со получим р (( ) = О. Кроме того, р (() ) 0 при (, == Г < („й = О, 1, ..., так что р (() ) 0 при (,(1<1„следовательно, Г, — искомое время. 3) Найдется номер ги такой, что Г„,~Т<1,„, р((;)) )О при всех ю'=О, ги. Тогда (, )Т. Таким образом, метод (34) позволяет выяснить, будет ли 1„принадлежать [(и Т1, и если это так, то найти приближение (» — („. В том случае, если выяснится, что 1 )Т, аналогично можно продолжать поиск 1„на отрезках )Т, ТД, ~Тм Т,~, ...; ири этом в методе (34) вао каждый раз нужно использовать ту константу Е, которая соответствует рассматриваемому отрезку.
6. Описанный метод (34) предполагает, что значения функции р (1) известны точно. Перейдем к изложению метода (34), свободного от этого недостатка. А именно, предположим, что при каждом 1, 1,(1( Т, вместо точного значения р(() может быть найдено его приближение рл(1) такое, что ~р(1) — рх (1) )($м рл(1) О, 1о =1==Т, (35) где (фл)-~0. Пусть известно (У вЂ” 1)-е приближение (х „ У)1 (начальное приближение 1, нам задано). Для определения следующего приближения (л рассмотрим итерационный процесс 1ньы=1мь+рл(1мл)71- й=О 1 " !мо=Го (36) аналогичный процессу (34). Из-за погрешностей равенство рл (1) = 0 не влечет за собой равенства р (1) = 0 и наоборот, поэтому процесс (36) следует прекращать не по критерию рл((х~)=0, как было выше, а по условию вида рх (1дъ)(Ех, где (Вм) стремится к Нулю при М-~-оо согласованно с погрешностью Д,ч) из (35).
В рассматриваемом случае, когда функция р(1) удовлетворяет условию (33), такое согласование означает, что Ех'--2Цх. У=1, 2, ... (37) Имеются две возможности: 1) (м,(Т при всех lг=О, 1, ... Так как ((„а) монотонна, то существует !пп (,чд(Т. Тогда из (36) следует, что !нп рл ((л'а) = О 2) Существует номер з такой, что (х, ( Т ((л,+,. Отсюда следует, что за конечное число шагов процесса (36) найдется номер В=А такой, что будет выполнено одно из двух следующих условий: р (1„,) ~ Ел, й = О, т р (1 ) е, 1„.
Т, (36) или Рл(Ь„,))ах, й=О, тл — 1, 1и -т~Т<Лл . (39) Положим (л = пнп (1лг , 'Т1, ,'40) Метод описан, ЗЕ1 Т е о р е м а 3. Пусть выполнены условия (33), (35), (37). Тогда последовательность ((н), определяемая методом (36), (38) — (40), сходится к ш)п (Юь; Т).