CC BY
20
Доклады БГУИР
2011 №4 (58)
УДК 681.3
МЕТОДИКА ОЦЕНКИ ПРОПУСКНОЙ СПОСОБНОСТИ ПОТОКОВОЙ СЕТИ НА ОСНОВЕ ВЕРОЯТНОСТНО-АЛГЕБРАИЧЕСКОГО МОДЕЛИРОВАНИЯ
Е.И. СУКАЧ
Гомельский государственный университет имени Ф. Скорины Советская, 104, Гомель, 246019, Беларусь
Поступила в редакцию 24 февраля 2011
Рассматривается подход к оценке динамических характеристик пропускной способности потоковой сети, основанный на методе вероятностно-алгебраического моделирования, позволяющий определить вероятностные значения характеристик функционирования исследуемых систем по аналогичным характеристикам составляющих их компонентов.
Ключевые слова: вероятностно-алгебраическое моделирование, потоковые сети, пропускная способность, цепи Маркова.
Введение
Объектом исследования являются потоковые системы. К ним относятся производственные, вычислительные и транспортные системы, которые можно представить в виде графовых структур и функционирование которых подчиняется следующим условиям: потоки на входе переходника суммируются; поток на выходе переходника равен входному потоку, деленному на число выходов; суммарный поток не может превышать на входе/выходе 100%.
Эксплуатация потоковых систем предполагает, как правило, различные варианты организации выполнения предписанных им функций. Выбор лучшего варианта функционирования системы в смысле исследуемого свойства должен основываться на математически обоснованном сравнении вариантов, отличающихся как структурной организацией, так и параметрами элементарных участков, составляющих систему.
С позиций анализа вероятностных значений пропускной способности потоковых систем применимы как статические, так и динамические модели. Для статического исследования характеристик потоковых систем используются аналитические методы и средства их автоматизации, такие как АРБИТР (ПК АСМ СЗМА) [1] - программный комплекс автоматизированного структурно-логического моделирования и Relex Reliability Studio [2] - программная среда, включающая различные методы анализа надежности и реализующая разнообразные формы задания моделей (графы, деревья отказов). Применение указанных методов ограничено числом участков, составляющих систему, учитывающих, как правило, не более двух возможных значений пропускной способности. Увеличение количества участков и рассмотрение множества их вероятностных значений пропускной способности приводит к росту пространства состояний модели и усложнению связей между состояниями.
К динамическим моделям можно отнести статистическое имитационное моделирование, предполагающее рассмотрение различных траекторий функционирования исследуемых систем во времени с последующим усреднением полученных результатов моделирования. В частности, с целью определения максимальной пропускной способности сети, функционирующей в условиях случайных воздействий, разработан метод, основанный на сочетании алгоритма Форда-Фалкерсона и метода Монте-Карло [3]. При этом изменение параметров функционирования исследуемой системы требует проведения очередной серии имитационных экспериментов и последующую обработку результатов.
В статье рассматривается методика определения пропускной способности потоковой системы на основе вероятностно-алгебраического моделирования (ВАЛМ) [4], позволяющая учесть множество возможных изменяющихся во времени значений пропускной способности любого числа участков, составляющих систему и е е структурную организацию.
Общая характеристика метода вероятностно-алгебраического моделирования
Вероятностно-алгебраическое моделирование реализует процесс формирования вектора вероятностей состояний системы по векторам вероятностей состояний составляющих систему компонентов с учетом установленных между ними связей [4].
Система формализуется в виде множества компонентов К = {К1}, ' = 1, т, согласованное взаимодействие которых обеспечивает ее функционирование. Компоненты могут находиться в одном из множества состояний £ = ^ ] = 1, п , характеризующих исследуемое свойство системы (надежность, пропускную способность и др.). Вероятности нахождения компонентов системы в каждом из состояний задаются векторами вероятностей:
р'=(р, р2,..., рп), £р• =1, i=1т. (1)
]=1
Предполагается, что компоненты системы независимы и между ними могут быть установлены функциональные связи с учетом целей исследования, которые описываются набором детерминированных и/или вероятностных функций F={Fz}.
Элементы вектора вероятностей результирующего вектора Р3, полученного в результате вероятностно-алгебраического умножения векторов Р1 и Р2, определяются по формуле:
п п
р1 = £ £ 4р;р;2 , где ', к = 1п . (2)
1=1 г=1
Коэффициенты ак называются коэффициентами вероятностно-алгебраического моделирования и удовлетворяют следующим требованиям:
V', 1,к аг]к > 0 и £^аг]к = 1. (3)
к=1
В случае детерминированных связей между компонентами исследуемой системы коэффициенты вероятностно-алгебраического моделирования определяются соотношениями:
\а'к =1 еслик = F (', j), (4)
I ак = 0, если к Ф F(г', 1).
Основу метода ВАЛМ составляет алгебраический аппарат, который является удобным для исследования вероятностных свойств сложных систем и позволяет распространить общие свойства алгебраических структур на исследуемую предметную область.
Методика оценки пропускной способности потоковой сети
Определение вероятностных характеристик пропускной способности потоковой системы реализуется с использованием системы вероятностно-алгебраического моделирования PALS (Probability-Algebraic Simulation) [5] следующей последовательностью шагов.
На шаге 1 выделяется множество элементарных участков потоковой системы и связей между ними. В диалоговом режиме с использованием стандартных компонентов K={Ki} и библиотеки функций F={Fz} формируется графическая схема G(F,K) исследуемой системы, которая представляет собой дерево. Листьями дерева являются компоненты исследуемой системы, а узлами - функции, определяющие связи между компонентами системы.
На шаге 2 на основе натурных экспериментов с прототипом потоковой системы или методом экспертных оценок реализуется подготовка параметров моделирования.
На шаге 3 для каждого компонента задается количество состояний S = ^. |, . = 1, п и
начальные вектора вероятностей, характеризующие эти состояния:
P'0 = (Р0, Р0,..., p"), £ pj = 1, i = 1,
m
(5)
j=1
Для этой цели используются соответствующие средства ввода информации PALS, которые автоматически заносят в базу данных исходные данные.
На шаге 4 осуществляется задание времени моделирования и выбор способа формирования значений векторов вероятностей
pu = (pi, р?,..., РП ), £ pj = i, t=i, т,
j=i
(6)
характеризующих изменение пропускной способности компонентов во времени.
Одним из способов, определяющих динамическое изменение векторов вероятностей (6), является первичное моделирование с использованием марковской модели с дискретным временем и дискретными состояниями, переходная матрица которой имеет вид:
PK =
Pi Р2
0 £ Р,
i=1
0 0
Рз Р4
Рз Р4
з
£ Рi Р4
•• Рп-1 Рп
■. Рп~1 Рп
\
" Р,
п-1 Fn
£ р i Рп
0
1
(7)
Возможен переход к рассмотрению процессов вероятностных изменений участков потоковой системы с восстановлением. В зависимости от стратегий восстановления компонентов выбирается одна из имеющихся параметризованных моделей Маркова. Параметры моделей, описывающих изменение вероятностных характеристик компонентов потоковой системы, определяются в результате алгоритмической обработки данных, полученных при проведении натурных экспериментов.
Альтернативным способом задания динамического изменения векторов вероятностей является описание элементов исходных векторов вероятностей (5) параметрическими функциями. Использование в качестве параметра функций времени позволяет управлять скоростью изменения вероятностных характеристик компонентов. Проведение статического моделирования, реализуемого одной итерацией, при параметрическом задании исходных векторов (5) позволяет сформировать результирующий вектор, характеризующий пропускную способность потоковой системы в символьном виде.
На шаге 5 в соответствии с заданным временем моделирования проводится динамическое вероятностно-алгебраическое моделирование, которое реализует компьютерные вычисления на каждой итерации моделирования с учетом вероятностного изменения состояний компонентов. При этом неоднократно реализуется статическое моделирование, осуществляющее свертку исходных векторов вероятностей состояний компонентов по формуле (2) с учетом уровня вложенности функций и коэффициентов вероятностно-алгебраического моделирования
(4).
На шаге 6 по желанию пользователя результаты моделирования графически отображаются в виде временных диаграмм, представляющих изменение вероятностных характеристик пропускной способности, как отдельных участков, так и всей потоковой системы.
На шаге 7 для случая исследования функционирующей потоковой системы осуществляется проверка адекватности построенной вероятностно-алгебраической модели реальному объекту. В PALS она проводится автоматически путем проверки близости средних значений откликов модели соответствующим характеристикам реальной потоковой сети. В случае отрицательных результатов осуществляется переход на шаг 4.
На шаге 8 определяется влияние вероятностных характеристик пропускной способности участков потоковой сети на значение компонентов вектора откликов всей системы при е е фиксированной структурной организации. С этой целью организуются модельные эксперименты, в которых, во-первых, варьируются значения векторов (5), позволяющие оценить влияние начального состояния участков системы на результирующее состояние всей сети. Во-вторых, для рассмотрения возможных режимов эксплуатации потоковой системы изменяются параметры марковских моделей участков, описывающие вероятностное изменение пропускной способности во времени.
На шаге 9 исследуется влияние структурной организации потоковой системы на векторы откликов моделирования при неизменных вероятностных значениях параметров компонентов. При этом могут быть рассмотрены случаи: введения новых участков, увеличивающих пропускную способность участков - «узких мест» сети, обнаруженных при исследовании исходного структурного варианта сети; исключения участков, пропускная способность которых обращается в ноль в результате аварий (профилактических ремонтов); рассмотрения альтернативных вариантов структурной организации потоковой сети. Сравнение результирующих векторов пропускной способности сети для различных вариантов ее структурной организации позволит выбрать лучший из них, а также оценить выигрыш при резервировании отдельных участков и проигрыш при рассмотрении усеченных вариантов.
Пример определения вероятностных характеристик пропускной способности
потоковой системы
Для демонстрации методики вероятностной оценки пропускной способности была выбрана транспортная система сообщения (ТСС), включающая участки дорог, соединяющие населенные пункты. Варианты структурной организации ТСС представлены на рис. 1,а,б.
Участки дорог исследуемой ТСС представлены компонентами K = {Kj|, i = 1,28, которые описываются однотипным образом и характеризуются множеством состояний S = {S j I, j = 1,10, соответствующих определенному уровню пропускной способности, вероятностно изменяющемуся во времени. Пропускная способность участка (PR) указывает на максимально возможное количество единиц транспорта, которое он способен пропустить за выбранную единицу времени. Состоянию Si соответствует максимальная пропускная способность PRi. Состоянию S10 соответствует пропускная способность PR10=0, при которой движение по участку дороги прекращается. Остальные состояния характеризуют промежуточные уровни пропускной способности участков, значения которых устанавливаются в зависимости от множества случайных факторов, одним из которых является износ дорог, возникающий в процессе обслуживания сетью транспортных потоков.
Начальные значения вероятностей состояний, характеризующих пропускную способность участков, задаются векторами:
10 _
Pi 0 = (p0, p20,..., К), £ pj = 1, i = 1,28,
j=1
которые изменяются во времени.
Предполагается, что вероятностное изменение состояний участков дорог описывается марковскими моделями с дискретными состояниями и дискретным временем, параметры которых задаются матрицами вида (7). Матрицы переходных вероятностей PK.,[10х 10],i = 1,28 для всех компонентов различны и формируются в результате статистического анализа данных, характеризующих износ участков.
Рис. 1. Графовая структура ТСС: а - вариант 1; б - вариант 2
В результате первичного моделирования формируются значения векторов вероятностей, описывающие состояния участков сети на заданном промежутке времени:
Рй = (рй, р2,..., К), £ р. = 1, i = 1,28, X = 1,100 .
з =1
а
Полученные вектора являются исходными данными для вероятностно-алгебраического моделирования исследуемой системы. Ставится задача определения вектора вероятностей, характеризующего состояние пропускной способности сети, по вероятностным характеристикам участков сети.
В расчетных вероятностных моделях используется функция F1 (/, у) = mm(i + у -1, п) для параллельных участков сети и функция F2(i, у) = max(i, у) для последовательных участков, что имеет естественную интерпретацию в случае исследования свойства пропускной способности ТСС. Результирующая пропускная способность параллельно расположенных участков является суммой их пропускных способностей. Пропускная способность последовательно расположенных участков определяется минимальным значением пропускных способностей участков.
Результатом моделирования являются динамические вероятностные характеристики пропускной способности, как отдельных участков, так и ТСС в целом для двух вариантов ее структурной организации, представленные на рис. 2, 3.
На рис. 2 представлена динамика состояний участка К10, отображающая переход от вероятностного распределения состояний, соответствующих максимальной пропускной способности, к более равномерному распределению вероятностей возможных состояний. На рис. 3 сопоставляются вероятностные характеристики пропускной способности двух вариантов сети (25 участков и 28 участков).
Дааампка тиенения вероятногтных характеристик Вероятности СОСТОЯНИЙ пропускной способности ТСС
Рис. 2. Динамика вероятностных состояний Рис. 3. Вероятностные характеристики
участка К10 пропускной способности двух вариантов сети
Сравнительный анализ данных позволяет сделать вывод о том, что изменение структурной организации ТСС заметно отражается на изменении вероятностных характеристик пропускной способности. Введение дополнительных участков (вариант 2) позволяет уменьшить снижение пропускной способности сети во времени, возникающее в результате естественного износа эксплуатируемых дорог ТСС.
Заключение
Применение методики оценки пропускной способности потоковой системы с использованием ВАЛМ позволяет учесть ее вероятностную природу и решить следующие задачи:
- статически и в динамике провести сравнительный анализ различных структурных вариантов системы;
- подобрать параметры компонентов, обеспечивающих необходимый уровень пропускной способности системы;
- оценить влияние параметров пропускной способности участков и их групп на пропускную способности всей сети;
- получить в символьном виде вектор вероятностей пропускной способности сети.
THE METHOD OF EVALUATION OF THE CAPACITY OF STREAMING NETWORKS BASED ON THE PROBABILISTIC-ALGEBRAIC MODELING
E.I. SUKACH
Abstract
An approach to estimating the dynamic characteristics of the capacity of streaming network, based on the method of probability-algebraic modeling is considered, which allows to determine the probability values of the characteristics of the systems' functioning using the same characteristics of their components, is consideted.
Литература
1. Нозик А.А., Можаев А.С // Монтаж и наладка средств автоматизации и связи. 2007. №2. С. 32-40.
2. [Электронный ресурс] Режим доступа: http://www.relex.com/products.
3. Сукач Е.И. // Реестращя, збертання i обробка даних (Data Recording, Storage & Processing). 2008. Т.10, №3. С. 37-45.
4. Сукач Е.И. // Математичш машини i системи (Mathematical Machines and Systems). 2010. №3. С. 116123.
5. Сукач Е.И., Ратобыльская Д.В. // VI Международная конференция-форум «Информационные системы и технологии», Академия управления при Президенте Республики Беларусь, 24-25 ноября 2010.
