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

ВЕСТНИК ТОМСКОГО ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА 2018 Управление, вычислительная техника и информатика № 44

УДК 519.872.1

DOI: 10.17223/19988605/44/7

Е.В. Мокров, К.Е. Самуйлов

СРЕДА СОВМЕСТНОГО ЛИЦЕНЗИРОВАННОГО ДОСТУПА В БЕСПРОВОДНОЙ СЕТИ КАК СИСТЕМА МАССОВОГО ОБСЛУЖИВАНИЯ С НЕНАДЕЖНЫМ РЕСУРСОМ

Публикация подготовлена при поддержке Министерства образования и науки Российской Федерации (проект № 2.882.2017/4.6).

В терминах теории массового обслуживания рассматривается математическая модель среды совместного лицензионного доступа (LSA, Licensed Shared Access), применяемой в беспроводных сетях LTE (Long Term Evolution). Модель представлена в виде системы массового обслуживания (СМО) с заявками случайного объема и с ненадежным ресурсом. Получены уравнения для общего случая распределения случайного объема заявок. Для дискретного распределения объема заявок с переразыгрыванием выделяемого им ресурса предложен метод анализа вероятностных характеристик.

Ключевые слова: система массового обслуживания; марковский процесс; ненадежный ресурс; заявка случайного объема.

Большая часть работ, посвященных анализу показателей эффективности системы совместного лицензионного доступа LSA, концентрируется либо на экспериментах с измерениями на лабораторном оборудовании [1-3], либо на оптимизации механизмов доступа к лицензируемой полосе частот [4-6]. В работах, исследующих LSA с помощью моделей систем массового обслуживания, как правило, рассматриваются либо различные законы поступления и обслуживания заявок, либо алгоритмы расчета вероятностных характеристик системы [7-11]. Однако в этих работах применяются классические СМО, в которых заявка при обслуживании занимает ровно один прибор. В отличие от указанных выше работ, в данной статье для обслуживания заявки в системе кроме прибора выделяется ресурс, объем которого является случайной величиной с заданным распределением. В отличие от работ [12-15], рассматривается многолинейная система, в которой имеется ресурс, частично выходящий из строя и восстанавливающийся в соответствии с марковским процессом. В этом случае два состояния, в которых может прибывать ресурс, моделируются при помощи надежной и ненадежной составляющих ресурса. Надежную составляющую ресурса можно рассматривать как, например, доступный оператору сети спектр частот, а ненадежную - как арендуемый им спектр, поскольку при использовании системы LSA часть частот оператора арендуется у третьей стороны и может изыматься в случайные моменты на случайное время.

В данной статье получены уравнения для стационарного распределения наличия в системе, находящейся в состоянии s, n заявок, каждая из которых занимает ri, i = l,...n, ресурсов для случая системы с непрерывным ресурсом, а для системы с дискретным ресурсом предложен метод анализа и расчета вероятности эвакуации и вероятности прерывания обслуживания.

1. Построение модели

Рассматривается С-линейная СМО, на которую поступает пуассоновский поток заявок с интенсивностью X. Обслуживание заявок происходит по экспоненциальному закону с параметром р (рис. 1). Система может находиться в одном из двух состояний s, переход между которыми происходит

в случайные моменты времени, распределенные согласно экспоненциальному закону, 5 е {0,1} . С интенсивностью а осуществляется переход из состояния 5 = 1 (состояние приема заявок) в состояние 5 = 0 и с интенсивностью Р - переход из состояния 5 = 0 в состояние 5 = 1 (состояние потери заявок). В нормальных условиях функционирования системы (5 = 1) заявкам доступен ресурс объема Щ, а при отказе (5 = 0 ) его объем уменьшается до величины Щ < Щ.

Рис. 1. Схема СМО с ненадежным ресурсом

В момент поступления заявки в систему для ее обслуживания должен быть выделен один прибор и случайный объем ресурса в соответствии с функцией распределения (ФР) Р(х), причем Р(х) = 0, х < 0. Если для обслуживания поступившей в нормальных условиях функционирования системы заявки нет свободных приборов и / или отсутствует в достаточном объеме свободный ресурс, заявка теряется.

Кроме того, при смене состояния системы с 5 = 1 на 5 = 0 из-за уменьшения объема доступного ресурса могут быть потеряны уже принятые в систему заявки. Это происходит, если в момент перехода в состояние 5 = 0 в системе имеется п заявок, каждая из которых занимает ресурс объема г , а суммар-

п

ный объем занятого ресурса г. = ^г > Щ . Ниже предложены два алгоритма выбора заявок, которые

г=1

покинут систему до окончания обслуживания (будут потеряны).

Алгоритм 1 (А1). Потеря последней поступившей заявки. Алгоритм предлагает в случае нехватки ресурса терять заявки, поступившие в систему последними. Таким образом, будут теряться те заявки, которые провели в системе меньше всего времени, и сохраняться те заявки, которые провели в системе больше времени. Алгоритм начинает работу при смене состояния системы с 5 = 1 на 5 = 0 . Шаги алгоритма А1 приведены ниже.

Шаг 1. Проверить суммарный объем занятого ресурса; если г. < Я0 закончить работу, иначе перейти к шагу 2.

Шаг 2. Найти заявку с минимальным временем пребывания в системе. Если таких заявок несколько, выбрать любую из них.

Шаг 3. Освободить занимаемый выбранной заявкой ресурс и прибор, удалить заявку из системы, перейти к шагу 1.

Алгоритм 2 (А2). Потеря заявки максимального объема. Алгоритм предлагает в случае нехватки ресурса терять заявки, занимающие максимальный объем ресурса. Таким образом, алгоритм минимизирует число потерянных заявок, поскольку каждый раз при потере заявки достигается максимальный выигрыш в освобождаемом ресурсе. Алгоритм также начинает работу при смене состояния системы с 5 = 1 на 5 = 0 и работает аналогично предыдущему.

Шаг 1. Проверить суммарный объем занятого ресурса; если г. < Щ закончить работу, иначе перейти к шагу 2.

Шаг 2. Найти заявку с максимальным объемом занимаемых ресурсов. Если таких заявок несколько, выбрать любую из них.

Шаг 3. Освободить занимаемый выбранной заявкой ресурс и прибор, удалить заявку из системы, перейти к шагу 1.

Таким образом, на шаге 2 алгоритма А1 будет теряться заявка с номером п в соответствии с порядком поступления заявок в систему. Этот же простой принцип реализации можно получить и для алгоритма А2, если в момент принятия в систему новой заявки перенумеровывать заявки в порядке возрастания объема занимаемого ресурса. Однако ввиду того, что подобная перенумерация достаточно сложна в случае неэкспоненциального распределения функции требований к ресурсам [16], далее в этой статье мы будем исследовать систему при использовании алгоритма А1.

2. Уравнения равновесия и вероятностные характеристики системы

- число заявок в системе,

Введем случайный процесс X(¿) = (и(г(7),-,Г (0),5(0), где п(7)

(г (0, •••, Г (0) - вектор объемов ресурсов, занимаемых заявками, 5 (¿)е{0,1} - состояние, в котором находится система в момент 7. По построению процесс X(7) является марковским на пространстве

состояний

Введем следующие обозначения

Рп (Гп) = р{п(7) = п(г (7),-,Гп (7)) = (?!,•..,гп),5(7) = 1} ,

Чп (Г1, — Гп ) = р{п (7) = п (Г1 (7^ Гп (7)) = (Г1, — Гп ), 5 (7) = 0} .

В случае, когда ФР Г(х) непрерывна с плотностью f (х) , система уравнений равновесия (СУР) для системы с алгоритмом потери заявок А1 может быть записана в следующем виде:

Чо ( ^ (") + Р ) = Ц | Ч\ (х) Лх + а

2 | Рь (Лхк) + Ро

к=1 х .<"

"1

, ро (^(") + а) = Ц|р (х)Лх + РЧо,

Чп ( Г1,-, Гп )( ^ ("о - х.) + + Р ) = 4,-1 (Г1 Гп-1 ) f (Гп ) +

"„-г.

+^2 | Чп+1 (Г1 Г'

х, Г+1 Гп ) Лх -

г=о о

+а

с

2 | Рк (Г1 Гп А ,-,Лхк-п ) + Рп ОЬ-, Гп )

к=п+1 х <"-г.

V Г+ х1й "о у

п = 1,-, С -1, г. < ",

Рп (Гп )(^ ("1 - х.) + + а) = ^Рп-1 (Гп-1) f (Гп ) +

п ) ",-г .

+^2 I Рп+1 (Г1Г ' х' Г+1 Гп ) Лх + РЧп (Г1,-, Гп ) ,

г=о о

п = 1,-,С -1, г. < ", Чс (Г1,—Гс)(СЦ+Р) = ^Чс-1 (Г1,—Гс-1) f (Гс) + аРк (Г1,-,Гс) Г• <"о. Рп (^ • ^Гс )(сЦ + а) = ^Рс-1 (Г1,•■■, Гс-1) f (Гс )+ РЧс (г1,-5Гс ), Г• < "о.

о

о

В качестве основной вероятностной характеристики системы рассмотрим вероятность успешной эвакуации Р{Х(г) = (п,(г ,...,г ),0) | X(г - А) = (п,(г,...,г ),1)} как вероятность того, что при изъятии ресурса ни одна заявка не будет потеряна:

а

Р =

с

21

а

+ АР(Щ - г.) + пц

Рп ( г )

п=0 г .<Щ,

и вероятность прерывания Р{Х (г ) = ( п, ( г,..., гп) ,0) | X (г - А) = ( т, ( г,..., гт ) ,1), т > п} как вероятность того, что при изъятии ресурса будет потеряна хотя бы одна заявка

а

с

е=221

п=1 г >Щ

а + АР (Щ - г.) + пц

Рп(г )-

3. Случай дискретного ресурса

В связи с тем что получить аналитическое решение и даже метод вычисления характеристик в общем случае произвольного распределения ресурса весьма затруднительно, далее рассмотрим случай дискретного распределения ресурса. Для дискретного случая в предположении, что р (х) есть вероятность вновь пришедшей заявке занять ресурс объема х, аналогично вышеизложенному получаем СУР следующего вида:

( \

^0

( АР (Щ ) + в ) = ц^! (-)

+ а

22 Рк (r!,., гк ) + Р0

к=1 г .<Щ

Р0 ((Щ) + а) = (7) + ,

г=0

Яп ( Г1,., гп )( (- г-)+ пЦ + в ) = АЯп-1 ( г1гп-1 ) Р ( гп ) +

«0-2г

п+1

+ц2 2 Яп+1 (гг-гп)+

+а 2 2 Рк (rl,., гп, rn+1,., гк ) + Рп (Г1,., гп )

к=п+1 г.<Я1

п

Щ0-гп+1 <2г <Я0 7=1

И = 1,., С — 1,

Рп (г1 , . * * , гп )( (Щ - г-) + пЦ + а )= АРп-1 (г1 , . * *, гп-1) Р (гп ) +

щ -2о

п+1 j &

+ц2 2 Яп+1 (г1Т гп )+ РЯп (г1гп ) , I=1 г = 0

п = 1,.,С -1,

Яс (^ ■■ - ,гс)(СЦ + Р) = ^с-1 (rl,.,гс-1)Р(гс) + «Рс (rl,гс), рс (гх,..гс)(Сц+ а) = (гх,..,гс_х)р(гс) + Рдя (гх,..,гп), Для упрощения вычислений и ввиду большой размерности пространства состояний X введем марковский процесс 7(?) = (я (?), г (?),.$•(?)) с пространством состояний У = {(я,г,я): 0 <п < С,5 е {0,1},г .

к

Заметим, что здесь учитывается только суммарный объем занятого ресурса, а не объем, занятый каждой заявкой. В момент окончания обслуживания заявки, а также в момент изменения состояния системы с 5 = 1 на 5 = 0, приводящего к потере заявки, освобождаемый ресурс переразыгрывается согласно формуле условной вероятности. Недостатком такого упрощения является то, что при его применении можно наблюдать разницу между занимаемым и освобождаемым заявкой ресурсом, однако, как показано в [16], для надежной системы данное упрощение позволяет достаточно точно оценить поведение исходной системы.

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

Рп (г) = Р{п(0 = п,(Г^)) = г,5^) = 1}, дп (г) = р{п(*) = п, (Г (*)) = Г,5^) = 0}, и тогда СУР для процесса У (I) может быть записана в следующем виде:

+а

^0

до (0)(ХЕ(Я) + в) = ц(1) +

1 =о

2 2 р(1) 2 р () р; 1(0) + ро(0)

¿=1 1 = Яо+1 5=Я0+1 р (1)

Р0 (0)(ХЕ(Я) + а) = ц2Р1 (1) + вд0 (0),

1 =0

дп (г) (ХЕ(Я - г) + пц + в) = Х2Р( 1 )дп-1 (г -1) ■

1 =0

Я„-г

+^2 1Щ- '-+■ (г+1)+

(

+а

1=0 Р (г +1

1 -г

2 2 р(л 2 р(г>р(5>р(г-5> ++

¿=п+11 =Я0 +1 5 =Я0 г+1 р (1)

п = 1,..., С — 1,

г

рЛг К - г)+пц+а)=х2р(1) р»- (г -1)+

1=0

+ц2 р. +(г+1 )+Рд. ( г ),

Л

1=0 ¿=1

р ( г +1)

(п+1)

п = 1,.,С -1,

г

дс (г) (пц + в) = Х2р(1)дс-1 (г -1) + арс (г)'

1=0

рс (г)(пц+а)=х2р( 1 рс-1 (г -1)+вдп (г).

1=0

Здесь р (г)(п) «-кратная свертка, определяющая вероятность того, что п заявок займут ресурс

(1)

(2)

(3)

(4)

(5)

(6)

объема г.

4. Вычисление матрицы интенсивностей

Полученная в предыдущем разделе СУР (1)-(6) может также быть записана в матричном виде. Для этого представим пространство состояний в виде непересекающихся подмножеств

г

С

= где э*=иэ£» *<е{0,1}, и у.2 =

¡=0

= {(У, Л +1,1), ..., (У, Л ,1), У = 1, ..., С} . В этом случае матрица интенсивностей переходов процесса У ( ^) будет иметь блочный вид

где

Л00 Л01 Л02

Л = Л10 Л11 Л12

V Л20 Л21 Л22 у

Л00 =

Го00 и00 О 0

V 0 0 Ьос

о ^ о

ип

Б

0С у

Лп =

(о и

Б00 и 00

1^01 О01

Л01 = <1а§ (Р)

Л02 =( 0), Лю = <Иа§ (а) О 0 ^ О

и,

'ОС-1

ч О О Ьос Бос ,

Г и.

Л12 -

20 0

о 0 ^ о

0 0 0

Л20 = (% у),(г ,п)) ) =

1 -г

р (г)(п)

а-— 2 Р(5)Р(У " г "

-1)

р (У)

$=Е0-г+1

г = 0,..., Л ,п = 0,..., у = Л +1,., Л„ ■ = 0,., С,

Л21 =

Ь21 о о

о о ь

2С 0у

А -

, Л22 =

Г Ою ию Ьп Вп

Ьм

О

о о

V

и01 =(}р(г - У))у=о,.,Цо,

'=Л. Rо,

и11 =(Ар(г - У))}=Ra,...,R1, и21 =(АР(■ - У)) У=о,-Л,,

■ =Л0 +1,.,Л1,

. Л°-Тр(У) (1 - р(У))(п+1-г) '

2--—

иначе О 0 ^ О

'•• и^

^1С О1С у

У=о

р (г + У )

(п+1)

р ( У ) (1 - р (У ))(п+1-г

2

у=о

р (г + У )

(п+1)

У=0,...Ло ■=0,...,У,

У = Л0,.,Л1, ■=Л0,.,У,

0

Li =

ip 2

p ( j ) (1 - p ( j ))'

( и+1-

j=0

p ( r + j )

Л(и+1)

-■R> +1.....

j=0,...,R„

D0i = diag(-AF(R - r) + ри + p), Djj = diag(-AF(Rj - r) + ри + a), D2i = diag(-ÀF (R - r) + ри + a).

В матричном виде СУР (1)-(6) может быть решена известными, например итерационными, методами.

5. Численный пример

Рассмотрим пример численного анализа вероятностных характеристик СМО с дискретным ресурсом. Исходные данные для расчетов взяты из работы [11]. В примере под ресурсом понимается емкость (скорость передачи данных) радиоканала беспроводной сети. Базовая полоса пропускания оператора сети составляет C0 = 10 МГц, а арендуемая им LSA полоса составляет C - C = 5 МГц. Таким образом, пропускные способности основной и арендуемой полос можно рассчитать по формуле R = C ' V , где V - спектральная эффективность, составляющая для сетей LTE величину 4 бит/с/Гц [17]. Тогда объем надежной составляющей ресурса равен R = 10 МГц - 4 бит/с/Гц = 40 Мбит/c, а объем ненадежной составляющей равен R - R = 5 МГц - 4 бит/с/Гц = 20 Мбит/c . Следовательно, общий объем ресурса равен R = 60 Мбит/c . Для рассматриваемого сценария будем считать, что пользо-

[0, х < 2

ватели порождают трафик 2 Мбит/c.

, т.е. F ( х ) = |0

. Среднее время обслуживания р 1 = 15 с, изъ-

х > 2

ятие ресурса происходит в среднем каждые а 1 = 2 мин , а длительность периода изъятия составляет в среднем р"1 = 1 мин .

0,01 0,0001 0,000001 1E-08 1E-10 1E-12

вероятность эвакуации, P вероятность прерывания, Q

50

100

150

Рис. 2 Вероятности эвакуации и прерывания для СМО с ненадежным ресурсом

Характеристики системы при различных значениях нагрузки для приведенных данных показаны на рис. 2. Выбор логарифмического масштаба оси ординат объясняется тем, что при малой нагрузке вероятность эвакуации Р стремится к вероятности изъятия ресурса, а при росте нагрузки резко убывает.

1

0

Р

Заключение

В статье рассмотрена СМО с ненадежным ресурсом и предложены два алгоритма потери заявок. Подробно исследован алгоритм А1, для которого получена СУР в общем случае и проведен анализ упрощенного процесса для получения метода расчета основных вероятностных характеристик. Для дискретного случая и упрощенного процесса получена СУР в матричном виде, а также приведен численный пример для близких к реальным исходных данных.

В дальнейшем планируется исследовать модель с очередью, другие варианты функции распределения ресурсов, а также алгоритм выбора заявок А2, как более эффективный с точки зрения потерь заявок. Кроме того, интересно исследовать модель с более сложным входящим потоком, например MAP-потоком.

В заключение авторы приносят благодарность профессору Валерию Арсентьевичу Наумову и доценту Эдуарду Сергеевичу Сопину за ценные замечания и внимание к исследованиям.

ЛИТЕРАТУРА

1. Masek P., Mokrov E., Pyattaev A., Zeman K., Ponomarenko-Timofeev A., Samuylov A., Sopin E., Hosek J., Gudkova I., Andreev S.,

Novotny V., Koucheryavy Y., Samouylov K. Experimental Evaluation of Dynamic Licensed Shared Access Operation in Live 3GPP LTE System // Global Communications Conference (GLOBECOM). IEEE. 2016.

2. Palola M., Rautio T., Matinmikko M., Prokkola J., Mustonen M., Heikkila M., Kippola T., Yrjola S., Hartikainen V., Tudose L.,

Kivinen A., Paavola J., Okkonen J., Makelainen M., Hanninen T., Kokkinen H. Licensed Shared Access (LSA) trial demonstration using real LTE network // International Conference on Cognitive Radio Oriented Wireless Networks (CROWNCOM). 2014.

3. Palola M., Matinmikko M., Prokkola J., Heiska K. Live field trial of Licensed Shared Access (LSA) concept using LTE network in

2.3 GHz band // IEEE International Symposium on Dynamic Spectrum Access. 2014.

4. Mokrov E., Ponomarenko-Timofeev A., Gudkova I., Masek P., Hosek J., Andreev S., Koucheryavy Y., Gaidamaka Yu. Modeling

Transmit Power Reduction for a Typical Cell with Licensed Shared Access Capabilities // IEEE Transactions on Vehicular Technology. 2017.

5. Ponomarenko-Timofeev A., Pyattaev A., Andreev S., Koucheryavy Y., Mueck M., Karls I. Highly dynamic spectrum management

within licensed shared access regulatory framework // IEEE Communications Magazine. 2016. V. 54, No. 3. P. 100-109.

6. Yrjola S., Kokkinen H. Licensed Shared Access evolution enables early access to 5G spectrum and novel use cases // EAI Endorsed

Transactions on Wireless Spectrum. 2017. V. 3, No. 10.

7. Suliman I., Lehtomaki J. Queueing Analysis of Opportunistic Access in Cognitive Radios // Cognitive Radio and Advanced Spec-

trum Management (CogART). 2009.

8. Giorgetti A., Varrella M., Chiani M. Analysis and performance comparison of different cognitive radio algorithms // Cognitive Radio

and Advanced Spectrum Management (CogART). 2009.

9. Ivanova D., Karnauhova E., Markova E., Gudkova I. Analyzing of Licensed Shared Access Scheme Model with Service Bit Rate

Degradation in 3GPP Network // Information Technologies and Mathematical Modelling Queueing Theory and Applications. 2017. P. 231-242.

10. Samouylov K.E., Naumov V., Sopin E., Gudkova I., Shorgin S. Sojourn Time Analysis for Processor Sharing Loss System with Unreliable Server // Analytical & Stochastic Modelling Techniques & Applications ASMTA, 2016. P. 284-297.

11. Острикова Д.Ю. Анализ вероятностно-временных характеристик схем доступа с прерыванием обслуживания в телекоммуникационных беспроводных сетях : дис. ... канд. физ.-мат. наук / Рос. ун-т дружбы народов. М., 2015. 124 с.

12. Тихоненко О.М., Климович К.Г. Анализ систем обслуживания требований случайной длины при ограниченном суммарном объеме // Проблемы передачи информации. 2001. Т. 37, вып. 1. С. 70-79.

13. Бородин Н.Н., Малинковский Ю.В. Характеристики суммарного объема требований в СеМО Джексона // Вестник Томского государственного университета. Управление, вычислительная техника и информатика. 2010. № 1 (10). C. 253-285.

14. Лисовская Е.Ю., Моисеева С.П. Асимптотический анализ немарковской бесконечнолинейной системы обслуживания требований случайного объема с входящим рекуррентным потоком // Вестник Томского государственного университета. Управление, вычислительная техника и информатика. 2017. № 39. С. 30-38.

15. Назаров А.А., Бронер В.И. Система управления запасами с гиперэкспоненциальным распределением объемов потребления ресурсов // Вестник Томского государственного университета. Управление, вычислительная техника и информатика. 2016. № 1 (34). С. 43-49.

16. Naumov V., Samouylov K., Yarkina N., Sopin E., Andreev S., Samuylov A. LTE Performance Analysis Using Queuing Systems with Finite Resources and Random Requirements // 7th International Congress on Ultra Modern Telecommunications and Control Systems and Workshops (ICUMT). 2015.

17. 3GPP TS 36.300: Evolved Universal Terrestrial Radio Access (E-UTRA) and Evolved Universal Terrestrial Radio Access Network (E-UTRAN); Overall description; Stage 2: Release 12. 2015. 177 p.

Поступила в редакцию 17 марта 2018 г.

Mokrov E.V., Samouylov K.E. (2018) LSA FRAMEWORK IN LTE WIRELESS NETWORK ENVIRONMENT AS A QUEUEING SYSTEM WITH UNRELIABLE RESOURCE. Vestnik Tomskogo gosudarstvennogo universiteta. Upravlenie vychislitelnaja tehnika i informatika [Tomsk State University Journal of Control and Computer Science]. 44. pp. 62-71

DOI: 10.17223/19988605/44/7

In this paper, we consider an LTE cell with LSA framework as a queueing system with unreliable resources. Here the resources can be assumed to represent e.g. frequency spectrum, since the LSA framework implies that part of the spectrum is leased from the third party and can be seized at a random time for a random duration.

A system with C servers and Rs resources is studied. The system can be found in one of s e {0,1} states. The state transitions are

distributed according to the exponential law with parameters an P , notably r0 < r1. The customers arrive according to the Poisson process with parameter X , and are served according to the exponential law with parameter p. When a new customer arrives it requires a server and a random resource volume with CDF F(x) . A customer is considered lost if there is either no enough available servers in the system or there are no enough available resources to accommodate it. Thus in the considered system a customer can be lost either upon arrival or while it is being served upon the system state change.

First consider a general case of continuous resources. Here we introduce a Markov process X(t) = (n(t),(r1 (t),...,rn (t)),s(t))

X = ),i) : 0 <n<C,se {o,l},r. = ^ j

with state space |vM1>--->'"./>v ■ u - n - e tu>1/>r- / ' J ■ However, since, unfortunately, we failed to obtain a closed

form solution for this case further we consider a discrete case. But, as it is shown in the paper, applying the same process for the discrete case means that the state space is growing exponentially against the maximum number of customers in the system. And can be obtained

n \

-I i №+Rs)

This equation considers the scenario when a new customer can require the resources in the range from

i+

'—> n!

j=1

0 to Rs. As one can note from the equation the state space grows extremely rapidly with the number of servers and volume of available

resources, which highly affects the matrix techniques performance.

Taking all this into account, we propose an approximation, where the resources allocated to the customer are reallocated when the customer leaves the system. Thus we only have to trace the volume of the allocated resources in total and not the resources allocated to each customer separately, thus greatly reducing computation difficulty. In this case the process can be obtained as

Y(t) = (n (t), r(t),s (t)) . Since we consider the Poisson arrival flow and exponential service time this approximation gives us a reasonable approximation of the initial process.

For the considered system evacuation and interruption probabilities are obtained as the main characteristics. Evacuation probability is a probability that when the frequency is being seized, no customer was interrupted while being served are the interruption probability is a probability that when the frequency is being seized at least one customer was interrupted while being served. As one can note, the sum of these probabilities is a constant and equals to the probability of a frequency being seized. In the last section we give a numerical example.

Keywords: queuing system; Markov process; unreliable resources; cumulative distribution function.

MOKROVEvgenii Vladimirovich (Assistant Professor, Peoples' Friendship University of Russia, Moscow, Russian Federation). E-mail: melkor77@yandex.ru

SAMOUYLOV Konstantin Evgenievich (Doctor of Technical Sciences, Professor, Head of Department of Applied Informatics and Probability Theory, Peoples' Friendship University of Russia, Moscow, Russian Federation). E-mail: samuylov_ke@pfur.ru

REFERENCES

1. Masek, P., Mokrov, E., Pyattaev, A., Zeman, K., Ponomarenko-Timofeev, A., Samuylov, A., Sopin, E., Hosek, J., Gudkova, I.,

Andreev, S., Novotny, V., Koucheryavy, Y., & Samouylov, K. (2016) Experimental Evaluation of Dynamic Licensed Shared Access Operation in Live 3GPP LTE System. Global Communications Conference (GLOBECOM), IEEE. DOI: 10.1109/GLOCOM.2016.7841826

2. Palola, M., Rautio, T., Matinmikko, M., Prokkola, J., Mustonen, M., Heikkila, M., Kippola, T., Yrjola, S., Hartikainen, V., Tudose, L.,

Kivinen, A., Paavola, J., Okkonen, J., Makelainen, M., Hanninen, T. & Kokkinen, H. (2014) Licensed Shared Access (LSA) trial demonstration using real LTE network. International Conference on Cognitive Radio Oriented Wireless Networks (CROWNCOM). DOI: 10.4108/icst. crowncom.2014.255724

3. Palola, M., Matinmikko, M., Prokkola, J. & Heiska, K. (2014) Live field trial of Licensed Shared Access (LSA) concept using LTE

network in 2.3 GHz band. IEEE International Symposium on Dynamic Spectrum Access. DOI: 10.1109/DySPAN.2014.6817778

4. Mokrov, E., Ponomarenko-Timofeev, A., Gudkova, I., Masek, P., Hosek, J., Andreev, S., Koucheryavy, Y. & Gaidamaka, Yu.

(2017) Modeling Transmit Power Reduction for a Typical Cell with Licensed Shared Access Capabilities. IEEE Transactions on Vehicular Technology. pp. 1-1. DOI: 10.1109/TVT.2018.2799141

5. Ponomarenko-Timofeev, A., Pyattaev, A., Andreev, S., Koucheryavy ,Y., Mueck, M. & Karls, I. (2016) Highly dynamic spectrum

management within licensed shared access regulatory framework. IEEE Communications Magazine. 54(3). pp.100-109. DOI: 10.1109/MCOM.2016.7432155

6. Yrjola, S. & Kokkinen, H. (2017) Licensed Shared Access evolution enables early access to 5G spectrum and novel use cases. EAI

Endorsed Transactions on Wireless Spectrum. 3(10). DOI: 10.4108/eai.12-12-2017.153463

7. Suliman, I. & Lehtomaki, J. (2009) Queueing Analysis of Opportunistic Access in Cognitive Radios. Cognitive Radio and Advanced

Spectrum Management (CogART). pp. 153 - 157. DOI: 10.1109/COGART.2009.5167252

8. Giorgetti, A., Varrella, M. & Chiani, M. (2009) Analysis and performance comparison of different cognitive radio algorithms.

Cognitive Radio and Advanced Spectrum Management (CogART). DOI: 10.1109/COGART.2009.5167247

9. Ivanova, D., Karnauhova, E., Markova, E. & Gudkova, I. (2017) Analyzing of Licensed Shared Access Scheme Model with Service

Bit Rate Degradation in 3GPP Network. Information Technologies and Mathematical Modelling. Queueing Theory and Applications. pp. 231-242. DOI: 10.1007/978-3-319-68069-9_19

10. Samouylov, K.E., Naumov, V., Sopin, E., Gudkova, I. & Shorgin, S. (2016). Sojourn Time Analysis for Processor Sharing Loss System with Unreliable Server. Analytical & Stochastic Modelling Techniques & Applications (ASMTA). pp. 284-297. DOI 10.1007/978-3-319-43 904-4_20

11. Ostrikova, D.Y. (2015) Analiz veroyatnostno-vremennykh kharakteristik skhem dostupa spreryvaniyem obsluzhivaniya v telekom-munikatsionnykh besprovodnykh setyakh [Analysis of probability-time characteristics of access schemes with service interruptions in wireless telecommunication networks]. Physics and Mathematics Cand. Diss. Moscow.

12. Tikhonenko, O.M. & Klimovich, K.G. (2001) Analysis of Queuing Systems for Random-Length Arrivals with Limited Cumulative Volume. Problems of Information Transmission. 37(1). pp. 70-79. (In Russian). DOI: 10.1023/A:1010451827648.

13. Borodin, N.N. & Malinovsky, Y.V. (2010) The characteristics of cumulative volume of demands in Jackson networks. Vestnik Tomskogo gosudarstvennogo universiteta. Upravleniye, vychislitel'naya tekhnika i informatika - Tomsk State University Journal of Control and Computer Science. 1(10). pp. 253-285. (In Russian).

14. Lisovskaya, E.Y. & Moiseeva, S.P. (2017) Asymptotical analysis of a non-Markovian queueing system with renewal input process and random capacity of customers. Vestnik Tomskogo gosudarstvennogo universiteta. Upravleniye, vychislitel'naya tekhnika i informatika - Tomsk State University Journal of Control and Computer Science. 39. pp. 30-38. (In Russian). DOI: 10.17223/19988605/39/5.

15. Nazarov, A.A. & Broner, V.I. (2016) Inventory model with hyperexponential distribution of demand's batch size. Vestnik Tomskogo gosudarstvennogo universiteta. Upravleniye, vychislitel'naya tekhnika i informatika - Tomsk State University Journal of Control and Computer Science. 1(34). pp. 43-49. (In Russian). DOI: 10.17223/19988605/34/5

16. Naumov, V., Samouylov, K., Yarkina, N., Sopin, E., Andreev, S. & Samuylov, A. (2015) LTE Performance Analysis Using Queuing Systems with Finite Resources and Random Requirements. 7th International Congress on Ultra Modern Telecommunications and Control Systems and Workshops (ICUMT). DOI: 10.1109/ICUMT.2015.7382412.

17. 3GPP (2015) TS 36.300: Evolved Universal Terrestrial Radio Access (E-UTRA) and Evolved Universal Terrestrial Radio Access Network (E-UTRAN). Overall description; Stage 2: Release 12.