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

книги из ГПНТБ / Острем К.Ю. Введение в стохастическую теорию управления

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

200

 

 

 

Глава

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-\-k1. Показать, что

закон управления

 

 

 

 

u(t)=-qkG'{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-\-k1 и п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