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

книги из ГПНТБ / Егоров С.В. Элементы идентификации и оптимизации управляемых систем учеб. пособие

.pdf
Скачиваний:
5
Добавлен:
24.10.2023
Размер:
9.47 Mб
Скачать

деления градиента функции. Она имеет первостепенное зна­ чение при создании экстремальных систем регулирования.

Пункт 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

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