х МОДЕЛИРОВАНИЕ СИСТЕМ И ПРОЦЕССОВ
УДК 004.75
ОЦЕНКА ВЕРОЯТНОСТИ ЭРЛАНГОВСКОГО СТАРЕНИЯ ИНФОРМАЦИИ
Б. Я. Советов,
доктор техн. наук, профессор
Санкт-Петербургский государственный электротехнический университет «ЛЭТИ»
М. О. Колбанёв,
доктор техн. наук, профессор
Санкт-Петербургский государственный экономический университет Т. М. Татарникова,
доктор техн. наук, доцент
Санкт-Петербургский государственный университет аэрокосмического приборостроения
Рассматривается процесс предоставления информационной услуги, который регламентируется соглашениями о качестве обслуживания. Предлагается процедура вычисления вероятности своевременной обработки запроса пользователя как одной из основных характеристик качества обслуживания. Приводится методика оценки вероятности своевременного предоставления информационных услуг центрами обработки данных. Показано, что искомая вероятность может быть вычислена, если преобразование Лапласа — Стильтьеса функции распределения времени обработки запроса пользователя представляет собой произведение или частное двух функций, является суперпозицией двух функций или задано неявно.
Ключевые слова — качество обслуживания клиента, информационная услуга, показатели качества предоставления услуг, центр обработки данных, вероятность своевременного предоставления информационной услуги.
Введение
В последние годы государственными органами и бизнесом все шире используются центры обработки данных (ЦОД), обеспечивающие надежное хранение данных, круглосуточную поддержку бизнес-процессов, качественное предоставление других информационных услуг [1, 2].
Процесс предоставления услуг регламентируется соглашениями о качестве обслуживания (Service Level Agreement — SLA). SLA — это часть договора между провайдером и клиентом, в котором приводится система показателей для контроля качества предоставляемых услуг, а также штрафы за нарушение соглашения и вознаграждения в случае успешного его выполнения.
Вероятность эрланговского старения информации как характеристика качества обслуживания
SLA не регламентирует содержательную сторону процесса информационного взаимо-
действия, а относится только к процессу предоставления информационной услуги. Качество этого процесса оценивается вероятностью своевременной обработки запроса пользователя, под которой понимается вероятность того, что ответ на запрос будет получен за допустимое время [3]. При моделировании сложных информационных систем и сетей, подобных ЦОД, эту характеристику обычно называют вероятностью старения информации или вероятностью своевременного предоставления информационной услуги (Р ). В общем случае Русл определяется выражением
Русл = |фдоп(*) dT(x),
(1)
0
где фдоп(х) — функция распределения (ФР) допустимого времени реализации моделируемого процесса предоставления информационной услуги, устанавливаемая SLA; T(x) — ФР времени реализации этого процесса.
При нормировании вероятности Русл, как правило, выбирается постоянное ограничение на время предоставления услуги фдоп(х). Приме-
ром такой нормы может служить рекомендация Q.543ITU-T, в которой один из показателей качества обслуживания вводится как 95%-квантили распределений времени своевременной обработки запросов пользователей. В ряде работ получены результаты для так называемого экспоненциального старения, когда функция Фдоп(х) имеет показательное распределение.
Наиболее общим случаем является использование в качестве ограничения обобщенного распределения Эрланга, т. е. функции вида
кг -1
ф(х) = 1 УгХ £
( Угх )
(2)
г =1
• =0
^ 1 П
= 1 - =£Чг/Уг У г=1
г=1
где 1/у — среднее допустимое время предоставления услуги. Эту характеристику можно назвать вероятностью эрланговского старения информации ЦОД. Путем соответствующего выбора значений параметров п, цг, кг и уг функцией (2) можно достаточно точно аппроксимировать любую реальную ФР допустимого времени предоставления услуги. Так, например, при п = 1 и кг ^ да функция (2) стремится к постоянному ограничению, а при п = 1 и кг = 1 превращается в экспоненциальное распределение.
Подставляя (2) в (1), после ряда преобразований получаем
к„ -1
уел
= 1 £ (-1)г^^0),
(3)
где в0 = йгуг; ¿(1)(в0) — значение 1-й производной преобразования Лапласа — Стильтьеса (ПЛС) функции распределения времени предоставления услуги в точке во-
Процедура вычисления вероятности Русл
Рассмотрим лежащую в основе расчета вероятностей Русл по формуле (3) процедуру вычисления производных от ПЛС ФР ¿(в) в точке во, предполагая, что функция ¿(в) может представлять собой произведение или частное двух функций, может быть сложной функцией (суперпозицией двух функций) или может быть задана неявно.
Пусть х(в) и у(в) п раз дифференцируемы
в точке в = в0. Тогда для п-й производной
0 «0 функции Z(в) = х(в)у(в) в точке в0 справедлива
формула Лейбница
„(и) ■чт^ • (п-I) (•)
Я ' = £ с„Х У ' 1=0
(4)
' тя тг '
где хУз и у« — значения производных функций х(в) и у(в) 0в точке в0 и предполагается, что
х(0) = х («0 ) у<°0) = У («0) и 4°0)= Я («0);
• п!
с1п = — биномиальный коэффициент.
• !(п - •)!
Если функция £(в) = х(в)/у(в), то п-я производная Я(п) может быть вычислена рекуррент-
но. Умножим левую и правую части последнего уравнения на у(в) и возьмем п-ю производную от левой и правой частей полученного равенства. Тогда, применяя (4), получим
( п) ■ЧТ^ • (п-• )„(•) •=0
(5)
откуда следует рекуррентное соотношение
(п) ^Т1 • (п-1)г,(1) г(п) = 1=0
У («0)
(6)
(п)
которое позволяет вычислить Я через производные £(в) в точке в0 более низких порядков.
Если функция £(в) = у[х(в)] и известны значения у()( = 1,п производных по х функции у(х)
х0 V /
в точке Х0 = х(в0), то ) может быть вычислена следующим образом: 0
Я1) = У(1)х(1)
Ух •
(7)
Переходя к производным порядка п и применяя (4), получаем
Я(п) =
«0
п-1
= £ Сп-! •=0
1(1)
(п-1-1)
X
(1+1)
(8)
где
— ¿-я производная по в функции
dy(х)/dх в точке в0.
Формула (8) позволяет вычислить Я , если
(п)
известны производные х
(1)
(п)
и У
(1)
(п-1)
(1)
■V ' х0
причем последние так-
же зависят от производных функции х(в) порядка не выше (п - 1) и могут быть вычислены рекур-рентно:
п
26 У ИНФОРМАЦИОННО-УПРАВЛЯЮЩИЕ СИСТЕМЫ
"7 №6, 2013
>)
(1)
(т+1) (1) -7
= ух 'х , т = 1, п -1;
(к) к-1 () = ц ск-1
во у=0
У(т+1)
Ух
(-1-У) (v+1) 1-к
х , т =1, п -к;
где
уХ1
уХт)
(п-1)
п-2
= 1 сП-2
у=0
(п)'
УХ ;
(п-2-у)
х
.(v+1)
(9)
(к)
— к-я производная по в функции
^у(х)^хт в точке в0
Пусть теперь х(в) = у[х(в)]+в. Тогда после дифференцирования левой и правой частей и при-
ведения подобных для первой производной х
(1)
получаем
х (1)= 1/11 -
(1 - у^),
(10)
а для производных более высокого порядка непосредственно из (8) после замены £(в) на х(в) получаем следующую рекуррентную формулу:
п-2
X 4-1
х(п) = _г=о_
Л
(п-1-1)
.(£+1)
1-у
,(1)
(11)
хо
где производные
(к)
( = 1, П - 1 могут
быть рассчитаны по формуле (9).
Таким образом, для того чтобы вычислить п-ю производную (п> 1) функции х(в) = у[х(в)]+в, необходимо вычислить:
х^ по формуле (10);
У(1)'
,(2)
(1)
,(п-1)'
(1)
по формуле (9) при
х^ по формуле (11);
у(1)
(2)
(п-2) Ух '
(2)
по формуле (9) при
к = 1 и т = 1, 2, ..., п - 2
(3)
х^ по формуле (11) и так далее, пока не будут получены значения
(п-1)
и х
(п)
Приложение вычисления вероятности эрланговского старения информации для ЦОД
Проиллюстрируем изложенную процедуру вычисления вероятности эрланговского старения на примере ЦОД, который содержит серверы, обменивающиеся блоками данных взаимодействия (БДВ). Предположим, что обмен осуществляется через среду множественного доступа, а при обработке данных серверами используется дисциплина обслуживания с относительными приоритетами. Структура современного ЦОД подробно рассматривается в работе [4].
В этом случае для ПЛС ФР времени th(в) процесса предоставления услуги в ответ на БДВ с относительным приоритетом h имеем
Ыв) = <бр № («Коот (в),
(12)
где Ю£6Р (в) и рЫбр (в) — соответственно ПЛС ФР времени ожидания начала обработки и времени обработки в сервере БДВ с приоритетом Ь,
дОСТ/-\ _
(в) — ПЛС ФР времени доставки БДВ к серверу через среду множественного доступа. Для ю°бр (в) справедливы выражения
(1 -р)Ок (в) + X Ху [1 -Ру (н (в))]
^0бр (в) = -
\-.hi >1г
в - X Ху [1 -Ру (Ян (в))] Он (в) = в + X Ху [1 -Ру (Ян (в))],
У.1; <1
(13)
(14)
где р — загрузка сервера; Ху — интенсивность поступления БДВ, имеющих приоритет у.
Преобразование юдост(в) определяется протоколами доступа управляющих элементов к среде передачи и управления логическим каналом и параметрами физического уровня и среды передачи [5]:
(в) =
в (1 -ру) а{евТ) в -Ху +\ё(еаТ)'
(15)
где ру и Ху — соответственно загрузка и интенсивность потока БДВ на среду множественного доступа; Т — временной параметр среды множественного доступа; g(eвT) — ПЛС ФР интервала обслуживания БДВ для выбранной совокупности указанных протоколов.
о
о
к = 1 и т = 1, 2, ..., п - 1
Методика вычисления вероятности своевременного предоставления информационной услуги
В соответствии с разработанным методом вычисление вероятности Русл осуществляется в следующем порядке.
1. Найти значения производных Щ' ^) от ПЛС Р;(в) в точке в0 (• = 1, кг -. Если, например, время обработки БДУ с приоритетом ] распределено по показательному закону со средним значением tj, то
Ру (s ) = 1/ (l + tjs) и искомые производные
(16)
j )(s ) = (-1)Ч !—^
(j)
(1 + jo)'
i+1
(17)
2. Найти значения производных Р(-1)(С0) функций Р;(вл(в)) по Сл(в) в точке в0 = Сл(в0). Решая (14) методом последовательных приближений при в = в0, находим ©0, а затем — искомые
значения Р( )(С0) по формуле (17) путем замены
в0 на °0. обр(;)
3. Найти значения производных юй (в0)
от преобразования ю)(в) в точке в0. Для этого обозначим числитель и знаменатель формулы
(13) (s) и ^ (s)
обр
соответственно.
Очевидно, что производные
„обР
обр
р (s0) равны
бР (so)
<f )(so ) = (!-p)Gi )(so )-
- X a-у [Ру (((s))^;
j:hj >h
1 + X aj [Ру (h (s))]s0) при i =1;
j:hj =h
[Ру (Gh (s))]s0) при i > 1.
o6p(i) 3h,2 =
(18)
(19)
j:hj =h
В свою очередь производные )(в0) функции (14) могут быть вычислены из (9)—(11), если положить х(в) = Сл(в), а производные по в
[^ (ь («)) ) сложной функции Ру(вл(в)) — по
(8) и (9), если положить £(в) = Р;(вл(в)), х(в) = Он\
Значения производных ю^6^^^) могут быть найдены по формуле (6) с заменой функций Z(в), х(в), у(в) и их производных на функции ю°]г"р ), юЬб1Р («), юЬбр (в) и их производные соответственно.
4. Аналогично п. 1—3 найти (кг - 1)-ю производную от (15).
5. Найти производные т^' ^) от ПЛС ФР (12), применяя последовательно (4) к следующим
парам функций: Рл(в) и ю£бр (я), [ю£брИРЛ(а)] и
дост
(s), вычисляя таким образом производные
произведения трех функций.
После вычисления производных )(й0), (• = 1, кг -1 искомая вероятность эрланговского старения Русл может быть найдена по формуле (3).
Заключение
Рассмотрен подход к вычислению вероятности старения информации, основанный на рекуррентных процедурах расчета производных большого порядка от ПЛС функции распределения времени формирования информационных услуг. Предложенный метод может использоваться для вычисления вероятности старения информации для центров обработки данных и других информационных систем при формировании и контроле выполнения SLA.
Литература
1. Советов Б. Я., Цехановский В. В. Информационные технологии. - М.: Юрайт, 2012. - 263 с.
2. Советов Б. Я., Цехановский В. В., Чертовской В. Д.
Интеллектуальные системы и технологии. - М.: Академия, 2013. - 320 с.
3. Колбанев М. О., Яковлев С. А. Модели и методы оценки характеристик обработки информации в интеллектуальных сетях связи. - СПб.: СПбГУ, 2002. - 230 с.
4. Колбанев М. О., Татарникова Т. М., Воробьев А. И.
Модель балансировки нагрузки в вычислительном кластере центра обработки данных // Информационно-управляющие системы. 2012. № 3(58). С. 37-41.
5. Колбанев М. О., Татарникова Т. М., Воробьев А. И. Модель обработки клиентских запросов // Телекоммуникации. 2013. № 9. С. 42-48.
и
28 Ж ИНФОРМАЦИОННО-УПРАВЛЯЮЩИЕ СИСТЕМЫ Ж №6, 2013