Рассматривается численная модель открытой системы массового обслуживания, степень влияния соотношения входных параметров системы на погрешность вычисляемых динамических значений характеристик.
Ключевые слова:система массового обслуживания, характеристики системы, моделирование, проверка динамических значений
Первоначальная оценка погрешностей на этапе отладки алгоритма [3] производилась вручную путем создания нескольких прогонок (прогонкой здесь и далее называется запуск модели СМО от начала до достижения критерия остановки) модели одного типа и изменения входных параметров системы. Но больший интерес представляет дальнейшее исследование влияния входных параметров системы на динамические значения характеристик. Оценка погрешности и установление оптимальных входных параметров производились по нескольким критериям:
1. Определение входных параметров, при которых максимальная погрешность определённой характеристики при некотором числе испытаний минимальна. Данную оценку можно применить к группам характеристик, объединенных общим признаком, например, число требований в системе, под обслуживанием и т. п.
2. В некоторых случаях результат каждой из этих оценок необходимо проверить, вручную запустив прогонку модели с полученными входными параметрами и оценить поведение характеристик. Например, относительная погрешность величины, равная , может говорить о том, что статическая величина, вычисляемая по формуле стремится к нулю, тогда как динамическая величина уже стала равной
, а относительная погрешность равная
, означает, что статическая величина чрезвычайно мала, что не позволяет хранить её значение в переменной, в результате чего, она также становится равной
. Однако такие случаи редки и довольно легко могут быть определены.
3. Повторение испытания с увеличением максимального числа требований в системе, а также числа прогонок.
4. Определение входных параметров, при которых максимальная погрешность испытаний по всем характеристикам созданной системы минимальна.
Погрешность испытаний также в значительной степени зависит от системы ГСЧ, используемой в созданной модели, поэтому этот параметр неизменен на протяжении всего испытания и выделен в отдельную характеристику при подсчёте погрешностей. По этой же причине, а также чтобы минимизировать влияние ГСЧ на результаты испытаний, оцениваться будут максимальные погрешности, а не средние или минимальные. Таким образом, оценка в целом будет заключаться в поиске такого сочетания входных параметров, при котором максимальная погрешность по всем испытаниям определенной характеристики будет минимальна. Поэтапная оценка алгоритма, моделирующего СМО открытого типа с очередью ограниченной длины [1], представлена ниже.
1. Первым испытанием станет «поверхностная» проверка многоканальной СМО с неограниченной очередью с использованием встроенного псевдослучайного ГСЧ среды Visual Studio 2010 [2] и с ограничением максимального числа требований в системе равным . Цель — выявление наилучшего порядка величин
и
. Также необходимо ограничить величину
интервалом
, т. к. сходимость динамической системы в граничных условиях не может быть гарантирована. При указанных ограничениях и числе испытаний каждого соотношения равном
, необходимо будет совершить
прогонок. Наилучшим сочетанием входных параметров в целом, на основе определения минимальной погрешности наборов и максимальной погрешности характеристик в каждом наборе по всем прогонкам, для данного испытания является следующий набор, показанный в табл. 1.
Таблица 1
Наилучшее сочетание входных параметров в целом
Тип характеристик |
|
Погрешность (%) |
Моменты 1-го порядка |
|
|
Момента порядка выше 1 |
|
|
Из-за различной скорости сходимости моментов разного порядка, данная оценка была разбита на 2 части.
2. Последующая оценка погрешности СМО с очередью ограниченной длины производится аналогично, однако порядок величин и
изменяться не будет. Ограничение максимальной длины очереди разумно ограничить слева
, справа средней длиной реальной очереди модели без ограничений, умноженной на
. Список характеристик, наилучших сочетаний входных параметров и максимальных погрешностей данного испытания показаны в табл. 2, где в столбце «Огр». указано ограничение максимального числа требований в очереди.
Таблица 2
Сводная таблица испытания № 2
Характеристика |
|
Огр. |
Погр. (%) |
Интенсивность потока заявок |
|
|
|
Интенсивность потока обслуживания |
|
|
|
Приведённая интенсивность потока заявок |
|
|
|
Вероятность обслуживания |
|
|
|
Вероятность простоя |
|
|
|
Вероятность ожидания |
|
|
|
Вероятность отказа |
|
|
|
Коэффициент загрузки |
|
|
|
Коэффициент простоя |
|
|
|
Среднее число требований под обсл. |
|
|
|
2 ЦМ числа требований под обсл. |
|
|
|
Число требований в очереди |
|
|
|
2 ЦМ числа требований в очереди |
|
|
|
Число требований в реальной очереди |
|
|
|
2 ЦМ числа требований в р. очереди |
|
|
|
Число требований в СМО |
|
|
|
2 ЦМ числа требований в СМО |
|
|
|
Время обслуживания |
|
|
|
2 ЦМ времени обслуживания |
|
|
|
Время ожидания в очереди |
|
|
|
2 ЦМ времени ожидания в очереди |
|
|
|
Время ожидания в реальной очереди |
|
|
|
2 ЦМ времени ожидания в р. очереди |
|
|
|
Время пребывания заявки в СМО |
|
|
|
2 ЦМ времени пребывания |
|
|
|
Ковариация числа заявок под обсл. и в оч. |
|
|
|
Коэффициент корреляции |
|
|
|
Примечание: из результатов исключена нулевая погрешность средней длины реальной очереди при ограничении максимальной длины очереди равной или
, для исключения влияния этого вырожденного случая на оценку погрешности. На основании данных результатов можно сделать вывод, что сходимость моментов порядка выше 1 медленнее и требует увеличить число требований в системе.
3. Основываясь на предыдущей оценке можно выбрать 5 наилучших сочетаний входных параметров для каждого момента порядка выше 1, что в результате даст различных уникальных сочетаний, в виду повторения некоторых из них. Во втором испытании также будет использован встроенный псевдослучайный ГСЧ Visual Studio 2010, число прогонок для каждого набора увеличится до
, максимальное число требований же возрастет до
, и оцениваться будут только моменты порядка выше 1. Список характеристик, наилучших сочетаний входных параметров и максимальных погрешностей данного испытания показаны в табл. 3.
Таблица 3
Сводная таблица испытания № 3
Характеристика |
|
Огр. |
Макс. (%) |
2 ЦМ числа требований под обсл. |
|
|
|
2 ЦМ числа требований в очереди |
|
|
|
2 ЦМ числа требований в р. очереди |
|
|
|
2 ЦМ числа требований в СМО |
|
|
|
2 ЦМ времени обслуживания |
|
|
|
2 ЦМ времени ожидания в очереди |
|
|
|
2 ЦМ времени ож. в р. очереди |
|
|
|
2 ЦМ времени пребывания |
|
|
|
Ковариация числа заявок |
|
|
|
Коэффициент корреляции |
|
|
|
Наилучшие сочетания входных параметров в целом показаны в табл. 4.
Таблица 4
Наилучшее сочетание входных параметров в целом
Тип характеристик |
|
Погрешность (%) |
Моменты 1-го порядка |
|
|
Момента порядка выше 1 |
|
|
Модель многоканальной СМО с отказами следует изучить отдельно.
Оценка погрешности данной системы производится аналогично, однако порядок величин и
изменяться не будет, т. к. наилучший порядок уже был найден в испытании первой модели. Список характеристик, наилучших сочетаний входных параметров и максимальных погрешностей данного испытания показаны в табл. 5.
Таблица 5
Сводная таблица испытания № 4
Характеристика |
|
Погр. (%) |
Интенсивность потока заявок |
|
|
Интенсивность потока обслуживания |
|
|
Приведённая интенсивность потока заявок |
|
|
Вероятность обслуживания |
|
|
Вероятность простоя |
|
|
Вероятность отказа |
|
|
Коэффициент загрузки |
|
|
Коэффициент простоя |
|
|
Среднее число требований под обслуживанием |
|
|
2 ЦМ числа требований под обслуживанием |
|
|
Число требований в СМО |
|
|
Время обслуживания |
|
|
2 ЦМ времени обслуживания |
|
|
Производить дальнейшие оценки данной модели, равно как и увеличивать максимальное число требований в системе, нет необходимости, динамические величины сходятся к значениям аналитических формул с приемлемой точностью. Наилучшие сочетания входных параметров в целом показаны в табл. 6.
Таблица 6
Наилучшее сочетание входных параметров в целом
Тип характеристик |
|
Погрешность (%) |
Моменты 1-го порядка |
|
|
Момента порядка выше 1 |
|
|
Исследовано поведение модели [3] при различных наборах входных параметров, что позволило выявить наиболее удачные из этих сочетаний для созданных моделей СМО и отдельных характеристик, использование которых позволит снизить влияние ГСЧ на результаты прогонки. Результаты данного исследования будут полезны при поиске неизвестных пока аналитических формул характеристик СМО открытого типа с ограничениями.
Литература:
1. А. П. Кирпичников, Методы прикладной теории массового обслуживания. Казань, Изд-во Казанского университета, 2011. 200 с.
2. Microsoft Developer Network, Руководство по программированию на C#. https://msdn.microsoft.com/, 2014.
3. Шемахин Е. Ю., Кирпичников А. П. Моделирование многоканальных открытых систем массового обслуживания с ограничениями в среде Visual Studio 2010, Вестник Казанского технологического университета. 2015. Т.18, № 3.