ЧЕБЫШЕВСКИЙ СБОРНИК
Том 21. Выпуск 3.
УДК 519.217.1 DOI 10.22405/2226-8383-2020-21-3-129-141
Математическая модель цифровой системы управления с контроллерами фон Неймановского типа сложными многоконтурными объектами1
Е. В. Ларкин, А. Н. Привалов, Т. А. Акименко, И. Н. Лариошкин
Евгений Васильевич Ларкин — доктор технических наук, профессор, Тульский государственный университет (г. Тула). e-mail: elarkin@mail.ru
Александр Николаевич Привалов — доктор технических наук, профессор, Тульский государственный педагогический университет им. Л. Н. Толстого (г. Тула). e-mail: privalov.61@mail.ru
Татьяна Алексеевна Акименко — кандидат технических наук, доцент, Тульский государственный университет (г. Тула). e-mail: tantan72@mail.ru
Иван Николаевич Лариошкин — аспирант, Тульский государственный университет (г. Тула).
e-mail: tantan72@mail.ru
Аннотация
В работе построена математическая модель цифрового управления многоконтурными объектами, учитывающая реальные характеристики цифрового контроллера, как элемента системы управления. Сформулирована проблема, заключающаяся в том, что методы моделирования цифровых систем управления известны и широко применяются в практике инженерной деятельности, однако в подавляющем большинстве они предполагают формирование моделей, не учитывающих наличие временных интервалов между транзакциями в ЭВМ фон Неймановского типа.
Для решения задачи разработана типовая структурная схема сложных многоконтурных систем управления с цифровыми контроллерами фон Неймановского типа, которая учитывает случайный характер обрабатываемых данных и реальные временных задержки между транзакциями.
Предложено с учетом случайности временного интервала между транзакциями и стохастического характера переключения в сопряженные операторы считать адекватной моделью алгоритма функционирования цифровых систем управления считать полумарковский процесс.
На основе полумарковских процессов предложен метод оценки параметров временных интервалов между транзакциями в циклических алгоритмах управления, который позволяет оценить характеристики системы на этапе ее проектирования, а следовательно является ключом к рациональному проектированию цифровых систем управления многоконтурными объектами с алгоритмами управления практически любой сложности. Представлен пример математического моделирования двухконтурной системы с цифровым управлением.
Ключевые слова: полумарковский процесс, время блуждания, цифровая система управления, контроллер, алгоритм управления цифровой системой управления, транзакция.
Библиография: 21 название.
1 Исследование выполнено при финансовой поддержке РФФИ по проекту 19-47-710004_р_а.
Для цитирования:
Е. В. Ларкин, А. Н. Привалов, Т. А. Акименко, И. Н. Лариошкин. Математическая модель цифровой системы управления с контроллерами фон Неймановского типа сложными многоконтурными объектами // Чебышевский сборник, 2020, т. 21, вып. 3, с. 129-141.
CHEBYSHEVSKII SBORNIK Vol. 21. No. 3.
UDC 519.217.1 DOI 10.22405/2226-8383-2020-21-3-129-141
Mathematical model of a digital control system with background controllers of the Neuman type for complex multicirculated objects
E. V. Larkin, A. N. Privalov, T. A. Akimenko, I. N. Larioshkin
Evgenii Vasil'evich Larkin — Doctor of technical sciences, Professor, Tula State University (Tula).
e-mail: elarkin@mail.ru
Alexander Nikolaevich Privalov — Doctor of technical sciences, Professor, Tula State Pedagogical University (Tula). e-mail: privalov.61@mail.ru
Tatyana Alekseevna Akimenko — Candidate of Technical Sciences, Associate Professor, Tula State University (Tula). e-mail: tantan72@mail.ru
Ivan Nikolaevich Larioshkin — postgraduate student, Tula State University (Tula). e-mail: tantan72@mail.ru
Abstract
In the work, a mathematical model of digital control of multi-circuit objects is built, taking into account the real characteristics of a digital controller as an element of a control system. The problem is formulated that the methods of modeling digital control systems are known and are widely used in engineering practice, however, in the overwhelming majority, they involve the formation of models that do not take into account the presence of time intervals between transactions in a Von Neumann type computer.
To solve the problem, a typical block diagram of complex multi-loop control systems with digital controllers of the Von Neumann type has been developed, which takes into account the random nature of the processed data and real time delays between transactions.
It is proposed, taking into account the randomness of the time interval between transactions and the stochastic nature of switching to conjugate operators, to consider a semi-Markov process as an adequate model of the algorithm for the functioning of digital control systems.
On the basis of semi-Markov processes, a method is proposed for estimating the parameters of time intervals between transactions in cyclic control algorithms, which makes it possible to evaluate the characteristics of the system at the design stage, and therefore is the key to the rational design of digital control systems for multi-circuit objects with control algorithms of almost any complexity. An example of mathematical modeling of a two-circuit system with digital control is presented.
Keywords: semi-Markov process, wandering time, digital control system, controller, digital control system control algorithm, transaction.
Bibliography: 21 titles. For citation:
E. V. Larkin, A. N. Privalov, T. A. Akimenko, I. N. Larioshkin, 2020, "Mathematical model of a digital control system with background controllers of the Neuman type for complex multicirculated objects", Chebyshevskii sbornik, vol. 21, no. 3, pp. 129-141.
1. Введение
Основной тенденцией развития современной промышленности является усложнение технологических процессов и связанное с этим широкое внедрение цифровой управляющей техники систем управления для обеспечения требуемых показателей качества производимой продукции [1, 2]. Переход при управлении объектами на цифровые системы связан с включением в их структуру контроллеров фон Неймановского типа, т.е. приборов, которые обладают новыми, по сравнению с аналоговыми контроллерами, свойствами [3, 4, 5]. Свойства обусловлены последовательной интерпретацией операторов алгоритма управления, развивающейся в реальном физическом времени, за счет чего управляющая ЭВМ, кроме того, что алгоритмически реализует заданный закон управления, вносит временные задержки в процесс управления, в связи с чем, временной фактор напрямую влияет на качественные характеристики системы в целом. Алгоритм обрабатывает случайные данные, формируемые на выходах сенсоров, и содержит операторы принятия решения в местах ветвления, поэтому в результате его интерпретации генерируются случайные временные интервалы между транзакциями, как при вводе-выводе данных (перекос данных), так и при расчете управляющего воздействия (чистое запаздывание) [5, 6, 7]. В то же время промышленные, военные и иные системы при переходе на цифровые технологии должны обеспечивать требуемую точность функционирования в широких диапазонах изменений параметров объекта управления, независимо от случайных временных задержек в контурах управления [8, 9, 10]. Поэтому возникает научная проблема создания моделей управления сложными объектами, учитывающих реальные физические характеристики управляющих ЭВМ, которые могли бы быть положены в основу методологии синтеза цифровых регуляторов.
Методы моделирования цифровых систем управления известны [11, 12, 13] и широко применяются в практике инженерной деятельности, однако в подавляющем большинстве они предполагают формирование моделей, не учитывающих наличие временных интервалов между транзакциями в ЭВМ фон Неймановского типа. Это, в свою очередь, неоправданно увеличивает сроки и объемы экспериментальных исследований при отработке систем управления и, в частности, приводит к необходимости дополнительных экспериментальных исследований систем при переводе отработанных алгоритмов управления на другую элементную базу, например на контроллеры с другим быстродействием, или контроллеры с другой системой команд.
Учет временного фактора при проектировании цифровых систем управления развит недостаточно, что обусловливает необходимость и актуальность исследований в данной области.
2. Структура цифровой системы
Структурная схема цифровой системы управления приведена на рис. 1. Блоки структурной схемы Фц (s), ..., Фи (s), ..., Фкк (s) описывают динамические процессы, протекающие в объекте управления, через передаточные функции, Перекрестные межблочные связи, когда k = I, отражают реальное взаимовлияние контуров в объекте. Блоки Ф01 (s), ..., Фо& (s), ..., Фок (s) описывают динамику сенсоров системы.
Система функционирует следующим образом. В цифровой контроллер (Digital controller) вводится вектор сигналов F (s) = [Fi (s) , ..., Fk (s) , ..., Fk (s)j, формирующих цель управления. Указанные сигналы могут также генерироваться самим контроллером, что показано на рис. 1 жирной штриховой линией. На выходах контроллера программно генерируются последовательности кодов е (s) = [e1 (s), ..., е^ (s), ..., £к (s)j, подаваемые на входы объекта управления, который на физическом уровне преобразует сигналы вектора е (s) в вектор параметров X (s) = [Xi (s), ..., Х^ (s), ..., Хк (s)j, описывающих состояние объ-
Рис. 1: Структурная схема цифровой системы управления.
екта. Параметры измеряются сенсорами, и результаты измерения в виде цифровых кодов Хо (в) = [Хо\ (в), ..., Хок ($), ..., Хок (з)] вводятся в контроллер. Во всех вышеприведенных векторах я — оператор дифференцирования преобразования Лапласа [8, 14].
Выходная последовательность кодов генерируется в соответствии с законом управления и циклограммой, заложенной в алгоритмическое и программное обеспечение контроллера. В дальнейшем будем считать, что временные интервалы в циклограмме отсчитываются от момента ввода первого элемента вектора Р (в), Р\ (в). Остальные элементы вводятся/выводятся с задержками по времени обозначенными следующим образом: элемент ^ (в) вводится с задержкой т/^, Tf 0 = 01 ^ к ^ К; элемент Х0^ (в) вводится с задержкой т0к, 1 ^ к ^ К; элемент ей (я) выводится с задержкой 1 ^ к ^ К. В соответствии с теоремой о смещении во временной области [14]
Цф - т)] = ехр(-тв)Ф(в), т > 0, (1)
где т — величина смещения; £ — время; р — импульсный отклик; Ь [...] — прямое преобразование Лапласа.
Из (1) следует, что
^ («) = ^ («) ■ Nf (8); (2)
Х8Н (з) = Хо (8) ■ N0 (8); (3)
(«) = £ («) ■ Х£ («), где Nf (,в), N0 (в), Хе (в) — матрицы смещений, равные
( 1 ... 0
(4)
М; (з) =
0
ехр (-т^к
... 0 ... 0
(5)
\ 0 ... 0 ... ехр (—т^к8) )
( ехр (—Г01«) ... 0 ... 0 \
N0 (8) =
0
ехр (—Т0к «)
V
/ ехр (—Те1£)
м£ (з) =
ехр (-Тек«)
V
0 ...
0
0
... ехр (-Т0К«) )
\
... ехр (-т£к в) )
(6)
(7)
В том случае, если цифровой контроллер реализует линейный закон управления Фс (з), на исполнительные приводы поступает сигнал
£ (5) = [Р (з) ■ ^ (з) - X (з) ■ Ф0 (5) ■ N0 (5)] ■ Фс (з) ■ Хе (з),
(8)
где Фс (в) = [Фск1 (8)] , 1 ^ к,1 ^ К — матрица линейных передаточных функций, реализуемых в цифровом контроллере.
Вектор X (в) определяется как
X (5) = £ (5) ■ Ф (в) = [^ (5) ■ Щ (з) - X (з) ■ Ф0 (5) ■ N0 (5)] ■ Фс (5) ■ М£ (з) ■ Ф (в)
(9)
где Ф (в) = [Фы (з)] , 1 ^ к,1 ^ К — матрица линейных передаточных функций, реализуемых в объекте управления.
Разрешение (9) относительно X (в) дает
X (5) = [Е - Ф0 (5) ■ Фс (5) ■ М0е (в)]-1 ^ (в) ■ Фс (в) ■ (з) ■ Ф (в)
(10)
Где е — единичная диагональная матрица;
( ехр (-Те13)
(з) = Щ (з) ■ Хе (з) = 0
... 0
\
V
0 ... 0 ...
ехр [- (т{к + Тек) «] ... 0
0 ... ехр [- (т1ь + теК) 8] /
(11)
Х0е (8) = N0 (з) ■ М£ (з) = / ехр [- (Г01 + т£1) 8] ... 0
0 ... 0 ...
... 0
ехр [- (Т0к + Тек) «] ... 0
0 ... ехр [- (Т0К + ТеК) «] /
0
0
0
0
0
0
Члены, характеризующие чистое запаздывание, появляются как в числителе, Nfs (s), так и в знаменателе, Ще (s), (10). Запаздывание в числителе определяет т.н. «перекос данных» (data screw) и общее запаздывание выполнения объектом внешних команд. Запаздывание в знаменателе изменяет значения полюсов характеристического многочлена, а следовательно меняет качественно характеристики переходных процессов.
3. Оценка временных интервалов
Для оценки временных интервалов задержки должны быть построена модель интерпретации алгоритма управления контроллером фон Неймановского типа. Для простоты, без нарушения общности, модель может быть построена только из операторов транзакций. Поскольку управление сводится к периодическому опросу элементов векторов F (s) = [Fi (s), ..., Fk (s), ..., FK (s)j, Хо (s) = [Xoi (s), ..., Хо'к (s), Хок (s)j, и e (s) = [ei (s), ..., ek (s) , ..., £k (s)j, алгоритм является циклическим, но в нем отсутствует эффект зацикливания. Кроме того, поскольку для расчета управляющего воздействия е (s) все элементы векторов F (s^ Хо (s) должны быть использованы, на структуру алгоритма накладываются условия сильной связности [15]. Подобными свойствами в самом общем случае обладает структура, представленная полным графом с петлями (рис. 2).
С учетом случайности временного интервала между транзакциями и стохастического характера переключения в сопряженные операторы, адекватной моделью алгоритма является полумарковский процесс [16, 17, 18, 19], представленный полумарковской матрицей
h (t) = [hktl (i)j = [gktl (i)j <g> [pk>l (i)j, (13)
где pk,i (t) — вероятности прямого переключения из к-го оператора транзакции в 1-й; gk,i (t) — плотность распределения времени прямого переключения из к-го оператора транзакции в 1-й; ® — знак прямого произведения матриц; t — физическое время.
В силу названных свойств моделируемого алгоритма, полумарковский процесс (13) является эргодическим, и в нем нет поглощающих и частично поглощающих состояний.
Рис. 2: Общая структура полумарковского процесса.
В силу эргодичности полумарковского процесса на плотности распределения gkj (i) и вероятности pk i (t) накладываются следующие ограничения:
0 < < arg [дк>1 (i)j < Тк*Г < 1 < k,l < 3К; (14)
3К
= 1; (15)
Timax
k=1
9k,i (í) dt = 1, (16)
где 3 К — общее количество операторов транзакций; Т™" и T^i* — нижняя и верхняя границы области определения плотности распределения gk,i (t).
При оценке временных интервалов между транзакциями в алгоритме не имеет значения, каким образом полумарковский процесс (13) попал из к-го состояния в 1-е. Определяющим в данном случае является то ,что это попадание является первым, а не вторым, третьим и т.п. Для оценки исходный полумарковский процесс (13) должен быть преобразован в процесс со структурой, приведенной на рис. 2 Ь, в котором к-е состояние является стартовым, а 1-е состояние — поглощающим. Для этого:
1-я строка и к-й столбец матрицы h (t) обнуляются;
вероятности pk,i (t) во всех строках, исключая 1-ю, и во всех столбцах, исключая к-т, пе-ресчитываются по зависимости:
р'т,п = 3К Рт'П-, 1 < т,п < 3К, т = к, п = I. (17)
п = 1 Рт,п П = I
Таким образом
h (t) ^ h (t) = [9k,i (t) ■ p'^]. (i8)
После пересчета вероятностей по зависимости (17) частично поглощающие состояния, т.е. состояния, для которых не выполняется свойство (15), уничтожаются, и события достижения 1-го состояния из к-е по одной из возможных траекторий блуждания начинают составлять полную группу несовместных событий [20]. Таким образом, плотность распределения времени достижения 1-го состояния из к-го определяется как
9ы (í) = 1'к ■ L
-1
3=1
(19)
где L [...^ L-1 [...] — прямое и обратное преобразование Лапласа; I' — вектор-
строка, к-й элемент которого равен единице, а остальные элементы равны нулю; знак транспонирования; If — вектор-столбец, 1-й элемент которого равен единице, а остальные элементы равны нулю [21].
Для плотности распределения (19) могут быть рассчитаны математическое ожидание и дисперсия по зависимостям [20]:
/■те
Tli = * ■ 9ы (t)dt; (20)
uk,l
f те 2
/ tt-T^) -oh (t)dt. (21)
J 0
Состояния полумарковского процесса (13) могут быть распределены между транзакциями, осуществляемые в соответствии со схемой, приведенной на рис. 1, следующим образом:
состояния с первого по К-е являются математическими аналогами операторов ввода элементов вектора Р (в)-,
состояния с (К+1)-го по (2К)-е моделируют ввод элементов вектора Хо («);
Ч ЦШП 1k,l
If
0
состояния с (2/Г+1)-го по (3К)-е моделируют вывод элементов вектора е («).
Оценка временного интервала может быть оценена либо по математическому ожиданию
(20), либо по максимальному значению [20] с использованием соотношения
^ = + (22)
Соответственно, для матриц (5), (6), (7) временные задержки будут определяться, как это показано в таблице 1:
Вид транзакции Номер вводимого элемента Среднее значение времени Максимальное значение времени
Ввод ^ («) 1 V1 =0 V1 = 0
2 < к < К ^к = т1к, Ук = + у
Ввод Х0 («) 1 < к < К ток = , ток = ^1к+к + 3\! Е^к+к
Вывод £ («) 1 < к < К т£ к = Т^2К+к, Т£к = Т^к+к + 3\/ В\,2К+к
Таблица 1: Оценки временных интервалов.
При учете временных задержек, приведенных в третьей колонке таблицы 1, оцениваются показатели качества управления в среднем. При учете временных задержек, приведенных в четвертой колонке, оцениваются предельные показатели качества управления.
4. Пример анализа системы управления
В качестве примера исследована простейшая двухконтурная система. В системе реализована единичная обратная связь по обоим контурам, а в цифровом контроллере реализован П-закон управления.
Передаточные функции, определяющие динамику объекта управления, равны:
Ф11 («) = Ф22 (5) = Ф12 (5) = Ф21 (в) = ^. (23)
Между транзакциями обеспечиваются следующие временные задержки:
1 = 0; 2 = 0, 05; 701 = 0,1; то2 = 0,15те1 = 0, 2; те2 = 0, 25. (24)
Характеристическое уравнение, определяющее характер переходных процессов, без учета погрешностей, вносимых цифровым контроллером, имеет вид:
+ 38в3 + 480з2 + 22728 + 3885 = 0. (25)
Характеристическое уравнение, определяющее характер переходных процессов, с учетом погрешностей, вносимых цифровым контроллером, имеет вид:
+ 5 + 10 ехр (-0, 38)] ■ + 3 + 12 ехр (-0, 4в)| ■ (¿2 + 88 + 16) -- exp (-0, 78) ■ (в2 + + 15) = 0.
(26)
Структура простейшей двухконтурной системы приведена на рис. 3 а, с цифровым контроллером, в котором реализован алгоритм со структурой, показанной на рис. 3 Ь и ее графическое решение Зс.
Р^) ^^ Г2(Б) 1 )-—Н 2 )--к 3
0 0
Х2к)
н-( 4
-Х1(5)|
V.
2 3
9 10
Рис. 3: Двухконтурная система (а); циклограмма управления (Ь); графическое представление (с).
Решение характеристических уравнений (25) и (26) дает следующие значения полюсов: без учета погрешностей
= -14 в2 = -16 8з = —4 + 0, 091г; 84 = —4 — 0, 091г, где г — мнимая единица, с учетом погрешностей
«1 = 23, 5; 82 = —21, 7; 8з = 0, 79 + 21, 5г; 84 = 0, 79 — 21, 5г, где г — мнимая единица.
Наличие пары комплексно-сопряженных корней:
с положительной действительной частью указывает на то, что в замкнутой системе возникает расходящийся периодический переходный процесс; наличие действительного положительного корня, в дополнение к расходящемуся периодическому процессу свидетельствует о неустой чивости замкнутой системы;
а
Ь
с
— с отрицательной действительной частью указывает на то, что в замкнутой системе возникает затухающий периодический переходный процесс, который обусловлен исключительно наличием перекрестах связей. При отсутствии перекрестных связей в объекте управления система распадается на два независимых контура управления объектами, которые описываются звеньями первого порядка, т.е. затухающая периодическая составляющая в них отсутствует.
Таким образом, замыкание обратных связей устойчивой аналоговой системы через цифровой контроллер при той же логике работы приводит к качественным изменениям характера переходных процессов, что необходимо учитывать при проектировании цифровых систем управления.
5. Заключение
В результате построена математическая модель цифрового управления многоконтурными объектами, учитывающая реальные характеристики цифрового контроллера, как элемента системы управления. Показано, что контроллеры фон Неймановского типа с последовательной интерпретацией команд алгоритма управления являются источником временных задержек, которые приводят таким явлениям, как перекос данных и чистое запаздывание в контуре управления. В свою очередь, и перекос данных и чистое запаздывание приводят к ухудшению таких качественных параметров управления, как время переходного процесса, перерегулирование и колебательность. Предложенный метод расчета временных интервалов позволяет оценить характеристики системы на этапе ее проектирования, а следовательно является ключом к рациональному проектированию цифровых систем управления многоконтурными объектами с алгоритмами управления практически любой сложности.
Дальнейшие исследования в этой области могут быть направлены на разработку алгоритмов цифрового управления, оптимальных по соотношению показателей сложность/качество, и обеспечивающих требуемые характеристики функционирования цифровых систем.
СПИСОК ЦИТИРОВАННОЙ ЛИТЕРАТУРЫ
1. Malin Löfving \!.. Säfsten К.; Winroth М. Manufacturing strategy formulation, leadership style and organizational culture in small and medium-sized enterprises // IJMTM. 2016. Vol. 30. no. 5, pp. 306-325.
2. Zhou M.C. Modeling, Simulation, and Control of Flexible Manufacturing Systems: A Petri Net Approach // Wspc. 1999, 428 p.
3. Landau I.D., Zito G., Digital Control Systems, Design, Identification and Implementation // Springer. 2006, 484 p.
4. Aström J., Wittenmark В., Computer Controlled Systems: Theory and Design // Tsinghua University Press. Prentice Hall. 2002, 557 p.
5. Larkin E., Bogomolov A., Privalov A., Estimation of Events Flow Generated with Ergodic Semi-Markov Processes // 2nd International Ural Conference on Measurements (URALCON), South Ural State University, Chelyabinsk, Russia, 2017, pp. 124-129.
6. Fadali M.S., Visioli A., Digital control engineering: Analysis and design // Elsevier Inc. 2013, pp. 239 - 272.
7. Auslander D.M., Ridgelv J.R. Jones J.C., Real-time software for implementation of feedback control // The control handbook. Control system fundamentals. Ed. by W.S. Levine. CRC Press. Taylor and Francis Group, US. 2017, pp.16-32.
8. Karnopp D.C., Margolis D.L., Rosenberg R.C., System dynamics: Modeling, simulation and control of mechatronic systems // John Willev k, Sons. New Jersey. 2012, 636 p.
9. Tzafestas S.G. Introduction to Mobile Robot Control. Elsevier, 2014. 750 p.
10. Babishin V., Taghipour S. Optimal maintenance policy for multicomponent systems with periodic and opportunistic inspections and preventive replacements // Applied Mathematical Modelling 2016, V. 40. no. 23, pp. 10480-10505.
11. Arnold K. A. "Timing analysis in embedded systems", In embedded hardware by J. Ganssler, K. Arnold et all, USA, Elsevier Inc. 2008. pp. 239 - 272.
12. Balsamo S., Harrison P.G., Marin A., "Methodological construction of product-form stochastic Petri nets for performance evaluation" //In Journal of Systems and Software. Elsevier Inc., 2012, Vol. 85, no. 7. pp. 1520-1539.
13. Hamann A., Racu, R., Ernst R., "Multi-dimensional robustness optimization in heterogeneous distributed embedded systems" //In Proceedings of the 13th IEEE Real Time and Embedded Technology and Applications Symposium, RTAS '07, IEEE Computer Society, Washington, DC, USA.'2007, pp. 269-280.
14. Schiff J.L., The Laplace transform: Theory and applications, USA, NY Springer Verlag, 199, 233 p.
15. Larkin E.V., Bogomolov A.V., Privalov A.N., "A Method for Estimating the Time Intervals between Transactions in Speech-Compression Algorithms", In: Automatic Documentation and Mathematical Linguistics. 2017, Vol. 51. Iss. 5. pp. 214-219.
16. Bielecki T.R., Jakubowski J., Niew§glowski M., "Conditional Markov chains: Properties, construction and structured dependence", In Stochastic Processes and their Applications. 2017, V. 127, no. 4. pp. 1125-1170.
17. Ching W.K., Huang X., Ng M.K., Siu T.K., "Markov Chains: Models, Algorithms and Applications", In International Series in Operations Research k, Management Science, V. 189. Springer Science + Business Media NY, 2013, 241 p.
18. Howard R.A., Dynamic Probabilistic Systems. Vol. 1: Markov Models. Vol. II: Semi-Markov and Decision Processes // Courier Corporation. 2012, 205p.
19. Janssen J., Manca R., Applied Semi-Markov processes. Springer US. 2006, 310 p.
20. Kobavashi H., Marl B.L., Turin W., Probability, Random Processes and Statistical Analysis: Cambridge University Press. 2012, 812 p.
21. Petersen P.,Linear algebra. Springer-Verlag. NY. 2012, 427 p. REFERENCES
1. Malin Lofving M., Sáfsten K.; Winroth M., 2016, Manufacturing strategy formulation, leadership style and organizational culture in small and medium-sized enterprises, IJMTM. Vol. 30. no. 5, pp. 306 - 325.
2. Zhou M.C., 1999, Modeling, Simulation, and Control of Flexible Manufacturing Systems: A Petri Net Approach. Wspc, 428 p.
3. Landau I.D., Zito G., 2006, Digital Control Systems, Design, Identification and Implementation. Springer, 484 p.
4. Aström J., Wittenmark B., 2002, Computer Controlled Systems: Theory and Design. Tsinghua University Press. Prentice Hall, 557 p.
5. Larkin E., Bogomolov A., Privalov A., 2017, Estimation of Events Flow Generated with Ergodic Semi-Markov Processes // 2nd International Ural Conference on Measurements (URALCON), South Ural State University, Chelyabinsk, Russia, pp. 124-129.
6. Fadali M.S., Visioli A., 2013, Digital control engineering: Analysis and design. Elsevier Inc., pp. 239 - 272.
7. Auslander D.M., Ridgelv J.R. Jones J.C., 2017, Real-time software for implementation of feedback control // The control handbook. Control system fundamentals. Ed. by WT.S. Levine. CRC Press. Taylor and Francis Group, US. pp.16-32.
8. Karnopp D.C., Margolis D.L., Rosenberg R.C., 2012, System dynamics: Modeling, simulation and control of mechatronic systems. John WTillev k, Sons. New Jersey, 636 p.
9. Tzafestas S.G. Introduction to Mobile Robot Control. Elsevier, 2014. 750 p.
10. Babishin V., Taghipour S.. 2016, Optimal maintenance policy for multicomponent systems with periodic and opportunistic inspections and preventive replacements // Applied Mathematical Modelling, V. 40. no. 23, pp. 10480 -10505.
11. Arnold K. A. "Timing analysis in embedded systems", In embedded hardware by J. Ganssler, K. Arnold et all, USA, Elsevier Inc. 2008. pp. 239 - 272.
12. Balsamo S., Harrison P.G., Marin A., 2012, "Methodological construction of product-form stochastic Petri nets for performance evaluation", In Journal of Systems and Software. Elsevier Inc., Vol. 85, no. 7. pp. 1520-1539.
13. Hamann A., Racu, R., Ernst R., 2007, "Multi-dimensional robustness optimization in heterogeneous distributed embedded systems", In Proceedings of the 13th IEEE Real Time and Embedded Technology and Applications Symposium, RTAS '07, IEEE Computer Society, Washington, DC, USA. pp. 269-280.
14. Schiff J.L., 199, The Laplace transform: Theory and applications, USA, NY Springer Verlag, 233 p.
15. Larkin E.V., Bogomolov A.V., Privalov A.N., 2017, "A Method for Estimating the Time Intervals between Transactions in Speech-Compression Algorithms", In: Automatic Documentation and Mathematical Linguistics, Vol. 51. Iss. 5. pp. 214-219.
16. Bielecki T.R., Jakubowski J., Niew§glowski M., 2017, "Conditional Markov chains: Properties, construction and structured dependence", In Stochastic Processes and their Applications, V. 127, no. 4. pp. 1125-1170.
17. Ching WT.K., Huang X., Ng M.K., Siu T.K., 2013, "Markov Chains: Models, Algorithms and Applications", In International Series in Operations Research k, Management Science, V. 189. Springer Science + Business Media NY,. 241 p.
18. Howard R.A., 2012, Dynamic Probabilistic Systems. Vol. 1: Markov Models. Vol. II: Semi-Markov and Decision Processes, Courier Corporation, 205p.
19. Janssen J., Manca R., 2006, Applied Semi-Markov processes. Springer US, 310 p.
20. Kobavashi H., Marl B.L., Turin W., 2012, Probability, Random Processes and Statistical Analysis: Cambridge University Press.. 812 p.
21. Petersen P., 2012, Linear algebra. Springer-Verlag. NY. 427 p.
Получено 12.06.2020 г. Принято в печать 22.10.2020 г.