Научная статья на тему 'Анализ модели схемы доступа к радиоресурсам беспроводной сети с прерыванием обслуживания пользователей'

Анализ модели схемы доступа к радиоресурсам беспроводной сети с прерыванием обслуживания пользователей Текст научной статьи по специальности «Электротехника, электронная техника, информационные технологии»

CC BY
107
23
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
СМО / СЛУЧАЙНАЯ СРЕДА / ПРЕРЫВАНИЕ ОБСЛУЖИВАНИЯ / РЕКУРРЕНТНЫЙ АЛГОРИТМ / ПРИБЛИЖЕННЫЙ МЕТОД / QUEUING SYSTEM / RANDOM ENVIRONMENT / SERVICE INTERRUPTION / RECURSIVE ALGORITHM / APPROXIMATE METHOD

Аннотация научной статьи по электротехнике, электронной технике, информационным технологиям, автор научной работы — Аду Кпангни Ивь, Маркова Екатерина Викторовна, Гудкова Ирина Андреевна

В связи с развитием телекоммуникационных беспроводных сетей последующих поколений, в частности сетей на базе технологии LTE, повышается спрос на высокоскоростные услуги, а также увеличивается число их пользователей. Это приводит к экспоненциальному росту объема передаваемого в мобильных сетях трафика. Согласно последним данным ведущих участников телекоммуникационного рынка, таких как, например, компания Cisco Systems, объем трафика, генерируемого в современных мобильных сетях, увеличится на 46 процентов в период с 2016 по 2021 года. В то же время повышаются требования к уровню качества обслуживания пользователей. Стремительный рост трафика в условиях предоставления пользователям мобильных сетей широкого спектра высококачественных мультимедийных услуг приводит к нехватке доступного частотного диапазона. Данная проблема становится главной в развитии современных технологий беспроводной связи. В связи с этим операторы сотовых сетей предлагают различные методы ее решения, например, применение системы LAA (англ. Licensed Assisted-Access) или LSA (англ. Licensed Shared Access). Эти системы позволяют более эффективно использовать имеющиеся частотные ресурсы. Помимо этого, в рамках реализации различных технологий для обеспечения требуемого уровня качества обслуживания пользователей могут быть применены такие механизмы управления радиоресурсами, как снижение мощности передачи данных и прерывание обслуживания пользователей. Модели реализующие такие механизмы могут быть описаны в виде систем массового обслуживания (СМО) с ненадежными приборами. В работе построена модель соты беспроводной сети, представленная в виде СМО, приборы которой находятся в случайной среде. Состояние случайной среды может изменяться, что с технической точки зрения эквивалентно уменьшению или увеличению мощности передачи сигнала. Снижение мощности может привести к прерыванию обслуживания пользователей. Предложены точный и приближенный методы для расчета стационарного распределения вероятностей состояний модели.

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по электротехнике, электронной технике, информационным технологиям , автор научной работы — Аду Кпангни Ивь, Маркова Екатерина Викторовна, Гудкова Ирина Андреевна

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

TO PERFORMANCE ANALYSIS OF A WIRELESS NETWORK WITH DYNAMIC POWER CONTROL POLICY

In connection with the development of next-generation wireless networks based on LTE technology, demands for high-speed services are increasing leading to an exponential increase in the volume of traffic generated in mobile networks. According to the latest data from Cisco Systems, this volume will increase by 46 percent by 2021. At the same time, requirements for the level of quality of service (QoS) are increasing. The rapid growth of traffic in the conditions of providing users with a wide range of high-quality multimedia services leads to a lack of available frequency resources. This problem becomes the main one in the development of modern wireless technologies. In this regard, cellular network operators offer various methods for its solution, for example, using the LAA (Licensed Assisted-Access) or LSA (Licensed Shared Access) systems allowing more efficient use of available frequency resources. In addition, in the framework of the implementation of various technologies to ensure the required QoS, radio resource management mechanisms (RRM) such as reducing data transmission power and interrupting users’ service can be applied. Models implementing such mechanisms can be described as queuing systems (QS) with unreliable servers. In the paper, a model of a wireless network cell is constructed, presented in the form of a QS, servers of which are in a random environment (RE). RE state can change, which is technically equivalent to a decrease or increase in the transmission power. Reducing power may interrupt users’ service. Exact and approximate methods are proposed for calculating the stationary probability distribution of the model.

Текст научной работы на тему «Анализ модели схемы доступа к радиоресурсам беспроводной сети с прерыванием обслуживания пользователей»

ТЕОРЕТИЧЕСКИЕ ВОПРОСЫ ИНФОРМАТИКИ, ПРИКЛАДНОЙ МАТЕМАТИКИ, КОМПЬЮТЕРНЫХ НАУК И КОГНИТИВНО-ИНФОРМАЦИОННЫХ ТЕХНОЛОГИЙ

УДК 519.872.1

DOI: 10.25559/SITITO.15.201903.563-571

Анализ модели схемы доступа к радиоресурсам беспроводной сети с прерыванием обслуживания пользователей

К. И. Аду1, Е. В. Маркова1, И. А. Гудкова1,2

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

* markova_ev@rudn.ru

2 Федеральный исследовательский центр «Информатика и управление» Российской академии наук, г. Москва, Россия

119333, Россия, г. Москва, ул. Вавилова, д. 44-2 Аннотация

В связи с развитием телекоммуникационных беспроводных сетей последующих поколений, в частности сетей на базе технологии LTE, повышается спрос на высокоскоростные услуги, а также увеличивается число их пользователей. Это приводит к экспоненциальному росту объема передаваемого в мобильных сетях трафика. Согласно последним данным ведущих участников телекоммуникационного рынка, таких как, например, компания Cisco Systems, объем трафика, генерируемого в современных мобильных сетях, увеличится на 46 процентов в период с 2016 по 2021 года. В то же время повышаются требования к уровню качества обслуживания пользователей. Стремительный рост трафика в условиях предоставления пользователям мобильных сетей широкого спектра высококачественных мультимедийных услуг приводит к нехватке доступного частотного диапазона. Данная проблема становится главной в развитии современных технологий беспроводной связи. В связи с этим операторы сотовых сетей предлагают различные методы ее решения, например, применение системы LAA (англ. Licensed Assisted-Access) или LSA (англ. Licensed Shared Access). Эти системы позволяют более эффективно использовать имеющиеся частотные ресурсы. Помимо этого, в рамках реализации различных технологий для обеспечения требуемого уровня качества обслуживания пользователей могут быть применены такие механизмы управления радиоресурсами, как снижение мощности передачи данных и прерывание обслуживания пользователей. Модели реализующие такие механизмы могут быть описаны в виде систем массового обслуживания (СМО) с ненадежными приборами. В работе построена модель соты беспроводной сети, представленная в виде СМО, приборы которой находятся в случайной среде. Состояние случайной среды может изменяться, что с технической точки зрения эквивалентно уменьшению или увеличению мощности передачи сигнала. Снижение мощности может привести к прерыванию обслуживания пользователей. Предложены точный и приближенный методы для расчета стационарного распределения вероятностей состояний модели.

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

Финансирование: публикация подготовлена при финансовой поддержке Российского фонда фундаментальных исследований в рамках научного проекта № 18-37-00231 мол_а «Разработка моделей и алгоритмов для анализа схем совместного использования радиочастот в беспроводных мультисервисных сетях с произвольным расположением устройств в пространстве».

Контент доступен под лицензией Creative Commons Attribution 4.0 License. The content is available under Creative Commons Attribution 4.0 License.

Vol. 15, No. 3. 2019 ISSN 2411-1473 sitito.cs.msu.ru

Modern Information Technologies and IT-Education

ТЕОРЕТИЧЕСКИЕ ВОПРОСЫ ИНФОРМАТИКИ, ПРИКЛАДНОЙ МАТЕМАТИКИ, КОМПЬЮТЕРНЫХ НАУК И КОГНИТИВНО-ИНФОРМАЦИОННЫХ ТЕХНОЛОГИЙ

Для цитирования: Аду К. И., Маркова Е. В., Гудкова И. А. Анализ модели схемы доступа к радиоресурсам беспроводной сети с прерыванием обслуживания пользователей // Современные информационные технологии и ИТ-образование. 2019. Т. 15, № 3. С. 563-571. DOI: 10.25559/ SITITO.15.201903.563-571

Современные информационные технологии и ИТ-образование

Том 15, № 3. 2019 ^ 2411-1473 sitito.cs.msu.ru

THEORETICAL QUESTIONS OF COMPUTER SCIENCE, COMPUTATIONAL MATHEMATICS, COMPUTER SCIENCE AND COGNITIVE INFORMATION TECHNOLOGIES

To Performance Analysis of a Wireless Network with Dynamic Power Control Policy

K. Y. Adou1, E. V. Markova1, I. A. Gudkova12

1 Peoples' Friendship University of Russia, Moscow, Russia 6 Miklukho-Maklaya Str., Moscow 117198, Russia

* markova_ev@rudn.ru

2 Federal Research Center "Computer Science and Control" of the Russian Academy of Sciences, Moscow, Russia

44-2 Vavilov Str., Moscow 119333, Russia

In connection with the development of next-generation wireless networks based on LTE technology, demands for high-speed services are increasing leading to an exponential increase in the volume of traffic generated in mobile networks. According to the latest data from Cisco Systems, this volume will increase by 46 percent by 2021. At the same time, requirements for the level of quality of service (QoS) are increasing. The rapid growth of traffic in the conditions of providing users with a wide range of high-quality multimedia services leads to a lack of available frequency resources. This problem becomes the main one in the development of modern wireless technologies. In this regard, cellular network operators offer various methods for its solution, for example, using the LAA (Licensed Assisted-Access) or LSA (Licensed Shared Access) systems allowing more efficient use of available frequency resources. In addition, in the framework of the implementation of various technologies to ensure the required QoS, radio resource management mechanisms (RRM) such as reducing data transmission power and interrupting users' service can be applied. Models implementing such mechanisms can be described as queuing systems (QS) with unreliable servers. In the paper, a model of a wireless network cell is constructed, presented in the form of a QS, servers of which are in a random environment (RE). RE state can change, which is technically equivalent to a decrease or increase in the transmission power. Reducing power may interrupt users' service. Exact and approximate methods are proposed for calculating the stationary probability distribution of the model.

Keywords: Queuing system, random environment, service interruption, recursive algorithm, approximate method.

Funding: The publication was prepared with the financial support of the Russian Foundation for Basic Research as a part of scientific project No. 18-37-00231 MO^_a «Development of Models and Algorithms for the Analysis of Radio Frequency Sharing Schemes in Wireless Multiservice Networks with an Arbitrary Location of Devices».

For citation: Adou K.Y., Markova E.V., Gudkova I.A. To Performance Analysis of a Wireless Network with Dynamic Power Control Policy. Sovremennye informacionnye tehnologii i IT-obrazovanie = Modern Information Technologies and IT-Education. 2019; 15(3):563-571. DOI: 10.25559/SITI-T0.15.201903.563-571

Abstract

Modern Information Technologies and IT-Education

ТЕОРЕТИЧЕСКИЕ ВОПРОСЫ ИНФОРМАТИКИ, ПРИКЛАДНОЙ МАТЕМАТИКИ, К. И. Аду, Е. В. Марк0ва,

КОМПЬЮТЕРНЫХ НАУК И КОГНИТИВНО-ИНФОРМАЦИОННЫХ ТЕХНОЛОГИЙ ' И. А. Гудкова

Введение

В современных мобильных сетях возникает проблема нехватки радиоресурсов [1,2], необходимых для предоставления пользователям услуг с требуемым уровнем качества обслуживания (англ. Quality of Service, QoS) [3, 4]. На данный момент существуют несколько способов решения этой технической проблемы, например, мультиплексирование с ортогональным частотным разделением каналов (англ. Orthogonal frequency-division multiplexing, OFDM) [5, 6, 7], многоантенная передача данных (англ. Multiple Input Multiple Output, MIMO) [8, 9, 10], мультимедийная, широковещательная и многоадресная передача (англ. Multimedia Broadcast Multicast Service, MBMS) [11, 12, 13, 14], совместное использование лицензированного спектра радиочастот при помощи системы LSA (англ. Licensed Shared Access) [15, 16, 17] или нелицензированного спектра частот при помощи системы LAA (англ. Licensed AssistedAccess) [18, 19]. При этом для более эффективного управления сетевыми ресурсами в рамках каждого из этих решений могут быть использованы различные механизмы доступа, например, снижение мощности передачи сигнала [20, 21] и прерывание обслуживания пользователей [22, 23]. Модели с такими механизмами могут быть описаны в виде систем массового обслуживания (СМО) с ненадежными приборами.

Одним из типов таких СМО являются системы, работающие в случайной среде [24, 25, 26]. Случайная среда может значительно повлиять на работу телекоммуникационных систем. В частности, она оказывает существенное влияние на такие показатели эффективности системы, как средняя скорость передачи данных, среднее число обслуживаемых запросов и т.д.

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

Построение математической модели

Рассмотрим некоторую случайную среду, оказывающую влияние на показатели качества обслуживания пользователей в соте беспроводной сети радиуса R . Введем упрощающее предположение: все пользователи находятся на одинаковом расстоянии d от базовой станции (БС). Пусть в систему поступают запросы пользователей на предоставление услуг, генерирующих потоковый трафик с гарантированной скоростью передачи r0.

Определим состояния случайной среды множеством {0,...,S}. Управление доступом к радиоресурсам системы реализовано таким образом, что под действием внешней среды

происходит снижение мощности передачи сигнала со значения р1 до значения ра-1, при этом среда с интенсивностью аз переходит из состояния я в состояние я-1, а с интенсивностью происходит обратный переход из состояния я-1 в состояние я и восстановление мощности до уровня р1.

Таким образом под влиянием случайной среды мощность передачи сигнала может меняться, а это, в свою очередь, отражается на максимально-достижимой скорости передачи данных г1 и на максимальном числе N одновременно обслуживаемых запросов. Это число может быть рассчитано как отношение достижимой скорости г1, определяемой по формуле Шеннона, к гарантированной скорости г0 следующим образом: ,

N,

ffllnl —— +1

d KK

= 0,.

(1)

где о - ширина полосы частот канала, G - коэффициент затухания сигнала, к - степень затухания сигнала, K0 - мощность шума. Из формулы (1) следует, что Ns > Ns_j >... > N0. При переходе среды из состояния я в состояние я -1 происходит снижение числа обслуживаемых запросов пользователей с Ns до Ns-1 за счет прерывания обслуживания Ns - Ns-1 запросов. В таблице 1 представлены основные обозначения, используемые в работе.

Таблица 1. Основные обозначения Table 1. Principal Symbols

Обозна-

Описание

R

радиус соты, м

® ширина полосы частот канала, Гц

интенсивность перехода в состояние я -1, я ф 0, с в- интенсивность перехода в состояние я +1,я ф $, с N число обслуживаемых запросов в зависимости от состо-_яния среды У = Г), , .V

Ps

максимальная мощность сигнала, Гц

p минимальная мощность сигнала, Гц

текущая мощность сигнала s = 0,...,S, Гц

d расстояние до БС, м

гарантированная скорость передачи данных, Мбит/с достижимая скорость передачи данных я = 0,..., 3,

Мбит/г

X интенсивность входящего потока запросов, с интенсивность обслуживания одного запроса, с р = Я/^ интенсивность предложенной нагрузки К мощность шума, Гц

О коэффициент затухания сигнала, константа к степень затухания сигнала

Поведение системы описывает марковский процесс (МП) X(/) = (N(t),S(t),t > 0) над пространством состояний:

X = {(п,s): 0 < п < Ns,s = 0,...,^, (2)

где п - число обслуживаемых запросов пользователей, я - состояние случайной среды.

Предположим, что входящий поток запросов является пу-ассоновским с параметром X, а время обслуживания запросов экспоненциальным со средним значением ^т1. Тогда диаграмма интенсивностей переходов, имеет вид, представленный на рис. 1.

r

r

Современные информационные технологии и ИТ-образование

Том 15, № 3. 2019 ISSN 2411-1473 sitito.cs.msu.ru

K. Y. Adou, E. V. Markova, I. A. Gudkova

THEORETICAL QUESTIONS OF COMPUTER SCIENCE, COMPUTATIONAL MATHEMATICS, COMPUTER SCIENCE AND COGNITIVE INFORMATION TECHNOLOGIES

N. i.s'-l s * 0; n > N i

n-1. 5 s < S; п> 0

-«s+l

n.s+1 s*S-, n < Ns

л .v-1 I.-D;u< JV, ,

!L_S

(n, s) с X

-A—1

-(П + 11/.Г

и a

/¡+1. .V

s s S-, n < N,

Р и с. 1. Диаграмма переходов для СМО в случайной среде F i g. 1. Transition diagram for QS in a random environment

Как видно из диаграммы, СУР для стационарного распределения p(n,s) имеет вид:

pin, s)■

(3)

'XI {s < S, n < Ns} + npI {s < S, n > 0} + +aj {s Ф 0, n < Ns } + PsI {s Ф S -1, n < Ns = p(n -1, s) -XI {n > 0} + p(n +1, s) ■ (n + 1)^I {n < Ns} + +p(n, s +1) ■ as+jI {s < S} + p(n, s -1) ■ в_1I {s > 0}, (n, s) e X, где I - функция индикатор.

Вследствие реализации механизма прерывания обслуживаемых запросов, МП X(t) не является обратимым, поэтому не следует ожидать решения СУР (3) в замкнутом аналитическом виде. В связи с этим стационарное распределение p = (p(n,s))( } х может быть получено путем численного решения СУР (3), для чего необходимо вычислить инфинитези-мальный оператор A МП X(t). Нетрудно убедиться, что элементы a((n,s)(n',s')) матрицы A вычисляются по формуле:

as, если n' = n, s' = s -1, n < Ns 1, s Ф 0,

или n'= Ns 1,s'= s — 1, n > Ns 1,s Ф 0,

Ps, если n'= n, s' = s +1, n < Ns, s Ф S,

X, если n' = n +1, s' = s, n < Ns, s < S,

пц, если n' = n — 1,s' = s,n > 0,s < S,

Ф, если n' = n, s' = s,

0, в противном случае

где

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

p = -(aj{s # 0} + PJ{s Ф S,n < Ns } + Xl{s < S, n < Ns} + n^l{n > 0}).

Из постановки задачи следует, что анализу подлежат следующие вероятностные характеристики модели - вероятность блокировки запросов пользователей B , среднее число обслуживаемых в системе запросов N, средняя скорость передачи данных r, вероятность прерывания обслуживания запросов П , определяемые соответственно по формулам (5) - (8).

в = É>( N, s); (5)

(4)

N=ZZ np(n'

s=0 n=1

r = Z ÍP(", s)/1l-Z P(0, s)

п = X

а,

n - Ns-l

=N-l +1аs + ßs + nß + X n

аs Ns - Ns-i

p(n, s) +

Ns -l

+ X

+ ßs + Nsß Ns

а n - Ns -i

P( Ns, s)

(6)

(7)

(8)

;-+1а£ + nß + X n а s Ns - Ns -i

p(n, S ) +

аs + Ns ß Ns

p(Ns, S ).

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

Модель с тремя состояниями случайной среды

На практике представляет интерес случай модели с тремя состояниями случайной среды, т.е. S = 2. Такая модель, как показано в [16], применима, например, к анализу системы LSA.

Рассматриваемая система описывается МП Y(t) = (N(t),S(t), t > 0) над пространством состояний:

Y = {(n,s): 0 < n < Ns,s = 0,1,2}, (9)

В данном случае для расчета стационарного распределения p(n,s) можно получить рекуррентный алгоритм. Для этого воспользуемся СУР (3), и введем вспомогательные функции, которые вычисляются по формулам:

9(0,0) = 1, ^(0,1) = х, q(0,2) = y, (10)

q(n,s) = Yns + Sns ■ x + 0ns ■ y. (11)

Заметим, что стационарное распределение вероятностей p(n,s) определятся по формуле:

p(n,s) = q(n,ü)j^(i )6Yqd, j),{n,s) e Y. (12)

Из формул (10) и (11), несложно получить рекуррентные соотношения для коэффициентов yn¡, 5n¡ и Qs, необходимых для расчета вспомогательных функций q(n, s):

Y 00 = 1,500 = °Ао = 0;

Yoi = 0Ai = 1,^0, = 0;

y 02 = 0,sm = 0,002 = 1;

Г10 = (^+Po)l нАо =-aJ HAo = 0; Y11 =-Po/^,5ii = (A+«i +Р1У vAi =-аг!

Y,2 = =-PJ HAI =(+«2

Yno = [(Г10 +(n - !)) ■ Yn-1,0 + Sio ■ Yn-1,1 - P ■ Yn-2,0 ]/n,n = 'í,...,No -1 S„0 = [(/10 +(n - !)) A-1,0 + S10 ■ 5„-1,1 - P ■ Sn-2,0 ]/" > n = ^..v - 1 8„0 = [(10 + (n-1))■ ^-,,0 + 5,0 ■ в„_1Л - p ■ en_2fi]/n,n = 2,..., Wo -1; r,1 =\_(p+(n-S12 )-Yn-1,1 -P-Yn-2,1 +011 -r,-1,2 +Yl1 ■ Yn-1,0 ]/» > n = 2,...>

Snl = \_(p+(n-S12 )'Sn-l,l -P-S„-2,1 +0Ц ■ S„-1,2 +Г11 -5,-1,0 ]/« > n = 2,..., N0, 0,1 =[(p + (n-1)-sio -S12)9,-1,1 -P■ 0,-2,1 +0ii ■ -i,2 +Yn ■ 0,-i,o]/я,n = 2,...,N„;

Y n2 =[(p+(n -1)-ви ) /„-1,2 -P- Yn-2,2 +5,2 ' /„-1,1 ]/« = ..., W„ 5„2 =[(p + (n-i)-0u)■ 5„-i,2 -P■ sn_2a + 5,2 ■ 8,-,,,]/«,n = 2,...,W„ e„2 =[(p+{" -1) - y e„-,,2 - P ■ в,-2,2 + 5,2 ■ en-,,, ]/«, n = 2,..., N

Ynl =[(p + (» -1) - ¿10 - 5,2) • 7,-1,1 - P ■ 7,-2.1 +в11 ■ Yn-1,2 ]/»>« = +-1>

5,1 =[(p + (» -!)- 5i0 - S12) ■ ¿,-1,1 - p ■ 8,-21 +en ■ Sn-12 ]/ n, n = Ж0 +1,..., -1, 0,1 =[(p + (n -1) - ¿10 -S12) • 0,-1,1 - p ■ 0,-2,1 +011 • 0,-1,2 ]/n, n = Ж0 +1,..., Ж, -1;

Modern Information Technologies and IT-Education

N-1

+

568

ТЕОРЕТИЧЕСКИЕ ВОПРОСЫ ИНФОРМАТИКИ, ПРИКЛАДНОЙ МАТЕМАТИКИ, К. И. Аду, Е. В. Марк0ва>

КОМПЬЮТЕРНЫХ НАУК И КОГНИТИВНО-ИНФОРМАЦИОННЫХ ТЕХНОЛОГИЙ ' И. А. Гудкова

Уп2 =[{р + (п - 1) - 01 ) • У„-1,2 - Р ■ Уп-2,2 ]/«, П = N + 1,..., N - 1,

8Л2 = \^(р + (п - 1)-0п )-8л-11 - р-8Л_2,2]/л,л = N +1,...,N2 -1,

вп 2 = \[р + (п -1) - в11) • вп-1,2 - р • вп-2,2 ]/п, п = N1 +1,..., N. -1.

Используя СУР (3), определим значения функций х и у, необходимых для расчета q(n, s):

х = (Р ■ Q - R ■ L)/(Р • Т - и ■ L),

(13)

у = ( • Т - и • Q)/(Р • Т - и • L),

(14)

где:

р=(а+в ) • ^ д-а • Е !=N в-2+я •-1,1;

2 = Уы2-1,2 - («2 + ЯгР)- Уы2,2~;

к = а2 ■ Е "=м у "2 + У* -1,1- («1 + ) • Ум, ,1;

I = (а + ^иУв^л - я •-1,2;

Т = ( + - Я • -1,2 ;

^ = («1 + N ^ + в )-8м„1 -а2 • Е П= N 8п 2 +х'5*1 -1,1.

Рассчитав стационарное распределение вероятностей с помощью рекуррентного алгоритма, вычислим основные показатели эффективности системы по формулам (5) - (8). Для сокращения сложности проводимых расчетов в следующем разделе рассмотрим приближенный метод анализа модели с тремя состояниями случайной среды.

Приближенный метод анализа вероятностных характеристик

Для получения решения в мультипликативном виде упростим модель. Предположим, что снижение мощности под влиянием внешней среды происходит только в том случае, когда число обслуживаемых запросов п, соответствующее состоянию среды 5 становится равным максимальному числу запросов, обслуживаемых в состоянии 5 -1, т.е. когда п = Ns-1 . Таким образом, в данном случае не происходит прерывания обслуживания. Рассмотрим такое упрощение на примере частного случая модели, представленного в предыдущем разделе.

Диаграмма интенсивностей переходов для модели без прерывания будет иметь вид, представленный на рис. 2.

I 0.2 I ' 'I 1.2 I ' 'I 2,2 -

I-1—5 I-—I I I 3/1 л

1 I л- /; I к- /■ I I /; I к-

1-1(Л/0+1)д ßl II"2 ßi 1 ■ 1 ,V1'' II

I I [«>1№ Nin 1-1

Р и с. 2. Диаграмма интенсивностей переходов для модели без прерывания F i g. 2. Transition intensities diagram for the model without interruption

С учетом введенного упрощения СУР представима следующим образом:

(Ä + ß0 )• />(0,0) = и- p(l,0) + а • p(0,l);

(Я + а, + ßi) • p(0,1) = И • p(1, l) + а2 • p(0,2) + ßQ • p(0,0); (Я + а2)• p(0,2) = p(1,2) + ßi • p(0,1); (Я + n^ + ß0 )• p(n,0) = Я^ p(n -1,0) + (n +1)^ p(n +1,0) + +а1 • p(n, 1),n = 1,...,N0 -1;

(Я + n^ + a1 + ß1 )• p(n, 1) = !• p(n -1,1) + (n +1)p(n +1,1) + +а2 • p(n,2) + ß0 • p(n,0),n = 1,...,N0; (Я + n^ + a2 )• p(n,2) = &• p(n -1,2) + (n + 1),u- p(n +1,2) + +ß1 • p(n, 1), n = 1,..., N1;

(Я + n^ + ß1 )• p(n, 1) = Я- p(n -1,1) + (n +p(n +1,1) + +а2 • p(n,2),n = N0 +1,...,N -1; (Я + n^)-p(n,2) = Я- p(n -1,2) + (n +1)^ p(n +1,2), n = N +1,...,N2 -1;

((V + ß0) • p( N0,0) = я • p( N0 -1,0)+а • /(N0,1);

(N^ + ß1) • p(N1,1) = Я • p(N1 -1,1) + а2 • p(N1, 2); N2iu • p(N2,2) = Я • p(N2 -1,2),(n, s) e Y.

(15)

В данном случае МП У^) = (N(0,S(0, t > 0), описывающий поведение системы, является обратимым. Стационарное распределение вероятностей может быть представлено в мультипликативном виде:

р(П, ^ = ^ — • Г( + • — • Г(П = V \ /(П О) (16)

где:

X TT-f(j) +1| f(s), e Y \ {0,0)

(i,j)eY\{0,0)} t ! I n!

ßk-1

f (z ) = |iinf" «k

11,

5Д 8 z > 0; 5Д 8 z = 0.

Также стационарное распределение вероятностей p(n, s) может быть найдено и при помощи рекуррентного алгоритма, предложенного в предыдущем разделе.

Численный анализ, сравнение точного и приближенного методов

Проведем оценку эффективности точного и приближенного методов на примере анализа модели, работающей в случайной среде с тремя состояниями.

Проиллюстрируем поведение показателей эффективности системы (5) - (8) в зависимости от расстояния d от пользователя до БС (рис. 3, 4) и максимальной мощности p2 передачи сигнала (рис. 5, 6). Также проанализируем зависимость исследуемых характеристик от интенсивности as, s = 1,2 снижения мощности сигнала. Сведем исходные данные в две таблицы. Таблица 2 содержит постоянные значения входных параметров, т.е. значения, которые одинаковы для всех проводимых численных экспериментов, таблица 3 - динамически меняющиеся, т.е. зависящие от рассматриваемого сценария.

Таблица 2. Постоянные значения входных параметров Table 2. Constant values of input parameters

ßß- ß-1 (О Я и'1 G к

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

30 с 60 с -60 дБм 10 Мц 5 0.1 с 1 Мбит/с 197.43 5

Таблица 3. Динамически меняющиеся значения входных параметров Table 3. Dynamically changing input parameters

а1 1 а21 d Р- pi Р0

Рис. 3 1800, 1200 с 200 - 300 м 42 дБм р2•2/3 р2 v3

Рис. 4 200 м 33 - 42 дБм

Современные информационные технологии и ИТ-образование

Том 15, № 3. 2019

ISSN 2411-1473

sitito.cs.msu.ru

K. Y. Adou, E. V. Markova, I. A. Gudkova

THEORETICAL QUESTIONS OF COMPUTER SCIENCE, COMPUTATIONAL MATHEMATICS, COMPUTER SCIENCE AND COGNITIVE INFORMATION TECHNOLOGIES

На рисунках 3 и 4 видно, что чем больше расстояние между пользователем и БС, тем хуже рассматриваемые показатели эффективности модели. Это связано с тем, что в соответствии с формулой (1) при фиксированной мощности передачи сигнала увеличение расстояния до БС приводит к уменьшению максимального количества запросов пользователей, которые могут быть обслужены с гарантированной скоростью передачи данных г0, и, следовательно, к увеличению вероятности блокировки (рис. 3) и вероятности прерывания (рис. 4).

d [м]

Р и с. 3. Вероятность блокировки B при изменении расстояния d до БС F i g. 3. The probability of blocking B when changing the distance d to the БС

3.5E-05 3.0E-05 2.5E-05 2.0E-05 1.5E-05 1.0E-05 5.0E-06 0.0E+00

Hl = a2 = l/1200c ai =a2 = l/1800c

= 1/120ÛC и of2 — l/1800c ai = 1/1800C n, = 1/1200C

12.5 14 p2 [W]

200 220 240 260 280 d Гм]

Р и с. 4. Вероятность прерывания П при изменении расстояния d до БС F i g. 4. The probability of interruption П when changing the distance d to the БС

И наоборот, при фиксированном расстоянии пользователей от БС и увеличении максимальной мощности передачи сигнала показатели эффективности модели улучшаются (рис. 5, 6), что также объясняется формулой (1). По обоим рисункам также видно, что чем реже происходит снижение мощности сигнала, тем лучше показатели эффективности системы.

2 3.5 5 6.5 8 9.5 11 12.5 14 Рг [W] Р и с. 6. Вероятность прерывания П при изменении максимальной мощности p2 передачи сигнала F i g. 6. The probability of interruption П when changing the maximum transmit power p2

В заключение проведем численный эксперимент, иллюстрирующий точность предложенного приближенного метода за счет оценки относительной погрешности S(B) расчета вероятности блокировки B

S(B ) = ((- b|/B ) х 100%, (17)

где B* - значение вероятности блокировки, полученное при использовании приближенного метода.

Из таблицы 4 видно, что чем реже происходит снижение мощности передачи сигнала, т.е. чем выше значения параметра а-, тем точнее работает приближенный метод. Это подтверждается и данными таблицы 5, показывающей, что чем чаще мощность передачи сигнала увеличивается, т.е. чем меньше значение показателя Д-1, тем точнее приближение.

Таблица 4. Относительная погрешность S (B) расчета вероятности блокировки при изменении интенсивностей а ,5 = 1,2 снижения мощности сигнала, d = 400 14 рг = 42 дБм, Д"1 = 30 с и Д"1 = 60 с Table 4. Relative error S(B) in calculating the probability of blocking when changing intensities in signal power decrease, d = 400 щр1 = 42 дБм, Д"1 = 30 сиД"1 =60 с

O-1 ) 1200 1800 2400

1200 1.451358E-02 1.45% 1.045163E-02 1.05% 8.161150E-03 0.82%

1800 1.304435E-02 9.388005E-03 7.327966E-03

1.30% 0.94% 0.73%

2400 1.230458E-02 1.23% 8.853554E-03 0.89% 6.909772E-03 0.69%

Таблица 5. Относительная погрешность S (B) расчета вероятности блокировки при изменении интенсивностей /} i = 0,1 увеличения мощности сигнала, d = 400 м,рг = 42 дБмио^"1 = сс^1 = 2400 с Table 5. Relative error S(B) in calculating the probability of blocking when changing the intensities of the increase in signal power, d =400i^^2 = 42 дБиис^"1 = д"1 = 2400 с

ßß- ) 30 60 90

30 3.673345E-03 0.37% 6.909772E-03 0.69% 9.770259E-03 0.98%

60 4.344384E-03 0.43% 8.163822E-03 0.82% 1.153455E-02 1.15%

90 5.012724E-03 0.50% 9.408299E-03 0.94% 1.327962E-02 1.33%

Р и с. 5. Вероятность блокировки B при изменении максимальной мощности p2 передачи сигнала F i g. 5. The probability of interruption B when changing the maximum transmit power p2

На основе результатов, полученных из таблиц 4 и 5, были взяты исходные данные для построения таблиц 6 и 7. Результаты, представленные в таблицах 6 и 7 показывают, что при выбранных исходных данных относительная погрешность 8(В) расчета вероятности блокировки В составляет менее 10%. Следует отметить, что относительная погрешность достигает минимального значения при рассмотрении максимального расстояния d от пользователя до БС, а также максимального значения мощности передачи сигнала р2.

Vol. 15, No. 3. 2019 ISSN 2411-1473 sitito.cs.msu.ru

Modern Information Technologies and IT-Education

570 ТЕОРЕТИЧЕСКИЕ ВОПРОСЫ ИНФОРМАТИКИ, ПРИКЛАДНОМ МАТЕМАТИКИ, К. И. Аду, Е. В. Марк0ва,

5/0 КОМПЬЮТЕРНЫХ НАУК И КОГНИТИВНО-ИНФОРМАЦИОННЫХ ТЕХНОЛОГИЙ ' И. А. Гудкова

Таблица 6. Относительная погрешность 5(B) расчета вероятности блокировки при изменении расстояния d от пользователя до БС, рг = 42 дБчffij"1 = 1200 с, af = 1800 с, Д1 = 30 с и Д"1 = 60 с. Table 6. Relative error 5 (B ) in calculating the probability of blocking when the [5] distance from the user to the БС changes.

p} = 42 дБч ffij"1 = 1200 c, a.-1 = 1800 с, Д1 = 30 с и Д"1 = 60 с.

d 300 350 400 450

S(B) 1.34417E-02 1.344% 5.84308E-06 0.001% 1.04516E-02 1.045% 5.34262E-04 0.053%

Таблица 7. Относительная погрешность 5(B) расчета вероятности блокировки при изменении максимальной мощности рг передачи [6]

сигнала, А = 300 м, ctf1 = аг-1 = 2400 с, Д"1 = 30 с и Д"1 = 60 с. Table 7. Relative error 5 (B ) in calculating the probability of blocking when the maximum signal transmission power changes, A = 300 м, of1 = a"1 = 2400 с, Д1 = 30 с и Д"1 = 60 с. [7]

p 33 36 39 42

7.16301E-02 S(B) 7.16% 5.34262E-04 0.05% 3.44599E-02 3.45% 3.20718E-03 0.32%

[8]

Заключение

В работе построена и исследована математическая модель схемы доступа к радиоресурсам беспроводной сети, реа- [9] лизующая механизмы снижения мощности и прерывания обслуживания пользователей. Модель описана с помощью СМО, работающей в случайной среде с произвольным количеством состояний. Рассмотрен один из возможных случаев реализации модели при помощи ограничения количества состояний [10] случайной среды тремя возможными. Предложен точный метод расчета (рекуррентный алгоритм) стационарного распределения вероятностей модели с тремя состояниями случайной среды, а также приближенный метод, основанный на [11] введении упрощающего предположения. Проведена оценка показателей эффективности модели - вероятности блокировки и вероятности прерывания - точным и приближенным методами. Сравнительный анализ расчета вероятности блокировки при помощи точного и приближенного методов показал, приближение максимально при высокой мощности передачи [12] данных и большом расстоянии от БС до пользователя.

References

[1] Ericsson mobility report, Q2 Update. August 2019. Ericsson, 2019. Available at: https://www.ericsson. [13] com/4a517b/assets/local/mobility-report/docu-ments/2019/emr-q2-2019-update.pdf (accessed 19.08.2019). (In Eng.)

[2] Cisco Visual Networking Index: Global Mobile Data Traffic Forecast Update, 2017-2022. White Paper. Cisco 2019. [14] Available at: https://www.cisco.com/c/en/us/solutions/ collateral/service-provider/visual-networking-index-vni/ white-paper-c11-738429.pdf (accessed 19.08.2019). (In Eng.)

[3] 3GPP TS 22.105: Services and service capabilities: Release [15] 15. 3GPP. 2018-07. ETSI, 2018. Available at: https://www. etsi.org/deliver/etsi_ts/122100_122199/122105/15.00.0 0_60/ts_122105v150000p.pdf (accessed 19.08.2019). (In Eng.)

[4] Basaure A., Sridhar V., Hammainen H. Adoption of dynamic [16]

spectrum access technologies: a system dynamics approach. Telecommunication Systems. 2016; 63(2):169-190. (In Eng.) DOI: 10.1007/s11235-015-0113-7 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

14. 3GPP. 2019-10. ETSI, 2019. Available at: https://www. etsi.org/deliver/etsi_ts/136300_136399/136300/14.11.0 0_60/ts_136300v141100p.pdf (accessed 19.08.2019). (In Eng.)

Acar Yu., Aldirmaz £olak S., Bajar E. Channel estimation for OFDM-IM systems. Turkish Journal of Electrical Engineering & Computer Sciences. 2019; 27(3):1908-1921. (In Eng.) DOI: 10.3906/elk-1803-101

Li J., Dang S., Wen M., Jiang X., Peng Y., Hai H. Layered Orthogonal Frequency Division Multiplexing With Index Modulation. IEEE Systems Journal. 2019; 13(4):3793-3802. (In Eng.) DOI: 10.1109/JSYST.2019.2918068 Ghallab R., Shokair M., Abou El-Azm A., Sakr A., Saad W., Na-guib A. Performance enhancement using multiple-input multiple-output (MIMO) electronic relay in massive MIMO cellular networks. IET Networks. 2019; 8(5):299-306. (In Eng.) DOI: 10.1049/iet-net.2018.5023 Garcia-Rodriguez A., Geraci G., Giordano L. G., Bonfante A., Ding M., Lopez-Perez D. Massive MIMO Unlicensed: A New Approach to Dynamic Spectrum Access. IEEE Communications Magazine. 2018; 56(6):186-192. (In Eng.) DOI: 10.1109/MCOM.2017.1700533

Ouyang F. Massive MIMO for dynamic spectrum access. In: 2017 IEEE International Conference on Consumer Electronics (ICCE). Las Vegas, NV, 2017, pp. 9-12. (In Eng.) DOI: 10.1109/1CCE.2017.7889210

3GPP TS 23.246: Multimedia Broadcast/Multicast Service (MBMS); Architecture and functional description: Release

15. 3GPP. 2019-10. ETSI, 2019. Available at: https://www. etsi.org/deliver/etsi_ts/123200_123299/123246/15.01.0 0_60/ts_123246v150100p.pdf (accessed 19.08.2019). (1n Eng.)

Kiji N., Sato T., Shinkuma R., Oki E. Virtual Network Function Placement and Routing Model for Multicast Service Chaining Based on Merging Multiple Service Paths. In: 2019 IEEE 20th International Conference on High Performance Switching and Routing (HPSR), Xi'An, China, 2019, pp. 1-6. (In Eng.) DO1: 10.1109/HPSR.2019.8807998

Ren Y., Chen J., Chin J., Tseng Y. Design and Analysis of the Key Management Mechanism in Evolved Multimedia Broadcast/Multicast Service. IEEE Transactions on Wireless Communications. 2016; 15(12):8463-8476. (1n Eng.) DO1: 10.1109/TWC.2016.2615605

Huang J., Zhong Z., Ding, J. An Adaptive Power Control Scheme for Multicast Service in Green Cellular Railway Communication Network. Mobile Networks and Applications. 2016; 21(6):920-929. (1n Eng.) DO1: 10.1007/ s11036-016-0712-x

Onidare S.O., Navaie K., Ni Q. On the Spectrum and Energy Efficiency in Dynamic Licensed Shared Access Systems: A Multiobjective Optimization Approach. IEEE Access. 2019; 7:164517-164532. (In Eng.) DOI: 10.1109/AC-CESS.2019.2952686

Markova E., Gudkova 1., Ometov A., Dzantiev 1., Andreev S.,

Современные информационные технологии и ИТ-образование

Том 15, № 3. 2019 ISSN 2411-1473 sitito.cs.msu.ru

K. Y. Adou, E. V. Markova, THEORETICAL QUESTIONS OF COMPUTER SCIENCE, COMPUTATIONAL MATHEMATICS,

I. A. Gudkova COMPUTER SCIENCE AND COGNITIVE INFORMATION TECHNOLOGIES

Koucheryavy Ye., Samouylov K. Flexible Spectrum Management in a Smart City within Licensed Shared Access Framework. IEEE Access. 2017; 5:22252-22261. (In Eng.) DOI: 10.1109/ACCESS.2017.2758840

[17] Borodakiy V.Y., Samouylov K.E., Gudkova I.A., Ostrikova D.Y., Ponomarenko A.A., Turlikov A.M., and Andreev S.D. Modeling unreliable LSA operation in 3GPP LTE cellular networks. In: 2014 6th International Congress on Ultra Modern Telecommunications and Control Systems and Workshops (ICUMT), St. Petersburg, 2014, pp. 390-396. (In Eng.) DOI: 10.1109/ICUMT. 2014.7002133

[18] Maule M., Moltchanov D., Kustarev P., Komarov M., Andreev S., Koucheryavy Y. Delivering Fairness and QoS Guarantees for LTE/Wi-Fi Coexistence Under LAA Operation. IEEE Access. 2018; 6:7359-7373. (In Eng.) DOI: 10.1109/AC-CESS.2018.2793941

[19] Markova E., Moltchanov D., Gudkova I., Samouylov K., Koucharyavy Y. Performance Assessment of QoS-Aware LTE Sessions Offloading Onto LAA/WiFi Systems. IEEE Access. 2019; 7:36300-36311. (In Eng.) DOI: 10.1109/AC-CESS.2019.2905035

[20] Borodakiy V., Gudkova I., Markova E., Samouylov K. Modelling and performance analysis of pre-emption based radio admission control scheme for video conferencing over LTE. In: Proceedings of the 2014 ITU kaleidoscope academic conference: Living in a converged world - Impossible without standards?, St. Petersburg, 2014, pp. 53-59. (In Eng.) DOI: 10.1109/Kaleidoscope.2014.6858480

[21] Basharin G.P., Samouylov K.E., Yarkina N.V., Gudkova I.A. A new stage in mathematical teletraffic theory. Automation and Remote Control. 2009; 70(12):1954-1964. (In Eng.) DOI: 10.1134/S0005117909120030

[22] Rykov V. Multidimensional Alternative Processes Reliability Models. In: Dudin A., Klimenok V. I., Tsarenkov G., Dudin S. Modern Probabilistic Methods for Analysis of Telecommunication Networks. BWWQT2013. Communications in Computer and Information Science, vol. 356. Springer-Verlag Berlin Heidelberg, 2013, pp. 147-156. (In Eng.) DOI: 10.1007/978-3-642-35980-4_17

[23] Ahmadian A., Galinina. O., Gudkova I., Andreev S., Shorgin S., and Samouylov K. On Capturing Spatial Diversity of Joint M2M/H2H Dynamic Uplink Transmissions in 3GPP LTE Cellular System. In: Balandin S., Andreev S., Koucheryavy Y. (eds) Internet of Things, Smart Spaces, and Next Generation Networks and Systems. ruSMART 2015, NEW2AN 2015. Lecture Notes in Computer Science, vol. 9247. Springer, Cham, 2015, pp. 407-421. (In Eng.) DOI: 10.1007/978-3-319-23126-6_36

[24] Ali A., Shah G.A., Arshad J. Energy Efficient Resource Allocation for M2M Devices in 5G. Sensors. 2019; 19(8):1830. (In Eng.) DOI: 10.3390/s19081830

[25] Su J., Xu H., Xin N., Cao G., Zhou X. Resource Allocation in Wireless Powered IoT System: A Mean Field Stackelberg Game-Based Approach. Sensors. 2018; 18(10):3173. (In Eng.) DOI: 10.3390/s18103173

[26] Chen S., Ma R., Chen H., Zhang H., Meng W., Liu J. Machine-to-Machine Communications in Ultra-Dense Networks - A Survey. IEEE Communications Surveys & Tutorials. 2017; 19(3):1478-1503, thirdquarter 2017. (In Eng.) DOI: 10.1109/œMST.2017.2678518

Поступила 19.08.2019; принята в печать 10.09.2018; опубликована онлайн 30.09.2019. Submitted 19.08.2019; revised 10.09.2019; published online 30.09.2019.

|об авторах:|

Аду Кпангни Ивь, магистр кафедры прикладной информатики и теории вероятностей, факультет физико-математических и естественных наук, Российский университет дружбы народов (117198, Россия, г. Москва, ул. Миклухо-Маклая, д. 6), ORCID: http://orcid.org/0000-0003-4669-0898, 1032135491@ rudn.ru

Маркова Екатерина Викторовна, доцент кафедры прикладной информатики и теории вероятностей, факультет физико-математических и естественных наук, Российский университет дружбы народов (117198, Россия, г. Москва, ул. Миклухо-Маклая, д. 6), кандидат физико-математических наук, ORCID: http://orcid.org/0000-0002-7876-2801, markova_ ev@rudn.ru

Гудкова Ирина Андреевна, доцент кафедры прикладной информатики и теории вероятностей, факультет физико-математических и естественных наук, Российский университет дружбы народов (117198, Россия, г. Москва, ул. Миклухо-Маклая, д. 6); старший научный сотрудник, Федеральный исследовательский центр «Информатика и управление» Российской академии наук (119333, Россия, г. Москва, ул. Вавилова, д. 442), кандидат физико-математических наук, ORCID: http://orcid. org/ 0000-0002-1594-427X, gudkova_ia@rudn.ru

Все авторы прочитали и одобрили окончательный вариант рукописи.

AbouttheaSS

Yves Adou, Master Student of the Department of Applied Probability and Informatics, Faculty of Science, Peoples' Friendship University of Russia (6 Miklukho-Maklaya Str., Moscow 117198, Russia), ORCID: http://orcid.org/0000-0003-4669-0898, 1032135491@ rudn.ru

Ekaterina V. Markova, Associate Professor of the Department of Applied Probability and Informatics, Faculty of Science, Peoples' Friendship University of Russia (6 Miklukho-Maklaya Str., Moscow 117198, Russia), Ph.D. (Phys.-Math.), ORCID: http://orcid.org/ 0000-0002-7876-2801, markova_ev@rudn.ru Irina A. Gudkova, Associate Professor of the Department of Applied Probability and Informatics, Faculty of Science, Peoples' Friendship University of Russia (6 Miklukho-Maklaya Str., Moscow 117198, Russia); Senior Scientist, Federal Research Center "Computer Science and Control" of the Russian Academy of Sciences (442 Vavilov Str., Moscow 119333, Russia), Ph.D. (Phys.-Math.), ORCID: http://orcid.org/0000-0002-1594-427X, gudkova_ia@rudn.ru

All authors have read and approved the final manuscript.

Modern Information Technologies and IT-Education

i Надоели баннеры? Вы всегда можете отключить рекламу.