CC BY
29
УДК 681.5 (519.95)
ПАРАМЕТРЫ ПОТОКА ТРАНЗАКЦИЙ, ГЕНЕРИРУЕМЫХ ПО ПРИНЦИПУ ПОЛЛИНГА
А. А. Аршакян, Е.В. Ларкин
Исследуется поллинг как способ опроса периферийных устройств. Определен типовой вид функции распределения времени между двумя транзакциями при поллинге. Показано, что типовую функцию распределения целесообразно аппроксимировать функцией, имеющей более компактный вид. С использованием метода неопределенных множителей Лагранжа получена система уравнений для расчета параметров аппроксимирующей функции. Показано, что если в структуре алгоритма имеется хотя бы один цикл, то функция распределения может быть аппроксимирована композицией экспоненциальных законов. С помощью компьютерного эксперимента, проведенного по методу Монте-Карло, показано, что в этом случае аппроксимирующей функцией может быть Гамма-распределение.
Ключевые слова: алгоритм, поллинг, поток транзакций, функция распределения, аппроксимация, гамма-распределение.
Поллинг как способ опроса периферийных устройств достаточно часто встречается в системах цифрового управления объектами и в информационных системах [1, 2, 3]. Алгоритмы, реализующие опрос периферийных устройств, обладают целым рядом особенностей, которые достаточно полно отражены в известных работах [4, 5, 6, 7]:
алгоритмы являются циклическими, т.е. они имеют оператор начала, но не имеют оператора окончания, окончание вычислительного процесса производится через внешнее прерывание;
опрос периферийных устройств производится за счет включения в алгоритм специальных операторов управления транзакциями, количество которых в алгоритме может быть произвольным;
выбор ветви продолжения вычислительного процесса в местах ветвления алгоритма является случайным и определяется условиями, включенными в операторы принятия решения и законами распределения обрабатываемых данных;
время выполнения операторов алгоритма является случайным, функции распределения времени выполнения операторов также определяются законами распределения обрабатываемых данных.
В подобном аспекте алгоритмы информационно-управляющих систем рассматривались в [8], где процесс интерпретации детерминированного алгоритма для внешнего, по отношению к системе, наблюдателя рассматривался как полумарковский процесс с непрерывным временем. Операторы алгоритма рассматривается как состояния процесса, а интерпрета-
ции алгоритма в целом рассматривается как блуждание по состояниям. Поток транзакций формируется когда полумарковский процесс через случайные моменты времени попадает в состояния, моделирующие обращения к периферийному устройству.
В [9, 10] показано, как можно определить время между двумя транзакциями при поллинге и получены зависимости для оценки названных временных интервалов, которые в общем случае имеют вид
f (t )= JPjfj (t), (1)
j=1
где P j - вероятность того, что реализацией блужданий по полумарковской
цепи между двумя состояниями, моделирующими один и тот же, или пару операторов транзакций будет j-я конкретная траектория; fj (t) - плотность
распределения времени блуждания по j-й конкретной траектории.
Использование f (t) в виде (1) для анализа информационно-управляющих систем вызывает известные трудности, и поэтому целесообразно аппроксимировать указанное выражение законом, часто используемым при анализе, например систем массового обслуживания [11]. Представим (1) в виде функции распределения
t J J
F (t )= J I Pjfj (t)dt = I PjFj (t), (2)
0 j=1 j=1
t
где Fj (t) = J fj (t)dt - функция распределения времени до наступления j-го 0
события.
Очевидно, что F (t) - функция неубывающая, так как является линейной композицией неубывающих функций Fj (t), F (0) = 0; lim F (t ) = 1,
t
J
т к. I pj = 1 и lim Fj (t ) = 1.
j=1 t®¥
Аппроксимируем функцию F (t) законом распределения g (t, A), где A = (ü1,..., an,..., а^) - вектор параметров. В общем случае ошибка аппроксимации определяется выражением
e= J
0
IPjFj (t)- G(t, A) j=1
dt, (3)
t
где G(t, A) = J g(t, A)dt. 0
2
СЮ
Отметим, что критерий (3) имеет вид суммы квадратов разности, а не взвешенной суммы квадратов разности, которая используется в задачах выравнивания статистических рядов [12]. Это происходит вследствие того, что плотность распределения вида (1) уже отражает только существенные особенности статистического материала, по которому были получены плотности fj (г) и вероятности pj. Случайности, вызывающие отклонение
статистик, связанные с недостаточным объемом экспериментальных данных, были учтены при получении казанных плотностей и вероятностей [3, 5, 7, 8].
Для аппроксимирующего закона g(г, А) должны выполняться ограничения:
I dG (г, А) = 1; г=о
I tdG (г, А) = Т; г=о
(4)
I(г - Т)kdG(t, А) = тк , 2 < к < К, г=о
3
где Т = |г £ PjdFj (г) - математическое ожидание времени между двумя о j=1
¥ 2 3
транзакциями при поллинге; тк = I (г - Т) £ PjdFj (г) - к-й центральный
0 j=l
момент плотности распределения (1), в частности, при к = 2 - это дисперсия плотности распределения времени между двумя транзакциями.
Для поиска оптимальных значений параметров можно воспользоваться методом неопределенных множителей Лагранжа. Функция Лагран-жа будет иметь вид [13]
Ф( А, В) = I
3
£ Р]¥] (г)- в(г, А)
j=1
dt +
+ Ьг
1 - I dG(г, А)
+ Ь1
Т - I tdG (г, А)
К
+ £ Ьк к=2
тк -1 (г - Т)kdG(г, А)
, (5)
где В = (Ьо, Ь1, Ь2,..., Ьк) - вектор неопределенных множителей Лагранжа.
В соответствии с методом неопределенных множителей Лагранжа необходимые условия существования экстремума будут иметь вид:
сю
сю
сю
2
сю
о
сю
о
о
о
ЭФ
Э а„
= 2|
J
7=1
Щи А) да,,
Ж +
о
Ь0+Ь^+ | Ък{1-Т)к к=2
¿дСМ = 0. ! < „ < дг да,,
(6)
Общая система (4) и (6) включает ТУ + + 1 уравнения и должна решаться относительно неизвестных элементов векторов Л = и Л = (Ь0,¿>2---
Отметим, что из распределений, имеющих область распределения / > О, экспоненциальное и вырожденное распределение включает один независимый варьируемый параметр, распределения равномерное, логнор-мальное, бета-, гамма-, Вейбулла, имеют по два независимых варьируемых параметра, ряд Грама-Шарлье имеет произвольное количество независимых варьируемых параметров [14].
Если плотности распределения fj(t) представляют собой вырожденные законы распределения, т.е. /у (/)= 8(/ -7у), то (1) представляет собой взвешенную сумму единичных функций Хевисайда г|(/-7у):
3
7=1
(7)
т.е. является ступенчатой функцией.
Предположим, что внутри алгоритма, реализующего поток транзакций, имеется хотя бы один фрагмент со структурой г = Гд о)' описывае"
О О
мый полумарковской матрицей /*(?) =
(рис. 1).
"2
<12 8(*-Т2)
Рис. 1. Внутренний цикл генератора транзакций
В полумарковском процессе, приведенном на рис. 1, время достижения состояния С2 из состояния с\ меняется по линейному закону в зависимости от количества возвратов в состояние с\. Вероятность достижения состояния <?2 в зависимости от количества возвратов в состояние С\ убы-
43
вает по экспоненте. Следовательно, ступенчатая функция распределения (7) для рассматриваемого случая достаточно хорошо ложится на экспоненту, что показано на рис. 2.
Рис. 2. Функция распределения временных интервалов
Проведем прямой компьютерный эксперимент с применением метода Монте-Карло по исследованию параметров полумарковского генера-
, приведенного на рис. 3, у которого
1 1Л
•5 1-5 ~ -5
тора транзакций М = [Л^, г
5Й= {1,2,3}; Зг =
1/з-/(?) 1/з-/(0 1/з-/(0' 1/3-/М 1/з-/(0 1/3-/Й 1/3-/(0 1/з-/(0 1/3-/(0
5/3, если 0,7 < Г < 1,3; 0, если / < 0,7, или ? > 1,3.
Математическое ожидание и дисперсия плотности распределения
^г) равны, соответственно Ттп = 1 с, Отп = о,о3 с ,1 < т, п < 3. Возврат
генератора в состояние 1 означает генерацию транзакции. Компьютерный эксперимент будем проводить по следующей методике [15, 16]:
1) обнуление счетчика числа реализаций, I = о.
2) присвоение статуса текущего первому состоянию щ = 1.
3) обнуление счетчика времени, г1 = о.
4) запуск генератора случайных чисел и получение случайного числа 1 < р < 1 с равновероятным распределением.
5) если о < р < о,333, то 01 = 1; если о,333 < р < о,666, то 01 = 2; если о,666 < р < 1, то 01 = 3.
6) запуск генератора случайных чисел и получение числа р.
7) Определение приращения времени по зависимости:
Аг = о,6р + о,7 ^
8) г1 := ц + А г;
9) если П1 ф 1, то 4).
10) разгрузка г1 в массив результатов.
11) I := I +1. Если I < Ь, то 3).
12) конец эксперимента.
В результате формируется гистограмма, показанная на рис. 4 штрихами.
Оценка математического ожидания и дисперсии плотности распределения времени между транзакциями по гистограмме дают ^ 1 = 3с;
1Л11 = 6,123 с2.
Гистограмму рекомендуется аппроксимировать с помощью Г-распределения, имеющего вид
f (г) =
о при г < о,
1 га-1 ехр(-1г) при г < о, (8)
Г(а)
где а > о - параметр формы; 1 > о - параметр масштаба;
¥
Г(а) = I ха-1 ехр(-1х)сХ. о
Параметры 1 и а распределения (8) определяются по зависимостям
1Т11
1=т^; (9)
Ч1
(10)
Vi,!
0,5 0,4 0,3 0,2 ОД
Рис. 4. Гистограмма, времени между транзакциями, сформированная по результатам эксперимента, и аппроксимирующее гамма-распределение
Экспериментальные значения а и X равны соответственно а = 1,47, А, = 0,49 с1.
Таким образом, в данном конкретном случае параметры потока транзакций, генерируемых в результате поллинга, достаточно точно аппроксимируется с помощью Г-распределения.
Список литературы
1. Ивутин А.Н., Ларкин Е.В. Обобщенная полумарковская модель алгоритма управления цифровыми устройствами // Известия Тульского государственного университета. Технические науки. Вып. 1. Тула: Изд-во ТулГУ, 2013. С. 221 -228.
2. Larkin E.V., Ivutin A.N. Estimation of Latency in Embedded Real-Time Systems // 3-rd Meditteranean Conference on Embedded Computing (MECO-2014). 2014. June 15-19. Budva, Montenegro, 2014. P. 236 - 239.
3. Ларкин E.B., Ивутин A.H. Определение временных интервалов в алгоритмах управления // Известия Томского политехнического университета. № 5. Т. 324. 2014. С. 6-12.
4. Клинцов Т.Н., Ларкин Е.В. Управление системами тоннелепро-ходческого комплекса при помощи программного пакета SC ADA // Известия Тульского государственного университета. Технические науки. Вып. 9. Ч. 1. Тула: Изд-во ТулГУ, 2014. С. 156 - 162.
5. Ивутин А.Н., Ларкин Е.В. Временные и вероятностные характеристики транзакций в цифровых системах управления // Известия Тульского государственного университета. Технические науки. Вып. 1. Тула: Изд-во ТулГУ, 2013. С. 252 - 258.
6. Ивутин А.Н., Ларкин Е.В. Прогнозирование времени выполнения алгоритма // Известия Тульского государственного университета. Технические науки. Вып. 3. Тула: Изд-во ТулГУ, 2013. С. 301 - 315.
7. Ларкин Е.В., Котов В.В., Котова Н.А. Оценка эффективности программного обеспечения робота с использованием сетей Петри-Маркова // Известия Тульского государственного университета. Технические науки. Вып. 9. Ч. 2. Тула: Изд-во ТулГУ, 2013. С. 156 - 163.
8. Клинцов Г.Н., Ларкин Е.В. Полумарковое моделирование алгоритмов выполнения программы на ЭВМ // Известия Тульского государственного университета. Технические науки. Вып. 9. Ч. 1. Тула: Изд-во ТулГУ, 2015. С. 143 - 149.
9. Ларкин Е.В., Привалов А.Н., Клепиков А.К. Алгоритм последовательных упрощений при оценке временной сложности алгоритмов с применением полумарковских моделей // Известия Тульского государственного университета. Технические науки. Вып. 9. Ч. 2. Тула: Изд-во ТулГУ,
2014. С. 83 - 98.
10. Larkin E.V., Lutskov Yu.I., Ivutin A.N., Novikov A.S. Simulation of concurrent process with Petri-Markov nets // Life Science Journal. 2014. N. 11. P. 506 - 511.
11. Ларкин Е.В., Котов В.В., Титов С.В. Аппроксимация взвешенной суммы плотностей распределения вероятностным законом // Известия ТулГУ. Серия: Проблемы специального машиностроения. Вып. 3. Ч. 1. Тула: ТулГУ, 2000. С. 389 - 393.
12. Румшинский Л.З. Математическая обработка результатов эксперимента: Справочное руководство. М.: Наука, 1971. 192 с.
13. Аоки М. Введение в методы оптимизации. М.: Наука, 1977.
314 с.
14. Хастингс Н., Пиккок Дж. Справочник по статистическим распределениям. М.: Статистика, 1980. 96 с.
15. Ivutin A.N, Larkin E.V. Simulation of Concurrent Games // Bulletin of the South Ural State University. Series: Mathematical Modelling, Programming and Computer Software. Chelyabinsk. 2015. Vol. 8. № 2. P. 43 - 54.
16. Котов В.В., Котова Н.А., Ларкин Е.В. Метод имитационного моделирования систем с использованием сетей Петри-Маркова // Известия Тульского государственного университета. Технические науки. Вып. 9.
2015. С. 164 - 170.
Аршакян Александр Агабегович, канд. техн. наук, докторант, elarkinamail. ru, Россия, Тула, Тульский государственный университет,
Ларкин Евгений Васильевич, зав. кафедрой, докт. техн. наук, проф., elarkin@,mail.ru, Россия, Тула, Тульский государственный университет
PARAMETERS OF FLOW OF TRANSACTIONS, WHICH ARE GENERATED ON POLLING PRINCIPLE
A,A. Arshakyan, E. V. Larkin
Polling as the technique of scanning of peripherals is investigated. Typical form of distribution function of time between neighbor transactions, when polling process, is defined. It is shown that the typical distribution function it is expediently to approximate with a compact function. With use Lagrange multipliers method the system of equations for calculation of parameters of approximated function is obtained. It is shown that when in a structure of algorithm there is at least one cycle, distribution function may be the exponential one. With use of Monte-Carlo method it is shown that in this case as the approximated function may be choose Gamma-function
Key words: algorithm, polling, flow of transactions, distribution function, approximation, Gamma-function.
Arshakyan Alexander Agabegovich, candidate of technical science, postgraduate, elarkin@,mail. ru, Russia, Tula, Tula State University,
Larkin Eugene Vasilyevich, head of chair, doctor of technical science, professor, elarkin@,mail. ru, Russia, Tula, Tula State University
УДК 623.320
АЛГОРИТМ УСТРАНЕНИЯ ИНТЕРФЕРЕНЦИИ ПОСЛЕ ЛИНЕЙНОГО СИНТЕЗА АПЕРТУРЫ АНТЕННЫ РАДИОЛОКАЦИОННОГО ДАТЧИКА
Ю.И. Вареница, Ю.И. Мамон, К.А. Хомяков
Предложен алгоритм устранения интерференции сигнала радиолокационного датчика без отказа от использования монохроматического зондирующего сигнала. Показано, что для уменьшения интерференции в лепестках восьмёрки передаточных функций точек траектории объекта, нужно просуммировать или перемножить значения взаимно корреляционных функций в точках пересечения корреляционных лучей. По наборам параллельных корреляционных лучей под разными углами можно восстановить приемлемое изображение объекта.
Ключевые слова: радиолокационный датчик, опорный траекторный сигнал, интерференция, синтезированная апертура, антенна, передаточная функция, корреляционная функция, матрица изображения.
Известен теоретический аппарат синтеза апертуры антенны радиолокационного датчика на базе применения взаимных корреляционных функций объектного траекторного сигнала с набором опорных траекторных сигналов
48
