Непараметрическое моделирование и управление многосвязными системами Текст научной статьи по специальности «Математика»

УДК 62-501

А. В. Медведев, А. Р. Низамеев

НЕПАРАМЕТРИЧЕСКОЕ МОДЕЛИРОВАНИЕ И УПРАВЛЕНИЕ МНОГОСВЯЗНЫМИ СИСТЕМАМИ

Приведена схема управления многосвязным технологическим процессом. Даны общая и математическая постановки задачи управления многосвязным процессом. Построены непараметрические модели идентификации входных и выходных значений управляемого процесса. Приведены некоторые результаты численных исследований построенных моделей.

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

Общая постановка задачи. Приведем схему управления процессом (рис. 1).

Рис. 1

Примем следующие обозначения: О - объект; АС -адаптивная система; ut - управляющее воздействие; ц? - контролируемое неуправляемое воздействие; ^ - случайное воздействие; х,, 2, - выходные переменные; Ни, кц , кх, кг - помехи в каналах измерения; иД цк, хД 2^ - измеренные значения соответствующих переменных; х* и 2*—заданные значения выходных переменных.

Следует заметить, что блок АС представляет собой достаточно сложную структуру, детализацию которой в данной статье приводить не будем, лишь укажем, что в ее состав входят как модель исследуемого процесса, так и соответствующая процессу иерархия блоков управления.

Целью системы управления является поддержание заданного значения х* и 2,*. Отметим существенное значение и отличие выходных переменных xtи 2,: переменная xtконтролируется через достаточно малые интервалы времени А?, как и переменные и,, ц?, а переменная 2, - через существенно большие интервалы времени АТ (АТ >> Аí).

С технологической точки зрения, для всего технологического процесса наиболее важным является контроль переменной 2,. Например, выходная переменная х, контролируется с помощью различного рода индукционных, емкостных и других датчиков, а 2, измеряется в результате химического анализа, физико-механических испытаний и др. Этим и обусловлено существенное отличие диск-

ретности контроля выходных переменных х, и 2,. Если дискретность измерения х, иц? - это секунды и минуты, то дискретность 2составляет смену, сутки, недели и более. Последнее обусловлено технологией проведения самого контроля, который обычно регламентируется государственными стандартами. Приведем сравнительно простую схему технологического процесса (рис. 2).

Рис. 2

На этой схеме опущены помехи, действующие в каналах измерения, чтобы не загромождать рисунок. Но предполагается, что они, конечно же, действуют (см. рис. 1). Содержание переменных соответствует описанному выше. Отметим лишь, что роль неуправляемых переменных в этой схеме играют как ц?, так и 2,1, 2,2, контроль которых производится через различные, но значительные интервалы времени. Блоки Р^ Р2 и Р3 - это регуляторы соответствующих объектов 01, О2 и О3,1 - вся имеющаяся информация о процессе.

Математическая постановка задачи. Пусть

ц = { ц/,..-, ц *} е як, и= {и,1,..., и™}е Я™, х, = {хД ...,

х"} е Я", 2, = {2Д ., 2Д е Я1. Характеристики объектов О!, О2 и О3 взаимооднозначны по вектору управляющих переменных, т. е. одному значению и, соответствует только одно значение ц?. Имеется обучающая выборка { ц?, и,, х, 21-Т, ? = \, ^ }. Здесь т указывает на запаздывание (смена, сутки и т. д.).

Сформулируем критерий оптимальности:

К(2) = М,ЛМ;(2 - 2)1 \ , ц} = тт, (1)

2

где 2 - оценка вектора выхода 2.

Используя необходимое условие минимума, т. е. (приравняв производную функции Я по искомой величине 2 к нулю), получим

2орр = М(2 \ и,ц}. (2)

Непараметрическая оценка выходной переменной 2{ для системы, представленной на рис. 1, имеет вид

z*+, =

Xz* Пф

i=1 j=1

«s+1j - «/ n П ф íцp+l -цp 1

Cj Cp

V ' { ц s

ХПФ

i=1 j=1

*s+1

C]

ц s

\ /

Пф

p=1

Ґцp+1 -цр Л

, (3)

R = X (*,■- (C))2 ^ min,

' 7*f zf 1 Z s+1 Li іПф |^p+1 -цp 1

Cf Cp

V У 1 ц s J

ХПФ

і=1 f=1

* f f * s+1-

Cf.

Пф

p=1

"цР+1 -цр4 ср.

S т

Х--“Пф

i=1 j=1

«1Sj+1- «lj

Ci

Пф

'ц*+1 -ц*4

C*1s

ХПф

i=1 j=1

M‘S+1 - Mlj

X *2*Пф

I=1 j=1

Cj1

u s

* = 1, N,

N

П ф

Пф

ц*+1 -ц*

C*

^1s J

(9)

Гіф fp=1 | 1p 1p 1 Xp+1- Xp

cp xs J

ХПф

i=1 j=1

Cj2

Пф

Пф

p=1

Ci1

X *3*Пф

i=1 j=1

CJ,

* = 1, Q,

Пф

(10)

Пф

p=1

ср_

ср

Ц і

-і

где ^ - размер обучающей выборки; іі+1 - оценка ^й компоненты вектора выхода, і = 1,1; - вектор входа

из экзаменующей выборки; Ф - ядро функции, выбирается по следующим условиям: | Ф(Ъ)аЪ = 1,

Ф(Ь1) <Ф(Ъ2) для | Ъ2|<|Ъ1|, | ЪФф)Л , здесь

р = 2, 3,..., Ь - аргумент функции Ф; Си - параметр размытости ядра Ф, удовлетворяет следующим условиям:

С ^ 0 , пС (п) ^ ~ •

и п^ п^~

При выполнении перечисленных выше условий оценка і является асимптотически несмещенной, состоятельной, асимптотически нормально распределенной.

Параметр размытости С выбирается по условию х*+1 выполнения минимума критерия

\2

ХПф

і=1 j=1

Пф

Пф

p=1

* = 1, A . (11)

Непараметрические оценки выходных переменных x t', xД x^записываются следующим образом:

X15 =

s+1

S т

Хх?Пф

i=1 j=1

«4+1 -

Cj

Пф

'ц *+1 -ц*4 Cis

ХПф

і=1 j=1

«І1 - ulj

CK

У i«s J

* = Гн,

Пф

'ц*+1 -ц* Л

C*

ц s

\ р /

(12)

X х“Пф

i=1 j=1

Cj2

Пф

Пф

p=1

Ci1

ХПф

i=1 j=1

Пф

C‘,

(5)

* = 1, F,

где s - размер обучающей выборки

Критерий оптимальности для оценки управляющего воздействия выглядит так:

R(u) = M {Mu (и - и)2 | д, z} = min, (6)

и

отсюда u°pt = M(и | д, z = z*}. (7)

Формула оценки j-го компонента вектора и в момент времени s + 1 записывается следующим образом:

X « Пф

i=1 f=1

Хх3*Пф

м* i=1 j=1

X„ + =-

Пф

2* 2*

Пф

=1

Q

Пф

cp

(13)

s Q

X П ф

I =1 j=1

Пф

2* 2 *

Пф

* = 1, T . (14)

Непараметрические оценки выходных переменных x t', xД x^записываются следующим образом:

X «1*Пф

i=1 j=1

Пф

ц*+1 -ц* с*,

Пф

p=1

сі1

ХПф

i=1 j=1

Cj.

Пф

b b ц+1 -ц*

Пф

p=1

CJ1

] = 1, т. (8)

Рассмотрим теперь схему (см. рис. 2). Пусть Ж, О , Л -размерности векторов гД zí2 и zí3, Ж, О и Л - размерности векторов ', г2 и г(3; Н, ^ и Т- размерности векторов ', хД х(3; ^, Г и Q - размерности векторов иД и(2 и и(3; D - размерность вектора .

Запишем модели прогноза выходных значении ', г(2, г (3 и х 1', х2, хД а также модели расчета управляющих воздействий и', иД и(3 для объектов 01, 02 и О3.

Непараметрические оценки выходных переменных г 1', г2, г(3 будут выглядеть следующим образом:

* = 1, т,

(15)

X «2*Пф

-2* i=1 j=1

| _2 j" 2j \ | ,,2 *" ,,2 * 1

z j+1- z і П ф X i+1 X і

C\ C*2

У z s x xs J

s Q

^ Пф

i=1 j=1

z2j* - z2j

z s+1 z i

Cj

z s

* = 1r,

X «3*Пф

-35 i=1 j=1

Cj3

Пф

1

I T ft ф

x2*" - x2* '

s+1 л I

с*

(16)

C*

s A

X П ф

i=1 j=1

z3j" - z3j

Z s+1 Z і

Cj

Пф

/ " - x3* 4 A s+1 A і

с*

x А

V x s У

5 = 1, Q • (17)

Численные исследования непараметрических моделей. Приведем пример численных исследований непараметрических моделей прогноза выходных значений объектов последовательного процесса, блок-схема которого представлена на рис. 3. Там применяются следующие обозначения: Ор 02, Оп- объекты; и, Ж, V- векторы входов; Z - вектор выхода; К - элемент запаздывания; X = ДЦ) + ^; 7 = ДХ, Ж) + ^; Z = /(У, V) + ^; h1, h3 -

помехи.

i=1

На вход первого объекта подается вектор и = {и ип}.

Выход первого объекта - векторX = {х хт} поступает

на вход объекта 02. Также на вход 02 поступает вектор Ж = {■м1, ..., wk}. Векторы Y = {у1,.,ур} и V = ..., уД

поступают на вход объекта 03. Выходом третьего объекта является скаляр Z. Ключ К указывает на то, что истинная величина выхода объекта будет известна только через определенный промежуток времени.

“ Ч - - т-------[ ІГ 1-^—

Рис. 3

=х Xе пф

¡=1 }=\

Сх

' !Пф

і =1 і =1

Схк

, (19)

где у,+1 = {у,+і1>->у,+і ,^+11>->^}; d=2;/ = 3 (размер

вектора V); г = d +/.

В модели второго вида процесс рассматривается как единый объект. Входами являются векторы и, Ж, V, выходом - вектор Z. Модель для данной постановки выглядит так:

2 = / (и, Ж ,У). (22)

Компоненты вектора Z находятся по формуле

і=1

і .

где +1 = {и,

Сх

5

/ 1Пф

і=1 і=1

... 4

Сх

, (23)

, +13}; г=12.

Пусть имеется исходная выборка и, Wi, V,, X., Y., Z, I = 1, ^ , содержащая информацию о входах и выходах всех объектов. Известно, что при измерении входных и выходных величин в каналах измерения действуют помехи с нулевым математическим ожиданием и ограниченной дисперсией.

На основании этой выборки требуется построить модель, которая по известному значению входа и первого объекта, а также известным значениям W и V могла давать прогноз выхода всего процесса.

Будем считать, что объекты относятся к классу статических с запаздыванием, а истинные зависимости - к классу непрерывных дифференцируемых функций. Также предположим, что связь выходов объектов с их входами является однозначной.

Дополнительно отметим, что в численных исследованиях для простоты исходная выборка входных-выходных значений взята не из реального процесса, а из процесса, сымитированного при помощи средств программирования.

Численные исследования представлены для следующих размерностей входов и выходов объектов: и = {иР ...,и5};х = X ...,х3}; W = ..., w4}; Y = {y1,у2};

V = ^х, ..., v3}; Z- скаляр.

В зависимости от требований реальной задачи возможно построение двух видов модели процесса.

В модели первого вида весь процесс рассматривается как последовательность объектов. Требуется вычислить выходы на каждом из объектов.

Для данной постановки модель процесса имеет вид 2 = / (7 (X (V ),Ж ),¥). (18)

Запишем формулу оценки выхода первого объекта:

В исследованиях для всех моделей использовались колоколообразные ядра вида

, , [0,335-0,067(В)2,і^В)2 < 5,

ф(в )=-! ’ 2 (24)

[ 0, if (В) > 5.

Приведем некоторые результаты численных исследований построенных непараметрических моделей.

Пусть истинные зависимости имеют линейный вид: х1 = 1,6 и1 + 0,3 и3; х2 = и2 + 0,7 и4; х3 = 0,5 и5;

у1 = 0,5 х3 + 0,8 w1 + 0,8 w3; у2 = 0,6 х1 + 0,8 х2 + 0,6 w4 +

+ 0^3; г = 0,7 у1 + 0,6 у2 + v1 + v2 + v3.

Ниже приведены результаты расчета выхода х модели (18) (ряд 2) в сравнении с результатами работы модели (22), (ряд 3) и реальными значениями Z (ряд 1) (рис. 4...9).

Объем обучающей выборки 5 = 200; и выбирались случайным образом из [0; 5]; Ж и V - аналогично из [0; 1]. Экзаменующая выборка 5 = 50; h = 3 %; среднеквадратичные ошибки прогнозирования модели (18) є = 0,38, модели (22) є = 0,37 (рис. 4).

где х5+1к - оценка ^й компоненты вектора выхода, к = 1,3 ; I = 5; 5 - размер обучающей выборки; и++1 - вектор входа из экзаменующей выборки.

Для второго объекта модель записывается следующим образом:

у.=±у; П ф

,■=1 і=1

где п = 1,2 ; +1 = {х,+1 ,...,х,+1

(размер вектора Ж); h = е + t.

Для третьего объекта модель строится аналогично:

Рис. 4

Часть зависимостей нелинейная.

Пусть истинные зависимости имеют вид х1 = sin(5м1) / (12+м1) + .6м1; х2 = 2+1^т(-.2и2)+.8и3;

х3=23т(0.3и4)+13т(.3и5)+ .6м1; у1 = w1exp(.7x3) + 2w2) / 5;

у2 =.6х1 + 1.5х2 + 1.7w2 + 2.5 w4; г = .1.3sin(.2y1) +.7у2 +

+ 1.5v1 + 1.3 v2 + v3.

Объем обучающей выборки 5 = 200; и выбирались случайным образом из [0;5]; W и V -аналогично из [0;1]. Экзаменующая выборка 5 = 50; помехи h = 3 %; среднеквадратичные ошибки прогнозирования для модели (18) £ = 1,27, для модели (22) е = 1,03 (рис. 5).

=ЁхПф

¡=1 і=1

у,+1 - у і

Сх

/ !Пф

¡=1 і=1

Сх

Рис. 5

Обучающая выборка содержит всего 20 точек; hm¡x= 3 %; среднеквадратичные ошибки прогнозирования для модели (18) е = 2,49, для модели (22) е = 1,87 (рис. 6).

Объем обучающей выборки 200 точек; и выбирались из [0;5], W - из [0;3] и V - из [0;4], = 3 %; среднеквад-

ратичные ошибки прогнозирования для модели (18) £ = 8,50, для модели (22) е = 3,20 (рис. 7).

График зависимости квадратичного критерия R от С для модели (18) (рис. 8) показывает, что С =1,28 - это наилучший параметр, определенный с помощью данного критерия ^(1,28)=1,27). На самом же деле наилучший параметр лежит левее найденного ^(0,93)=1,03). Данный факт объясняется тем, что мы имеем дело со стохастическими данными.

На основании графиков (см. рис. 4.9) можно сделать вывод, что модель (22), в которой процесс рассматривается как единый объект, в среднем дает лучший прогноз выхода всего процесса, чем модель (18), в которой процесс рассматривается как последовательность объектов. А если влияние входов W и V значительно, то модель (18) работает плохо.

Анализ численных исследований. Для исследований моделировались обучающие выборки малых объемов, при наличии помех и «дыр», что максимально приближено к реальным условиям.

В результате сравнительного анализа, как мы уже отмечали выше, было выяснено, что модель (22) дает более объективный прогноз выхода всего процесса, чем модель (18). Это можно объяснить тем, что во втором случае на вход модели, кроме измеренных величин, подаются оцененные величины. Во-первых, повышается размерность задачи, а во-вторых, поскольку оцененные величины являются не точными, то они искажают конечную оценку.

Кроме того, численные исследования модели последовательного процесса (18) показали, что если влияние неуправляемых входов значительно, то данная модель работает плохо. Причина заключается в выборе параметра размытости, поскольку для модели (18) требуется более сложный механизм настройки параметров размытости, чем для модели (22).

Таким образом, произведенные исследования показали, что непараметрические модели последовательных многомерных процессов дают хорошие прогнозы при обучении на достаточно малых обучающих выборках. Предложенные в статье обучающиеся непараметрические алгоритмы управления могут быть модифицированы и использованы при разработке реальных компьютерных систем управления многосвязными процессами.

Возникает случай, когда в некоторых областях наблюдения отсутствуют обучающие точки, т. е. случай с «дырами» в обучающей выборке: ио6уч е (0; 1) и (2; 4) и (4,5; 5), Ш е (0; 1), V е (0; 1), hmx¡=І %; среднеквадратичные ошибки прогнозирования для модели (18) е = 1,71, для модели (22) е = 1,31 (рис. 9).

Библиографический список

1. Медведев, А. В. Актуальные проблемы информатики, прикладной математики и механики. Ч. 3. Информатика / А. В. Медведев. Новосибирск ; Красноярск : Изд-во СО РАН, 1996. 170 с.

2. Medvedev, A. V. Nonparametric theory of control system design / A. V. Medvedev // Вест. Сиб. гос. аэрокос-мич. ун-та им. акад. М. Ф. Решетнева / под ред. проф. Г. П. Белякова ; Сиб. гос. аэрокосмич. ун-т. Вып. 3. Красноярск, 2002. С. 44-55.

Рис. 9

A. V. Medvedev, A. R. Nizameev

NONPARAMETRIC MODELING AND CONTROLING OF MULTIPLY CONNECTED SYSTEMS

Offered controlling scheme of multiply connected technological process. Given general and mathematical statements of multiply connected processes controlling problem. Constructed nonparametric identification models of input and output values of the process.