Россия
Россия
с 01.01.2018 по настоящее время
Государственный аграрный университет Северного Зауралья
с 01.01.2020 по настоящее время
Описывается концепция нового вида абстрактных автоматов «BP-схема» (Business Process), эффективно распределяющих управленческие команды между сущностями при их перманентном росте в парадигме сложно моделируемых систем. Исследована предлагаемая модель реализации BP-схемы, представляющая собой искусственную нейронную сеть, функционирование которой основано на правилах нечеткой логики с применением алгоритма Такаги-Сугено, для вывода правил в форме функциональных зависимостей. Модель использует QPU при вычислениях поведения системы, что предполагает быстрый отклик на изменения сущностей в системе. Использованный с учетом выработанных правил нечеткой нейронной сети позволяет совершенствовать сценарии и прогностические модели поведения сложной системы.
моделирование систем, система управления, искусственная нейронная сеть, квантовые технологии, абстрактные автоматы, BP-схема
Введение
Технические системы идут по пути усложнения, всего за несколько десятилетий был осуществлен переход от моно- к мульти (много) агентности. В виртуальной среде в эффективном тандеме друг с другом «работают» искусственные нейронные сети, в физическом мире разрабатываются и реализовываются концепции взаимодействия между собою таких сложных систем, как технологические линии, цеха и заводы. Причем, уже на этапе разработки закладываются возможности быстрой «перенастройки» подобных сложных систем с одного вида продукции на другой. Подобные вещи невозможно было сделать без возможности предварительного моделирования сложных систем [1, 2].
Методики моделирования систем
Моделирование систем, как научная область, в настоящее время обладает разнообразным методическим инструментарием, на рис. 1 показана его базовая классификация.
Рис. 1. Классификация методик моделирования систем
Fig. 1. Classification of system modeling techniques
Фундаментальный математический аппарат данной научной области, на прикладном уровне представлен следующими абстрактными схемами:
1. D-схемы, которые представляют собой непрерывно-детерминированные модели, отражающие моделируемую систему в динамике (ее изменение, в зависимости от не статичного параметра, которым в большинстве случаев является время).
2. F-схемы, которые реализовываются в виде математической абстракции конечного автомата со множеством регулируемых или учитываемых внутренних состояний, а также входных и выходных сигналов, которые в большинстве случаев представляются в виде конечного множества.
3. Р-схемы, которые реализовываются в виде вероятностных автоматов, через дискретно-стохастический подход, что позволяет сымитировать условный преобразователь информации с возможностью запоминания и, как следствие, корректировки последующей своей работы на основе предыдущего функционирования.
4. Q-схемы, являющиеся имитацией очерёдности элементов в целом и систем массового обслуживания в частности. При этом схема учитывает вероятность появления случайного «возмущения» в случайный момент времени, т.е. реализуется возможность учета стохастического характера процесса.
5. N-схемы, имитирующие трудноформализуемые системы, т.к. учитывают описания структуры и взаимодействия параллельных систем и процессов, с возможностью запоминания и обучения причинно-следственных связей в них.
6. А-схемы, которые являются всевозможными комбинациями вышеперечисленных схем.
Прикладная реализация данных схем приводит к усовершенствованию текущих технологических процессов и, как следствие, приводит к снижению логистических издержек, затрат ресурсов, повышению точности управления и иных сопутствующих явлений. Однако даже подобные подходы со временем станут неэффективными, поскольку представляемые производством задачи опережают потенциальные возможности абстрактной схемы.
Чаще всего подобные схемы оцениваются/сравниваются по критерию оптимальности затрат ресурсов на создание управленческих команд. Логично, что чем больше связей между сущностями, тем больше ресурсов затрачивается на создание подобных сигналов. Это в свою очередь означает необходимость создания технологии, которая способна осуществлять минимальное создание управленческих команд при увеличении сущностей.
В связи с этим предлагается концепция нового вида абстрактных автоматов «BP-схема» (Business Process), эффективно распределяющих управленческие команды между сущностями при их перманентном росте в парадигме сложно моделируемых систем.
Концептуальные особенности нового вида абстрактных автоматов
Постановка задачи. Описание и принцип действия BP-схемы рассматриваются применительно к сложной системе, которая состоит из n-сущностей, функционал которых заключается в реагировании на поступающие триггеры в диапазоне
Входные данные:
– Сложная система в двумерном пространстве, состоящая из связанных друг с другом n-сущностей.
– n-сущности, каждая из которых может прибывать в состояниях
– Триггеры n-сущностей в диапазоне
Требуется: разработать концептуальную модель BP-схемы; описать принцип действия BP-схемы в виде конечного автомата; сформулировать критерий оптимальности; провести модельные эксперименты и оценить эффективность предложенных решений, по сравнению с другими видами конечных автоматов.
Описание моделей и алгоритмов. В соответствии с постановкой задачи система рассматривается как сложная, а значит она состоит из множества взаимодействующих друг с другом сущностей. Взаимодействие сущностей приводит к появлению и приобретению системой новых свойств, которые отсутствуют на локальном инфраструктурном уровне. Это означает, что динамику развития подобной системы сложно моделировать ввиду перманентного изменения зависимостей/отношений, между сущностями и, как следствие, появления таких свойств, как нелинейность, эмерджентность, спонтанный порядок, адаптация и петли обратной связи. При визуализации подобных систем их изображают в виде сети, где узлы представляют компоненты и связи между ними. Также необходимо учитывать, что ряд сущностей может получать триггеры параллельно или одна сущность может получать одновременно несколько управленческих команд.
Для описания информационного потока триггеров и управленческих команд определим интервал времени, между их взаимодействиями с сущностями системы. Интервал времени определяется, как:
где
Логично, что данная формула расчетов будет адекватна при рассмотрении сущностей, которые находятся на одном маршруте информационного потока системы, иными словами, когда сущности системы линейно взаимосвязаны друг с другом. Это, в свою очередь, приводит к суждению о том, что сложную систему, состоящую из различных сущностей, взаимосвязь между которыми может вероятностно измениться в каждый момент времени, можно подвергнуть процедуре аппроксимации и представить ее в следующем виде:
Рис. 2. Сложная система до аппроксимации (слева) и после (справа)
Fig. 2. Complex system before approximation (left) and after (right)
Подобное упрощение системы, а именно понимание какие сущности взаимодействуют друг с другом в каждый момент времени дает возможность расчета
Учет вышеперечисленных моментов приводит к необходимости описать плотность распределения как:
И через нее выводим интенсивность как:
Временной интервал взаимодействия триггеров и управленческих команд с сущностью в системе определяется интервалом времени от момента начала взаимодействия до его завершения. При моделировании подобной системы необходимо учитывать приоритетность подобных информационных потоков для сущности. Количество подобных взаимодействий может быть неограниченно. Сущности не могут взаимодействовать с одним информационным потоком группой, но каждая сущность может создать на основе триггера или управленческой команды ряд других триггеров или управленческих команд на другие сущности, в этом случае сущность переходит в иную позицию. Всего у сущностей может быть три позиции:
1) исполнитель – когда сущность принимает триггеры или управленческие команды;
2) заказчик – когда сущность направляет триггеры или управленческие команды;
3) исполнитель/заказчик – комбинация из двух вышеперечисленных.
В итоге получается следующая схема:
а) б) в)
Рис. 3. Виды позиций сущностей в моделируемой системе:
а – сущность в позиции заказчик; б – сущность в позиции исполнитель; в – сущность в позиции исполнитель/заказчик
Fig. 3. Types of entity positions in the simulated system:
a – entity in the position customer; b – entity in the position of performer; c – entity in the position performer/customer
В данной системе информационный поток будет рассматриваться как однородный в связи с тем, что разный объем триггеров или управленческих команд, может быть обслужен любой сущностью.
В связи со всем вышеизложенным, предложена следующая схема и временная диаграмма моделируемой системы (рис. 4).
а) б)
Рис. 4. Схема моделируемой системы (а) и временная диаграмма (б)
Fig. 4. Diagram of the simulated system (a) and the timing diagram (b)
В соответствии с временной диаграммой, представленной на рис. 4, приходящие на сущности триггеры или управленческие команды могут выполняться параллельно. В силу того, что триггеры или управленческие команды различаются по степени важности (ранжируются), то длительность взаимодействия с ними и их место в очереди для сущности может изменяться.
Таким образом, требуют решения следующие задачи:
1) разработать методику корреткного ранжирования триггеров/управленческих команд;
2) определить алгоритм создания параллельных очередей триггеров/управленческих команд на одну сущность;
3) минимизировать время нахождения триггеров/ управленческих команд в параллельных очередях на одну сущность;
4) создать процедуру определения позиции одной сущности относительно других.
Для оценки эффективности решения поставленных задач сформированы следующие критерии:
1) количество обрабатываемых одной сущностью триггеров/управленческих команд за определенный временной интервал;
2) максимальная длинна очереди на одну сущность;
3) максимальное количество параллельных очередей на одну сущность;
4) равномерность распределения рангов триггеров/управленческих команд в очередях на одну сущность;
5) соотношение распределения состояний сущностей (исполнитель, заказчик, исполнитель/заказчик) внутри системы.
О пятом критерии поговорим подробнее, поскольку представляемая концепция нового вида автоматов «BP-схема», на прикладном уровне представляется авторами, как модель экономической системы взаимодействия друг с другом организаций-игроков, то необходимо обозначит следующее:
– При перенасыщении системы сущностями, прибывающими в состоянии «исполнитель» и «исполнитель/заказчик» система становиться негативно нестабильной.
– При перенасыщении системы сущностями, прибывающими в состоянии «заказчик» система становиться стабильно развивающейся.
В соответствии с вышеизложенными теоретическими положениями предлагается следующая формальная модель:
где
Опираясь на разработанную модель системы (1) мы можем представить BP-схему в виде графа (рис. 5).
Рис. 5. Граф автомата BP-схемы:
Fig. 5. The graph of the BP scheme automaton:
a, b – internal states; x – transition conditions; n – triggers/control commands; 1 – oracle
Если указанные выше элементы легенды графа классически, либо интуитивно понятны, «оракул» является уникальным для описания элементом. Естественно, необходимо объяснить данный инструмент, который является ключевым, для функционирования BP-схемы.
«Оракул» BP-схемы представляет собой искусственную нейронную сеть, функционирование которой основано на правилах нечеткой логики. Приходящие на «оракул» информационные потоки в виде триггеров/управленческих команд могут быть выражены в виде цифрового сигнала и формировать динамическую базу данных для обучения нейронной сети и выработки уникальных правил для каждой ситуации.
Нечёткая нейронная сеть (ННС) «оракула» для аппроксимации входящих данных по информационному потоку использует алгоритм Такаги-Сугено для вывода правил в форме функциональных зависимостей:
где
Структура ННС представлена на рис. 6.
Рис. 6. Структура ННС «оракула»:
Fig. 6. The structure of a neural network «oracle»:
Функциональные особенности каждого слоя ННС следующие:
1. Распределение входных данных по потокам триггеров/управленческих команд, их преобразование в нечёткие переменные и параллельное сравнение с имеющимися данными.
2. Определение метрик истинности для правил, уже имеющихся в базе данных.
3. Обучение ННС, результатом которого является выражение входных данных в виде одной объединенной функции принадлежности.
4. Составление списка выработанных правил.
5. Переход от функции принадлежности к ее измеряемому значению.
Итогом работы ННС является сценарий распределения потоков триггеров/управленческих команд на сущности, входящие в систему для поддержки корректного соотношения состояний сущности с целью сохранения системы в режиме стабильно развивающейся.
Как можно заключить, «оракул» является ключевым элементом BP-схемы и описание алгоритма его работы, а также последующее моделирование является превалирующим для логики данной работы.
Поведение сложной системы, которую BP-схема должна поддерживать в стабильно развивающемся состоянии функционирует в логике, что плотность потока триггеров/управленческих команд зависит от степени свободы сущности. Такими образом
Ввиду того, что триггеры/управленческие команды между сущностями должны распределяться в порядке параллельных очередей и не выстраиваться в одну последовательную очередь, то классические алгоритмы теории массового обслуживания в данном случае будут не эффективны. Поэтому для сохранения высоких метрик
Обобщённая схема данного алгоритма представлена на рис. 7.
Рис. 7. Обобщенная форма алгоритма Гровера для произвольного числа кубитов
Fig. 7. The generalized form of Grover's algorithm for an arbitrary number of qubits
При поступлении потока триггеров/«управленческих импульсов» в систему, «оракулу» необходимо осуществить расчет
В подобной постановке задачи оценка «оракулом» критерия
где
Аппроксимируем сложное пространство сущностей до уровня отрезка-взаимодействия между двумя сущностями. Подробно рассмотрим одно взаимодействие между двумя сущностями с потоком триггеров/управленческих команд. Как можно понять существуют следующие сценарии:
1) одна сущность в позиции «заказчик», другая в позиции «исполнитель»;
2) одна сущность в позиции «заказчик/исполнитель», другая в позиции «исполнитель»;
3) обе сущности в позиции «заказчик/исполнитель».
Поэтому на уровне полноценной системы, которая состоит из подобных отрезков-взаимодействий, возникает задача нахождения оптимального пути для информационного потока триггеров/управленческих команд от одной сущности к другой через ряд произвольных сущностей, пока информационный поток не будет исполнен.
Логично, что подобную задачу можно рассмотреть, как поисковую задачу на графе, а именно:
где
Сложную систему взаимодействия сущностей, можно выразить в виде сетки и, как следствие, в виде графа, где его ребра – это сущности, а маршруты информационных потоков между ними – это вершины. Отсюда следует, что для рассмотрения каждого маршрута информационного потока необходимо найти некую последовательность ребер графа. Очевидно, что существует m маршрутов от стартовой до финальной сущности и чем дальше они будут находиться друг от друга, тем больше данных маршрутов может быть.
Рассмотрим процедуру поиска оптимального маршрута из m маршрутов с использованием инструментов и терминов комбинаторики. Систему и пролагаемый по ней маршрут информационного потока от стартовой сущности до финальной представлена на рис. 8.
Рис. 8. Сложная система в виде пространства графа
Fig. 8. Complex system represented as a graph space
В этом случае стартовая позиция информационного потока (стартовая сущность, которая создаёт поток) и его конечная позиция (финальная сущность, которая поток выполняет) обязательны, а значит могут быть выражены, как
Примем, что мощность
Учитывая, что маршрут от стартовой сущности к конечной по пространству графа может быть длительным, но не бесконечным (т.к. количество сущностей ограничено и маршрут не может проходить более одного раза по одной сущности), то для поиска наиболее короткого маршрута информационного потока триггеров/управленческих команд необходимо:
1) найти все возможные комбинации упорядоченных установок без повторений;
2) вычислить возможную длину информационного потока от стартовой до конечной сущности;
3) разработать инструмент выбора маршрута с наименьшим количеством сущностей в промежутки от стартовой до конченой.
Наименьший маршрут находится по формуле:
где
Отсюда следует, что с помощью данного выражения можно осуществить нахождение наименьшего маршрута, превратив выражение (8) в систему уравнений:
Затем выразить в матричном виде и учесть величину
Итоги расчета (10) покажут, что одна из величин больше других и, как следствие, имеет больший маршрут/сроки выполнения сущностей триггера/управленческой команды.
Алгоритм распределения информационного потока представляется в виде следующей последовательности:
1. Ожидание триггера/управленческой команды. Если действие получено, то осуществляется переход на следующий этап.
2. Сущностью выполняется, либо не выполняется условие перехода
3. Оракулом производится расчет
4. Производится выбор наиболее короткого маршрута.
5. Происходит перерасчет режима системы. Если обнаружены отклонения от стабильно развивающегося режима, то «оркаул» производит процедуру корректировки состояний сущностей.
6. Осуществляется вывод правил в форме функциональных зависимостей и информация по ним поступает на «оракула»
7. После выполнения триггера/управленческой команды неактивные сущности переходят в режим ожидания.
Во время работы «оракулом» используется квантовый алгоритм Гровера (для сохранения высоких метрик
Результаты экспериментов
Для моделирования работы BP-схемы разработана программа для QPU (Quantum Processor Unit). Схематичное изображение модели, выраженное с помощью QPU представлено на рис. 9.
Рис. 9. BP-схема, выраженная с помощью средств QPU
Fig. 9. BP scheme expressed using QPU means
Моделировался процесс поиска (с помощью алгоритма Гровера) нахождения определённой (наиболее приоритетной) сущности, которая бы удовлетворяла вышеописанным критериям в неструктурированном пространстве сложной системы (неструктурированной базе данных).
Алгоритм работы данной модели выглядит следующим образом:
1. Инициализация кубитов: вначале появления триггера/управленческой команды, инициализируется пять кубитов, четыре из которых хранят информацию о критериях, по которым происходит поиск, а оставшийся – используется как «оракул» для фазового сдвига.
2. Переход кубитов в состояние «суперпозиции»: все кубиты критериев и оракула переводятся в суперпозицию с помощью преобразования Адамара. Кубит «оракула» также получает дополнительный фазовый сдвиг на 180 градусов (
3. Работа «оракула»: «оракул» производит изменение фазы состояния кубитов, отвечающих за критерии потока триггера/«управленческого импульса» с наивысшим приоритетом. При программировании данного алгоритма, оракул идентифицирует состояние |1101> как наиболее приоритетное и прикладывает к нему фазовый сдвиг.
4. Работа оператора диффузии: данный оператор усиливает вероятность измерения состояния, фаза которого была инвертирована оракулом, и уменьшает вероятность всех остальных. Это достигается путем преобразования Адамара, применения операции NOT ко всем кубитам, применения фазового сдвига и повторения преобразования Адамара и операции NOT.
5. Выполнение алгоритма Гровера: поисковый алгоритм повторяет комбинацию применения оракула и оператора диффузии, а именно
6. Измерение: после выполнения необходимого числа итераций алгоритма, кубиты измеряются для получения конечного результата.
7. Вывод результата: измеренное значение показывает состояние с наивысшим приоритетом. Результат измерения выводится в консоль. В двоичном виде это будет представлять собой сочетание критериев, которое «оракул» определил, как приоритетное (рис. 10).
Рис. 10. Визуализация работы BP-схема, по определению оптимальной сущности
Fig. 1. Visualization of the BP-scheme operation, defining the optimal entity
На данной визуализации мы видим, что изменились состояния кубитов №3 и №19. Это означает, что данные кубиты находятся в равной суперпозиции (это видно по площади заполнения круга) состояний ∣0〉∣0〉 и ∣1〉∣1〉. Положение стрелок на данных кубитах указывает сдвиг на 00 и 1800 соответственно. В контексте алгоритма Гровера, стрелка вниз указывает на кубиты, чья фаза была инвертирована оракулом — это соответствует маркировке правильного ответа. В данном случае это означает, что кубиты №3 и №19 помечены, как соответствующие условию поиска. В результате выполнения алгоритма получено число «11», это означает, что после измерения кубиты критериев находились в состоянии, которое соответствует двоичному представлению числа 3 (в двоичной системе «11» равно десятичному числу 3). Это означает, что некий третий условный элемент списка/сценария, по которому должна перестроиться BP-схема, как имитационная модель сложной системы является приоритетным, согласно правилам ННС, на основе которых выведены критерии для «оракула».
Заключение
Представленная в данной работе BP-схема может быть использована при имитационном моделировании сложных систем в области экономики или гибкого производства. Предлагаемое решение отличается использованием QPU в вычислениях поведения системы, а значит быстрым откликом на изменения состояния сущностей в системе. Помимо этого, использование поискового алгоритма Гровера с учетом выработанных правил ННС позволяет совершенствовать сценарии и прогностические модели поведения сложной системы. В перспективе авторы научного коллектива планируют исследовать поведение «оракула» BP-схем в сложных системах и определить дальнейшие способы его совершенствования.
1. Петров А.М., Попов А.Н. Разработка метода математического моделирования термодинамических процессов однофазных потоков наружных сетей теплоснабжения // Строительство и техногенная безопасность. - 2022. - № 26(78). - С. 59-63.
2. Петров А.М., Попов А.Н. Разработка интеллектуальной системы поддержки принятия решений по оценке состояния объектов системы теплоснабжения // Автоматизация и информатизация ТЭК. - 2023. - № 6(599). - С. 15-21.
3. Markov L., Saeedi M. Constant-optimized quantum circuits for modular multiplication and exponentiation. Quantum Information and Computation 12, 0361-0394, 2012.
4. Ambainis A.A nearly optimal discrete query quantum algorithm for evaluating NAND formulas, 2017.
5. Маслов И.О., Бачило А.О., Черкесова Л.В. Повышение быстродействия квантового алгоритма Гровера путем применения инверсии вокруг среднего // Молодой исследователь Дона. - 2019. - № 2(17). - С. 35-41.
6. Ulyanov S., Reshetnikov A., Tyatyushkina O. Modelling of Grover's quantum search algorithms: implementations of simple quantum simulators on classical computers // Системный анализ в науке и образовании. - 2020. - № 3. - С. 65-128.
7. A Fixed-Phase Quantum Search Algorithm with More Flexible Behavior / Х Li et al. // J. of Quantum Information Science. - 2012. - No 2. - Pp. 28-34.
8. Johansson N., Larsson J.А. Quantum simulation logic, oracles, and the quantum advantage // Entropy. - 2019. - Vol. 21, No. 8.
9. Оценка возможностей классических компьютеров при реализации симуляторов квантовых алгоритмов / С. В. Ульянов, Н. В. Рябов, П. В. Зрелов и др. // Программные продукты и системы. - 2022. - № 4. - С. 618-630.
10. Quantum Simulators: Architectures and Opportunities / E. Altman, I. Siddiqi, K. R. Brown et al. // PRX Quantum. - 2021. - Vol. 2, No. 1. - P. 017003.
11. Quantum Computer Systems for Scientific Discovery / Y. Alexeev, D. Bacon, K. R. Brown et al. // PRX Quantum. - 2021. - Vol. 2, No. 1. - P. 017001.
12. Gyongyosi L., Imre S. Circuit Depth Reduction for Gate-Model Quantum Computers // Scientific Reports. - 2020. - Vol. 10, No. 1. - P. 11229.
13. A Modified Quantum Search Algorithm / H. Mehri-Dehnavi, H. Dashtianeh, H. Ya. Kuchaksaraei et al. // International Journal of Theoretical Physics. - 2018. - Vol. 57, No. 12. - P. 3668-3681.
14. Химено-Сеговиа М., Хэрриган Н., Джонстон Э. Программирование квантовых компьютеров. Базовые алгоритмы и примеры кода. - СПб.: Питер, 2021. - 336 с.