книги из ГПНТБ / Острем К.Ю. Введение в стохастическую теорию управления
.pdf200 |
|
|
|
Глава |
6 |
|
|
|
|
|
|
|
единичного круга, а В2— полином степени п2, |
все нули |
которого |
||||||||||
находятся вне единичного |
круга. |
|
|
|
|
|
|
|
|
|||
Для |
решения уравнения |
(5.5) |
наложим |
дополнительное ог |
||||||||
раничение на G(z), |
состоящее в том, что G(z) содержит |
B2(z) |
||||||||||
в качестве множителя, т. е. вместо выражения |
(5.5) использу |
|||||||||||
ем равенство |
|
|
|
|
|
|
|
|
|
|
|
|
|
^ k - ^ |
= |
Д |
р' |
+ |
£ 2 |
^ ^ |
|
|
( 5 |
g ) |
|
Рассуждая так же, как и при доказательстве |
теоремы |
4.1, на |
||||||||||
ходим закон управления |
|
|
|
|
|
|
|
|
|
|
||
|
|
u(t) = |
|
q k G ' i q ) |
|
у {t). |
|
|
|
|
(5.10) |
|
При этом ошибка управления равна |
|
|
|
|
|
|
|
|||||
у (0 = X {е (0 + f ,е (/ - 1) + • • • + fk_,e |
(t - k + 1) - f |
|
||||||||||
|
+ fke(t-k)+...+fk+rh_le{t-k-n2^\)). |
|
|
|
|
|
(5.11) |
|||||
Закон |
управления |
(5.10), |
который |
не является |
оптимальным, |
|||||||
дает ошибку с дисперсией |
|
|
|
|
|
|
|
|
|
|
||
|
var у = min (var у) + W {/f t 2 |
+ |
• • • + |
|
|
) . |
|
(5.12) |
Закон управления (5.10) не очень чувствителен к изменениям параметров системы. Для доказательства этого предположим опять, что система задается моделью (А0 , 5°, С°, Х°), а закон управления определяется исходя из модели (А, В, С, X) с не много другими параметрами. Уравнение, описывающее управ ляемую систему, приобретает вид
[Л» (<7) Вх (q) F (q) + В0 (q) G' (q)] у (t) = X°C° (q) BL (q)F (q) e(t). (5.13)
Когда параметры модели равны истинным параметрам, ха рактеристическое уравнение системы имеет вид
|
z „ 2 + f e - i B o ( 2 ) C o ( z ) = |
0 . |
|
( 5 . 1 4 ) |
|
Из |
определения В\ и допущения, сделанного относительно С0, |
||||
следует, что все колебания устойчивы |
при равенстве расчетных |
||||
и |
истинных параметров. Устойчивость |
при небольших |
измене |
||
ниях параметров следует теперь из непрерывности. |
|
|
|||
|
Сравнение выражений (5.7) и (5.14) |
показывает, |
что при |
||
законе управления (5.10) нули полинома В2 переносятся |
в на |
||||
чало координат. |
|
|
|
|
|
В2 |
При других квазиоптимальных стратегиях управления |
нули |
|||
переносятся не в начало координат, |
а в произвольные |
точки |
внутри единичного круга. Следовательно, всякий раз, когда по лином 5° имеет нули вне или на границе единичного круга, оп тимальная система будет настолько чувствительной к измене-
Стратегии управления, минимизирующие дисперсию |
201 |
ниям параметров, что становится непригодной. Бесконечно малые изменения параметров могут дать бесконечную диспер сию. Даже если В° или С° имеют нули внутри единичного кру га, оптимальная система может все же быть чувствительной к изменениям параметров. Проиллюстрируем сказанное на при мере.
Пример |
|
|
|
|
Рассмотрим систему, описываемую |
уравнением |
|
||
у (/) + ahj (t — 1) = и (t — 1) + |
b°u (t — 2) |
+ е (t) + с°е (t— |
1), (5.15) |
|
где а ° = 0 , 7 , й°=0,99 и |
с°=0,95. Для |
получения стратегии уп |
||
равления, минимизирующей |
дисперсию, используем |
равенство |
||
(4.13): |
|
|
|
|
(l+c0q-l) |
= |
{l+a°q-l)+q-lgo. |
|
Следовательно, g"o=c°—а°. Стратегия управления, минимизи рующая дисперсию, имеет вид
и(0 = = |
C°~a\y(t) |
= |
l-^—1y(t). |
(5-16) |
' |
1+6°<7 |
|
1 + 0,99?—1 w |
' |
Допустим, что закон управления определяется исходя из си стемы с другим значением коэффициента при u(t—2). Тогда по лучим следующий закон управления:
u{t) = |
h&—y{f). |
При этом законе управления выходная переменная равна
= |
0> + |
а{Ь + |
с«)+ЬсО |
|
<72 + |
Я (Ь + |
с°) + |
(6 — 6°) а° + с°6° |
w |
Ее дисперсия определяется выражением
V(b) = Еу* = ±§Н(г)Н |
(г"1 ) -f- . |
На рис. 6.5 показано, как дисперсия выходного сигнала за висит от параметра Ъ. Из рис. 6.5 видно, что оптимальная стра тегия более чувствительна к изменениям параметра Ъ, чем ква зиоптимальная.
Для получения закона управления, не являющегося опти мальным, но менее чувствительного по сравнению с оптималь ным к изменениям параметра Ь, можно использовать закон управления (5.10). Из равенства (5.9) получаем
(1 + с У 1 ) = (1 + а У ) (1 + У ) + |
H+bYl)8- |
202 Глава 6
Следовательно, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
f = |
^ (со |
аруфо |
_ |
а о) = |
0,966, |
|
|
|
||||
|
|
g = — а? (с0 — с0 ) (6° — а0 ) = |
0,684, |
|
||||||||||
и закон управления |
(5.10) |
принимает вид |
|
|
|
|
|
|||||||
|
U (Л = |
|
|
ё - т - y(t)= |
|
|
0,684 |
, |
• у (0 |
(5.17) |
||||
|
W |
1 + |
|
w |
|
|
Ц - О . Э б б ? - 1 |
у |
w |
|
||||
|
з.о |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2,5 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Квазиоптимальная |
|
|
|
|
|
|
|||
|
1,5 |
|
|
|
стратегия |
|
|
|
|
|
|
|
||
|
с |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
=4 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1.0О |
|
i |
I |
i |
i |
1 |
I |
I |
|
I |
1,0 |
|
|
|
|
|
|
|
0,5 |
|
|
|
|
|
|
|||
Р и с . |
6.5. Зависимость |
дисперсии |
выходного |
сигнала |
от |
параметра Ь. |
||||||||
Из |
выражения |
(5.11) |
следует, |
что |
стратегия |
управления |
||||||||
(5.17) дает следующую ошибку управления: |
|
|
|
|
|
|||||||||
|
у (t) = |
е (t) + |
fe (i — 1) = |
е (t) + |
0,966е (t — 1). |
|
||||||||
Таким образом, |
стратегия |
управления |
(5.17) |
дает |
диспер |
|||||||||
сию vary=l,93 . (Минимальная дисперсия равна |
1.) |
|
|
|||||||||||
Допуская, что стратегия управления |
определена исходя из |
|||||||||||||
модели с параметром Ь, отличным от истинного значения |
Ь° для |
|||||||||||||
системы, находим, что ошибка |
управления |
равна |
|
|
||||||||||
|
УК> |
|
q2 + |
q(a° |
+ f + g) + |
a°f + |
b°g |
|
|
|
|
|
Вычисляя дисперсию выходной переменной для различных значений параметра Ь, получаем результат, приведенный на
Стратегии управления, минимизирующие дисперсию |
203 |
рис. 6.5. Видно, что закон управления (5.17) гораздо менее чув ствителен к изменениям параметра Ь, чем оптимальный закон управления (5.16).
5 |
Время, t |
|
|
|
Р и с . 6.6. Моделирование вы.ходкого и управляющего сигналов системы |
(5.15) |
||
при оптимальной стратегии управления (5.16). |
|
||
На рис. 6.6 и 6.7 показаны модели выходных и |
управляющих |
||
сигналов, когда система управляется |
по оптимальному |
(5.16) |
|
и квазиоптимальному законам управления (5.17). |
|
|
|
Упражнения |
|
|
|
1. Для системы |
|
|
|
у (t) + 0,640 (* — 1) + 0,22// (/ - |
2) = 6,4м (/ - |
3) + |
|
+ 19,2и (t — 4) + Я [е (t) — 0,82е |
0,21е (t — 2)] |
|
204 |
Глава 6 |
Р и с . 6.7. Моделирование выходного и управляющего сигналов системы (5.15) при квазиоптимальной стратегии управления (5.17).
определить |
закон |
управления, |
минимизирующий |
дисперсию, |
|
и исследовать |
его чувствительность. Построить |
закон управле |
|||
ния, который |
менее |
чувствителен |
к изменениям |
параметров. |
|
2. Сравнить дисперсию ошибки управления |
при |
законе уп |
равления (5.10) с дисперсией ошибки {k—п2) -шагового упредителя.
3. Заменить равенство (5.9) следующим равенством:
H(q)C(q) |
= A(q)F(q) |
+ |
Ba(q)G'(q), |
|
где Я — произвольный |
полином |
степени n2-\-k—1. Показать, что |
||
закон управления |
|
|
|
|
u(t)=-qk—G'{q) |
|
|
y{t) |
|
W |
4 В, |
(9) |
F' (q) |
* W |
приводит к ошибке управления |
|
|
|
H{q)y{t) = -kF(q)e(t),
а характеристическое уравнение замкнутой системы имеет вид
H(z) Вг (z)C(z) = 0.
4. Для системы, заданной уравнением (4.8), показать, что стратегия управления
Стратегии управления, минимизирующие дисперсию 205
"(0 = • y(t),
где полиномы Fi и G\ степеней n-\-k—1 и п—1 соответственно определяются равенством
дает замкнутую систему с характеристическим уравнением 2п + * - , С(г) = 0,
которая не очень чувствительна к изменениям параметров. Вы числить дисперсию ошибки управления и сравнить с минималь ной дисперсией.
5. Для системы
y(t) + ay(t—l) = u(i—l) + 2,5u(t-2) +
+ u(t—3) + e(t) + ce(t— 1)
определить стратегию управления, минимизирующую дисперсию и соответствующую ошибку управления. Показать, что опти мальная стратегия чрезвычайно чувствительна к изменениям па раметров, и вывести с помощью методов, рассмотренных в раз деле и в упражнении 4, квазиоптимальные стратегии управле ния, которые не очень чувствительны к изменениям параметров. Определить ошибку управления для двух квазиоптимальных стратегий управления.
6. ПРОМЫШЛЕННОЕ ПРИМЕНЕНИЕ
Введение
В данном разделе рассмотрен |
пример промышленного при |
менения теории, изложенный в разд. 4 и 5. Материал взят из ра |
|
бот, выполненных в скандинавской |
лаборатории фирмы I B M в |
Стокгольме в связи с установкой управляющей вычислительной
машины |
на |
бумагоделательной |
фабрике |
Billerud |
в |
Гравене |
||||
(Швеция). Большая часть материалов взята из журнала |
фирмы |
|||||||||
I B M 1 . Автор |
выражает глубокую |
признательность |
фирме |
IBM |
||||||
за разрешение использовать эти материалы. |
|
|
|
|
|
|||||
Стохастическая теория |
управления |
была |
использована |
для |
||||||
решения |
некоторых задач |
управления |
качеством, |
а также |
уп |
|||||
равления |
процессом в связи с |
исследованиями |
на |
фабрике |
Billerud. Мы опишем лишь приложение к одной типичной зада че — к задаче управления весом бумаги. Рассмотрим примени-
1 Computer Control of a Paper |
Machine — an |
Application of Linear |
Sto |
chastic Control Theory, IBM Journal |
of Research |
and Development, July |
1967. |
2 0 6 |
Глава S |
мость линейной |
стохастической теории управления и дадим ма |
тематическую постановку задачи. Кратко обсудим также зада чу получения математических моделей возмущений и динамики процесса по экспериментальным данным. Естественно, что это очень важная задача для практического приложения. Наконец,
приведем |
некоторые эксперименты |
по управлению |
в |
контуре |
с помощью оптимальных стратегий. |
|
|
||
Перед установкой вычислительной машины для управления |
||||
ситемой были исследованы флуктуации веса бумаги. |
Флуктуа |
|||
ции имели стандартное отклонение, равное 1,3 гс/м2 . |
Плановая |
|||
величина |
стандартного отклонения |
веса бумаги при |
управлении |
|
с помощью вычислительной машины была установлена |
равной |
|||
0,7 гс/м2 . |
|
|
|
|
Управляющая вычислительная машина была установлена в |
||||
декабре |
1964 г. Два эксперимента |
для определения |
динамики |
процесса были проведены в марте 1965 г., а первая работа с ис пользованием машины в контуре управления весом бумаги была выполнена в апреле 1965 г. Эксперимент длился 10 ч. С тех пор выполнено большое количество экспериментов, а с начала 1966 г. непрерывно действует замкнутая система управления ве сом бумаги. Теперь при работе системы можно постоянно дости гать стандартного отклонения 0,5 гс/м2 для веса во влажном со стоянии и 0,3 гс/м2 для веса бумаги в сухом состоянии.
Применимость стохастической теории управления
Теория, описанная в разд. 4, базируется на следующих до пущениях:
1. Динамика процесса описывается линейными дифференци альными уравнениями с постоянными коэффициентами. Имеет ся одна выходная и одна входная переменная.
2. Возмущения являются стационарными гауссовыми процес
сами |
с дробно-рациональными спектральными |
плотностями. |
3. |
Критерий оптимальности — минимизация |
дисперсии вы |
ходного сигнала. |
|
В начальной стадии проекта были выполнены эксперименты, которые показали, что возмущения, встречающиеся в процессе нормальной эксплуатации, настолько малы, что систему можно описать линейными уравнениями. Результаты, полученные в этих экспериментах, были впоследствии проверены в экспери ментах над системой, управляемой с помощью вычислительной машины.
Обычный спектральный анализ показал, что флуктуации мож но представить стационарными процессами. Было показано так же, что возмущения имеют почти нормальное распределение.
Иногда возникают «нарушения», которые приводят к боль-
Стратегии управления, минимизирующие дисперсию 207
шим отклонениям. Причинами возникновения нарушений мо гут быть, например, неисправности оборудования. Для устране ния нарушений требуются особые корректирующие действия. Маловероятно, что эти типы возмущений можно описать вероят ностными моделями. В дальнейшем эти нарушения не будут учитываться.
Вес является важной качественной переменной для крафтбумаги. При продаже бумаги обычно указываются допуски пе ременных качества. Контрольные процедуры потребителей обыч но осуществляются таким образом, что бумага принимается, ес ли переменные качества для контрольной выборки лежат внутри контрольных пределов с заданной вероятностью. Так как в про цессе нормальной эксплуатации всегда существуют флуктуации качества, производители бумаги выбирают задание для регуля тора веса бумаги значительно выше нижнего контрольного пре дела, чтобы быть уверенным в том, что их продукция удовлет ворит требованиям потребителей. Уменьшение колебаний веса бумаги позволяет сдвинуть задание для регулятора ближе к приемлемому пределу и при этом не изменить вероятность попа дания внутрь контрольных пределов (рис. 6.8). Это дает воз-
Р и с . 6.8. Распределение плотности вероятности для двух процессов.
208 Глава 6
можность увеличить прибыль, которая может быть выражена в экономии сырья пли в увеличении выпуска продукции. В нашем частном случае эта прибыль составила существенную часть при были, полученной от внедрения управления с помощью вычис лительной машины.
Производство бумаги с небольшими срлуктуациями перемен ных качества имеет также другие преимущества, например уп рощение последующей обработки бумаги покупателем. Однако эти преимущества очень трудно оценить объективно.
Так как флуктуации веса бумаги нормально распределены, задача установления задания как можно ближе к допустимому пределу сводится к такому управлению процессом, при котором дисперсия выходного сигнала была бы наименьшей.
Подведем итоги. Были найдены веские доводы, что теорию, описанную в разд. 4, можно применить к задаче управления ве сом бумаги. И ничего нет странного в том, что эта теория была действительно применена. Доводы о применимости теории были проверены в том смысле, что результаты, предсказанные стоха стической теорией управления, были получены на практике.
Математическая постановка задачи управления весом бумаги
На рис. 6.9 показана упрощенная схема той части бумагоде лательной машины, которая представляет интерес для управле
ния весом бумаги. Густая волокнистая |
масса, т. е. смесь фибры |
с водой (концентрация фибры ~ 3 % ) , |
поступает из машинной |
камеры. Волокнистая масса перемешивается с чистой водой, так что концентрация фибры в напорном ящике уменьшается до 0,2—0,5%. На сетке фибра отделяется от воды, и образуется по лотно бумаги, из которго на прессах удаляется вода. После это го бумага обезвоживается на сушильных цилиндрах сушильной части машины.
В данном случае на вес бумаги можно влиять путем измене
ния потока густой волокнистой массы |
и (или) |
ее плотности |
(т. е. концентрации фибры в густой волокнистой |
массе). Обе |
|
эти переменные непосредственно влияют |
на количество фибры, |
вытекающей из напорного ящика, а следовательно, и на вес бу маги.
Воздействие на управляющие переменные осуществляется путем выбора задания аналоговых регуляторов, которые управ ляют задвижками, регулирующими поток и плотность густой во локнистой массы.
Вес бумаги определяется измерителем с бета-лучами. Выход ной сигнал с бета-измерителя пропорционален массе фибры и воды на единицу площади, т. е. весу влажной бумаги, так как ко-
|
Стратегии |
управления, минимизирующие дисперсию |
209 |
|||
эффициеиты |
абсорбции |
бе |
|
|||
та-лучей в фибре и воде при |
|
|||||
близительно одинаковы. Для |
|
|||||
получения веса |
сухой бума |
|
||||
ги, т. е. массы фибры на еди |
|
|||||
ницу |
площади, |
показания |
|
|||
бета -из м ер ителя |
|
дол ж н ы |
|
|||
быть |
скорректированы |
на |
|
|||
вел ич и ну, |
уч итыв а ющую |
|
||||
влажность бумаги. |
Измери |
|
||||
тель |
влажности |
емкостного |
|
|||
типа в нашем случае может |
|
|||||
перемещаться вдоль полотна |
|
|||||
бумаги, хотя обычно он ус |
|
|||||
танавливается в |
фиксиро- |
|
||||
в а и ном полож еи ии. |
|
|
||||
Бета-измеритель |
уста |
|
||||
навливается |
также |
и перед |
|
|||
сушильной |
частью. |
Вес |
бу |
|
||
маги измеряется |
оператором |
|
ив контрольных лаборато риях. Когда рулон бумаги готов, определяют его вес и размер, что дает очень точ ную величину среднего веса рулона. Эта информация ис пользуется для градуировки других измерителей. Анализ информации о весе бумаги показал, что:
1.Информация о размере
ивесе рулона бумаги может быть использована для кор рекции измерений при дрей фе бета-измерителя.
2.Высокочастотные флук туации в измерителе влаж ности и бета-измерителе имеют сходные характери стики, а хорошей оценкой веса сухой бумаги служит величина
у = WSP (1 — MSP), (6.1)
где WSP — градуированный сигнал с бета-измерителя, а
14—403