Научная статья на тему 'НЕПАРАМЕТРИЧЕСКИЕ АЛГОРИТМЫ ИДЕНТИФИКАЦИИ И УПРАВЛЕНИЯ МНОГОМЕРНЫМИ БЕЗЫНЕРЦИОННЫМИ ПРОЦЕССАМИ'

НЕПАРАМЕТРИЧЕСКИЕ АЛГОРИТМЫ ИДЕНТИФИКАЦИИ И УПРАВЛЕНИЯ МНОГОМЕРНЫМИ БЕЗЫНЕРЦИОННЫМИ ПРОЦЕССАМИ Текст научной статьи по специальности «Математика»

CC BY
64
11
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
НЕПАРАМЕТРИЧЕСКОЕ МОДЕЛИРОВАНИЕ / УПРАВЛЕНИЕ / БЕЗЫНЕРЦИОННЫЕ СИСТЕМЫ / Т-ПРОЦЕССЫ / Т-МОДЕЛИ / ЦЕПОЧКА АЛГОРИТМОВ / NONPARAMETRIC MODELING / CONTROL / INERTIALESS SYSTEMS / T-PROCESSES / T-MODELS / ALGORITHM CHAIN

Аннотация научной статьи по математике, автор научной работы — Медведев Александр Васильевич, Ярещенко Дарья Игоревна

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

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

NONPARAMETRIC IDENTIFICATION AND CONTROL ALGORITHMS FOR MULTIDIMENSIONAL INERTIALESS PROCESSES

Problems of identification and control of multidimensional inertialess systems with delay in conditions of nonparametric uncertainty are considered, i.e. in conditions when the form of parametric equations for various channels of the object is absent due to lack of a priori information.

Текст научной работы на тему «НЕПАРАМЕТРИЧЕСКИЕ АЛГОРИТМЫ ИДЕНТИФИКАЦИИ И УПРАВЛЕНИЯ МНОГОМЕРНЫМИ БЕЗЫНЕРЦИОННЫМИ ПРОЦЕССАМИ»

ВЕСТНИК ТОМСКОГО ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА

2020 Управление, вычислительная техника и информатика № 53

УДК 519.876.5

DOI: 10.17223/19988605/53/7

А.В. Медведев, Д.И. Ярещенко

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

Рассматриваются задачи идентификации и управления многомерными безынерционными системами с запаздыванием в условиях непараметрической неопределенности, т.е. в условиях, когда вид параметрических уравнений по различным каналам объекта отсутствует из-за недостатка априорной информации. Ключевые слова: непараметрическое моделирование; управление; безынерционные системы; Т-процессы; Т-модели; цепочка алгоритмов.

Идентификация и управление многомерными безынерционными процессами (объектами) с запаздыванием продолжают оставаться довольно актуальными проблемами в настоящее время. При этом следует учитывать, что по различным каналам многомерного объекта процессы чаще всего могут быть и динамическими, но контроль переменных осуществляется через дискретные интервалы времени. В качестве примера можно привести самые различные технологические, производственные и активные процессы, в частности в стройиндустрии при сухом измельчении клинкера (клинкер - это продукт обжига), который в последующем подлежит измельчению, а измельчение приводит к получению цемента. Основным параметром, определяющим его качество, является активность цемента, т.е. его прочность при сжатии. Но важными технологическими показателями являются также тонкость помола, удельная поверхность, расплыв конуса и др. Их измерение осуществляется через несколько часов: два, три и более, в то время как постоянная времени по различным каналам многомерного объекта составляет 3-5 мин, следовательно, переходный процесс длится в течение 20-25 мин. Это обстоятельство приводит к тому, что мы вынуждены рассматривать те или иные каналы как статические с запаздыванием, а также к зависимости выходных переменных; в частности, активность цемента зависит от тонкости помола и др. В этом случае объект описывается в виде некоторых неявных функций. Таким образом, можно сделать вывод, что многомерный объект в общем виде описывается в виде системы неявных функций. Особенностью настоящей задачи, о которой пойдет речь в статье, является то, что вид такой системы функций априори оказывается неизвестным. Тем не менее исследователь стоит перед необходимостью решить эту систему функций относительно компонент выходных переменных x(t) объекта при известных входных u(t). Подобные многомерные процессы были названы Т-процессами, а объекты, соответственно, Т-объектами [1].

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

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

1. Т-процессы

Т-процесс - это многомерный дискретно-непрерывный процесс, имеющий стохастическую зависимость выходных переменных х(V) = (х1 (V),х2 (V)хп (V)), 7 = 1,п [1]. Обозначим вектор входных

компонент - и (V) = (м1 (V),и2 (V),..,ит (V)), £ = 1,т . Это приводит к тому, что математическое описание

объекта представляется в виде некоторой системы неявных функций вида ^ (и (V),х (V)) = 0, 7 = 1,п . Задача идентификации сводится к задаче решения системы нелинейных уравнений

(и (г),х(г)) = 0, 7=1п, (1)

относительно компонент вектора х (V) = (х1 (V),х2 (V),...,хп (V)), при известных значениях и (V ) = (u](t), и2 (V),..., ит ()). В общем виде исследуемая многомерная система, реализующая Т-процесс, может быть представлена следующим образом (рис. 1).

Рис. 1. Многомерный объект Fig. 1. Multidimensional object

На рис. 1 приняты следующие обозначения: U(t) = (u (t),u2 (t),...,um (t)) - га-мерный вектор

входных переменных; x(t) = (x1 (t),x2 (t),..., xn (t)) - «-мерный вектор выходных переменных, который напоминает мультиколлинеарный вектор с той лишь разницей, что компоненты вектора x(t) могут зависеть и нелинейно; t) - случайные помехи, действующие на объект; вертикальные стрелки показывают стохастическую зависимость выходных переменных. По различным каналам исследуемого объекта зависимость j-й компоненты вектора х может быть представлена в виде некоторой зависимости от тех или иных компонент вектора u: x<j > = f (u<j>) j = 1, n. Каждый j-й канал зависит от нескольких компонент вектора u, например x<5> =(u,U,U), где x<5> - составной вектор. Составной вектор - это вектор, составленный из компонент векторов входных и выходных переменных [3]. В этом случае система уравнений (1) примет вид:

Fj(u<j>(t),X<j>(t))= 0, j = 1n, (2)

где U<j>(t), X<j>(t) - составные векторы. В общем виде зависимость той или иной компоненты вектора выхода может быть от всех компонент вектора входа. Одним словом, составные векторы выписываются исследователем на основании имеющейся априорной информации.

Заметим, что вид уравнений Fj (•), j = 1, n продолжает оставаться неизвестным и не может интерпретироваться как модель исследуемого процесса. Задача состоит в моделировании подобных процессов, т.е. Т-процессов.

2. Т-модели

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

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

*>,) = <>> 7 = 1Я (3)

где ХЛ.. г/л - временные векторы (набор данных, поступивший к 5-му моменту времени), но и в этом случае Р-(•), , = 1, п продолжают оставаться неизвестными. В теории идентификации подобные задачи не только не рассматриваются, но и не ставятся. Обычно идут по пути выбора параметрической структуры (2), но, к сожалению, преодоление этого этапа затруднено из-за недостатка априорной информации. И требуется длительное время для определения параметрической структуры, т.е. представления модели в виде:

¥} (Г^), Г^), а) = 0, , = Щ, (4)

где а - вектор параметров. Далее следуют процедура оценки параметров по элементам обучающей выборки , х^, г = 1,5 , и последующее решение системы нелинейных взаимосвязанных соотношений (4). Успех построения модели в данном случае будет зависеть от качественной параметризации системы (4).

Рассмотрим задачу построения Т-моделей в условиях непараметрической неопределенности, т.е. в условиях, когда система (4) неизвестна с точностью до параметров. Таким образом, задача моделирования Т-процессов сводится к прогнозу значений выходных переменных X (г ) = X (г), х2 (г),..., хп (г)) при известных входных и (г )= ),и2 (г),..., ит (г)).

В результате измерений, проводимых на объекте, формируется обучающая выборка хг, иг, г = 1,5 . В этом случае оценка компонент вектора выходных переменных х(?) = (х(?), х2(?),..., хи (?)) при известных значениях входных и()=(их(г),и2(г),...,ит(г)), как уже было отмечено выше, приводит к необходимости решать систему уравнений (3).

В итоге задача идентификации сводится к тому, что при заданном значении вектора входных переменных и = и', необходимо решить систему (3) относительно вектора выходных переменных х (г ) = (хх (г) , х 2 (г),..., хп (г)).

По этому поводу следует сделать некоторые специальные замечания, которые и будут содержать сущность отличия проблематики идентификации Т-процессов от общепринятой схемы идентификации [3, 4]. Она состоит в том, что алгоритм идентификации не может быть представлен так, как обычно он принят в теории идентификации:

х (г ) = ^ (и (г), а) (5)

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

еи 7=1,«, (6)

где в у (й<]> ,х<]> | находятся с использованием непараметрической оценки функции регрес-

сии Надарая-Ватсона [5]:

^ < т>( ик - и, тЛ ^ х, [г ] П Ф к - к

в, (г) = Ф, (и<*>, х, (г)) = х, (г) - --к=' ^ С"* ^, , = 1п, (7)

5 <т> и'к- ик[г]

X П Ф

г =1 к=1

С.Шк

где <т> - размерность составного вектора щ, < т > < т, в дальнейшем это обозначение используется и для других переменных.

Колоколообразные функции Ф ющим условиям [5]:

( ' г-Л

ик- ик [7]

и параметр размытости си удовлетворяют следу-

ик

Ф (•)<»;

I Ф (с-1 (ик - иы ))ш =1;

П(и)

с-1ф (с-1 (ик - ик7)) = 5 (% - ик7);

11т ^ С = 0;

11т *с7 =

Следующий шаг состоит в оценивании условного математического ожидания:

xj = М[х] |и<7>,8 = о}, 7 = 1,п .

В качестве оценки (13) примем непараметрическую оценку регрессии [Там же]:

*7 =

I ^ [7] •П Ф

7=1 к1 =1

<т> К - иК [7] [7] ^

П Ф

V яи У к2 =

4 <т> (иК - иК [7] [7] ^

7 = 1, п

1П Ф

%=1 к =1

V яи У

П Ф

к2 =1

V С4В у

где колоколообразные функции Ф (•) примем в виде треугольного ядра:

- ик1[7]Л

Ф

V яи у

1 К - ик1[7]\ Iик1 - ик1[7]\ Л

0 К- ик1[7]| ^

(8) (9)

(10)

(11) (12)

(13)

(14)

(15)

Ф

V С4В У

1 |0 - Вк2 [7 ]| |0 -Вк2[7]\ ^

о, , 1.

(16)

Осуществляя эту процедуру, получаем значения выходных переменных х при входных воздействиях на объект и = и', а в этом и состоит основное назначение искомой модели, которая в дальнейшем мжет быть использована в различных системах управления [6], в том числе в организационных [7, 8].

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

Ех. - х (и.)

| 7 4 V 7 /|

5= -7=-

7

7 = 1, п,

(17)

где х7 - наблюдения на объекте, х4 (и7) - прогноз выхода объекта, х - среднее значение по каждой

компоненте вектора х.

4

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

С

яи

и

С

¿=1

Поясним механизм работы предложенной цепочки (7) и (14), составляющей алгоритм работы идентификации. Учитывая локальные свойства непараметрических оценок функции регрессии, можно увидеть и проанализировать технологию функционирования цепочки (7) и (14), которая и дает прогноз каждой компоненты вектора выхода х при известных значениях вектора входа ы. Действительно, если значения вектора входных переменных ы известны, то они (7), (14) тем самым локализуют некоторую подобласть в пространстве входных-выходных переменных (х, и). Далее, при оценивании (прогнозе) каждой компоненты вектора х в соответствии с (14), в силу локальности непараметрических оценок, выделяется та подобласть, для которых соответствующие е^- близки к нулю.

При других значениях входных переменных ы подобная ситуация повторяется.

3. Управление дискретно-непрерывным процессом

Рассмотрим задачу управления многомерным Т-объектом в условиях непараметрической неопределенности. При этом приведем следующую схему (рис. 2).

На рис. 2 введены следующие обозначения: и(г) = (щ (г), и2 (г),..., «ш (г)) - входные управляемые

переменные процесса; ) = (^(г),ц2(г),...,^р(г)) - входные неуправляемые, но контролируемые переменные процесса; х(г ) = (х (г), х2 (г),..., хп (г)) - выходные переменные процесса; х* = (х],...,х*)ео(х*)с Я" - задающее воздействие, £(г) - случайные стационарные помехи, действующие на объект.

Рис. 2. Схема непараметрической системы управления безынерционным объектом Fig. 2. Diagram of a nonparametric control system of an inertialess object

Переменная ^,(t) является входной контролируемой, но неуправляемой переменной, в частности это могут быть технологические параметры, для измерения которых используется физико-химическая технология. Например, клинкер может быть недообожженным или переобожженным. Такая переменная существенно влияет на объект и на значения выходных переменных x(t). Поэтому задача управления таким объектом усложняется, так как необходимо поддерживать на выходе объекта заданное значение x*(t) при известном значении ^(t).

Управление дискретно-непрерывным процессом рассматривается в условиях непараметрической неопределенности, т.е. в условиях, когда модель процесса с точностью до вектора параметров отсутствует полностью. В этом случае известные приемы неприменимы [4], и следует использовать другие подходы для решения задачи [9].

В случае отсутствия достаточной априорной информации об исследуемом объекте целесообразно использовать непараметрический алгоритм управления изложенный в [1]:

к 1=1

Z ик П Ф

j=i

!П ф

1=1 /=1

г,

- xj ) p П ф Jv=1 ( * i\

'x V c^ J

x1) xj p П Ф v=1 1 * i \

J c, J

к = 1, m,

(18)

где (и1, ц1,х , I = 1, ^) - обучающая выборка.

Однако в данной задаче естественно использовать следующую цепочку: входную переменную п\(7) берем произвольно из области О(м). Входная переменная И2(0 может быть определена в соответствии со следующим алгоритмом:

( Л (

1щ2Ф и и1- пФ

1 V си ) 1—1

щ — -

V Си1 J

\ г

п

п Ф

1—1

* 1 К* - х1 \ п ф ( * Л к-к

сх х1 V J

1 Л - х п ф У=1 V * 1 Л

с х1 J

(19)

Для входной переменной ыз(^ алгоритм управления будет выглядеть следующим образом:

г 'Л (

^^ ф и'

I и?Ф

V Си1 J

V Си2 J

1=1

I Ф

,. Л

и - и

Ф

V и1 J

, л

и^ щ

V Си2 J

п Ф

1=1

( * г х*- х \ пФ ( * 1 Л к-к

сх V 1 V J

* 1Л х* - х п ф У=1 V * Л

сх х1 J

(20)

И так далее для каждой компоненты входа нДО объекта. В общем виде для многомерной системы алгоритм управления будет выглядеть следующим образом:

ик =

* к -1

Ы п ф

1=1 к =1

II п Ф

и, - и.

ик Л п П Ф } 1=1 * 1 х* -х- Л * П Ф ^ ( * - Л к - к

сх V х1 с V J

1 Л (к п ( п ф 1=1 V * 1 ^ хз - * ( П Ф У=1 V *1 К "К

J J J

к = 1, га

(21)

В реальных задачах часто число компонент вектора и больше числа компонент вектора х. Если же размерность вектора и превышает размерность вектора х, т.е. т > п, то обычно поступают следующим образом: в число компонент вектора ц включают компоненты вектора и, с тем чтобы размерность векторов и и х сделать одинаковой [10].

Настраиваемыми параметрами будут параметры размытости с^ , сх и сц , для них будем использовать следующие формулы: с% = а\ик -, сх =р|х* -х'.| и с^ = у|ц*, -ц1|, где а, в и у - некоторые параметры, большие 1: а > 1, в > 1, у > 1. Следует заметить, что выбор см^, сх и сц осуществляется на каждом такте управления. При этом если сначала определен сщ , то определение сх и сц

осуществляется с учетом этого факта. Однако может быть и наоборот, например сначала определяется г или с , а потом остальные.

1=1

1=1

щ =

1=1

1=1 к=1

4. Экспериментальная часть

Для проведения вычислительного эксперимента был взят объект с пятью входными переменными и(?)=(и1 (?),и2(?),и3(?),щ(?),и5(?)), принимающими случайные значения в интервале и(?)е[0, 3], и тремя выходными переменными х(?) = (х ), X ), хз )), принимающими значения в интервалах: х1 ({)е[-2; 15], х2 ^)е[1; 10], х3 ({)е [6; 23]. Для данного объекта была сформирована выборка входных и выходных переменных и найдены прогнозные значения выходных переменных при известных входных. Для вычисления использовались объем выборки 5 = 2 000, параметр размытости с = 0,3, помеха, действующая на компоненты вектора выходных переменных, ^ = 0,07 . Описание объекта

с точностью до параметров было принято только для проведения компьютерного исследования и оставалось неизвестным для изложенной выше теории.

На рис. 3 и 4 по оси абсцисс представлены такты I, по оси ординат - значения выхода объекта и модели. На рисунках представлены 20 точек выборки из-за простоты представления результатов, т.е. каждая сотая точка выборки. «Точкой» обозначены значения выхода объекта, а «крестиком» - значения выхода модели. Как видно из рисю 3 и 4, прогноз значений выходных переменных многомерного объекта по известным входным переменным достаточно удовлетворителен.

Рис. 3. Прогноз выходной переменной xi(t) объекта, измеренный с равномерной помехой 7% Fig. 3. Forecast by the output variable xi(t) of the object, measured with uniform interference of 7%

Рис. 4. Прогноз выходной переменной xi(t) объекта, измеренный с равномерной помехой 7% Fig. 4. Forecast by the output variable X2(t) of the object, measured with uniform interference of 7%

Далее приводятся результаты вычислительных экспериментов для данного объекта при использовании алгоритма управления (21). В проведенном вычислительном эксперименте число компонент вектора u больше числа компонент вектора х. Если же размерность вектора u превышает размерность вектора х, т.е. m > n, то заменим и4 (t) = | (t), а щ (t) = ц2 (t), чтобы размерность векторов u и х сделать одинаковой. Так как входные переменные u(t) принимали случайные значения в интервале u(t)е[0, 3], то | (t) и ц2 (t) также принимают случайные значения в заданном интервале.

Обратим еще раз внимание на то, что исследователю неизвестен вид системы уравнений, описывающих управляемый объект. В качестве информации о последнем используются измерения входных и выходных переменных (u1, |1, x1, i = 1,5) .

Как видно из рис. 5, при управлении объектом выход объекта xi(t) близок к задающему воздействию x** (t) .

Далее в качестве задающего воздействия х* (г) будем принимать случайные воздействия в интервале выходной переменной х2 (г)е[1; 10].

Рис. 5. Управление при задающем воздействии x* (t) в виде ступенчатой функции Fig. 5. Control under the setting action x* (t) in the form of a step function

Рис. 6. Зависимость выхода объекта X2(t) от задающего воздействия x* (t) , носящего случайный характер Fig. 6. The dependence of the output of the object X2(t) on the driving influence x* (t) , which is random

Из рис. 6 можно увидеть, что выход объекта X2(t) также близок к задающему воздействию x2(t) .

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Заключение

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

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

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

ЛИТЕРАТУРА

1. Медведев А.В. Основы теории непараметрических систем / Сиб. гос. аэрокосмич. ун-т. Красноярск, 2018. 727 с.

2. Михов Е.Д. Определение наиболее существенных переменных в задачах моделирования и управления безынерционными

стохастическими процессами // Системы автоматизации в образовании, науке и производстве AS'2017 : тр. конф. Новокузнецк, 2017. С. 424-429.

3. Цыпкин Я.З. Основы информационной теории идентификации. М. : Наука, Гл. ред. физ.-мат. лит., 1984. 320 с.

4. Эйкхофф П. Основы идентификации систем управления // пер. с англ. В.А. Лотоцкого, А.С. Манделя. М. : Мир, 1975. 680 с.

5. Надарая Э.А. Непараметрическое оценивание плотности вероятностей и кривой регрессии. Тбилиси : Изд-во Тбилисского

ун-та, 1983. 194 с.

6. Антомонов Ю.Г., Харламов В.И. Кибернетика и жизнь. М. : Сов. Россия, 1968. 327 с.

7. Медведев А.В., Ярещенко Д.И. О моделировании процесса приобретения знаний студентами в университете // Высшее

образование сегодня. 2017. Вып. 1. С. 7-10.

8. Агафонов Е.Д., Медведев А.В., Орловская Н.Ф., Синюта В.Р., Ярещенко Д.И. Прогнозная модель процесса каталитиче-

ской гидродепарафинизации в условиях недостатка априорных сведений // Известия Тульского государственного университета. Технические науки. 2018. № 9. С. 456-468.

9. Банникова А.В., Корнеева А.А., Корнет М.Е. О непараметрическом дуальном управлении многомерным объектом с за-

паздыванием // Идентификация систем и задачи управления : тр. X Междунар. конф. М., 2015. С. 191-200.

10. Фельдбаум А.А. Основы теории оптимальных автоматических систем. М. : Гос. изд-во физ.-мат. лит., 1963. 553 с.

Поступила в редакцию 17 февраля 2020 г.

Medvedev A.V., Yareschenko D.I. (2020) NONPARAMETRIC IDENTIFICATION AND CONTROL ALGORITHMS FOR MULTIDIMENSIONAL INERTIALESS PROCESSES. Vestnik Tomskogo gosudarstvennogo universiteta. Upravlenie, vychislitelnaja tehnika i informatika [Tomsk State University Journal of Control and Computer Science] 53. pp. 72-81

DOI: 10.17223/19988605/53/7

The work is devoted to the problems of nonparametric identification and control of multidimensional discrete-continuous processes. Discrete-continuous processes occur continuously in space, but their variables are monitored at discrete time instants. Such processes are considered as inertialess with delay. This is explained by the fact that, through various channels of a multidimensional system, the measurement of output variables is carried out at different time intervals. Examples of such processes can be the processes of the mining or processing industries, in particular, in the construction industry (cement production process), metallurgy (steel smelting process), oil refining (diesel hydrotreatment process) and many others

If the output variables of a multidimensional object are somehow stochastically dependent, and this dependence is unknown, then such processes were called T-processes. Such processes require a special look at the identification problem, which is somewhat different from the generally accepted ones. The main thing here is that the identification of such objects should be carried out in a way that is not traditional for the existing theory of identification.

This emphasizes the importance of the identification problem for many real-life processes of discrete-continuous nature. A feature of such processes is that the vector of output variables x(i) = (л^ (t),xz (t),...,xn (i)), j = I,n , consisting of n components, is such that the components of this vector are stochastically dependent in an unknown manner. Denote the vector of input components by u(t) = («j (i),M2(i),...,Mm (i)), k = 1 ,m . Such a formulation of the question leads to the fact that the mathematical description of

the object can be represented in the form of some system of implicit functions of the form Fj (u (t),x(i)) = 0, j = 1 ,n . The main feature of this modeling problem is that the class of dependencies F(-) is a priori unknown, due to a lack of a priori information. That is, there is no parametric class of vector functions F^(u[t),x (i),a), j = 1 ,n , where a is a vector of parameters, which does not

allow the use of parametric identification methods, because the class of functions cannot be determined a priori and the known identification methods are not suitable. As a result of the above, the identification problem is reduced to the problem of solving a system of nonlinear equations j = l,n with respect to the components of the vector x(t) = (хг (t),x2(t),...,*„(i)) , for

known values of м(?) = (щ (i),M2(i),...,Mm(i)). A situation arises when it is necessary to solve a system of interrelated equations, but this system is not defined. In this case, it is possible to use a sequential chain of algorithms to find the values of the components of the vector of output variables from the known input й(/).

The task of controlling a discrete-continuous process is considered under conditions of nonparametric uncertainty, that is, when the object in question is not described up to the parameter vector a. In this case, it is advisable to use a chain of control algorithms to search for the corresponding control action u(t) at each step. If the dimension of the input vector u{t) exceeds the dimension of

the output vector x{t), then some of the components of the input variables u{t) can be interpreted as uncontrolled, but controlled.

This often corresponds to actual technological processes.

We have carried out numerous computational experiments on the identification and control of T-processes, which have shown a fairly high efficiency. The article presents some fragments of numerical studies. In this case, options were considered with different dimensionalities of objects, the level of interference, the volume of training samples, and others.

Keywords: nonparametric modeling; control; inertialess systems; T-processes; T-models; algorithm chain.

MEDVEDEV Alexsander Vasilevich (Doctor of Technical Sciences, Professor, Siberian Federal University, Krasnoyarsk, Russian Federation).

E-mail: [email protected]

YARESCHENKO Darya Igorevna (Senior Lecturer, Siberian Federal University, Krasnoyarsk, Russian Federation). E-mail: [email protected]

REFERENCES

1. Medvedev, A.V. (2018) Osnovy teorii neparametricheskikh sistem [Fundamentals of the theory of nonparametric systems].

Krasnoyarsk: Siberian State Aerospace University.

2. Mikhov, E.D. (2017) [Determination of the most significant variables in the problems of modeling and control of inertialess

stochastic processes]. Sistemy avtomatizatsii v obrazovanii, nauke i proizvodstve AS '2017 [Automation systems in education, science and production]. Proc. of the Conference. pp. 424-429. Novokuznetsk: [s.n.]. pp. 424-429.

3. Tsypkin, Ya.Z. (1984) Osnovy informatsionnoy teorii identifikatsii [Fundamentals of Information Theory of Identification].

Moscow: Nauka.

4. Eykkhoff, P. (1975) Osnovy identifikatsii sistem upravleniya [Fundamentals of Identification Management Systems]. Translated

from English by V.A. Lototsky, A.S. Mandel. Moscow: Mir.

5. Nadaraya, E.A. (1983) Neparametricheskoe otsenivanie plotnosti veroyatnostey i krivoy regressii [Nonparametric estimation of

probability density and regression curve]. Tbilisi: Tbilisi University Press.

6. Antomonov, Yu.G. & Kharlamov, V.I. (1968) Kibernetika i zhizn' [Cybernetics and Life]. Moscow: Sovetskaya Rossiya.

7. Medvedev, A.V. & Yareshchenko, D.I. (2017) O modelirovanii protsessa priobreteniya znaniy studentami v universitete

[On modeling the process of acquiring knowledge by students at a university]. Vysshee obrazovanie segodnya. 1. pp. 7-10.

8. Agafonov, E.D., Medvedev, A.V., Orlovskaya, N.F., Sinyuta, V.R. & Yareshchenko, D.I. (2018) Prognoznaya model' protsessa

kataliticheskoy gidrodeparafinizatsii v usloviyakh nedostatka apriornykh svedeniy [Predictive model of the process of catalytic hydrodewaxing in the absence of a priori information]. Izvestiya Tul'skogo gosudarstvennogo universiteta. Tekhnicheskie nauki -Izvestiya Tula State University. Technical Sciences. 9. pp. 456-468.

9. Bannikova, A.V., Korneeva, A.A.& Kornet, M.E. (2015) [On nonparametric dual control of a multidimensional object with

delay]. Identifikatsiya sistem i zadachi upravleniya [System Identification and Management Tasks]. Proc. of the Conference. Moscow. pp. 191-200.

10. Feldbaum, A.A. (1963) Osnovy teorii optimal'nykh avtomaticheskikh sistem [Fundamentals of the theory of optimal automatic systems]. Moscow: Gos. izd.-vo fiz.-mat. lit.

i Надоели баннеры? Вы всегда можете отключить рекламу.