CC BY
25МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ, КОМПЬЮТЕРНЫЕ ТЕХНОЛОГИИ В ТЕОРИИ УПРАВЛЕНИЯ СЛОЖНЫХ ПРОЦЕССОВ
ОБ ОДНОЙ ЗАДАЧЕ СИНТЕЗА ЗАМКНУТОЙ СИСТЕМЫ МАССОВОГО ОБСЛУЖИВАНИЯ
А.А. Таранцев, доктор технических наук; П.А. Манин;
О.В. Щербаков, доктор технических наук, профессор. Санкт-Петербургский университет ГПС МЧС России
Рассмотрена возможность нахождения числа каналов обслуживания и их производительности для замкнутой стандартной многоканальной системы массового обслуживания по требованиям к её характеристикам. Приведён способ уточнения параметров системы при ограниченной готовности каналов обслуживания.
Ключевые слова: система массового обслуживания, источник заявок, канал обслуживания
ABOUT THE PROBLEM OF SYNTHESIS OF THE CLOSED SERVICE SYSTEM TAKING INTO ACCOUNT READINESS OF CHANNELS
A.A. Tarantsev; P.A. Manin; O.V. Sherbakov.
Saint-Petersburg university of State fire service of EMERCOM of Russia
The possibility of location of number of service channels and their productivity for the closed standard multichannel service system under requirements to its characteristics is considered. The method of specification parameters of system is resulted at the limited readiness of service channels.
Key words: service system, source of orders, service channel
На практике широкое распространение получили замкнутые системы массового обслуживания (СМО), используемые, например, для моделирования процессов поддержания в готовности нескольких объектов - станков, автомобилей, серверов и т.п. [1, 2]. Такие объекты являются источниками заявок - сигналов, вызовов, сообщений и др., которые обслуживаются операторами, ремонтниками, аварийно-наладочными бригадами, именуемыми каналами обслуживания.
Применительно к СМО решаются следующие задачи:
- анализ (по параметрам СМО определение её характеристик);
- синтез (по заданным характеристикам СМО и некоторым параметрам определяются остальные параметры СМО);
- идентификация (по результатам наблюдения за входным и выходным потоками заявок определение некоторых параметров и характеристик СМО) [3].
Первая задача (анализ) не представляет особой трудности, третья (идентификация) встречается довольно редко, а второй задаче (синтезу) следует уделить особое внимание применительно к замкнутой СМО. Тем более, что для незамкнутой СМО такая задача в графическом виде была решена в работе [4].
Математическое описание замкнутой стандартной СМО
Под стандартной замкнутой СМО будем понимать систему, содержащую N источников одинаковых простейших потоков заявок, обслуживаемых п каналами обслуживания. Причём N>n, очередь заявок организована по принципу «раньше пришёл, раньше обслужился», время обслуживания случайное и подчинено экспоненциальному закону с параметром (скоростью) ц=об-1, взаимопомощь между каналами обслуживания отсутствует, заявки от каждого источника поступают с частотой у=п-1, а отношение среднего времени обслуживания заявок к среднему времени поступления Р=об/^п является приведённой нагрузкой.
Особый интерес представляет именно стандартная СМО, поскольку при тех же параметрах N п и в её основные характеристики (вероятности незанятости каналов обслуживания р0 и полного останова рп) не отличаются от характеристик замкнутой СМО с другими законами поступления-обслуживания (эрланговскими потоками заявок и времени обслуживания, подчинённому распределению Стьюдента [5]), а математическое описание значительно проще.
Стандартная замкнутая СМО с N источниками и п каналами обслуживания может пребывать в N+1 состояниях, вероятности которых для установившегося случая имеют вид [1]:
р. =1
С^РР^г е[1,п]
N !пп [п!( N - г) !]-1 (в / п)гр0, г е[п +1, Ы]' (1)
п N
Р0-1 = 1 + + X в'С^г!(ппп!)-1, (2)
г=1 г=п+1
где 1 = N! [г!(N - г)!] - число сочетаний из N по г.
Принято различать следующие основные характеристики замкнутой СМО:
- вероятность незанятости всех каналов (ни один источник не подал заявок, полная готовность СМО) - р0;
- вероятность полного останова СМО (все источники подали заявки, п обслуживается, N-n - ожидают):
Рп = Р^ (3)
- вероятность немедленного принятия заявки к обслуживанию (свободен хотя бы один канал обслуживания):
Рн = Р0+Р1+ ...+Рп-1; (4)
средняя величина очереди:
N
т
= П(г - п)Р.; (5)
оч
г=п+1
нагрузка на канал:
р = Щ/п.
(6)
Задача анализа ставится и решается следующим образом: по заданным параметрам СМО - N п и в требуется оценить характеристики СМО. Для этого сначала с использованием выражений (1) и (2) определяются вероятности состояний СМО, а затем по выражениям (3)-(6) - её характеристики.
Задача синтеза замкнутой стандартной СМО
Задача синтеза замкнутой СМО может иметь следующую постановку: известны число источников заявок N и частота V и заданы требования к характеристикам СМО:
ро > род (7)
Рп < рпд, (8)
Рн > рнд, (9)
т0ч< точд, (10)
Р < Рд (11)
д д д д д
где р0 , рп , рн , точ , р - допустимые значения соответствующих характеристик.
Требуется подобрать такие значения параметров п и р, чтобы выполнялись условия (7)-(11).
Задачу синтеза СМО, в отличие от задачи анализа, можно решать только численными методами, что не всегда удобно с точки зрения практики. Для преодоления этой трудности, как и в работе [4], целесообразно использовать номограммы.
Такая номограмма для нормативного значения рп=0,001 [6] приведена на рис. 1 (подобные номограммы могут быть построены и для других значений рп). На ней в координатах характеристик «вероятность незанятости каналов р0 - вероятность немедленного принятия заявки рн» нанесены линии равных значений параметров N и п, а также характеристик точ и р. В табл. даны значения, приведённой нагрузки в для той же вероятности рп=0,001.
Пользоваться номограммой для синтеза замкнутой стандартной СМО можно следующим образом:
- по заданным ограничениям на характеристики р0д, рнд, точд и рд (рис. 1) выделяется рабочая область, в пределах которой на пересечении линий N и п определяется рабочая точка А - соответствующее ей число каналов обслуживания п и является требуемым (синтезированным) параметром СМО;
- по величинам N и п по табл. определяется величина приведённой нагрузки в;
- по величинам параметров в и V находится другой параметр - требуемая величина скорости обслуживания:
и= у/в; (12)
- по величинам параметров в, N и п с использованием выражений (1)-(6) находятся величины характеристик СМО и производится проверка выполнения условий (7)—(11).
Это целесообразно пояснить на примере. Пусть некоторая техническая система содержит четыре агрегата (N=4), которые выходят из строя в среднем один раз в сутки ^=1 сут-1 -0,0417 час-1), причём поток отказов простейший. К характеристикам этой системы предъявляются следующие требования:
- вероятность того, что откажут все агрегаты (то есть полный выход из строя технической системы в целом), не должна превышать 0,1 % (рпд=0,001);
- вероятность полной готовности всех агрегатов должна быть не менее 50 % (р0д=0,5);
- вероятность того, что при отказе агрегата к нему немедленно прибудет наладчик для проведения ремонта, должна быть не менее 90 % (рнд=0,9);
- в очереди на проведение ремонта должно быть не более одного агрегата (даочд=1);
- наладчики, исходя из эргономических соображений [7], должны заниматься ремонтными работами не более 40 % продолжительности смены (рд=0,4).
Требуется определить минимальное число наладчиков в смене п и среднюю скорость проведения ремонта ц, если время, затрачиваемое на ремонт, имеет случайный характер и подчинено экспоненциальному закону.
Для решения такой задачи синтеза технической системы, являющейся замкнутой 4- канальной СМО, сначала используем номограмму (рис. 1). На ней выделим рабочую область, исходя из требований к характеристикам, и найдём рабочую точку А на пересечении линий N=4 и п=2 (рис. 2). Таким образом, для обслуживания агрегатов технической системы необходимо не менее двух наладчиков. Затем по табл. для N=4 и п=2 находим Р=0,1565, а из (12) - ц=0,0417/0,1565-0,2662 час-1. Другими словами, каждый наладчик должен суметь отремонтировать агрегат в среднем не менее, чем за ¿об=3,75 часа.
Проверка, проведённая по выражениям (1)-(6) для параметров N=4, п=2 и Р=0,1565 (она же - задача анализа), позволила получить следующие значения характеристик СМО: р0=0,5563; рп=0,001; рн=0,9045; даоч=0,0148; р=0,3130, находящиеся в пределах допусков, что подтверждает правильность решения задачи синтеза.
Следует отметить, что, во-первых, иногда требования к характеристикам СМО могут быть заданы так, что решение задачи синтеза становится невозможным (например, для рассмотренного случая, если бы р0д=0,6). Тогда необходимо проводить корректировку требований к характеристикам СМО.
Во-вторых, надёжность каналов в плане их готовности к обслуживанию заявок может быть не 100-процентной, например, поскольку наладчик может отлучаться с рабочего места. Это количественно может быть учтено через коэффициент готовности:
(13)
где ¿с, ¿от - продолжительность смены и общее время отсутствия наладчика.
Таблица. Значения приведённой нагрузки в для вероятности ^п=0,001
N п в N п в N п в N п в
2 1 0,0229 1 0,1233 1 0,1375 7 0,8012
2 0,0327 2 0,2352 2 0,2700 9 8 0,8455
1 0,0585 6 3 0,3286 3 0,3935 9 0,8662
3 2 3 0,0957 0,1111 4 5 0,3984 0,4425 8 4 5 0,5036 0,5952 1 2 0,1426 0,2834
1 0,0884 6 0,4625 6 0,6642 3 0,4206
4 2 3 0,1565 0,1983 0,2163 1 2 3 0,1320 0,2564 0,3676 7 8 0,7086 0,7292 4 5 0,5518 0,6735 0,7816
4 1 0,1407 10 6
5 1 2 3 4 5 0,1095 0,2029 0,2727 0,3162 0,3355 7 4 5 6 7 0,4601 0,5296 0,5739 0,5943 9 2 3 4 5 6 0,2784 0,4102 0,5329 0,6419 0,7326 7 8 9 10 0,8716 0,9397 0,9839 1,0048
Рис. 1. Номограмма «вероятность незанятости каналов р0 - вероятности немедленного принятия заявки к обслуживанию рн» с линиями равных значений числа источников заявок N числа каналов обслуживания п, величины очереди точ и нагрузки на канал р при вероятности
полного останова рп=0,001
Рис. 2. Пример нахождения числа каналов п при заданных требованиях к характеристикам замкнутой СМО
Учёт готовности каналов обслуживания
Если готовность каналов обслуживания ограниченная и может быть учтена через коэффициент &г<1, то задача анализа СМО решается следующим образом. Вычисляется действующее число каналов:
п
*=
кгп.
(14)
Затем при заданных параметрах N и в, полагая число каналов п переменной величиной в пределах пт„<п*<птах, с использованием выражений (1)-(6) вычисляются характеристики СМО как функции от п. После путём аппроксимации находятся значения характеристик СМО для величины п*, что и будет решением задачи анализа.
Рис. 3. Пример влияния числа каналов п на характеристики замкнутой СМО при N=4 и
в=0,1565 с учётом ограниченной готовности каналов (А - рабочая точка при п=2; А* - рабочая точка при п*=1,6; Б - рабочая точка при п=3;
Б* - рабочая точка при п*=2,4)
Рис. 4. Пример влияния числа каналов п на характеристики замкнутой СМО при N=4 и в=0,085 с учётом ограниченной готовности каналов (А - рабочая точка при п=2; А* - рабочая точка при п*=1,6)
Для условий рассмотренного примера при N=4 и в=0,1565 зависимости характеристик СМО от числа каналов приведены на рис. 3, рабочая точка А при п=2 обозначена пунктирной линией. Если готовность каналов ограничена, например, ^=0,8, то п*=1,6 и рабочая точка А* смещается влево. Тогда р0-0,55; рп-0,0016; рн-0,8; mоч-0,03; р-0,4 - две характеристики рп и рн вышли за пределы допуска.
Решение задачи синтеза СМО при ^<1 начинается с ранее рассмотренного подхода, когда по номограмме (рис. 1) определяются параметры п и в, как если бы готовность каналов была полной. Затем проводится проверка выполнения условий (7)—(11) при n*=kгn. Если характеристики СМО остаются в пределах допусков, задача синтеза считается решённой.
Если хотя бы одно из условий (7)—(11) не выполняется, приходится либо увеличить число каналов при той же скорости обслуживания р, либо увеличить величину р (то есть снизить приведённую нагрузку в) при прежнем числе каналов обслуживания.
В первом случае для условий рассмотренного примера число каналов нужно увеличить до трёх (п=3), оставив в=0,1565. Тогда число действующих КО п*=2,4 и, согласно рис. 3 (новая рабочая точка Б* обозначена штрих пунктирной линией), характеристики СМО войдут в пределы допусков: р0-0,55; рп-0,0008; рн-0,955; mоч-0,015; р-0,26.
Во втором случае, оставив п=2 (п*=1,6), нужно снизить приведённую нагрузку в. Поскольку из (14) следует, что при п*=1,6 критичной является характеристика рн, по выражениям (1), (2) и (4) подбираем новую величину в*=0,085. Из рис. 4 (новая рабочая точка А* показана штриховой линией) видим, что все характеристики вошли в пределы допусков: р0-0,61; рп-0,0008; рн-0,9; mоч-0,01; р-0,2. Но в этом случае скорость ремонта агрегатов каждым наладчиком должна быть увеличена: р*=0,0417/0,085-0,4902 час-1. Другими словами, каждый наладчик должен суметь отремонтировать агрегат в среднем не менее, чем за ^б=2,04 часа. Это может быть достигнуто путём оснащения наладчиков лучшим инструментом, повышением их квалификации и т.п.
Таким образом, приведён способ решения задачи синтеза стандартной замкнутой СМО с использованием номограммы, когда по заданным характеристикам могут быть определены такие параметры, как число каналов обслуживания и приведённая нагрузка, а по ней - необходимая скорость обслуживания заявки. Также показана возможность уточнения параметров СМО при ограниченной готовности её каналов.
Литература
1. Вентцель Е.С. Исследование операций. М.: Сов. радио, 1972.
2. Гнеденко Б.В., Коваленко И.Н. Введение в теорию массового обслуживания. М.: Наука, 1966.
3. Таранцев А. А. Инженерные методы теории массового обслуживания. 2-е изд., перераб. и доп. СПб.: Наука, 2007.
4. Таранцев А.А. О способе выбора параметров системы массового обслуживания с очередью // Автоматика и телемеханика. 1999. № 7. С. 172-176.
5. Таранцев А.А., Эрюжев М.В. Об аналитических соотношениях в замкнутых одноканальных системах массового обслуживания // Известия АН. Теория и системы управления. 2004. № 4. С. 84-91.
6. Городские и сельские телефонные сети. Нормы технологического проектирования. РД 45.120-2000 (НТП 112-2000).
7. Шаровар Ф.И. Автоматизированные системы управления и связь в пожарной охране. М.: Радио и связь, 1987.
