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

книги / Системы экстремального управления

..pdf
Скачиваний:
8
Добавлен:
19.11.2023
Размер:
33.28 Mб
Скачать

Для иллюстрации на рис. 5.2.4 показаны примеры про­ цессов поиска для различных значений q. Здесь же пока­ зано построение процессов по касательным к характери­ стике объекта (пробные шаги не показаны). Хорошо вид­ но изменение характера поиска с увеличением q.

В случае, когда величина Q* неизвестна, ее можно оцевить адаптивным образом. Основная идея здесь сводится к следующему. Если направления рабочих шагов изме­ няются, то при q 1 это означает, что выбранное значение Q* меньше истинного. Если же знаки двух последователь­ ных рабочих шагов Ах совпадают, то следует уменьшить Q*. Этот алгоритм можно записать в виде

Q*N+1 = QN — v sign (AxNAXn+1),

(5.2.18)

где v > 0 —постоянная. При 1 /2 < £ < Ч , как легко убе­ диться, следует действовать наоборот, т. е. изменить знак в (5.2.18) на обратный.

Любопытно, что уменьшение Q* по сравнению с истин­ ным значением позволяет повысить быстродействие поис­ ка, хотя и вносит некоторую «нервозность» в его работу.

Как видно, этот метод поиска особенно эффективен для кусочно-линейных объектов, однако он работает и для других объектов, в том числе и клювообразных, но при

Я> 1/2.

§5.3. Поиск с квадратичной экстраполяцией

Впредыдущем параграфе рассмотрен поиск, предпо­ лагающий достаточную линейность характеристики объ­ екта (в кусочном смысле). При этом было показано, что по мере отклонения характеристики объекта от линейной скорость сходимости уменьшается. Однако подобный экстраполяционный алгоритм может опираться не только на линейность, но и на другие свойства характеристики объекта. Так, предполагая, что характеристика объекта достаточно близка к квадратичной параболе, можно по­ строить поиск с квадратичной экстраполяцией.

Построим поиск с экстраполяцией в предположении о квадратичности экстремальной характеристики объекта

Q = к (х х*)ъ +Ç*-

(5.3.1)

Для удобства запишем это выражение в виде

Q = ахг

-\-Ьх -|- с.

(5.3.2)

Минимум этой функции,

как известно,

расположен в

точке

_Ь_

 

х*

(5.3.3)

2а

и равен

 

 

 

9 ’ =

« - -С - •

(5-3.4)

Задача экстраполяционного поиска заключается в опреде­ лении смещения

Да?! — х* х0,

(5.3.5)

которое сразу переводит систему из исходного состояния в экстремальное х* (в предположении о квадратичности характеристики объекта).

Для определения положения цели нужно определить коэффициенты параболы (5.3.2). Это делается при помощи трех замеров значений этой параболы. Пусть замеры про­ изводятся в районе исходного состояния в точках х0 g, х0 и а:0 + g и при этом получаются три значения показа­ теля качества, которые и составляют исходную информа­ цию для процесса экстраполяции:

<?i = Q (*о — g),

Qz =

Q (s„),

(5.3.6)

Qz =

Q (яо +

g)-

Подставляя эти значения в уравнение объекта (5.3.2), получаем в результате систему из трех линейных уравне­ ний с тремя неизвестными a, b и с:

л{х0

g)2 +

Ъ(х0g) -j- с =

Qi,

 

 

 

ÜX\ -j- ЬХу -J- c =

Q2,

(5.3.7)

a (x0 +

g)a +

b (x0-j- g) + c =

Q3.

 

Как известно, для того чтобы эта система имела решение, необходимо, чтобы ее определитель не был бы равен ну­ лю, т. е.

(®о — g f

(х0 — g)

1

 

XQ

XQ

1 Ф 0.

(5.3.8)

(s0 + g)2 (s0 + g) 1

Раскрывая этот определитель, записываем условие в виде - 2 ^ 0 , (5.3.9)

что выполняется всегда при g Ф 0.

Разрешая (5.3.7), получаем интересующие нас значения параметров объекта а, 6 и с:

Qi — 2Qa -b Qa

2g2

Ъ— "2^г ((?з — Qi)

2ж0(Çx — 2Qz +

Q$)], >(5.3.10)

c = Qa — (Q3- Q 1)2 Zo

Qi — 2Qa + Q3

j i

2g

2g5

X ° '

Подставляя эти значения в (5.3.3), получим окончательно величину рабочего смещения (5.3.5) в виде

Д д.__ 8_

Qa— Qi

(5.3.1 1

2

"Q i— 2Q a+ Qa

 

Любопытно, что полученное значение рабочего смещения не зависит от исходного состояния объекта ха. При этом ожидаемое значение экстремума (5.3.4) равно

(Qa — Q i)a

.

(5.3.12)

8 (Q i — 2 Q3 + Q3) ’

 

как видно, оно также не зависит от ж0.

 

После выполнения рабочего шага

Ахх следует прове­

рить, действительно ли найден экстремум. Для этого достаточно определить показатель качества в предпола­ гаемом экстремуме хх = х0 + Aæx и сопоставить его с оценкой (5.3.12). Если эти величины отличаются не более

чем на е > 0, то

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

можно

считать выполненной. Если же

 

|Д<?ж1 =

l« - e ( * i ) l > e ,

(5-3.13)

то это означает, что исходное предположение о квадратичности характеристики не выполняется и следует продол­ жить процесс поиска, т. е. выполнить следующий цикл, но уже эксперименты ставить в районе точки xv И так далее.

Блок-схема описанного алгоритма показана на рис. 5.3.1. Здесь, как видно, для функционирования алгоритма необходимо иметь два блока памяти для запоминания промежуточных значений показателя качества. Оператор

«Дж» производит определение рабочего смещения по формуле (5.3.11) на основе имеющейся информации. Оператор «AQ» обозначает вычис­

ление отклонения показателя ка­ чества, полученного в эксперимен­ те, от вычисленного по формуле (5.3.12) путем экстраполяции. Опе­

ратор

«Г» представляет собой вы­

 

держку системы в течение времени

 

Т, после чего снова начинается

 

поиск. Введение

этого оператора

 

связано

с возможностью

дрейфа

 

экстремума объекта в процессе его

 

эксплуатации.

Если

экстремум

 

уходит медленно, то значение вре­

 

мени выдержки Т выбирается боль­

 

шим. В противном случае вели­

 

чина

Т

выбирается малой. Если

 

же задача заключается в одно­

 

кратном определении экстремума,

 

то оператор «Г»

является

опера­

 

тором

полного

останова

поиска.

 

Функционирование такого ал­

 

горитма

обеспечивает

отыскание

 

экстремума за один цикл для объ­

 

ектов с квадратичной

характери­

 

стикой и весьма быстрое прибли­

 

жение к экстремуму для объек­

 

тов

с

характеристикой

более

 

высокого порядка.

Работа алго­

 

ритма показана

на рис. 5.3.2 на

 

кубическом

объекте.

 

 

Рис. 5.3.1. Блок-схема

На

рис.

5.3.3

для

приме­ программыпоиска с квад­

ра представлены результаты со­

ратичной экстраполя­

поставления

двух

 

экстраполя­

цией.

ционных

поисков

(линейного и

 

квадратичного) на объекте с кубической характеристикой

Q = 1*1* +<?*.

(5.3.14)

минимум которой расположен в начале координат. Нагляд­ но видно, что квадратичная экстраполяция для этого объекта более эффективна, чем линейная.

Однако следует помнить, что поиск с квадратичной экстраполяцией, как и всякий поиск, учитывающий априорные особенности объекта, имеет ограниченное при­ менение. Более того, существует целый класс объектов, для которых этот поиск неприемлем, так как на них он

Рис. 5.3.2. Поведение объекта с кубической характеристикой при поиске с квадратичной экстраполяцией.

становится неустойчивым. Так, при оптимизации объек­ тов с кусочно-линейной или клювообразной характеристи­ ками поиском с квадратичной экстраполяцией пользовать­ ся нельзя.

Рассмотрим теперь аппаратурную реализацию этого алгоритма. На рис. 5.3.4. показана блок-схема экстремаль­ ного регулятора. Здесь использованы те же элементы памяти, которые применяются и в предыдущих схемах экстремальных регуляторов. В результате последователь­ ной работы двух элементов памяти П получаем величины

à ?Q

=

AQ i — AÇg = Çi —

Qsi

 

AÇi

+

Д(?2 = Qs Qu

(5.3.15)

которые необходимы для определения

рабочего

шага

Д.г вычислительным устройством ВУ по (5.3.11).

В заключение этого параграфа отметим, что аналогично

могут быть построены улучшенные алгоритмы

поиска

при учете других априорных сведений об объекте. Очевид­ но, что каждый такой алгоритм будет очень хорошо работать на объектах того класса, для которых он был создан.

Рис. 5.3.3. Сравнение поведения объекта с кубической характери­ стикой при поиске с линейной (а) и квадратической (б) экстрапо­ ляцией.

Рис. 5.3.4. Блок-схема экстремального регулятора, работающего по алгоритму поиска с квадратичной экстраполяцией.

В качестве упражнения читателю рекомендуется ана­ логично изложенному выше построить и проанализировать алгоритм поиска с квадратичной экстраполяцией, но с известным экстремальным значением функции качества Q*. В этом случае для организации поиска можно ограничить­ ся двумя замерами функции качества в районе исходной точки.

§ 5.4. Обобщение экстраполяции

Рассмотрим теперь произвольный вид модели объекта в виде разложения по некоторой заданной системе функций

m

 

 

 

(5.4.1)

Здесь

фх (х),

., срт (х) —заданная система

функций,

кг,

кт —неизвестные коэффициенты.

Например,

если систему функций выбрать в виде

 

«Pi (я) = 1, ф2 (ж) = ж, ф3 (я) — ж2,

то получим известную квадратичную модель (3.1.9), где

кх Q* -J-кх*, кг = 2кх*, kz = к.

Экстремальное управление объектом, функция качества которого имеет вид (5.4.1), в этом случае будет сводиться к определению неизвестных коэффициентов кх, ., кт и к последующему расчету положения экстремума х*. Это и есть экстраполяция.

На первой стадии для определения k t (i = 1, ., т) необходимо сделать не менее чем тп замеров показателя качества в точках xi = х ± ig, i = 0, 1, ., где х — исходная точка, вокруг которой определяются значения показателя качества. Так, для m = 5 показатель качества определяется в точках х 2g, х g, х, х + g, х -f 2g, и для определения коэффициентов А* составляется система линейных уравнений:

5

Q i = Q (ж — 2g) = 2

kWi ix — 2£)>

i=l

 

5

 

< 2= Q (x g) = 2

S)>

i=l

 

5

(5.4.2)

 

i=l

 

6

 

Qi = Q(x + g) = S

+ e)i

i=l

 

5

 

Qe = Q - f 2g) = 2 (* + 2g),

где в левой части стоят конкретные значения показателя качества. Разрешая эту систему относительно к{ (г = = 1, .,5 ) получаем конкретное выражение для функ­ ции качества (5.4.1), из которого легко вычислить положе­ ние экстремума. Для этого следует решить очевидное уравнение

<iQ(др

2 h

(x)

= o.

(5.4.3)

dx =

dx

 

i=l

 

 

 

Решение xxявляется первым приближением. Если экстре­ мальная характеристика в точности аппроксимируется мо­

делью (5.4.1), то полученное решение хх совпадает с положением экстремума объекта, т. е. решает поставлен­ ную задачу:

Q (хх) = щ т Q (х).

(5.4.4)

X

 

Если же представление (5.4.1) имеет локальный характер, то придется сделать несколько шагов экстраполяции, где исходной точкой N -то этапа будет (N — 1)-е приближе­

ние XN-I .

Критерием остановки процесса будет выполнение условия

N

I Q (XN) — 2 MW)(Pi(tfîf) I < e,

(5.4.5)

i= l

 

где e 0 — заданная точность оценки экстремального

значения, k[N) —оценка i-ro коэффициента на N -м. шаге. Блок-схема описанной системы экстремального регу­ лирования с экстраполяцией показана на рис. 5.4.1. Здесь П —блок памяти, в котором запоминаются значе­

ния показателя качества Ç{N\ ., Qm* на каждом (N-м) этапе экстраполяции. Вычислительное устройство (ВУ) выполняет две функции: определяет значение коэффициен­

тов ki*\ • -, f f î на каждом шаге, т. е. решает систему уравнений типа (5.4.2), и затем определяет положение эк­ стремума на N-м шаге, т. е. решает уравнение (5.4.3). Исполнительный механизм (ИМ) отрабатывает полу­ ченное значение. Координирует работу системы блок

управления БУ. На него же возложена проверка выпол­ нения условий остановки (5.4.5).

Как видно, рассмотренная схема экстраполяционного поиска может эффективно решать задачи оптимизации

Рис. 5.4.1. Блок-схема

системы экстремального регулирования

с

экстраполяцией.

объектов с весьма сложными экстремальными характе­ ристиками. Однако следует отметить, что для этого необ­ ходимо иметь специализированные или универсальные вычислительные устройства, решающие задачи (5.4.2) и (5.4.3). Естественно, что введение такого устройства в схему экстремального регулятора весьма усложнит его и тем самым значительно снизит надежность эксплу­ атации. Именно этим обстоятельством следует объяснить отсутствие экстремальных регуляторов с экстраполяци­ онным поиском. Ио появление ВУ на интегральных схемах, по-видимому, даст возможность преодолеть ука­ занную трудность и реализовать экстраполяционные ал­ горитмы для экстремального регулирования.

НЕПРЕРЫВНЫЕ АЛГОРИТМЫ ПОИСКА

(БЕЗЫНЕРЦИОННЫЕ ОБЪЕКТЫ)

В предыдущих трех главах рассмотрены шаговые ал­ горитмы поиска. Отличительным свойством этих алгорит­ мов является дискретный характер их работы. Это ока­ зывается удобным для объектов, в которых функция ка­ чества образуется дискретным образом, например, когда показатель качества определяется путем расчетов для фиксированного состояния оптимизируемого параметра.

Однако имеется весьма широкий класс экстремальных объектов, для которых существует непрерывная динами­ ческая связь между входом и выходом. Использование этой связи дает возможность воспользоваться большей информацией об объекте и приводит к непрерывным алго­ ритмам поиска.

В простейшем случае объект может быть, например, безынерционным и тогда всякое изменение входа х мгно­ венно отражается на выходе объекта Q. Но часто связь между входом и выходом имеет инерционный характер. Это накладывает определенный отпечаток на поведение объекта в процессе оптимизации.

Однако разделение объектов на инерционные и безы­ нерционные, как указано выше (см. § 2.2), имеет относи­ тельный характер и зависит от скорости изменения управ­ ляемого параметра в процессе оптимизации. Если эта скорость мала, то динамическими свойствами объекта моясно пренебрегать и считать объект безынерционным. Но для повышения быстродействия необходимо увеличи­ вать интенсивность поиска путем повышения скорости изменения параметра х. В этом случае нельзя не учи­ тывать инерционных свойств объекта.

Рассмотрим отдельно поведение безынерционных и инерционных экстремальных систем в процессе поиска. Начнем с управления безынерционными системами, ко­ торым и посвящена эта глава.