К вычислению вероятностных характеристик СМО ограниченной ёмкости со случайными требованиями к ресурсам Текст научной статьи по специальности «Компьютерные и информационные науки»

Математическая теория телетрафика и сети телекоммуникаций

УДК 519.21;51-74

DOI: 10.22363/2312-9735-2017-25-3-209-216

К вычислению вероятностных характеристик СМО ограниченной ёмкости со случайными требованиями к

ресурсам О. Г. Вихрова

Кафедра прикладной информатики и теории вероятностей Российский университет дружбы народов ул. Миклухо-Маклая, д. 6, Москва, Россия, 117198

Для современных сетей связи характерен высокий уровень роста мобильного трафика данных. Устойчивые тенденции роста нагрузки в беспроводных сетях ускоряют развитие технологий и переход к сетям нового поколения (5G). Планируемые улучшения позволят в несколько раз увеличить пропускную способность каналов связи и позволят устройствам одновременно поддерживать как соединения сотой связи, так и, например, подключение к Wi-Fi сетям, возможна передача данных от устройства к устройству напрямую (device-to-device, D2D). В условиях вынужденной гетерогенности сетей связи предлагается отказаться от традиционной «парной» ассоциации восходящего (ВК) и нисходящего (НК) каналов и разделять их при условии гарантии необходимого уровня качества. Разделение ресурсов в современных гетерогенных сетях предлагается моделировать в виде системы массового обслуживания (СМО) со случайными требованиями. Подобные модели к анализу показателей качества в беспроводных сетях ранее не применялись. Исследуется многолинейная СМО с различными классами заявок, где каждой поступившей заявке выделяется некоторый вектор случайных требований к ресурсам. Было доказано, что при объединении потоков заявок различных классов в один поток со средневзвешенным требованием стационарные вероятности не зависят от порядка поступления заявок, а зависят от их общего числа в системе и объёма занимаемых ресурсов. Получен более простой вид формул для вероятности блокировки и среднего объёма занятых ресурсов, однако аналитические формулы требуют вычисления n—кратных свёрток для всех возможных наборов векторов занимаемых ресурсов, где п — количество заявок в системе. Был разработан эффективный алгоритм вычисления нормировочной константы, с помощью которой получены рекуррентные формулы для стационарных вероятностей и основных вероятностных характеристик СМО.

Ключевые слова: ресурсная СМО, ограниченная ёмкость, требования случайного объёма, нормировочная константа, вероятностные характеристики, рекуррентный алгоритм, гетерогенная сеть

1. Введение

Стремительный рост популярности мобильных приложений, по мнению представителей инфокоммуникационной отрасли, способствует активному развитию технологии 4G и уже к 2020 году приведёт к полномасштабному развёртыванию инфраструктур для технологии поколения 5G [1]. Изменения затронут физический уровень с целью повышения скорости передачи данных, а также коснутся сетевой архитектуры, современные мобильные устройства могут работать в двух режимах — поддерживать сотовое соединения и подключаться к частным сетям малых сот

Статья поступила в редакцию 13 апреля 2017 г.

Исследование выполнено при финансовой поддержке Минобрнауки России (соглашение № 02.а03.21.0008) и РФФИ в рамках проектов № 15-07-03051, 16-07-00766.

(микро-, пико-, фемтосот) [2]. Новый класс прямых соединений устройств Б2Б, которые могут осуществляться как внутри полосы пропускания, так и вне выделенной полосы, позволят производить обмен данными минуя базовые станции и не создавать дополнительную нагрузку на существующие каналы [3]. Все предложенные изменения подразумевают более эффективное использование и разделение ресурсов и позволят существенно улучшить пропускную способность беспроводных каналов передачи данных.

Принцип гетерогенности сетей, согласно которому происходит наложение зон покрытия макросот и малых сот, оказался самым доступным и бюджетным способом улучшить пропускную способность сетей ЬТЕ. Концепция разделения ВК и НК была предложена в [4] и, по мнению авторов, должна повысить эффективность использования радиоресурсов беспроводных сетей. Приведённые в [5,6] аналитические исследования и результаты имитационного моделирования свидетельствуют о существенном выигрыше в производительности сети при раздельной ассоциации.

Существующие модели разделения радиоресурсов в гетерогенных сетях позволяют определить параметры сети при фиксированном количестве пользовательских устройств и не принимают во внимание процессы установления новых соединений и завершения существующих. Чтобы иметь возможность наблюдать за показателями качества сети в динамике, в [7] предложена модель разделения ресурсов со следующими предположениями:

1) поступивший запрос на установление соединения с базовой станцией требует некоторый случайный объем ресурсов нескольких типов;

2) если для обслуживания соединения недостаточно ресурсов, то запрос будет отклонён;

3) устройство занимает требуемый объем ресурсов на все время обслуживания и освобождает его по завершении обслуживания.

Модель с требованиями случайного объёма была представлена в [8], её анализ в общем виде предложен автором в [9]. В [10] исследуется общая модель с ограниченным ресурсом и требованиями случайного объёма к различным типам ресурсов. Ресурсы занимаются на все время обслуживания, а по его окончании освобождается ровно тот объем ресурсов, который занимался. Случайный процесс, описывающий поведение системы, должен учитывать объёмы занятых ресурсов каждой из заявок, что существенно усложняет его анализ. Было предложено упрощение, согласно которому по завершении обслуживания заявки освобождается некий случайный вектор ресурсов, не зависящий от поведения системы в прошлом при заданной числе заявок и векторе занятых ресурсов.

В [11] и [12] с помощью имитационного моделирования были получены вероятностные характеристики исходной и упрощённой моделей, а также были найдены стационарные вероятности в предположении о пуассоновском входящем потоке заявок и экспоненциальном времени их обслуживания. Было доказано, что полученные значения стационарных вероятностей для обеих моделей эквиваленты. Представленные в работе результаты являются продолжением исследований в [13,14].

2. Математическая модель

Рассматривается многолинейная СМО ограниченной ёмкости. Число приборов для обслуживания в системе не превышает N, а доступный ресурс ограничен вектором К = (Дх,..., Км). В систему поступает поток заявок Ь классов, каждой из заявок класса I для обслуживания необходимо выделить вектор гI = (гц,... ,гм1) случайных требований к ресурсам М типов. Все значения векторов гI независимы от входящего потока и независимы в совокупности. Моменты поступления заявок 1-го класса подчинены пуассоновскому распределению с интенсивностью Аг, а времена обслуживания заявок независимы от процесса поступления заявок, независимы в совокупности и распределены экспоненциально с параметром ^.

Заявка 1-го класса займёт г i единиц ресурсов с вероятностью pi,ri. Если выра-

(к) Д (к-1) , „

жение р) r = Pi jР) \ является к-кратной свёрткой вероятностей pi r, тогда все j=o ' 'j

ki заявок класса I займут г единиц ресурсов с вероятностью р^. В [14] было доказано, что распределение стационарных вероятностей рассматриваемой СМО имеет вид (1) с начальным значением (2):

__ki къ

^ (r1,...,ri) = 9o £ « ...РЙ^ (1)

к\ +-----bk_L=n

/ N kt kL \-1

*> = ЕЕ E « (2)

V=0^ +-----hkL=n ri+-----hrL =R 1 L J

Введём функцию gk(г), просуммировав вероятности qk ¿(ri,..., г^) по всем векторам занятых ресурсов г = г1 + ... + г L:

9k (г)= Е Е П р^ ib. (3)

r!+...+r;=r kt +... + kL=k 1 = 1

Теорема 1. Распределение стационарных вероятностей СМО с объединённым

L L

входящим потоком со средневзвешенным требованием pr = Y1 тгViri, где р = Y1 Pi>

1=1 р 1=1

задаётся формулами:

nk fN R nk \-1 ,k(г) = .oP¥dk), Яo = ЕЕ ^rkM . (4)

' \k=o r=o ' )

Доказательство. Просуммируем обе части равенства (1) по всем векторам г i, г1 + ... + г l = г

Е С.^п,..,ГL)= Е Е £...%,,

ri+...+rL=r ri+...+r_L=r ^ +...+kL=k

Е ^(п,..г^ = ?o Е Е П«^

ri+...+r^=r ri+...+r;=r kl + ... + k_L=k г = 1

L Ь

Qk(г) = qo Е Е П PiS JTT, qk(г) = Qo§k(г).

ri+...+rL=r ki+...+kL=k г=1 l'

Покажем теперь, что gk(г) является k-кратной свёрткой вероятностей рг и gk(г) = рГь). Пусть zr = (zp), тогда производящая функция П(z) для §k(г) имеет вид

R R / l ki \

П (z) = Е 9k (г) zr = Е Е Е П ^ Ir

r=o r=o \ri+...+r _L=r k\+...+kL=k 1=1 '■/

zr

= £ £ Е П |г= Е П ^"?(■) =

( е Ап(■))")=й =й (Е^«У■

\к1+...+кь=к 1 Ь г=1 / \г=1 / \г=1 ^ /

а г (к )

где = ^ X] Р1г>7,Г1 также является производящей функцией. Произведение к

г=0 гг=0 ' 1

функций Ж1(z) является к-кратной свёрткой средневзвешенного требования рг, т.е.

/ , Л (к)

„к / Ь \ У 7 ^

к / Ь \ (к) к

»«=йТ (е >•) =

Теорема 1 позволяет получить аналитические формулы для вероятности блокировки В и среднего объёма занятых ресурсов Ь:

N-1 п Я N Л

В = 1 - <70 £ п£й("+1), Ь = ,0 £ п£^ (5)

п=0 j п=0 j

3. Рекуррентный алгоритм

Формулы (5) требуют вычисления й-кратных свёрток вероятностей рг для всех возможных векторов г^И. Чтобы понизить вычислительную сложность формул, был разработан алгоритм для вычисления нормировочной константы С (М, И) = 00-1, аналогичный алгоритму в [15]. С помощью нормировочной константы были получены рекуррентные формулы для нахождения вероятностных характеристик

СМО.

Введём обозначение:

с(п,г) = ££и(к)й-, п > о, г > о.

к=0j=0

Теорема 2. Функция С (п, г) удовлетворяет рекуррентному соотношению (6) с начальными условиями (7):

г

С (п, г) = С (п - 1, г) + пй £ Pj (С (п - 1, г - ^ -С (п - 2, г - j)), (6)

п' j=о

Г

С (0, г) = 1, г > 0, С (1, г) = 1 + р£Pj. (7)

j=о

Доказательство. Зададим шаг индукции:

N Я.

С (п, г) - С (п - 1, г) = £ £ А н п - £ £ л(и) п = п £ *(п) =

п=0 г=0 ' п=0 г=0 ' ' г=0

п £ £ -V-" = п £ * £ й-.(-1) = п £» ((Ь), £ 1) 1 =

3=0 1=0 1=0 1=0 ' j=0

рп г

- - - ■ (п - 1, г - 1) - С (п - :

(С (п - 1, г -1) -С (п - 2, г -1))

п

j=о

Следствие 1. Вероятность блокировки системы В вычисляется по формуле:

Я

В = 1 -С-1 И) £ (^ - 1, И -1). (8)

j=о

Доказательство. Получим из (4) рекуррентную формулу для вероятности блокировки СМО с учётом (6):

N-1 рп а N-1 я j

В = 1 - <70 £ п Е«(П+1) = 1 - С-1 (", И) Е п ^^й«-1(п) =

п=0 j=0 п=0 j=0 1=0

Я N-1 п Я-1 Я

= 1 - С-1 И) £й £ п Ей(п) = 1 - С-1 и) Е«С (^ - 1, И -1).

1=0 п=0 ' j=0 j=0

Следствие 2. Средний объем занятых ресурсов Ь вычисляется по формуле:

м нт

Ь = И-С-1 И) £ ет И - гего). (9)

т=1 г=1

Доказательство.

ЯМ Я N ЯМ Я N

Ь = ££ (Я = ££ е^ЕЕС-1 (*, И) п =

г=0 т=1 j=г п=0 г=0 т=1 j=r п=0

Я М / Я N п Я-г N п N п

= С-1 № К) £ £ е„, (е е ^ -££п*Т + Е^п

г=0 т=1 V j=0 п=0 j=r п=0 п=0

ЯМ / N рп \

= С-1 И)££ ет С И)-С И - г) + £ п^п =

г=0 т=1 \ п=0 /

м нт

= И-С-1 И) £ ет И - гет).

т=1 г=1

Следствие 3. Второй момент числа занятых ресурсов может быть найден по формуле:

Г Я

Ь + С-1 (Л", И) £ г (С (Л", И) - С (Л", г))

г=0

Ь(2) = 2

(10)

где Ь2 = (Ь12,..., Ьм2) • Тогда дисперсия объёма занятых ресурсов а2 = Ь(2) — Ь2. Доказательство.

я. я. N П П N

b(2) = 2£ rEE = ^ ^Е G-1 w R) ^

r=0 j=r j=0 r=0 j=r j=0

R / R N R-r N N

2G-1 (N. R) E r (EE ^ - E E ^ + E ^

r=0 \j=r j=0 j=0 j=0 j=0

R ( N oJ \

= 2G-1 (N, R) E M G (N, R) - G (N, R - J) + E ^\PJ =

r=0 V j=0 ' /

R

b + G-1 (N, R) E r (G (N, R) - G (N, r))

= 2

r (

r=0

4. Заключение

Многолинейную СМО ограниченной ёмкости с вектором случайных требований к ресурсам и несколькими классами заявок предлагается анализировать в виде СМО с объединённым входящим потоком средневзвешенных требований без потери точности. В качестве эффективного метода расчёта стационарных вероятностей и вероятностных характеристик упрощённой СМО был разработан рекуррентный алгоритм вычисления нормировочной константы и получены рекуррентные формулы для вероятности блокировки и среднего объёма и дисперсии занятых ресурсов.

Литература

1. Cisco Visual Networking Index: Global Mobile Data Traffic Forecast Update, 2016-2021. — 2017. — http://www.cisco.com/c/en/us/ solutions/collateral/service-provider/visual-networking-index-vni/ mobile-white-paper-c11-520862.html.

2. Femtocells: Past, Present, and Future / J. Andrews, H. Claussen, M. Dohler, S. Rangan // IEEE JSAC. — 2012. — Vol. 30, No 3. — Pp. 497-508.

3. Device-to-Device Communications for National Security and Public Safety / G. Fodor, S. Parkvall, S. Sorrentino, P. Wallentin, Q. Lu, N. Brahmi // IEEE Access. — 2014. — Vol. 2, No 1. — Pp. 1510-1520.

4. Why to Decouple the Uplink and Downlink in Cellular Networks and How To Do It / F. Boccardi, J. Andrews, H. Elshaer, M. Dohler, S. Parkvall, P. Popovski, S. Singh // IEEE Communications Magazine. — 2016. — Vol. 54, No 3. — Pp. 110-117.

5. Singh S., Zhang X, Andrews J. Joint Rate and SINR Coverage Analysis for Decoupled Up-link-Downlink Biased Cell Associations in HetNets // IEEE Transactions on Wireless Communication. — 2015. — Vol. 14, No 10. — Pp. 5360-5373.

6. Downlink and Uplink Decoupling: A Disruptive Architectural Design for 5G Networks / H. Elshaer, F. Boccardi, M. Dohler, R. Irmer // Global Communications Conference (GLOBECOM). — 2014. — Pp. 1798-1803.

7. LTE Performance Analysis Using Queuing Systems with Finite Resources and Random Requirements / V. Naumov, K. Samouylov, E. Sopin, N. Yarkina, S. Andreev, A. Samuylov // 7th International Congress on Ultra Modern Telecommunications and Control Systems and Workshops, 6-8 Oct. 2015. — IEEE, 2015. — Pp. 100-103.

8. Тихоненко О. М. Определение характеристик систем обслуживания с ограниченной памятью // Автоматика и телемеханика. — 1997. — Вып. 6. — С. 105-110.

9. Тихоненко O. M. Обобщенная задача Эрланга для систем обслуживания с ограниченным суммарным объемом // Проблемы передачи информации. — 2005. — Т. 41, вып. 3. — С. 64-75.

10. Наумов В. А., Самуйлов К. Е. О моделировании систем массового обслуживания с множественными ресурсами // Вестник РУДН. Серия: Математика. Информатика. Физика. — 2014. — № 3. — С. 58-62.

11. Наумов В. А, Самуйлов К. Е., Самуйлов А. К. О суммарном объеме ресурсов, занимаемых обслуживаемыми заявками // Актоматика и телемеханика. — 2016. — Вып. 8. — С. 125-135.

12. Two Approaches to Analysis of Queuing Systems with Limited Resources / V. Nau-mov, K. Samuoylov, E. Sopin, S. Andreev // Ultra Modern Telecommunications and Control Systems and Workshops, 6-8 Oct. 2014. — IEEE, 2014. — Pp. 485-488.

13. Samouylov K., Sopin E., Vikhrova O. Analyzing Blocking Probability in LTE Wireless Network via Queuing System with Finite Amount of Resources // Communications in Computer and Information Science. — 2015. — Vol. 564. — Pp. 393-403.

14. Самуйлов K., Сопин Э., Вихрова О. К разработке эффективных вычислительных алгоритмов нахождения вероятности блокировки для системы со случайными требованиями // Информационные технологии и математическое моделирование (ИТММ-2016): 15-я Международная конференция имени А. Ф. Терпу-гова. — Т. 1. — Изд-во Том. ун-та, 2016. — С. 192-197.

15. Buzen J. P. Computational Algorithms for Closed Queueing Networks with Exponential Servers // Communications of the ACM. — 1973. — Vol. 19, No 9. — Pp. 527-531.

UDC 519.21;51-74

DOI: 10.22363/2312-9735-2017-25-3-209-216

About Probability Characteristics Evaluation in Queuing System with Limited Resources and Random Requirements

O. G. Vikhrova

Department of Applied Probability and Informatics Peoples' Friendship University of Russia (RUDN University) 6 Miklukho-Maklaya St., Moscow, 117198, Russian Federation

Mobile data traffic increases on everyday basis for the last decade and is going to keep this trend in the near future. Exponential growth of data traffic in wireless networks accelerates the development of new technologies and the transition to 5G networks. The ongoing improvements will allow to increase the channels throughput and will allow devices to simultaneously support both cellular and Wi-Fi networks, and even allow direct device-to-device (D2D) connection without any base station involved. Evolving heterogeneous networks promise more efficient radio resources usage by a macro cell traffic offloading to the small cells and uplink and downlink decoupling (DUDe). A resource-sharing model in heterogeneous networks is for the first time proposed to be analyzed in terms of queuing system with random requirements. We suggest a multiserver queuing network with limited resources where each class of customers requires a random vector of resources to be served. It has been proved that stationary probabilities of the system with aggregated flow of customers with mean-weighted requirement are equal to the stationary probabilities of the suggested system. The analytical method for the key probability characteristics evaluation requires calculating all fc-fold convolutions for each set of vectors requirements. We propose a recurrent computation algorithm for normalization constant evaluation and efficient formulas for blocking probability, mean volume and variance of the occupied resources.

Key words and phrases: queuing system with limited resources, random requirements, recurrent algorithm, probability characteristics, normalization constant, heterogeneous networks

References

1. Cisco Visual Networking Index: Global Mobile Data Traffic Forecast Update, 2016-2021 (2017).

URL http://www.cisc0.c0m/c/en/us/s0luti0ns/c0llateral/

service-provider/visual-networking-index-vni/ m0bile-white-paper-c11-520862.html

2. J. Andrews, H. Claussen, M. Dohler, S. Rangan, Femtocells: Past, Present, and Future, IEEE JSAC 30 (3) (2012) 497-508.

3. G. Fodor, S. Parkvall, S. Sorrentino, P. Wallentin, Q. Lu, N. Brahmi, Device-to-Device Communications for National Security and Public Safety, IEEE Access 2 (1) (2014) 1510-1520.

4. F. Boccardi, J. Andrews, H. Elshaer, M. Dohler, S. Parkvall, P. Popovski, S. Singh, Why to Decouple the Uplink and Downlink in Cellular Networks and How To Do It, IEEE Communications Magazine 54 (3) (2016) 110-117.

5. S. Singh, X. Zhang, J. Andrews, Joint Rate and SINR Coverage Analysis for Decoupled Up-link-Downlink Biased Cell Associations in HetNets, IEEE Transactions on Wireless Communication 14 (10) (2015) 5360-5373.

6. H. Elshaer, F. Boccardi, M. Dohler, R. Irmer, Downlink and Uplink Decoupling: A Disruptive Architectural Design for 5G Networks, Global Communications Conference (GLOBECOM) (2014) 1798-1803.

7. V. Naumov, K. Samouylov, E. Sopin, N. Yarkina, S. Andreev, A. Samuylov, LTE Performance Analysis Using Queuing Systems with Finite Resources and Random Requirements, in: 7th International Congress on Ultra Modern Telecommunications and Control Systems and Workshops, 6-8 Oct. 2015, IEEE, 2015, pp. 100-103. doi:10.1109/ICUMT.2015.7382412.

8. O. M. Tikhonenko, Destricted Capacity Queueing Systems: Determination of their Characteristics , Automation and Remote Control 58 (1997) 969-973.

9. O. M. Tikhonenko, Generalized Erlang Problem for Service Systems with Finite Total Capacity, Problems of Information Transmission 41 (2005) 243-253.

10. V. A. Naumov, K. E. Samouylov, On the Modeling of Queuing Systems with Multiple Resources, Bulletin of Peoples' Friendship University of Russia. Series: Mathematics. Information Sciences. Physics (3) (2014) 58-62.

11. V. A. Naumov, K. E. Samouylov, A. K. Samouylov, On the Total Amount of Resources Occupied by Serviced Customers, Automation and Remote Control 77 (2016) 1419-1427.

12. V. Naumov, K. Samuoylov, E. Sopin, S. Andreev, Two Approaches to Analysis of Queuing Systems with Limited Resources, in: Ultra Modern Telecommunications and Control Systems and Workshops, 6-8 Oct. 2014, IEEE, 2014, pp. 485-488. doi:10.1109/ICUMT.2014.7002149.

13. K. Samouylov, E. Sopin, O. Vikhrova, Analyzing Blocking Probability in LTE Wireless Network via Queuing System with Finite Amount of Resources, Communications in Computer and Information Science 564 (2015) 393-403.

14. K. Samouylov, E. Sopin, O. Vikhrova, On Design of Efficient Algorithm for Blocking Probability Calculation in Queuing System with Random Requirements, in: Information technologies and mathematical modelling (ITMM-2016): XV International Scientific Conference named after A. F. Terpugov, Vol. 1, Tomsk State University, 2016, pp. 192-197.

15. J. P. Buzen, Computational Algorithms for Closed Queueing Networks with Exponential Servers, Communications of the ACM 19 (9) (1973) 527-531.

© Вихрова О. Г., 2017