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

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

CC BY
53
10
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ИДЕНТИФИКАЦИЯ / УПРАВЛЕНИЕ / МНОГОМЕРНЫЙ ОБЪЕКТ / СОСТАВНЫЕ ВЕКТОРА / НЕПАРАМЕТРИЧЕСКИЕ АЛГОРИТМЫ

Аннотация научной статьи по математике, автор научной работы — Ликсонова Д. И., Раскина А. В.

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

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

NON-PARAMETRIC IDENTIFICATION AND CONTROL ALGORITHMS FOR T-PROCESSES

In this paper, we consider nonparametric identification and control methods for multidimensional discrete-continuous processes with delay, which are typical for many real industries. Of course, such systems are typical for practice, including in the rocket and space industry, as well as in technological processes for the production of space technology. In multidimensional processes, we must take into account the relationships between input and output variables, as well as their relationship with each other. Moreover, these connections are not always known to the researcher. Taking into account the unknown connections of the input variables, the researcher will deal with tubular processes or H-models, and if the unknown connections of the output variables are taken into account, the model along one or another channel of the object will be analogs of implicit functions. In general, the model of a multidimensional object will be represented as a system of nonlinear implicit equations. In this case, the solution to the identification problem will be reduced to finding the forecast of the vector of output variables from the known values of the vector of input variables and can be obtained only as a result of solving the corresponding system of equations, which were called T-models, which will be discussed in this article. The solution of a system of nonlinear implicit equations by parametric identification methods will not lead to the desired result, due to the lack of sufficient a priori information, this is where the need to use nonparametric identification methods arises, as well as the necessary use of system analysis methods. A priori information in problems of nonparametric statistics is insufficient, which cannot be dealt with by generally accepted identification methods. When managing multidimensional processes, the dependencies of the output variables should be taken into account. Here another important feature arises, which consists in the fact that random values from the range of definition of output variables cannot be used as reference influences, they must be selected from their common intersection.

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

УДК 519.711.3

Doi: 10.31772/2712-8970-2021-22-4-600-612

Для цитирования: Ликсонова Д. И., Раскина А. В. Непараметрические алгоритмы идентификации и управления для Т-процессов // Сибирский аэрокосмический журнал. 2021. Т. 22, № 4. С. 600-612. Doi: 10.31772/2712-8970-2021-22-4-600-612.

For citation: Liksonova D. I., Raskina A. V. Non-parametric identification and control algorithms for T-processes. Siberian Aerospace Journal. 2021, Vol. 22, No. 4, P. 600-612. Doi: 10.31772/2712-8970-2021-22-4-600-612.

Непараметрические алгоритмы идентификации и управления для Т-процессов

Д. И. Ликсонова, А. В. Раскина

Сибирский федеральный университет Российская Федерация, 660074, ул. Академика Киренского, 26, к. 1 E-mail: [email protected]

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

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

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

Non-parametric identification and control algorithms

for T-processes

D. I. Liksonova, A. V. Raskina

Siberian Federal University 26 k., 1, Academician Kirensky St., 660074, Russian Federation E-mail: [email protected]

In this paper, we consider nonparametric identification and control methods for multidimensional discrete-continuous processes with delay, which are typical for many real industries. Of course, such systems are typical for practice, including in the rocket and space industry, as well as in technological processes for the production of space technology. In multidimensional processes, we must take into account the relationships between input and output variables, as well as their relationship with each other. Moreover, these connections are not always known to the researcher. Taking into account the unknown connections of the input variables, the researcher will deal with tubular processes or H-models, and if the unknown connections of the output variables are taken into account, the model along one or another channel of the object will be analogs of implicit functions. In general, the model of a multidimensional object will be represented as a system of nonlinear implicit equations. In this case, the solution to the identification problem will be reduced to finding the forecast of the vector of output variables from the known values of the vector of input variables and can be obtained only as a result of solving the corresponding system of equations, which were called T-models, which will be discussed in this article. The solution of a system of nonlinear implicit equations by parametric identification methods will not lead to the desired result, due to the lack of sufficient a priori information, this is where the need to use nonparametric identification methods arises, as well as the necessary use of system analysis methods. A priori information in problems of nonparametric statistics is insufficient, which cannot be dealt with by generally accepted identification methods.

When managing multidimensional processes, the dependencies of the output variables should be taken into account. Here another important feature arises, which consists in the fact that random values from the range of definition of output variables cannot be used as reference influences, they must be selected from their common intersection.

Keywords: identification, control, multidimensional object, composite vectors, nonparametric algorithms.

Введение

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

Подобные многомерные объекты (процессы) имеют чаще всего зависимости выходных переменных х = (,х2,...,xn), стохастически зависимых заранее неизвестным образом (Т-процессы). Это приводит к тому, что математическое описание многомерного объекта будет представлять-

ся в виде некоторого аналога системы неявных функций вида — (и, х) = 0, у = 1, п , где и = ( и1, и2,..., ит ) - входные переменные процесса. Основная особенность подобных объектов состоит в том, что класс зависимостей — (•) неизвестен с точностью до параметров. В этом случае классическая теория идентификации не применима. Задача идентификации сводится к задаче решения системы нелинейных уравнений — (и, х) = 0, у = 1,п относительно компонент

вектора х = (х1,х2,...,хп) при известных значениях и, х и наличии обучающей выборки (иг-, х{, 1 = 1,5) , которая может быть решена при использовании теории непараметрических систем [1]. Подобные задачи ранее не рассматривались [3; 4].

Воспользуемся конкретизацией понятия, связанного с термином «непараметрический», которое было использовано в работах М. Розенблатта и Э. Парзена [5; 6], а также в монографии Ф. П. Тарасенко [7].

«Непараметрическая задача - это статистическая задача, определенная на таких классах распределений, среди которых хотя бы один не сводится к параметрическому семейству функций».

«Непараметрическая задача оценивания неизвестных распределений - это задача нахождения процедуры, с помощью которой можно оценивать непараметризованные распределения из класса, например, всех непрерывных функций распределения или класса распределений, имеющих ряд производных и т. д.».

В настоящей же статье термин «непараметрический» означает отсутствие достаточной априорной информации для представления объекта с точность до параметров.

Исходя из данных определений, можно сказать, что непараметрические методы намного выигрывают в эффективности, если параметрическая модель не описывает адекватно наблюдаемые данные [8].

Задача идентификации

Рассмотрим многомерный процесс, представленный на рис. 1.

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

На рис. 1 приведены следующие обозначения: u (t) = (ux (t),u2 (t),...,uk (t),...,um (t)), k = 1,m - m-мерный вектор входных переменных; x(t) = (x1 (t),x2 (t),...,Xj (t),...,xn (t)), j = 1,n - n-мерный вектор выходных переменных, которые принадлежат соответствующим областям: Xj eQj (x) ; ^(t) - случайные помехи, действующие на объект; hu, hx - случайные помехи

измерении соответствующих переменных процесса; пунктирные линии свидетельствуют о наличии зависимости входных и выходных переменных; (V) - непрерывное время; и(, х1 - измерения входных и выходных переменных в дискретный момент времени V.

На основании вышеизложенного можно сказать, что многомерные процессы окружают нашу повседневную жизнь и представляют собой процессы, имеющие много входов и выходов, а также связи входных и выходных переменных между собой, связи только между входными переменными [9; 10] и только между выходными. Причем все эти связи не всегда известны исследователю [11], а может также присутствовать разнотипная априорная информация [12]. В вязи с тем, что выходные переменные многомерного процесса связаны между собой, следует отметить и составной (ситуационный) вектор, который был введен Я. З. Цыпкиным в [13]. Составной вектор - это вектор, составленный из входных и выходных переменных. Это может быть любой набор, например х<3> = (и2, и5, ц^ х4). Составной вектор известен исследователю

из априорной информации.

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

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

Рис. 2. Многомерный дискретно-непрерывный процесс Fig. 2. Multidimensional discrete-continuous process

На рис. 2 на вход рассматриваемого процесса поступает вектор входных управляемых переменных и (V) = (и (V),...ит (V)) е (и) с Ят и вектор входных неуправляемых, но контролируемых

переменных ) = (^1 (V),...Цр (V)}е(ц)с Яр, на выходе наблюдается вектор выходных переменных х(V) = (х1 (V),...хп (V))е(х)с Я", X(V) - выход модели. Входные и выходные перемен-

ные контролируются в дискретные моменты времени через интервал Аг. Средствами контроля Ни, НЦ, Нх получаем выборку наблюдений или обучающую выборку. В каналах измерения переменных действуют случайные помехи Ни, ИЦ, Нх.

Как было сказано выше, процессы, рассматриваемые в настоящей работе, имеют неизвестные зависимости компонент выходных переменных. Поэтому исследуемый процесс будет описываться системой неявных стохастических уравнений:

- (и (г), ц(г), х(г + т), (г)) = 0, у = 1Я (1)

где функции — (•) не известны, так как не известны зависимости выходных переменных;

т - известное запаздыванием по различным каналам изучаемого процесса.

Задача идентификации состоит в построении модели системы, которая представлена на рис. 2, при наличии выборки наблюдений над объектом: и{ =(ия, и{ 2, ..., ит),

Ц/ = (Ц/1, Ц/2, ..., Цр ), х/ =(хЯ, х/2, "., х/п ), 1 = 1,5

В этом случае Т-модель процесса с неизвестными зависимостями компонент вектора выходных переменных будет рассматриваться в виде системы:

- (и<у>, Ц<у>, х<у>, й,, Ц5

, ) = 0, у = 1,п,

(2)

где и<у>, ц<у>, х<у> - составные векторы, и5, Ц5, х5 - временные векторы (т. е. набор данных, которые поступили к 5-му моменту времени).

В результате измерения входных и выходных переменных может быть получена обучающая

выб°рка иг ^^ и/2, ит ), Ц = ( Цг2, Цр), х, = ( x/1, Xг2, - х,п ), 1 = 1 5 которЗЯ используется при построении модели многомерного объекта. Так как входные воздействия и заданы и известны, то решаем систему (2) и получаем оценки хсу- компонент вектора выходных

переменных х у при соответствующих значениях входных воздействий и . Здесь естественно воспользуемся методами непараметрического оценивания [14].

Для начала подставляем 1-е поступление входных переменных и} = (и

), I = 1,5, где

5 - объем обучающей выборки, в формулу (2). Далее подставляем выходные переменные из обучающей выборки х{ =(хг1, ..., хпп ) для определения невязок е1, / = 1,5. Невязки е1, / = 1,5 вычисляются по следующей формуле:

гу =

/у (и <1 > , Ц<1 > , х<1 > , и5 , Ц5 , х5 ) , у = 1,n,

(3)

где функции /у (•) принимаются в виде разницы между измеренными значениями компонент выхода и их оценками [15]:

I ху [/]П Ф

<т> т^ гЦ;-Цу[/]^

() = ху ()-

/=1

к=1

и'к - ик [/]

П Ф

"5ик У У=1

\ С5Цу

(и'к- ик[/] Ф [/] ^

1П Ф

/=1 к=1

5ик У

У=1

V С5Цу

у = 1, п,

(4)

где < т > - размерность составного вектора ик, < т >< т. Далее оцениваем условное математическое ожидание

Xj = М {Xj

< j > < j > ' , f , S =

0}, j = 1,n,

(5)

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

I xj и П Ф

<m>> „ f % - uki [i]) <p> ^ f ^ - [i]

x j =■

/=1

ki=i

1=1 v

П Ф

-su /Vl = 1

-V1 1"V1 I

V Csf

П Ф

k2=1

Sk2[i]

V css у

1П Ф

¿=1 k1=1

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

f uk1 - Uk1[i'] Y<P>

П Ф

-su /V1 = 1

V Csf

П Ф

k2=1

fSk2[i] ^

V Css

j = 1, n,

(6)

где в качестве колоколообразных функций Ф (•) были приняты функции различных типов и коэффициенты си, с, с^ удовлетворяли условиям сходимости [16].

Оценка (6) и является прогнозом выходных переменных х при известных входных и.

Задача управления

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

Рис. 3. Схема непараметрической системы управления объектом с запаздыванием Fig. 3. Scheme of a nonparametric control system for an object with delay

На рис. 3 приняты следующие обозначения: u (t) = (u1 (t), ..., um (t)) - управляемые входные переменные; f(t) = (f (t), ..., fp (t)) - неуправляемые, но контролируемые входные переменные; x(t + т) = (x1 (t + т), ..., xn (t + т))е Rn - выходные переменные процесса; x* (t + т) = (x* (t + т), ..., x* (t + т))е Rn - задающие воздействия; E,t, hf, h* - случайные стационарные помехи, действующие на объект и в каналах измерения входных и выходных пере-

менных; т - запаздывание по различным каналам многомерной системы. Запаздывание т известно по всем каналам многомерной системы, и в данном случае является одинаковым для каждой компоненты вектора выходных переменных.

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

Рис. 4. Схема непараметрической системы управления с дополнительным блоком Fig. 4. Diagram of a nonparametric control system with an additional block

Рис. 4 в отличие от рис. 3 дополнен блоком формирования задающих воздействий для их определения. На рис. 4 приняты следующие обозначения: x (t) - исходные значения задающих воздействий; x** (t) - задающие воздействия, которые необходимо найти из системы уравнений

F] (u<]> (t), ц<]> (t), /<]> (t)) = 0, j = Vn, (7)

где функции F] (•) продолжают оставаться неизвестными. Изложим процедуру управления, начиная с конкретного момента времени t.

Пусть имеется обучающая выборка (и', ц', x', i = 1,5). В конкретный момент времени t на

вход поступают неуправляемые, но контролируемые переменные ц, при этом управляющие воздействия ut и выходные переменные xt пока неизвестны. Далее из всей первоначальной

выборки (и', ц', x', i = 1,5) выбираются те строки, в которых значения ц наиболее близки

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

воздействия х** находятся из новой сформированной выборки х** е А(ху- ), а именно из решения системы (7). Решение системы (7) сводится к последовательности алгоритмов, изложенной ниже.

Для задающего воздействия хфФ берем произвольные значения из области ^ (хх) . Вторую переменную хфФ определяем с учетом выбранной компоненты хфФ из следующего выражения:

ФФ''

I х2 ф '=1 (х** х' Л< т> П Ф У к=1 ( ' ик - ик Л< р> П У У=1 ( ' Л Му -Му

V сх1 V % V с*м У

51 I Ф '=1 V х** х' л <т> ( П Ф к=1 V г Л ик- ик < р>( П У=1 V ' Л Му -Му

Сх1 У % У с*м У

(8)

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

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

Л<*>^-М1Л

ху = "

5у-1 у-1 I 4 П '=1 у=1 ( ФФ ' ху - ху Л<т> П Ф Ук=1 ( ик - и

Сх, \ У V С"к

5 у-1 у-1 ( IП '=1 у=1 V ФФ ' \ х - х ху Ху ст> ( П Ф к=1 V ик - ик

С ху У С"к

п

у=1

V У

<р» му-м1л

(9)

п

У=1

V У

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

области 0.(и1). Входная переменная и2 (t) может быть определена в соответствии со следующим алгоритмом:

I и2ф

'=1

и9 =-

и1 - и1

V Си1 У

' Л п (

п Ф

у =1

(.. ..< Л

'=1

IФ П Ф

V Си1 У

У=1

ху Л р П Ф >=1 ( ' Л Му-Му

Сц, V ^у У

' л у р П Ф У=1 V ' Л Му-Му

У СМу У

(10)

где (и', ц', х', ' = 1,5) - обучающая выборка; цУ - поступившие входные неуправляемые, но

контролируемые переменные.

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

I и3Ф

( ' Л

и, - и

'=1

и =-

1_"1

V и1 У

Ф

и2 - и2

V и2 У

у=1

'=1

(. .г Л

V и1 У

Ф

V и2 У

' Л п (

П Ф

у=1

( ФФ ' ху - ху Л р П Ф м ( ' Л Му -МУ

Сх, V у СМ V цу У

ФФ ' Л ху - ху р ( П Ф У=1 V ' Л Му -Му

Сх, У Сцу У

(11)

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

uk =

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

s k-1 z uk п ф i=1 k=1 f i uk - uk \ n п Ф ) j=1 f ** i xj- xj \ P п Ф )v=1 f i ^ Uv -Uv

I Cuk Cx, V j C u V uv )

s k-1 f Z п Ф i=1 k=1 V i \ uk uk n f п Ф j=1 I ** i ^ Xj - Xj p f п Ф v=1 V i > Uv -Uv

C uk cx. xj ) c Uv )

, k = 1, m .

(12)

В реальных задачах часто число компонент вектора входных переменных и больше числа компонент вектора выходных переменных х . В этом случае в число компонент вектора ц включают компоненты вектора и,, с тем, чтобы размерность вектора и и х сделать одинаковой. Настраиваемыми параметрами будут параметры размытости сщ

cx и cu , для них можно

л j H-V

использовать следующие формулы: сщ = а

u,, -u.

S =Р

Xj - Xj

и cuv =Y

Uv-Uv

где а, p

и у некоторые параметры большие 1, а > 1, в > 1, у > 1. Следует заметить, что выбор си

Uv

осуществляется на каждом такте управления. При этом если сначала определен сик, то оп-

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

Вычислительные эксперименты

Для примера был взят объект с 4-мя входными переменными и (t) = (и1 (t), и2 (t), и3 ^), и4 (t)) и тремя выходными переменными х (7) = (хх ^), х2 (t), х3 (t)) .

Подобные объекты характерны для реальных производств, например, для производства цемента, нефтепереработки, металлургии и т. д., где протекают только многомерные процессы. Для данного объекта была сформирована выборка входных и выходных переменных и найдены прогнозируемые значения выходных переменных при известных входных с использованием алгоритма (4) и (6). Для вычисления использовался объем выборки 5 = 1000, параметры размытости сш = 0,4; с5е = 0,2, помеха, действующая на компоненты вектора выходных переменных а = 0,1, запаздывание т = 2 такта (где 1 такт равен 1 минуте либо иной). Описание объекта с точностью до параметров было принято только для проведения компьютерного исследования и оставалось неизвестным для изложенной выше теории. Представим полученные результаты для третьей компоненты выхода многомерного объекта.

Рис. 5. Прогноз выходной переменой x3 (t) объекта, измеренный с равномерной помехой 10 %, 83 = 0,07 Fig. 5. The predicted output variable x3 (t) of the object, measured with a uniform noise of 10 %, 83 = 0,07

c

и

На рис. 5 представлены полученные прогнозируемые значения х3 (у) для третьей компоненты выхода х3 () . Погрешности определения переменных обусловлены наличием помех. Как видно из рис. 5, прогноз значений выходных переменных х3 (t) многомерного объекта по известным входным переменным и (t) с учетом 10 % помехи достаточно удовлетворителен

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

выходных переменных (и', ц', х', ' = 1,5).

Далее приводятся результаты вычислительных экспериментов для данного объекта для первой компоненты выхода при использовании алгоритма управления (12). В данном вычислительном эксперименте число компонент вектора и больше числа компонент вектора х. Если же размерность вектора и превышает размерность вектора х, т. е. т > п, то сделаем замену: и4 (t) ^ ^), чтобы размерность вектора и их сделать одинаковой. Далее значения задающих воздействий х" (), которые были найдены при помощи алгоритма (9), представим в виде ступенчатой функции, как показано на рис. 6.

Рис. 6. Зависимость выхода объекта xj (t) от задающего воздействия xj**(t) Fig. 6. Dependence of the output of the object xj (t) on the setting action xj**(t)

Как видно из рис. 6, выход объекта xj (t) достаточно близок к задающему воздействию xj** (t) при наличии 10 % помехи и не превышает 1,7 % от значений х** (t), что удовлетворяет

большинству практических задач.

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

Заключение

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

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

Еще основной особенностью является то, что и при задаче идентификации и задаче управления вводятся последовательности непараметрических алгоритмов. Сущность непараметрических алгоритмов идентификации и управления состоит в том, что для задачи идентификации прогноз выстраивается как условное математическое ожидание компонент вектора выходных переменных x (t) при известных компонентах входных u (t), а при управлении - как условное

математическое ожидание входных переменных u (t) при найденных задающих воздействиях

x"(t).

Исходя из полученных прогнозируемых значений в вычислительном эксперименте, в процессе использования непараметрических алгоритмов, получаются достаточно хорошие с практической точки зрения результаты, которые можно использовать и при решении реальных задач на производствах [17].

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

1. Медведев А. В. Основы теории непараметрических систем. Красноярск : СибГУ им. М. Ф. Решетнева, 2018. 727 с.

2. Васильев В. А., Добровидов А. В., Кошкин Г. М. Непараметрическое оценивание функционалов от распределений стационарных последовательностей. М. : Наука, 2004. 508 с.

3. Методы классической и современной теории автоматического управления: в 5 т. Т. 1: Математические модели, динамические характеристики и анализ систем управления / под. ред. К. А. Пупкова, Н. Д. Егупова. М. : МГТУ им. Н. Э. Баумана, 2004. 748 с.

4. Методы классической и современной теории автоматического управления: в 5 т. Т. 2: Статистическая динамика и идентификация систем автоматического управления / под. ред. К. А. Пупкова, Н. Д. Егупова. М. : МГТУ им. Н.Э. Баумана, 2004. 640 с.

5. Rozenblatt M. Remarks on some nonparametric estimates of density function // Ann. Math. Statist. 1956. Vol. 27. P. 832-837.

6. Parzen E. On estimation of probability density function and mode // Ann. Math. Statist. 1962. Vol. 33. № 3. P. 1065-1076.

7. Тарасенко Ф. П. Непараметрическая статистика. Томск : Изд-во Том. ун-та, 1976. 292 с.

8. Кошкин Г. М., Пивен И. Г. Непараметрическая идентификация стохастических объектов. Хабаровск : РАН Дальневосточное отделение, 2009. 336 с.

9. Медведев А. В. Н-модели для безынерционных систем с запаздыванием // Вестник Сиб-ГАУ. 2012. № 5 (45). С. 84-89.

10. Medvedev A. V., Mihov E. D., Nepomnyashchiy O.V . Mathematical modeling of H-processes // Журнал Сибирского федерального ун-та. Сер. «Математика и физика». 2016. Т. 9, № 3. С.338-346.

11. Медведев А. В., Ярещенко Д. И. Непараметрические алгоритмы идентификации и управления многомерными безынерционными процессами // Вестник Томского гос. ун-та. Управление, вычислительная техника и информатика. 2020. № 53. С. 72-81.

12. Медведев А. В. Теория непараметрических систем. К-модели // Вестник СибГАУ. 2011. № 3 (36). С. 57-62.

13. Цыпкин Я. З. Адаптация и обучение в автоматических системах. М. : Наука, 1968. 399 с.

14. Терешина А. В., Ярещенко Д. И. О непараметрическом моделировании безынерционных систем с запаздыванием // Сибирский журнал науки и технологий. 2018. Т. 19, № 3. С. 452-461.

15. Хардле В. Прикладная непараметрическая регрессия. М. : Мир, 1993. 349 с.

16. Надарая Э. А. Замечания о непараметрических оценках плотности вероятности и кривой регрессии // Теория вероятностей и ее применение. 1970. Т. 15, вып. 1. С. 139-142.

17. Прогнозная модель процесса каталитической гидродепарафинизации в условиях недостатка априорных сведений / Е. Д. Агафонов, А. В. Медведев, Н. Ф. Орловская и др. // Изв. ТулГУ. 2018. № 9. С. 456-468.

References

1. Medvedev A. V. Osnovy teorii neparametricheskikh sistem [Fundamentals of the theory of non-parametric systems]. Krasnoyarsk, SibGU im. M. F. Reshetneva Publ., 2018, 732 p.

2. Vasil'ev V. A., Dobrovidov A. V., Koshkin G. M. Neparametricheskoe otsenivanie funktsion-alov ot raspredeleniy statsionarnykhposledovatel'nostey [Nonparametric Estimation of Functionals of Distributions of Stationary Sequences]. Moscow, Nauka Publ., 2004, 508 p.

3. Pupkov K. A., Egupov N. D. Metody klassicheskoy i sovremennoy teorii avtomaticheskogo upravleniya: v 5 t. T. 1: Matematicheskie modeli, dinamicheskie kharakteristiki i analiz sistem uprav-leniya [Methods of classical and modern theory of automatic control: in 5 vol. Vol. 1: Mathematical models, dynamic characteristics and analysis of control systems]. Moscow, MGTU im. N. E. Baumana Publ., 2004. 748 p.

4. Pupkov K. A., Egupov N. D. Metody klassicheskoy i sovremennoy teorii avtomaticheskogo upravleniya: v 5 t. T. 2: Statisticheskaya dinamika i identifikatsiya sistem avtomaticheskogo upravleniya [Methods of classical and modern theory of automatic control: in 5 vol. Vol. 2: Statistical dynamics and identification of automatic control systems]. Moscow, MGTU im. N. E. Baumana Publ., 2004. 640 p.

5. Rozenblatt M. Remarks on some nonparametric estimates of density function. Ann. Math. Statist. 1956, Vol. 27, P. 832-837.

6. Parzen E. On estimation of probability density function and mode. Ann. Math. Statist. 1962, Vol. 33, No. 3, P. 1065-1076.

7. Tarasenko F. P. Neparametricheskaya statistika [Nonparametric statistics]. Tomsk, Izd-vo Tom. un-ta Publ., 1976, 292 p.

8. Koshkin G. M., Piven I. G. Neparametricheskaya identifikatsiya stokhasticheskikh ob"ektov [Nonparametric identification of stochastic objects]. Khabarovsk, RAN Dal'nevostochnoe otdelenie Publ., 2009, 336 p.

9. Medvedev A. V. [H-models for inertialess lag systems]. VestnikSibGAU. 2012, No. 5, P. 84-89 (In Russ.).

10. Medvedev A. V., Mihov E. D., Nepomnyashchiy O. V. [Mathematical modeling of H-processes]. Zhurnal Sibirskogo federal'nogo universiteta. Ser. Matematika i fizika. 2016, Vol. 9, No. 3, P. 338-346 (In Russ.).

11. Medvedev A. V., Yareshchenko D. I. [Nonparametric algorithms for identification and control of multidimensional inertialess processes]. Vestnik Tomskogo gosudarstvennogo universiteta. Uprav-lenie, vychislitel'naya tekhnika i informatika. 2020, No. 53, P. 72-81 (In Russ.).

12. Medvedev A. V. [Theory of nonparametric systems. K-models]. Vestnik SibGAU. 2011, No. 3, P. 57-62 (In Russ.).

13. Tsypkin Ya. Z. Adaptatsiya i obuchenie v avtomaticheskikh sistemakh [Adaptation and training in automatic systems]. Moscow, Nauka Publ., 1968, 399 p.

14. Tereshina A. V., Yareshchenko D. I. [On nonparametric modeling of inertialess systems with delay]. Sibirskiy zhurnal nauki i tekhnologiy. 2018, Vol. 19, No. 3, P. 452-461 (In Russ.).

15. Khardle V. Prikladnaya neparametricheskaya regressiya [Applied nonparametric regression]. Moscow, Mir Publ., 1993, 349 p.

16. Nadaraya E. A. [Notes on Nonparametric Probability Density Estimates and Regression Curves]. Teoriya veroyatnostey i eeprimenenie. 1970, Vol. 15, No. 1, p. 139-142 (In Russ.).

17. Agafonov E. D., Medvedev A. V., Orlovskaya N. F., Sinyuta V. R., Yareshchenko D. I. [Predictive model of the catalytic hydrodewaxing process under conditions of a lack of a priori information]. Izv. TulGU. 2018, No. 9, P. 456-468 (In Russ.).

@ ^HKCOHOBa fl. H., PacKHHa A. B., 2021

Ликсонова Дарья Игоревна - кандидат технических наук, доцент базовой кафедры интеллектуальных систем управления Института космических и информационных технологий, Сибирский федеральный университет. E-mail: [email protected].

Раскина Анастасия Владимировна - кандидат технических наук, доцент кафедры информационных систем Института космических и информационных технологий; Сибирский федеральный университет. E-mail: [email protected].

Liksonova Darya Igorevna - Cand. Sc., docent of the department Intelligent Control Systems of the Institute of Space and Information Technologies; Siberian Federal University. E-mail: [email protected].

Raskina Anastasia Vladimirovna - Cand. Sc., docent of Department of Information Systems of the Institute of Space and Information Technologies; Siberian Federal University. E-mail: [email protected].

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