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

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

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

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

Добавлен: 23.10.2024

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

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

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

Поскольку риск г* s |_т инвариантен относительно u[s—1], то

u*[s—1] имеет вид, аналогичный (1. 119), и оптимальное уп­ равление для любого момента времени определяется форму­ лой (1. 119). Конкретный вид формул для вычисления инфор­ мационных координат т и D зависит от статистических ха­ рактеристик p[s] и h[s]. Например, для независимых и нор­ мально распределенных /г [s] и p[s] = р,= const (как в подраз­ деле 1. 3. 1) имеем

 

 

S—1--X

 

 

ms_i[s]

 

1

q { /+т]

'

 

Н10+ S

«[/1

(1. 121)

 

 

 

 

/= 1

 

 

 

 

Dя—'2

 

 

 

D s—1

S

 

( 1. 122}

 

1+ ■Qs—а

 

 

 

 

 

 

 

°2/г

 

 

где а2Л ,

, т0— априорные дисперсии

и математическое

ожидание.

 

 

 

 

1. 4. ДУАЛЬНОЕ УПРАВЛЕНИЕ

МАРКОВСКИМ РАСПРЕДЕЛЕННЫМ ОБЪЕКТОМ

1.4. 1. Управление по одной выходной точке

Вданном подразделе мы рассмотрим задачу оптималь­ ного дуального управления, имеющую практическую значи­ мость для ряда процессов промышленной технологии непре­ рывных или полунепрерывных производств. Нормированные корреляционные функции р (v) возмущений в таких объектах часто аппроксимируют [1. 18, 1. 19] экспонентами вида

 

—«114 I

 

 

 

(1. 123)

 

P(v) = e

.

 

 

 

Пусть требуется

стабилизировать выходной

сигнал

q(t)=>

= <7д(/н>$ реактора (1. 75)

при условии,

что

/(м,р-)=«+р и

возмущение р.(/)

представляет собой

марковский центриро­

ванный гауссов

случайный

процесс

с

дисперсией

о3^ и

80


корреляционной функцией

в стационарном режиме

вида

(1. 123)-Запишем уравнения

v

 

объекта для i >

 

qR (x,i)=K0(x)

i

 

q{t)=K0[и(/—т„)+ р-(/-тн)1,

 

 

(1.

124)

y(t) = q(t)+h(t)

или в дискретном времени:

^[S] ==/<o(M[S— X] +{X[S—х]Х, t =

125)

y[s]=qls] + h[s].

Помеха /z[s]— центрированная гауссова с дисперсией о?Л. Функция потерь имеет вид

V^s=(^*[s] —^[s])a.

Плотность вероятности перехода для р. задана выражением

Л И Л /И /-!]):

(|a[z]—pix[z—I})2

 

2aV ('1 -Р‘

 

/ 2 * о -р*) ехрГ

 

 

(;1.

126)

где р=р( А /).

Рассматриваемая система является приводимой и нейтраль­

ной. Оптимальное управляющее

воздействие w*[s] находим

из условия минимума функции as |_T. Получаем

а+1

«*[s] = К ; я* [s+x] —р

(At) — S + т

5+т

где

A p-J[Q]-f-

6

2247

81

S— 1— x

+ 2 (h-U'I—PF-ti—ll)sAi-h

i=\

5—1

+ У^у[Л — Kau[i—x] —A^U’—*])2 ) ^[0]й[л[1]...^|х[5—x —1],

 

 

Д =

 

->

A i— ~r

 

 

 

 

 

 

 

 

 

-,X

 

eV (l-p * )

 

 

 

 

Функция a"s__

отличается

от a's ^_x

лишь

наличием

мно­

жителя (i[s—х—1] под знаком интеграла-

После

вычисле­

ний (подробные выкладки приведены в [В.

16])

находим

“' lSl

Л о ^ +

Р

W hPs- , - \ ,

 

^ -y[s— 1]— Ills—X— 1] +

 

Ло

 

 

 

 

 

+

 

 

 

-J-y[s—2] —w[s—х—2] ]H-

 

 

Ps—X—2

 

Ao

 

 

 

 

(1.

127)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

 

 

 

— y[s-3] —

 

 

 

 

 

 

 

Ps- * - 3 Lv^o

 

 

 

 

 

 

u[s—x— 3] ) + ...+

p -f yJ -yH -H l— ^[1]

 

 

 

где-

 

 

pAi

_ ______ p

 

 

 

(1.

128)

 

 

 

 

 

2 a \

(1—Pa)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P1=

- ^ - ( 1 +

Дх);

 

 

 

 

 

 

 

 

 

 

h

 

 

 

 

 

 

Pj

n i - ( l +

A x + p 'A i)--^ —

при y = 2 ,3 ,...,s -x -2 ;

 

za h

 

 

 

 

1 /-1

 

 

 

 

 

 

P s -

,

-

1

za Л

 

 

( 1 +

A

i ) -

- p —

 

 

 

 

 

 

 

r s —z—2

 

 

 

причем1 справедливо

 

 

 

f

 

 

 

 

неравенство q—< 1 при /= l,...,s —х—1.

Формулу (ll. 127)

можно

записать

в более

компактном

виде:

 

 

 

 

 

 

 

 

 

 

 

 

82-


5ф1 в ? ! | ^ _ Р(,11+д ,) £

(1. 130)

i= t+ l

Весовые коэффициенты g £- удовлетворяют условию

l> ^ - i> g 's - a > .- .> £ 't+ l >0.

■Их значения легко могут быть определены путем сопоставле­ ния выражений (1. 127) и (1. 130). Второй член в формуле (1. 130) представляет собой оптимальную оценку т возму-

•щения в (s+x)-fi момент времени. Это условное математи­

ческое ожидание

p [ s + t ]

при известных y[s—1] и «[s—1].

Как видно из (1.

127) и (1.

130), разности

у[г']—u[i—t] j

входят в выражение с различными весами, более старая ин­ формация деградируется и одновременно накапливается но­ вая. Появление множителя р(тн+ Д t) перед знаком суммы в (1. 130) означает, что в УУ должно быть прогнозирование в статистическом смысле будущих значений р. Действительно, для определения условного математического ожидания p[s+r] при известном p[s—1] необходимо умножить p[s—1] на ко­ эффициент корреляции р(тн + Д 0 ■ Оценка т является мар­ ковской достаточной статистикой возмущения. Другую доста­ точную статистику — апостериорную дисперсию—мы не при­ водим, поскольку она не используется при вычислении опти­

мальных

управляющих

воздействий.

 

 

С учетом сказанного преобразуем алгоритм управления

•{1. 130) и приведем его к марковскому виду:

 

 

 

И*[S] = ^(<7* [s+ т ]- b s-iq* [ s + T - 1 ]) +

 

 

+ 6s_itt[s—11—g s - l

y[s— 1] —u[s—T—1]

(1.

131)

 

 

 

Ко

 

 

где

 

£>s—i= P

2° \ P S

(1.

132)

 

 

 

 

 

1:

Ps

1 - -

P4 Д<

(1.

133)

 

 

:1+р2Д 1+ Д 1-

 

 

 

 

 

 

P s - , - 3

83


Переменные коэффициенты bs- 1, g s-i могут быть заранее вычислены и введены в память управляющей вычислительной машины.

Закон управления (1. 130) или (1. 131) получен в пред­ положении, что р — центрированный процесс. Если на объ­ ект действует возмущение с отличным от нуля средним значе­ нием X, то его можно представить в виде суммы ^+ p [s], где ц[5] по-прежнему центрировано. При известном априори К алгоритм остается прежним (1. 131), но непосредственно на объект вместо m*[s] следует подавать управляющие воздейст­ вия и* [s] = u*[s]—X. Если же параметр X не известен заранее и является случайной величиной, распределенной по нормаль­ ному закону

то описанная выше общая методика синтеза алгоритма УУ применима с непринципиальными изменениями. Закон уп­ равления линейный

(1. 134)

:+1

но весовые коэффициенты Р/ отличаются от полученных выше

(Р £/)•

Перейдем к рассмотрению непрерывной системы в ста­

ционарном режиме. При больших s

и малых Л/(Д/->0) имеем

a*[s]=a*[sA /] = —- q*(s А / + т Д /)—р(тн+ Д ()MiX

 

Л"

 

(1.

135)

Aj) —u (s A t —- A l —jA t)

Устремляя At к нулю, получим

 

 

 

t

 

 

 

_Z,vW _v)

—u(t—тн—v)

dv,

к'

о

А0

 

 

 

 

 

v = /A /,

 

(1.

136)

84


Рис. 1.8

где

/ - й , г +л

(1.

137)

 

;

 

 

ai°V

(1.

138)

 

1 «S^L + aJ

 

 

 

Блок-схема оптимального управляющего устройства

(УУ)

для i7* [s]= 7 * приведена на рис. 1. 8.

УУ включает модель М

объекта с передаточной функцией Кое

~р-.ч

и вычислительное

устройство, выполняющее операцию с выходными сигналами объекта и модели: усреднение — фильтром Ф и статистический

прогноз на время тн — блоком S2 (усилитель с коэффициен- —ai|"« I

том усиления^=р("н)=(? ). Фильтр Ф представляет со­ бой инерционное звено с постоянной времени Тф и коэффи­

циентом

усиления

Si

 

 

 

 

 

 

_1

Д1°У

2

 

 

 

 

 

(1.

139)

 

 

L

а 2!

5

 

 

 

 

 

 

М х

 

ацау

 

 

 

L

(

а

, /

Sh) V a ^

+iza1%Sh

 

 

 

 

 

 

(I.

НО)

85