Scientific journal
Modern high technologies
ISSN 1812-7320
"Перечень" ВАК
ИФ РИНЦ = 0,940

QUASI-STATIC NEIGHBORHOOD SYSTEMS

Shmyrin A.M. 1 Mishachev N.M. 1 Kanyugina A.S. 1
1 Lipetsk State Technical University
In this paper, we define static and dynamic neighborhood systems in the language of metagraphs. The control of production process near a nominal mode is usually carried out on the assumption that this mode corresponds to a stationary point of a hypothetical global dynamic model, and therefore static neighborhood systems in this case are a convenient means of local modeling. A typical situation is that when there are several nominal modes and the transitions between them are carried out discretely or continuously depending on the input data. We consider the problem of integration of local static neighborhood systems (models) into a combined quasi-static system with coefficients that depend on the input data. Two versions of this integration are proposed: discrete and continuous. Both are based on the clustering of the input data of the neighborhood system. In the discrete version, the choice of the local model is determined by the cluster center. In the continuous version, linear combinations of local models are used with weights forming the partition of unity associated with clustering. As an example of the application, the problem of integration of local neighborhood models for the stage of diffusion of sugar production, obtained by the parametric identification based on a sample of production data of JSC «Avrora», Borinsky Sugar Plant, is considered.
neighborhood structures
neighborhood systems
metagraphs
quasi-static systems
clustering
parametric identification
sugar production

Цели и задачи исследования

Окрестностные системы, или окрестностные модели – это статические или дискретные динамические системы уравнений на графе. Уравнения окрестностной системы могут соответствовать либо вершинам графа, либо его ребрам (см. [1]). В первом случае мы называем окрестностную систему вертексной, во втором – реляционной. Уравнения окрестностной системы структурно идентифицированы по наборам входящих в них переменных: этот уровень структурной идентификации задается графом. В задачах управления производственными процессами такие системы возникают на этапе перехода от структурной схемы производства (графа) к математической модели. Для моделирования вблизи номинальных режимов, как правило, достаточно рассматривать статические линейные или билинейные системы уравнений. При наличии нескольких номинальных режимов во многих случаях имеет смысл, вместо усложнения уравнений, перейти к квазистатической схеме. В данной статье мы обсуждаем общую теоретическую задачу объединения нескольких локальных статических окрестностных систем в квазистатическую систему и, далее, применяем предлагаемую схему объединения в задаче моделирования стадии диффузии производства сахара, рассматривавшейся в [2].

Окрестностные структуры и окрестностные системы

Окрестностные структуры (см. [1]) являются средством формализации связей между элементами моделируемой системы. На их основе удобно строить математические модели следующего уровня (системы уравнений). Окрестностной структурой мы называем ориентированный граф hmir01.wmf, содержащий вершины hmir02.wmf трех типов: входы U, узлы V и выходы W, при этом:

- каждые два узла hmir03.wmf, могут быть соединены между собой не более чем двумя противоположно ориентированными ребрами-связями e(v', v») и e(v», v');

- каждый узел hmir04.wmf имеет петлю e(v, v);

- каждый узел hmir05.wmf имеет входящие и выходящие ребра (помимо петель);

- каждый вход hmir06.wmf имеет только выходящие ребра e(u, v);

- каждый выход hmir07.wmf имеет только входящие ребра e(v, w).

Как обычно, источниками вершины называются все входящие в нее вершины и стоками – все исходящие. Все узлы (то есть вершины hmir08.wmf) из-за наличия петель являются своими стоками и источниками, все входы имеют только стоки и все выходы – только источники.

Окрестностные системы удобно описывать с помощью понятия метаграфа (см., например, [3]). Нам потребуется только определение метаграфа: метавершинами MV метаграфа hmir09.wmf называются подмножества некоторого конечного множества V, а метаребрами ME – пары метавершин. Или, на языке теории множеств: метавершины метаграфа – это элементы первого булеана hmir10.wmf, то есть hmir11.wmf, а метаребра метаграфа – это двухэлементные подмножества второго булеана hmir12.wmf, то есть hmir13.wmf. По аналогии с обычными графами можно определить ориентированные метаграфы и двудольные метаграфы. В окрестностной структуре hmir14.wmf каждый узел hmir15.wmf порождает метавершину источников (метаисточник) hmir16.wmf и метавершину стоков (метасток) hmir17.wmf, при этом hmir18.wmf. Обозначим через V+ и V– множества всех метаисточников и метастоков. Добавим еще к множеству метавершин все узлы hmir19.wmf. Каждый узел hmir20.wmf порождает метаребро (v+, v), соединяющее метаисточник узла с этим узлом, и метаребро (v+, v–), соединяющее метаисточник узла и его метасток. Таким образом, каждая окрестностная структура hmir21.wmf порождает ориентированные двудольные метаграфы hmir22.wmf и hmir23.wmf. Метаграф hmir24.wmf мы будем называть вертексным, а метаграф hmir25.wmf – реляционным. Эти метаграфы соответствуют вертексным и реляционным окрестностным системам над окрестностной структурой hmir26.wmf, или системам для узлов и для ребер (см. [1]). Опишем подробнее эти системы. Обозначения множеств вершин U, V, W и метаузлов v+, v– можно понимать и как обозначения соответствующих наборов чисел – номеров вершин. Каждому ребру e(i, k), включая петли, соответствует переменная hmir27.wmf входа-выхода из i-го узла в k-й. Переменные Y(i, k) = U(i, k) с ihmir28.wmf мы называем входами, переменные Y(i, k) = W(i, k) с hmir29.wmf – выходами, переменные Y(i, k) = V(i, k) с i, hmir30.wmf – внутренними переменными. Для петель e(i, i) положим Y(i, i) = X(i) и hmir31.wmf. Переменную Y(i, i) = X(i) мы называем состоянием узла. В вертексной модели вход-выход Y(i, k) не зависит от k, то есть узел hmir32.wmf передает по всем исходящим связям одну и ту же переменную своего состояния X(i) = Y(i, i), а вход hmir33.wmf – одну и ту же переменную входа hmir34.wmf. Далее, метаисточникам hmir35.wmf узлов hmir36.wmf соответствуют переменные состояния метаисточника hmir37.wmf, где hmir38.wmf – это произведение всех пространств hmir39.wmf c hmir40.wmf. В вертексном случае мы полагаем hmir41.wmf, где hmir42.wmf – это произведение всех пространств hmir43.wmf c hmir44.wmf. Метастокам hmir45.wmf узлов hmir46.wmf соответствуют переменные состояния метастока hmir47.wmf, где hmir48.wmf – это произведение всех пространств hmir49.wmf c hmir50.wmf. Вертексному метаграфу hmir51.wmf соответствует набор (мета)функций

hmir52.wmf, hmir53.wmf, (1)

каждая из которых преобразует состояние X+(i) метаисточника узла в состояние X(i) этого узла. Реляционному метаграфу hmir54.wmf соответствует набор (мета)операторов

hmir55.wmf, hmir56.wmf, (2)

каждый из которых преобразует состояние Y+(i) метаисточника узла i в состояние Y–(i) метастока этого узла. Уравнения статической вертексной системы имеют вид

hmir57.wmf, hmir58.wmf, (3)

здесь hmir59.wmf уравнений, а уравнения статической реляционной системы имеют вид

hmir60.wmf, hmir61.wmf, (4)

здесь hmir62.wmf операторных уравнений. Подробнее:

hmir63.wmf, hmir64.wmf, (5)

hmir65.wmf,

hmir66.wmf, hmir67.wmf. (6)

В последнем случае количество уравнений равно hmir68.wmf, где hmir69.wmf – это количество вершин в метастоке v–. Динамические версии систем (3) и (4) имеют соответственно вид

hmir70.wmf, hmir71.wmf, (7)

hmir72.wmf, hmir73.wmf. (8)

Квазистатические окрестностные системы

При моделировании производственного процесса вблизи номинального режима во многих случаях можно считать, что такому режиму соответствует стационарная точка гипотетической глобальной динамической модели процесса, и потому в качестве локальной модели можно использовать статическую линейную или билинейную окрестностную систему. Далее, достаточно типичной является ситуация, когда номинальных режимов несколько и диспетчер производства изменяет управляющие переменные в зависимости от входных данных. Построения более сложной динамической модели, описывающей переходные режимы (а не только номинальные) на основе данных наблюдений, то есть регрессионными методами, обычно не дает хороших результатов. Поэтому при моделировании процессов с несколькими номинальными режимами имеет смысл, вместо усложнения уравнений, перейти к следующей квазистатической схеме. Множество вершин-входов U можно считать состоящим из независимых внешних входов hmir74.wmf и зависимых (управляемых) входов hmir75.wmf. Соответственно, переменные U(i, k) бывают двух типов: внешние hmir76.wmf и управляемые hmir77.wmf. Обозначим через hmir78.wmf множество всех кортежей состояний системы (то есть значений входов, внутренних переменных и выходов) во все моменты наблюдения. Пусть NU – размерность пространства всех входов. Тогда hmir79.wmf, где hmir80.wmf – размерность пространства внешних входов, и hmir81.wmf. Можно предполагать, что множеству номинальных режимов соответствует кластеризация множества входных данных hmir82.wmf с центрами кластеров hmir83.wmf. Для определенности мы рассматриваем далее только случай реляционных систем (для вертексных все аналогично). Обозначим через hmir84.wmf регрессионные модели, то есть системы вида (4), построенные для каждого из номинальных режимов hmir85.wmf в результате параметрической идентификации по соответствующим кортежам данных («номинальные» модели). Здесь hmir86.wmf – покомпонентное действие вектора метаоператоров hmir87.wmf на вектор метаисточников hmir88.wmf (операторное произведение Адамара). Очевидным образом определены линейные комбинации hmir89.wmf этих моделей. Квазистатической окрестностной системой мы называем систему вида

hmir90.wmf (9)

c коэффициентами, зависящими hmir91.wmf, такими что hmir92.wmf, hmir93.wmf. Для hmir94.wmf обозначим через s(d) номер ближайшего к d центра кластера; если таких центров несколько, то s(d) – наименьший из номеров. Пусть hmir95.wmf – соответствующее разбиение пространства hmir96.wmf (разбиение Дирихле – Вороного) и χr(d) – характеристические функции множеств hmir97.wmf. Дискретной композицией номинальных моделей hmir98.wmf мы называем квазистатическую систему

hmir99.wmf. (10)

Для внешнего входа hmir100.wmf эта система совпадает с номинальной hmir101.wmf. Далее, для каждого из множеств hmir102.wmf вычислим ковариационную матрицу и обозначим через φ1,…,φS соответствующие hmir103.wmf мерные нормальные плотности с центрами hmir104.wmf. Положим

hmir105.wmf. (11)

Функции hmir106.wmf образуют разбиение единицы на пространстве hmir107.wmf. Непрерывной нормальной композицией номинальных моделей hmir108.wmf мы называем квазистатическую систему

hmir109.wmf. (12)

Замечание. Вместо многомерного нормального распределения можно, если это удобно по каким-либо соображениям, использовать любое другое унимодальное многомерное распределение. Описанную конструкцию можно, при желании, интерпретировать как нечеткую модель Такаги – Сугено [4, 5], но в данном случае язык нечетких моделей Такаги – Сугено, отягощенный логическими обозначениями, является излишним.

Квазистатическая реляционная модель стадии диффузии производства сахара

Сахарное производство (см. [6, 7]) является сложным многостадийным производственным процессом. В работе [2] была предложена реляционная окрестностная модель стадии диффузии производства сахара и были приведены результаты параметрической идентификации этой модели в предположении линейности и на основе месячной выборки данных производства АО АПО «Аврора» «Боринский сахарный завод» с временным интервалом 10 минут, объем выборки n = 3472. При локальной идентификации регрессионных коэффициентов по данным, относящимся как к разным дням, так и к разным интервалам времени в течение одного дня, наблюдались значимые отклонения этих коэффициентов от глобальных, вычисленных по всей выборке. Кроме того, в некоторых случаях наблюдались изменения наборов значимых предикторов. В то же время все полученные локальные линейные регрессионные зависимости соответствовали технологии сахароварения, а изменения коэффициентов достаточно хорошо объяснялись изменениями параметров свекловичной стружки. Поэтому, сохраняя структуру реляционной окрестностной модели, предложенной в [2], мы рассмотрели квазистатическую линейную модель, в которой коэффициенты зависят от входных параметров: дигестии (сахаристости) S и длины L свекловичной стружки (в [2] эти параметры обозначались соответственно, Y4 (0,1) и Y2 (0,1)). На рис. 1 приведены точки рассеяния двумерного вектора входных данных, где L – длина 100 г стружки (в метрах), S – дигестия сахарной стружки (в процентах). Гистограмммы каждой из координат имеют по две достаточно хорошо заметных моды, и потому при кластеризации по методу k-средних мы полагали k = 4. Средние кластеров указаны в таблице; на рис. 2 четыре кластера выделены разными цветами.

schmir1.tif

Рис. 1. Точки рассеяния двумерного вектора входных данных

Средние кластеров

 

Кластер 1

Кластер 2

Кластер 3

Кластер 4

L

6,68

6,20

6,71

5,99

S

18,07

18,40

17,25

17,13

schmir2.tif

Рис. 2. Кластеры входных данных

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

Для кластера 1:

hmir110.wmf;

Для кластера 2:

hmir111.wmf;

Для кластера 3:

hmir112.wmf;

Для кластера 4:

hmir113.wmf.

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

Выводы

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

Работа выполнена при финансовой поддержке Российского фонда фундаментальных исследований (проект № 16-07-00854).