Файл: Живоглядов, В. П. Адаптация в автоматизированных системах управления технологическими процессами.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 23.10.2024

Просмотров: 82

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

N — константа, не зависящая от я [я—х—1]. Минимизируя Ri по и[я—х—1], находим

и*[я—х— 1 ]= и[п —х—2] Н— 1 + g^ j (д*—т[п—х—2]—

 

 

—и[п—х—2}).

(1.

92)

Продолжая вычисления далее, устанавливаем,

что /?я_ т_ 5

и «[s] при любом

s имеют

вид, аналогичный (1- 91) и

(1. 92), т. е.

 

 

 

 

и*[«] = ф - 1 ] +

1_1: , г„1

- - • о1- ( 7 * - /и[^ -1]

- “[«—!]),

 

1+с[п—х—s]

(1.

93)

 

5=1,2,.--,Я —х

где

 

 

 

 

 

 

■с[п—х—s]= -

 

(1.

94)

 

1 +

1+с[я —X—S—1]

 

 

 

1+ ...

 

 

с[0]

—с

 

 

 

1+

с

 

 

 

 

Т+~с

 

 

 

2(n—x—s) раз

 

 

 

При (я—s)-+°о получаем

 

____

 

 

с[п—S—х ]^ с[я -х —S1 ]—Сое = __

 

(1.

95)

'Структура управляющего

устройства показана на

рис.

1.5,

где Д t— звено задержки

на

один такт,

Ф— дискретный

•фильтр, Т— усилитель с коэффициентом

---- •

Из

схе-

 

 

1

1^оо

 

 

мы (рис. 1.5) легко получить блок-схему системы управле­

ния инерционным объектом 6 первого порядка

с запазды­

ванием. Однако критерий оптимальности при

этом име­

ет вид

 

6!)



 

W={(f-q[s])2+c(q[s]—g[s— I])®,

(1. 96)

где '

c = i_[(2Cco+1)2- 1 ] .

 

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

^ = ( ? - ф ] - ^ 2+ + ( ф ] - + о м ) - 2

U- 97)

Введем обозначения

/Ч/

л*/

<7* = 9* _ н ном,

e[s] = -a[s]—ином.

При этом

Ws=*(q*—iф ] - н .) 2+ с 0и2И .

Пользуясь изложенной выше методикой, устанавливаем, что оптимальное управление в s-м такте равно

=

(т*— —1]).

(1. 98)

5=1,2,-..,/1—т

70

Таким образом, во всех рассмотренных случаях задача синтеза алгоритма управляющего устройства (УУ) распада­ ется на две: а) синтез УУ в предположении полной инфор­ мации о р и б) нахождение оптимальной оценки параметра р — информационной координаты т. Система управления яв­ ляется нейтральной с пассивным накоплением информации.

Оценки

т определяются формулами (1. 83), (1. 87) или

(1. 84),

(1. 85). При отсутствии априорных сведений о р дис­

персия

устремляется в бесконечность и переменные коэф­

фициенты в формулах оценок

принимают вид

ф ] = - .

s = l,2,...

Ф ) = у ,

*>0.

Вторая достаточная статистика — апостериорная дисперсия Ds (или D( ) — непосредственно в алгоритмы управления

(1. 83), (1. 84) не входит.

Посмотрим, какой получится алгоритм управляющего уст­ ройства, если применить при нахождении оценки т метод стохастической аппроксимации в том виде, как он изложен, в работе [В. 17]. Оценку т выберем такой, чтобы выражение

Я = М { Ш - Ф т = Щ Ш - К о “[8-']-Ко*п)*\

достигало минимума. В точке оптимума градиент AR равен нулю:

Текущая оценка возмущения р определяется следующей формулой:

0- 99)

Аналогично для непрерывного времени получаем

(1 . 1 0 0 )

71


Коэффициенты a[s] и a(t) должны удовлетворять следующим условиям [В. 27, В. 28]:

00

^ ] ф ] = ос',

s=0

СО

[ a{t)dt = °° ,

о

СО

 

 

 

a*[s]<oo,

(1.

101)

II

О

 

 

со

 

 

00

 

 

 

Гa2(/)d f< « -

(1-

102)

О

 

 

 

Алгоритмы (1. 101), (1. 102) обеспечивают сходимость оценок к истинным значениям ц с вероятностью, равной еди­ нице, и в среднеквадратичном, причем последовательность a [s] и функция a(t) могут выбираться в некотором отноше­ нии произвольно. Сравнивая алгоритмы (1. 99), (1. 100) и (1. 87), (1. 85), полученные в теории дуального управления,

замечаем, что они совпадают, если выбрать

a[s]

и a(t) в

(1. 99), (1. 100) в соответствии с формулами

(1. 83),

(1. 84).

Таким образом, оптимальные в смысле минимума полного риска коэффициенты a[s], a(t) даются теорией дуального управления.

1. 3. 2. Управление линейным объектом

при равномерно распределенных помехах

Следует отметить, что метод стохастической аппроксима­ ции далеко не всегда приводит к оптимальным алгоритмам. Пусть, например, действующие в цепи обратной связи помехи h и возмущение ц имеют не нормальное распределение, а рав­ номерное, причем —B^h*CB,— Л"0=1, функция потерь вида (1. 82). Значения помехи в различные моменты времени статистически независимы. Алгоритм, основанный на методе стохастической аппроксимации, остается прежним, ли­ нейным. Оптимальный же алгоритм приобретает вид

u-*s= Q*

2

1 us—1—х )max”b

us—1—x

)mln]>

где (y^_i—^

t )

max — максимальное значение

(1- 103)

из £ раз­

ностей

72.


(у9- и _ т ),...,(Ув-1 - « в - 1- т

Ш - А ) , a

(ys- , —«5- 1-х )min—

минимальное

значение из

(у0 —и _ .

),...,(ys- i —us—'l—z )>

{А—В). В непрерывной

системе оптимальное

управление

■определяется формулой

 

 

 

 

iu*(t) =<7*

2"[(у(0

ll{t

^иМтах-!- {У(0 u(t

'ш)}min] ■

.1. 3. 3. Применение распределенного контроля

Исследуем, измениться ли структура управляющего уст­ ройства при использовании распределенного контроля выход­ ного сигнала q% (х, t). Обращаясь вновь к методу информа­

ционных координат, убеждаемся, что введение дополнитель­ ных точек измерения состояния объекта (распределенного контроля) влияет лишь на определение апостериорной плот­ ности Я ^ р ). Алгоритмы управления остаются прежними, т. е. (1. 83), (1. 93) или (1. 84). Форма Ps (р) не меняется; меня­ ется лишь алгоритм вычисления информационных координат m[s—1], m(t). Пусть объект описывается уравнением (1. 71) с начальным условием q (х, о )= р и граничным условием q (0, i) = p-f-н (t). Функция состояния объекта при этом равна

7/? (*-0 = 11+ u(t —-^-), /> 0 , 1и> х > 0 .

В УУ поступает смесь сигнала и распределенной центрирован­ ной гауссовой случайной помехи h(x, t) типа белого шума

y(x,t) — q(x,t)+h(x,t).

Задача, как и выше, сводится к нахождению оптимальной

•оценки т возмущения р. Применяя для этой цели относитель­ но простой метод максимума правдоподобия, находим

 

 

U

т —с

j-B

de},

=«(*){J

 

(1- 104)

где

х

a{t) =

1_

 

7 ’

 

t '

73