книги из ГПНТБ / Егоров С.В. Элементы идентификации и оптимизации управляемых систем учеб. пособие
.pdfделения градиента функции. Она имеет первостепенное зна чение при создании экстремальных систем регулирования.
Пункт 3. Задача определения динамических характеристик
Эта задача является обобщением задачи п. 1, когда связь (1-3) между входными и выходными переменными является операторной, а статическая модель не устраивает исследо вателя, что обычно наблюдается при анализе и оптимизации переходных процессов в системе.
Рассмотрим для простоты определение динамической ха рактеристики только одного канала объекта — между точкой
ции
приложения воздействия v(t) и точкой измерения выходной переменной y (t). Действие помех измерения и остальных
воздействий учтем в приведенной к выходу аддитивной по мехе
Схему идентификации удобно представить как схему сравнения выходных переменных объекта и его неизвестной
А
модели (рис. 1-4), где до, до означают операторы соответст венно указанного канала объекта и его модели.
А
Требуется найти такой оператор до при заданных v(t), y(t) и некоторых дополнительных ограничениях, чтобы не
которая мера несоответствия между объектом и его моделью была минимальной. В качестве дополнительных ограничений принимают (Колмогоров А. Н., 1933 г., Винер Н., 1942 г.):
А
стационарность оператора ДО, бесконечность интервала на блюдения, эргодичность и некоррелированность между со-
бой оипналов v(t) и £(f). а также условие, что оператор до
принадлежит определенному классу, например, физически реализуемых линейных операторов (см. гл. 2).
21
'Если мерой точности идентификации (критерием) слу жит величина среднеквадратичной ошибки
М {е2} = М {{у —y)J} -* min, |
(1-11) |
W
то inp-и вышеуказанных условиях оптимальный по указанно му критерию оператор находится из решения относительно
А
до уравнения Винера—Хопфа
00
Ryv(t) = [w (x)R av{ t- x )d x |
(М 2) |
О |
|
при известных корреляционных функциях Ryv, Rvv, при этом
А
до (т) — импульсная характеристика указанного канала объек та. Однако .практическое решение уравнения (1-12) сталки вается с явлением нерегулярности (см. гл. 3), когда даже небольшие неизбежные погрешности в определении корре ляционных функций вызывают громадные ошибки в полу-
А
чаемой характеристике ДО(т). Кроме того, многие из указан ных дополнительных ограничений на практике не выпол няются.
Разработка методов решения указанной задачи определе ния динамических характеристик систем в настоящее время идет как по пути развития методов регуляризации решения (Тихонов А. Н.), так и по пути преодоления обременитель ных на практике дополнительных ограничений (Пуга чев В. С., Солодовников В. В., Заде Л. и Рагаццини Дж. и
ДР-)-
Обратим внимание, что рассмотренная задача, как и за дача п. 1, формулируется как некоторая оптимизационная задача.
Рассмотрим теперь типичные задачи оптимизации.
Пункт 4. Задача раскроя
Имеются заготовки, например, прямоугольные листы ма териала определенного размера, из которых .надо накроить детали различного профиля <в количествах, не меньших за данных величин b{(i—1, .... |§), при этом известно п спосо
бов раскроя одной заготовки.
Рассмотрим /-й .способ раскроя, когда из одной заготовки
получается |
aiS деталей /-го профиля. Обозначим |
щ —число |
заготовок, |
раскраиваемых по /-му способу (/ |
= 1, .... т ) - |
Тогда оптимальным решением задачи раскроя будет такое, при котором необходимое количество деталей будет получе но из минимума израсходованных заготовок, т. е. минимизи руется критерий оптимальности — линейная форма
_ |
/П |
|
|
V |
пип |
|
|
|
|
|
|
при линейных ограничениях |
|
|
|
^ ац Uj > |
b, |
i = 1, |
л , |
/=1 |
|
/' = 1, . . . , т . |
|
.щ > 0 |
Эта задача оптимизации относится к задачам линейного программирования. Она была первой практической задачей,
решенной в СССР методами математического программиро вания (Канторович Л. В., 1939 г.).
Пункт 5. Транспортная задача
Имеется п пунктов производства однородной продукции,
на складах которых сосредоточены объемы продукции в ко личествах а\, .... а„. Имеется m пунктов потребления с объе мами потребления Ь\, ..., Ьт. Известны величины cti — затра
ты на перевозку единицы продукции /-го пункта производ ства в /-й пункт потребления.
Требуется найти такой |
объем |
(план) |
перевозок u,j |
из |
|
t-х пунктов в /‘-е (/= 1, ..., п; |
/ = |
1, |
..., т ) , при котором были |
||
бы удовлетворены потребности |
в |
пунктах |
потребления, |
а |
затрать^на переврзки были минимальны. При этом предпо лагается, что производство и потребление сбалансировано, т. е.
пгп
i- 1 /■=- 1
Таким образом минимизируется критерий оптимально сти — линейная форма
п т |
|
L — V. 2 . Crf U[j -v |
mill |
«=i/=i |
!“«/} |
33
при линейных ограничениях:
т
a) |
U/j — (i{ i — 1< • • • > |
/=i
т. e. из каждого t-го пункта производства продукция выво зится полностью;
|
л |
|
б) |
2 w"' — hi |
/— i , . . . , |
|
t=i |
|
т. e. потребность в' /-м пункте потребления удовлетворяется полностью.
Кроме того, очевидно, что при отсутствии обратных пере возок все
Как и в предыдущем примере, имеем задачу Минимизации линейной формы при линейных ограничениях, причем время в явном виде нигде не фигурирует |(задача линейного про граммирования (АП). По имеющимся сведениям, примерно
85% всех решенных в США задач ЛП являются задачами транспортного типа.
Пункт в. Двойственная задача хозяйственного планирования
Рассмотрим задачу взаимоотношения предприятия с внешним рынком.
Пусть предприятие располагает ресурсами bI, ..., b„, ко-
торые надо использовать для |
выпуска т видов |
продукции, |
||
объем производства которой для каждого вида |
и}*»...... и',1’ |
|||
надо выбирать таким образом, чтобы: |
|
|
||
а) |
общая стоимость G| |
продукции на рынкебыла мак- |
||
симальной: |
|
|
|
|
|
т |
|
шах |
|
|
Ci и( 1) |
|
||
|
/= * |
|
ц( 1)1 |
|
б) расход ресурсов не превышал их наличия:
т |
|
2 ац и\Х) <■ b{ |
i =z 1 , . . . . /I. |
/=i |
|
При атом очевидно
4 ] > 0 / = 1, . . т .
24
Здесь Cj — установленная на рынке стоимость единицы
/'-го вида продукции;
a,j — технологическая норма расхода i-ro ресурса на
производство единицы /-го вида продукции. Допустим, что некоторая рыночная организация решила
купить все ресурсы, которыми располагает хозяйство. Тре
буется установить цены на эти ресурсы uf\ ..., |
. таким |
||
образом, чтобы: |
|
|
|
а) |
общая стоимость Gj |
покупаемых ресурсов |
была ми |
нимальной: |
|
|
|
|
0 t = |
-»■ min; |
|
|
* |
Ю |
|
б) при этом продажа ресурсов была выгодна и предприя тию, т. е. надо, уплатить предприятию неменее тех сумм, которые оно само может выручить при переработке ресурсов в готовую продукцию:
2 a</U*'2) > С! |
|
II| |
|
Та |
|
- |
V о |
-Б |
/=
II
>
1, • .., ш.
п.
Таким образом, здесь имеются две взаимосвязанные за дачи оптимизации, когда каждая из сторон (предприятие и рынок) стремится получить выгоду. По форме критериев оп тимальности и характеру ограничений эта задача относится к двойственным задачам линейного программирования.
Пункт 7. Задача размещения центра
При построении телемеханических систем с радиальной структурой, телевизионных станций для обслуживания ряда населенных пунктов, помещений КИП на крупных промыш ленных комбинатах и т. п. возникает задача нахождении координат пункта |(центра), суммарное расстояние которого до остальных пунктов обслуживания минимально. При этом на допустимую область размещения центра имеются огра ничения, связанные с рельефом местности и наличием вод ных поверхностей, с (необходимостью размещения в районе источника энергоснабжения и т. п. Некоторым обобщением этой задачи на многомерный случай является задача Штай нера, формулировка которой следующая.
25
В п-мерном евклидовом пространстве заданы N точек
С,= (с|„ .... c„i) и выпуклое множество £2ИТребуется отыс
кать во множестве Q„ такую |
точку U=(u\, .... |
и„), |
чтобы |
|
суммарное расстояние от нее до точек С( было |
минималь |
|||
ным: |
|
|
|
|
N |
N |
/ п |
|
|
q = У , р(и, < ? , ) = |
2 2 |
V 2 2 ( « / ■ - сп? - |
, n i " |
■ |
<=1 |
«-I |
/-1 |
(«/’ |
|
Поскольку функция р является выпуклой, а сумма выпук лых функций тоже выпукла, задача ее минимизации при вы пуклом ограничении £2„ является задачей выпуклого про граммирования. В более общем и на практике часто встре
чающемся случае ограничения не являются выпуклыми, более того — допустимая область размещения может быть неодносвязной. Это соответствует более общей задаче нели нейного программирования.
Пункт 8. Оптимизация режима дистилляционной колонны
Следуя Хэммонду и Рису, рассмотрим диотилляционную колонну (рис. 1-5), где приняты'обозначения:
J-----
Ч
- - - Н -*-f
|
|
Рис. 1-5. |
Дистилляционная |
|
|
колонна как объект управле |
|
|
|
|
ния |
г ], |
Zi — расход первичного продукта в колонну и его состав; |
||
г/i, |
г/2 — выход вторичного продукта и его состав; |
||
|
L,t — расход вторичного продукта, возвращаемого в колон- |
||
|
• |
«у; |
|
|
и1 |
— тепло, подводимое к кипятильнику; |
и2 —коэффициент рецикла (Mj= Ln-# r1).
25
Управлять режимом дистилляционной колонны можно, варьируя переменные и\, м2 в некоторой ограниченной обла
сти (рис. 1-6).
Критерий оптимальности является функцией от количе ства получаемого вторичного продукта, его состава и расхо да тепла:
Q= —Pit/i + М г/г*—уг\ +Р3И1-*- min
N ии иг ’
где <Pi, р2, р3 — заданные весовые коэффициенты, у2* — ка чественный состав стандарта.
На рис. 1-6 показаны экспериментально полученные за
висимости критерия |
оптимальности |
Q (щ, иг) |
(при Pi = 5, |
Р2 = 10, Рз= 0,1, Z| = |
18,3 ед., z 2 = 0,5 |
молярной |
фракции), |
для различных у\, у г*, при этом заштрихованная зона соот
ветствует неустойчивому режиму. Таким образом, перемен-
Рис. 1-6. Область изменения управлений и линии ранных значений критерия оптимальности для трех режимов
ные «|, |
Щ ограничены прямыми со штриховкой и должны, |
|
кроме |
того, удовлетворять очевидным условиям |
гц^О, |
«2^ 0. Линии уровня на рис. 1-6 проградуированы |
так, что |
|
Qi <Q2<Q3- |
|
27
Как видно из (рис. 1-6, при изменении режима работы ко лонны, оптимальные значения управлений изменяются весь ма значительно. При этом зависимость Q(ui, щ) при раз
личных возмущающих воздействиях имеет нелинейный ха рактер, .но математическая модель процесса отсутствует. Если пренебречь инерционностью колонны, то получаем за дачу минимизации нестационарной нелинейной функции, имеющей экстремум и задаваемой экспериментально, при ог раничениях (задача экстремального регулирования).
Пункт 9. Задач^ минимального расхода топлива
Такая задача возникает, например, при управлении лета тельными аппаратами (от самолетов до космических кораб лей), морскими судами и т. д. Поскольку топливо имеет су щественную стоимость, надо найти такое управление аппа ратом, чтобы, например, перелет .самолета из одного пункта в другой был произведен при минимальном расходе топлива. Задача экономии топлива весьма остро встает также в кос мических полетах на заданное или максимальное расстояние, когда дополнительно приходится учитывать ограниченность топливного ресурса.
Пусть уравнения движения управляемой системы имеют вид
X = F(X , U, t), |
X(t0)=Xo, |
(1-13) |
где Х=(х\, *2, Хз) — вектор |
пространственных |
координат |
системы; |
|
U — вектор тяговых сил двигателей.
Допустим, что управление системой осуществляется с по мощью m двигателей, потребляющих топливо и создающих
определенную тягу. Далее будем считать, что конструкция аппарата или двигателей такова, что позволяет получать тягу обеих знаков, соответствующую разгону и торможению. Если учесть, что топливо расходуется независимо от знака управления (и при разгоне, и при торможении, например, космического корабля, выполняющего мягкую посадку), то расход топлива /-м двигателем можно принять пропорцио нальным абсолютной величине управляющего воздействия, гак что суммарный мгновенный расход топлива «семи дви гателями равен
m
M = 2c;|U;|, (Cj > 0),
/= 1
28
при этом за время Т— /о будет израсходована масса топлива,
равная
т
У — J М dt. • min. |
( Ы 4 ) |
I.и
Очевидно, что величина топливного ресурса и тяговые усилия двигателей ограничены:
J ^ M 0,
I Щ \ |
макс. / = 1) •••» т - |
(1*15) |
Пусть задано целевое множество й т '(подвижное или не подвижное). Тогда задачу минимального расхода топлива можно сформулировать следующим образом.
Для заданной системы (1-13) найти такое управление, удовлетворяющее условиям (1-15), которое переводит систе му из начального положения Хо в йт при минимальной ве личине функционала (1-14). Время Т в приложениях может
быть фиксировано, ограничено или свободно.
Отметим, что коэффициенты уравнений движения (1-13),
.в которые неявно входит масса систему, являются перемен ными, поскольку в массу системы входит масса топлива, убы вающая во время полета. Однако в ряде случаев, например, при полетах самолетов, можно считать, что изменение мас сы системы за счет расхода топлива несущественно.
Рассмотренная задача является задачей динамической оп тимизации.
Пункт 10. Задача об оптимальном регуляторе
Многие управляемые системы (энергосистемы, непрерыв ные технологические процессы и др.) наиболее эффективно функционируют лишь в равновесном состоянии или вблизи него. Под действием возмущений или вследствие внутренней неустойчивости такие системы могут быть выведены из рав новесного состояния. Тогда возникает задача: найти такое допустимое управление, которое наилучшим образом воз вратило бы систему в состояние равновесия и удерживало бы ее в этом состоянии. Предполагается при этом, что функция управления выполняет автоматическое устройство (регулятор), которое формирует закон управления по ин формации о действительном состоящий системы и, если воз можно, о возмущениях. Такие задачи рассматривались Лето-
29
ным А. М., Красовским Н. Н., Калманом Р. и др. Пусть X*
является состоянием равновесия системы
X = F(X , U, Z), Х(/0)= Х 0 |
(1-16) |
(далее будем считать Х* = 0), а возмущения |
Z заданы как |
функции времени. Их действие может также проявляться в
произвольности начального состояния |
системы (1-16). |
Заданы ограничения ,на управления |
|
H (U )^ 0 . |
(1-17) |
Требуется перевести систему (1-16) с помощью управле |
|
ния U(X} /), удовлетворяющего условию |
(1-17), из началь |
ного состояния Ха в начало координат таким образом, чтобы
величина функционала 7 была минимальна. При атом за кри терий оптимальности могут быть взяты:
|
|
|
|
|
|
|
|
•• |
|
|
|
а) |
быстродействие, когда 1 J = |
J Gi (X)dt, |
|
||||||||
где |
G,(X) |
О |
при |
|
|| X |
|| < |
е, |
|
|
|
|
.1 |
при |
|
ИX |
|| > |
•; |
|
|
|
|
||
|
|
|
|
|
|
|
|||||
б) |
квадратичное отклонение, |
когда |
J |
Хт Q(t)Xdt, |
|||||||
где Q —матрица коэффициентов |
(пХп) |
положительно оп |
|||||||||
ределенной квадратичной формы; |
|
|
|
|
|||||||
в) |
квадратичное отклонение и энергия управления, когда |
||||||||||
|
|
•о |
|
|
|
|
|
|
|
|
|
|
J = |
)' (X* Q(t) ■X + |
U' |
R(i) |
U)d(, |
|
|||||
|
|
0 |
|
|
|
коэффициентов |
|
|
|||
где Q, R — матрицы |
|
(nXn), |
(tnXm) |
||||||||
положительно |
определенных форм, |
|
|
|
|||||||
г) |
быстродействие |
|
и |
энергия |
управления, |
когда |
|||||
|
|
J ~ |
J00 |
(0 ,(Х) | |
и * |
RU)dt, |
|
||||
|
|
|
о |
|
|
|
|
|
|
|
|
положительно определенных форм;
Заметим, что указанная регулярная постановка задачи не всегда правомерна, поскольку реальные возмущения часто являются случайными процессами и поэтому не могут быть заданы как заранее известные функции времени.
30