ВАК 05.13.01 Системный анализ, управление и обработка информации (по отраслям)
ВАК 05.13.06 Автоматизация и управление технологическими процессами и производствами (по отраслям)
ВАК 05.13.10 Управление в социальных и экономических системах
ВАК 05.13.18 Математическое моделирование, численные методы и комплексы программ
ВАК 05.13.19 Методы и системы защиты информации, информационная безопасность
ГРНТИ 20.01 Общие вопросы информатики
ГРНТИ 28.01 Общие вопросы кибернетики
ГРНТИ 49.01 Общие вопросы связи
ГРНТИ 50.01 Общие вопросы автоматики и вычислительной техники
ГРНТИ 82.01 Общие вопросы организации и управления
К рассмотрению предлагается разработка математической модели процесса функционирования системы управления информацией и событиями безопасности, известной как SIEM-система. Данная модель представляет собой формализованное аналитическое описание (в терминах марковской цепи в форме разностных стохастических уравнений) динамики смены состояний показателей качества, характеризующих существенные свойства процесса функционирования системы управления информацией и событиями безопасности, в пространстве состояний. Модель представляет собой традиционную для цепи Маркова, в форме конечных разностей, систему уравнений состояния и наблюдения. Научная задача состоит в усовершенствовании (модификации) алгоритмов преобразования шума возбуждения, используемого в модели. Предложен механизм определения значений приращения математического ожидания моделируемого процесса, полученных на основе априорных данных о цепи Маркова, по отношению к математическому ожиданию белого гауссовского шума, возбуждающего этот процесс. Этот механизм позволяет, опираясь на несложные вычисления, принимать решение о том, какие значения будут принимать элементы вектора компенсационных добавок в уравнении состояния вспомогательного вектора индикаторов данной модифицированной модели с учетом преобразования шума возбуждения. Это позволяет упростить модель, снизить ее вычислительную сложность без существенных потерь в точности (адекватности). Практическое применение усовершенствованной модели возможно как в рамках исследовательских работ, так и в системах автоматизированного контроля информационной безопасности.
математическое ожидание, система управления информацией и событиями безопасности, показатель качества, процесс функционирования, матрица, состояние.
Введение
Системы информационной безопасности компьютерных и телекоммуникационных сетей имеют сложную интегрированную структуру. Они включают в себя комплексы аппаратных и программных средств (серверов, баз данных, средств защиты, автоматизированных рабочих мест), а также персонал. Все эти компоненты организационно и технически объединены и предназначены для реализации процесса защиты информации. Цель этого процесса – минимизация потерь, которые могут быть вызваны нарушением доступности, целостности или конфиденциальности информации [1–4].
Одним из возможных и эффективных подходов к достижению данной цели является разработка и применение систем управления информацией и событиями безопасности (SIEM – Security Information and Event Management). Они предназначены для анализа в реальном времени событий (угроз, тревог) безопасности, исходящих от сетевых устройств и приложений, и позволяют реагировать на эти угрозы до наступления существенного ущерба [5–7].
Актуальность и практическая значимость исследований в рамках данной темы заключаются в формулировке новой модели для описания процесса функционирования SIEM-систем. Это позволит детально аналитически описать работу основного технического решения по обеспечению контроля информационной безопасности (ИБ) – системы управления событиями и инцидентами ИБ, SIEM-системы. Они анализируют информацию, поступающую от различных подсистем ИБ (управления доступом, антивирусной защиты, защиты межсетевого взаимодействия, анализа защищенности, систем обнаружения и предотвращения вторжений, контроля целостности и др.) и выявляют отклонения состояния ИБ по различным критериям. При выявлении отклонения в состоянии ИБ реализуется алгоритм, называемый «управление событиями безопасности» – в подсистеме управления событиями SIEM-системы создается «уведомление», и оповещаются заинтересованные лица. На основе этого оповещения реализуются технические или организационные управленческие решения по предотвращению угроз и противодействию им. Таким образом, SIEM-системы являются инструментом централизованного просмотра и обработки информации, регистрируемой на большом количестве источников за счет графического представления, фильтрации и группировки данных. Выявление отклонений позволяет своевременно и в автоматизированном режиме выявлять угрозы ИБ или предпосылки к их возникновению с учетом совокупности регистрируемых событий и собираемых данных.
Традиционным подходом к аналитическому описанию процессов функционирования сложных систем, например таких, как SIEM-система, являются попытки формально математически представить данные процессы через динамику линейного или нелинейного изменения состояния (пошаговых значений) показателей качества (ПК) SIEM-системы и процесса ее функционирования в пространстве состояний. Марковские модели процессов в этом случае являются незаменимым инструментом.
Под показателями качества SIEM-системы понимается количественная характеристика одного или нескольких свойств этой системы, обуславливающих ее качество. Эта количественная характеристика может иметь векторный вид (вектор ПК), обозначаться как вектор, например, . Здесь компоненты вектора ПК
представляют собой численные значения одного или нескольких параметров, описывающих конкретное свойство SIEM-системы, а также процессов функционирования различных ее элементов или системы в целом. При этом ПК как количественная характеристика одного или нескольких свойств системы, обуславливающих ее качество, рассматривается применительно к определенным условиям создания SIEM-системы и условиям реализации процесса ее функционирования.
Например, в качестве ПК, характеризующего своевременность выявления SIEM-системой отклонения в состоянии ИБ, может выступать среднее время выявления такого отклонения по конкретному критерию (появление вируса, нарушение прав доступа и т. д.). При этом под своевременностью выявления отклонения в состоянии ИБ понимается свойство SIEM-системы, характеризующее ее способность обеспечивать идентификацию отклонения и оповещение о нем в установленные сроки или в реальном масштабе времени.
В состав вектора ПК, характеризующего надежность SIEM-системы, например, могут быть включены параметры: среднее время наработки SIEM-системы на отказ , среднее время восстановления SIEM-системы
, коэффициент готовности SIEM-системы
и коэффициент технического использования
SIEM-системы. При этом под надежностью SIEM-системы подразумевается ее способность сохранять в установленных пределах времени значения всех ключевых параметров, определяющих выполнение ею требуемых функций.
Состояние ПК – понятие, обозначающее множество устойчивых значений этого показателя. Их может быть несколько либо большое множество. Состояние ПК характеризуется тем, что описывает переменные свойства этого показателя. При этом состояние какого-либо показателя качества процесса устойчиво до тех пор, пока процесс не реализуется, ведь сам процесс – это последовательная смена состояний показателей, его описывающих. Математическим языком описания процесса функционирования SIEM-системы как объекта исследования, на наш взгляд, может выступать язык теории множеств и функциональных пространств. В нашей статье рассмотрен подход к описанию математической модели процесса функционирования SIEM-системы как к формализованному описанию процесса смены состояний ПК (движения ПК) в пространстве состояний, причем данное движение является функцией времени. При этом предполагается, что аналитическая, вероятностно-временная модель смены состояний ПК процесса функционирования SIEM-системы может быть построена на основе марковских последовательностей.
В нашем случае в качестве примера могут быть рассмотрены три состояния ПК – три устойчивых значения j-го ПК функционирования SIEM-системы. Три состояния ПК – «отлично» (значения ПК удовлетворяют требованиям), «удовлетворительно» (значения ПК «на грани») и «плохо (авария)» (значения ПК не удовлетворяют требованиям).
Необходимыми базовыми исходными данными для аналитического описания процесса функционирования SIEM-системы в данной постановке являются вектор показателей качества (ВПК) SIEM-системы и требования к этим ПК. Тогда модель процесса функционирования SIEM-системы можно представить через аналитическое описание процесса смены состояний (значений) ее ВПК. При этом аналитическое описание ВПК процесса функционирования SIEM-системы, обозначим его х(k), аналитическая взаимосвязь отдельных ПК системы в общем случае динамического, вероятностного, нелинейного и нестационарного процесса ее функционирования связаны с необходимостью задания многомерных функций F(x, , t0, …, tk, …, tK) либо плотностей W(x,
, t0, …, tk, …, tK) распределения вероятностей значений параметров (показателей) процессов функционирования различных элементов SIEM-системы и процесса ее работы в целом на интервале функционирования этой системы (t0–tK). Математическое описание, с учетом расширения размерности (m – число переменных состояния процесса х(k); s – число параметров распределений; λ и k – количество отсчетов времени), в этом случае затруднительно.
Анализ релевантных работ
Анализ подобных ограничений приводит к необходимости поиска достаточно строгого аналитического описания процесса функционирования SIEM-системы в динамике. На наш взгляд, наиболее корректно данные ограничения могут быть учтены в рамках марковских моделей, рассматриваемых в ряде современных работ.
Работа [8] рассматривает математическое моделирование, основанное на общей теории систем. Здесь предложена модель в абстрактных терминах, что не позволяет полноценно использовать ее для моделирования конкретных угроз в рамках работы SIEM-системы.
В работах [9, 10] утверждается, что модель сложного управляемого процесса может быть реализована на основе стохастических дифференциальных уравнений. Но эти модели требуют значительных затрат на сбор статистики исходных данных для моделирования динамики изменений значений параметров и показателей качества SIEM-систем. Они не обеспечивают потребности пользователя по унификации и удобству (простоте) использования. Эти модели сложны с точки зрения описания динамики процессов. Так, в работах [11–13] предложены приложения моделей, построенных на стохастических дифференциальных уравнениях, но использованные в них традиционные методы определения параметров динамики изменений значений показателей качества объекта моделирования малоэффективны в силу необходимости рассмотрения множества вспомогательных параметров для полноценного описания всех свойств SIEM-систем, а это не всегда возможно.
Работы [14–19] посвящены возможностям математического описания процессов в динамике с использованием марковских цепей. Такие подходы будут рассмотрены и в нашей работе, это будут динамические модели процесса функционирования SIEM-систем. Они являются базой для моделирования, однако эти работы не рассматривают проблему расчета совместной плотности распределения вероятностей различных показателей качества процесса функционирования SIEM-систем. Частично свободна от этих недостатков модель, рассматриваемая в работе [20], однако, хотя она и опирается на цепи Маркова в форме разностных стохастических уравнений, но ориентирована на моделирование процессов обнаружения вредоносной информации. Наиболее близкой по сущности к идее, предлагаемой в данной статье, является работа [21], где предложена марковская модель процесса принятия решений по управлению информацией и событиями безопасности.
Таким образом, из анализа релевантных работ следует, что марковские модели процессов обладают высокой универсальностью, а сочетание теории марковских процессов с теорией переменных состояния открывает широкие возможности для исследования сложных систем, таких как SIEM-система. Иными словами, можно предположить теоретическую и практическую возможность построения аналитических, вероятностно-временных моделей смены состояний процесса функционирования SIEM-системы (как математических моделей смены состояний ПК этих систем) на основе марковских последовательностей, добиваясь требуемой степени адекватности вероятностно-временных свойств процесса функционирования систем такого класса при сокращении размерности их математического описания.
Теоретическая часть
Для построения аналитической вероятностно-временной модели смены состояний ВПК SIEM-системы, учитывающей динамический и вероятностный характер, нестационарность процесса ее функционирования и управления ею, воспользуемся аппаратом управляемых цепей Маркова (УЦМ), описываемых в форме разностных стохастических уравнений (РСУ). Анализ работ [8–21] позволяет говорить, что все известные классы марковских случайных процессов могут быть сведены к эквивалентным им (с точностью до допустимой ошибки моделирования по времени Δtдоп и по состоянию Δxдоп) цепям Маркова. Этот вид модели отличается от используемых ранее, поскольку пока не исследовались возможности применения УЦМ в форме РСУ (УЦМ-РСУ) для векторного анализа качества процесса функционирования SIEM-систем. С учетом результатов анализа можно утверждать, что от стохастических дифференциальных уравнений для разрывных процессов нетрудно перейти к РСУ, которые для дискретных моделей смены состояний ПК процесса функционирования SIEM-системы в форме УЦМ (Т = const) будут иметь вид
; (1)
; (2)
; (3)
, (4)
где выражение (1) – уравнение состояния процесса х на (k + 1)-м шаге (в нашем случае это состояние некоторого ПК процесса функционирования SIEM-системы), в котором – матрица-строка возможных значений отклонений ПК процесса функционирования SIEM-системы;
– вспомогательный вектор индикаторов состояния ПК процесса функционирования SIEM-системы, принимающий значения 0 или 1 и вводимый для удобства записи динамики переходов процесса функционирования системы такого класса из состояния в состояние
– число состояний).
Выражение (2) – уравнение состояния вспомогательного вектора индикаторов, в котором – матрица вероятностей перехода процесса, обуславливающего смену состояний ПК процесса функционирования SIEM-системы;
– вектор значений индикаторов состояния на предыдущем шаге;
– вектор компенсационных добавок, элементы которого предназначены для компенсации нецелочисленной части уравнения и получены в результате коррекции исходного шума возбуждения – белого гауссовского шума (БГШ) с математическим ожиданием и дисперсией, соответствующими начальному состоянию оцениваемого ПК процесса функционирования SIEM-системы
– уравнение (3).
Справедливость применения в роли шума возбуждения БГШ с адекватными реальному процессу параметрами обусловлена природой случайных значений ПК SIEM-системы; – диагональная блочная матрица компенсационных добавок, являющаяся ступенчатым мартингалом, элементы которого предназначены для компенсации нецелочисленной части в уравнении (2);
– вектор значений модифицированного возбуждающего шума, определяющий значения
на (k + 1)-м шаге (значения вспомогательных индикаторов процесса х (k + 1)).
Выражение (4) – уравнение наблюдения за процессом смены состояний ПК процесса функционирования SIEM-системы, где – вектор значений наблюдаемого ПК,
– матрица наблюдения, содержащая известные значения наблюдения за состоянием процесса х(k):
; (5)
– вектор шумов наблюдения за процессом смены состояний ПК SIEM-системы х(k) с нулевым средним и матрицей дисперсий
.
Аналитическая модель процесса смены состояний конкретного ПК, например среднего времени выявления отклонения в состоянии ИБ, характеризующего своевременность (оперативность) выявления SIEM-системой отклонения в состоянии ИБ, имеет вид
(6)
(7)
(8)
(9)
где выражение (6) – уравнение состояния показателя своевременности (оперативности) выявления отклонения в состоянии ИБ на (k + 1)-м шаге функционирования SIEM-системы, в котором – вектор-столбец значений показателя своевременности (оперативности) выявления отклонения в состоянии ИБ на (k + 1)-м шаге;
– транспонированная диагональная квадратная матрица (порядка m) возможных значений показателя своевременности (оперативности) выявления отклонения в состоянии ИБ на (k + 1)-м шаге функционирования
SIEM-системы, причем число m (строк и столбцов) зависит от выбранного числа состояний (глубины моделирования); – вспомогательный вектор-столбец индикаторов состояния показателя своевременности (оперативности), вводимый для удобства записи динамики перехода показателя своевременности (оперативности) из состояния в состояние. Выражения (7)–(9) имеют физический смысл, аналогичный выражениям (2), (3), и являются уравнением состояния вспомогательного вектора индикаторов, уравнением для определения приращений индикаторов состояния и уравнением наблюдения за процессом смены состояний среднего времени
выявления отклонения в состоянии ИБ соответственно.
Методологическая часть
Вектор значений модифицированного возбуждающего шума , определяющий значения
на (k + 1)-м шаге (значения вспомогательных индикаторов процесса х (k + 1)) формируются на основе сравнения с пороговыми значениями отклонений ПК процесса функционирования SIEM-системы выборочных значений БГШ
, преобразованного в возбуждающий шум
путем линейной процедуры
при
(10)
которая в матричной форме имеет вид
, (11)
где – приращения математического ожидания формируемого процесса, полученные на основе априорных данных о цепи Маркова, по отношению к матожиданию БГШ, возбуждающего этот процесс;
– коэффициент, предназначенный для коррекции дисперсии гауссовского шума на основе априорных данных о цепи Маркова для данного шага;
– матрица диффузии (возбуждения) процесса
;
– дисперсия процесса
;
– элементы матрицы интенсивностей перехода вспомогательных индикаторов состояния
процесса функционирования SIEM-системы; М – символ математического ожидания формируемого процесса на (k + 1)-м шаге; D – символ дисперсии формируемого процесса на (k + 1)-м шаге; Г (k + 1) – коэффициент, предназначенный для коррекции дисперсии на (k + 1)-м шаге; xпор m – граничное, «пороговое» значение конкретного m-го состояния ПК; Г2(k) – символ взятия в квадрат коэффициента, предназначенный для коррекции дисперсии гауссовского шума.
Рассмотрим детально шаги (этапы) методологии определения элементов вектора компенсационных добавок в уравнении состояния данной модифицированной модели системы управления информацией и событиями безопасности с учетом преобразования шума возбуждения.
Результаты анализа (на основе априорных данных о цепи Маркова), к какому состоянию «стремится» процесс на следующем шаге, дают нам решающее правило для выбора вектора компенсационных добавок и, в конечном итоге, для определения значений элементов вектора вспомогательных индикаторов состояния (k + 1) процесса функционирования SIEM-системы.
Вычисление значений , корректирующих вероятности состояния исходного шума возбуждения
, может быть осуществлено для случая однородной цепи априорно, на основе выражения
, (12)
где (–1) – символ взятия обратной функции; – интеграл вероятности;
– вероятность принятия процессом значения
, численно равная вероятности принятия процессом
значения
(k + 1) = 1:
, (13)
где М – количество состояний процесса; i – текущее i-е состояние процесса, где i = 1, … , m, … , M; pi(k –1) – вероятность того, что на предыдущем (k – 1)-м шаге процесс находился в i-ом состоянии; pim(k – 1/k) – вероятность перехода процесса из i-го состояния на предыдущем (k – 1)-м шаге в m-ое состояние на нынешнем k-ом шаге.
Вычисление значений коэффициента Гm(k), предназначенного для коррекции дисперсии случайных дискретных реализаций исходного БГШ , происходит следующим образом. Учитывая, что спектральная плотность мощности шума возбуждения
связана с его дисперсией выражением
, где Nϑ – коэффициент спектральной плотности мощности шума возбуждения; T – обозначение такта, отрезка времени, временного эквивалента одного шага
k процесса, за которое процесс может изменить свое состояние; а интенсивность перехода процесса Θ(k) из состояния в состояние – с соответствующими вероятностями – выражением , где pmm – элементы матрицы вероятностей перехода процесса Θ(k) из состояния в состояние, уравнение для коэффициента Гm(k) может быть записано в виде
. (14)
Если выбранное с учетом введенных коэффициентов значение модифицированного шума возбуждения попадает в границы допусков
(что показывает «стремление» Θm на (k + 1)-м шаге к 1), в уравнении (3) принимается решение о том, что элемент вектора компенсационных добавок
– положительный, а остальные элементы этого вектора – отрицательны, что при решении уравнения (2) дает нам Θm(k + 1) = 1. Формирование текущих значений матрицы компенсационных добавок ϑ(k) происходит в соответствии с правилом: если процесс описывается двумя состояниями Θ0 и Θ1 (т. е. m =
), то вектор последовательностей возбуждения для определения значения вспомогательного индикатора состояния Θ1 процесса функционирования SIEM-системы находится из матрицы
; (15)
если для описания процесса смены состояний ПК функционирования SIEM-системы используем три и более состояния Θm (т. е. m = 1, 2, …, M), то выбор матрицы компенсационных добавок для m-го состояния ϑm(k) (на первом этапе решения уравнения (3)) и вектора компенсационных добавок из этой матрицы (на втором этапе решения уравнения (3)) происходит из блочной диагональной матрицы вида
. (16)
Элементы матриц текущих компенсационных добавок ϑm(k) для различных состояний ПК процесса функционирования SIEM-системы находятся в соответствии с правилом
(17)
Количество состояний m определяется исходя из требуемой точности оценивания ПК процесса функционирования SIEM-системы. Оптимизация размерности пространств состояния и наблюдения процесса функционирования SIEM-системы с учетом требований к допустимой погрешности оценивания представляет собой отдельную теоретическую и практическую задачу. Окончательным результатом решения уравнения (3) является вектор , содержащий компенсирующие добавки для каждого рассматриваемого m-го состояния ПК SIEM-системы:
. (18)
Таким образом, показано, что предложенный методологический подход позволяет определять значения приращения математического ожидания моделируемого процесса, полученные на основе априорных данных о цепи Маркова, по отношению к математическому ожиданию БГШ, возбуждающего этот процесс. Этот подход сочетает вычисление значений приращения математического ожидания формируемого процесса и значений коэффициента, предназначенного для коррекции дисперсии случайных дискретных реализаций исходного БГШ, и позволяет принимать решение о том, какие значения будут принимать элементы вектора (18) в уравнении состояния данной модифицированной модели SIEM-системы с учетом преобразования шума возбуждения.
Заключение
Таким образом, выражениями (1)–(5) и (10)–(18) описан вариант модели процесса функционирования SIEM-системы. Эта модель представляет собой формализованное аналитическое описание, в терминах марковской цепи в форме разностных стохастических уравнений (1)–(4), динамики смены состояний показателей качества, характеризующих существенные свойства процесса функционирования SIEM-системы, в пространстве состояний.
В статье решена научная задача, заключавшаяся в усовершенствовании (модификации) алгоритмов преобразования шума возбуждения, используемого в модели. Описан новый механизм определения значений приращения математического ожидания моделируемого процесса, полученных на основе априорных данных о цепи Маркова, по отношению к математическому ожиданию белого гауссовского шума, возбуждающего этот процесс. Преобразование шума возбуждения служит для получения приращений математического ожидания моделируемого процесса, помогает определить эти приращения нецелочисленных значений индикаторов состояния. Полученные математические ожидания определяют «тренд», направление (ближе к 0 или 1) движения индикатора состояния ПК в пространстве состояний. В сущности, этот механизм позволяет, опираясь на несложные вычисления, принимать решение о том, какие значения будут принимать элементы вектора компенсационных добавок в уравнении состояния вспомогательного вектора индикаторов данной модифицированной модели с учетом преобразования шума возбуждения. Это позволяет упростить модель, снизить ее вычислительную сложность без существенных потерь в точности (адекватности).
Практическое применение усовершенствованной модели возможно как в рамках исследовательских работ, так и в системах автоматизированного контроля информационной безопасности. Направлением дальнейших исследований может быть разработка модели, учитывающей как количественный, так и качественный (например, описываемый с помощью лингвистических переменных, типичных для нечетких множеств) характер показателей качества.
1. Kizza J. M. Guide to Computer Network Security. New York: Springer, 2015. 545 p.
2. Sun Y. Research on Security Issues and Protection Strategy of Computer Network // The Open Automation and Control Systems Journal. 2015. N. 7. P. 2097-2101.
3. Dordal P. L. An Introduction to Computer Networks. Release 1.9.0. January 27, 2017, 745 p. URL: https://archive.org/details/academictorrents_958e2487d2db5f41f9c056bb35cf547edf38528f/mode/2up (дата обращения: 22.12.2019).
4. Pfleeger C. P., Pfleeger S. L. Security in Computing. New Jersey: Prentice Hall, 2015. 944 p.
5. Miller D. R., Harris S., Vandyke S. Security Information and Event Management (SIEM) Implementation. New York: McGrawHill, 2011. 430 p.
6. Kotenko I. V., Parashchuk I. B. Synthesis of controlled parameters of cyber-physical-social systems for monitoring of security incidents in conditions of uncertainty // Journal of Physics: Conference Series, IOP Publishing. 2018. N. 1069 (1): 012153. P. 1-6.
7. Kotenko I. V., Doynikova E. V. Countermeasure selection in SIEM systems based on the integrated complex of security metrics // Proceedings of the 23th Euromicro International Conference on Parallel, Distributed and Network-based Processing (PDP 2015). IEEE Computer Society, 2015. P. 567-574.
8. La Padula L. J. Secure Computer Systems: A Mathematical Model // MTR-2547. Massachusetts: The MITRE Corporation, Bedford, 1993. V. I. 33 p.
9. Higham D. J. An Algorithmic Introduction to Numerical Simulation of Stochastic Differential Equations // SIAM REVIEW. 2001. N. 43 (3). P. 525-546.
10. Iacus S. M. Simulation and Inference for Stochastic Differential Equations: With R Examples. New York: Springer Verlag, 2008. 214 p.
11. Yuksel S. Control of Stochastic Systems // Queen’s University Mathematics and Engineering and Mathematics and Statistics. 2017. 167 p.
12. Van Handel R. Stochastic Calculus, Filtering and Stochastic Control. New York, Springer, 2007. 261 p.
13. Oksendal B. Stochastic Differential Equations. An introduction with applications. Berlin Heidelberg, Springer Verlag, 2007. 311 p.
14. Stewart N. E., Thomas G. K. Markov processes: Characterization and Convergence. Wiley Series in Probability and Statistics. New York, John Wiley & Sons Inc., 1986. P. 214-234.
15. Bini D., Latouche G., Meini B. Numerical Methods for Structured Markov Chains. New York: Oxford University Press, 2005. 215 p.
16. Dobre T. Gh., Marcano J. G. S. Chemical Engineering: Modelling, Simulation and Similitude. Weinheim: Wiley-VCH, 2007. 568 p.
17. Oliver D., Kelliher T., Keegan J. Engineering Complex Systems With Models and Objects. New York: McGraw-Hill, 2007. 340 p.
18. Quarteroni A. Mathematical Models in Science and Engineering // Proc. of Notices of the AMS. 2009. V. 56. N. 1. P. 9-19.
19. Tweedie R. L. Drift conditions and invariant measures for Markov chains // Stochastic Processes and Their Applications. 2001. V. 1. P. 345-354.
20. Kotenko I. V., Parashchuk I. B. Determining the Parameters of the Mathematical Model of the Process of Searching for Harmful Information // Cyber-Physical Systems: Industry 4.0 Challenges. Studies in Systems, Decision and Control 260. A. G. Kravets et al. (eds.). Springer Nature Switzerland AG 2020, 2019. P. 225-236.
21. Kotenko I. V., Parashchuk I. B. An approach to modeling the decision support process of the security event and incident management based on Markov chains // 9th IFAC Conference on Manufacturing Modelling, Management and Control (MIM-2019) (Berlin, Germany, 28-30 August 2019). IFAC-PapersOnLine, 2019. V. 52. Iss. 13. P. 934-939.
