ТЕОРИЯ ВЕРОЯТНОСТЕЙ И МАТЕМАТИЧЕСКАЯ СТАТИСТИКА
УДК 519.214
DOI: 10.18101/2304-5728-2018-2-3-12
О ЧИСЛЕ ЕДИНИЦ В ОДНОЙ МУЛЬТИЦИКЛИЧЕСКОЙ ПОСЛЕДОВАТЕЛЬНОСТИ С ЗАВИСИМЫМИ ЗНАКАМИ
© Меженная Наталья Михайловна
кандидат физико-математических наук, доцент,
Московский государственный технический университет им. Н. Э. Баумана Россия, 105005, г. Москва, ул. 2-я Бауманская, 5 E-mail: [email protected]
В работе рассмотрено одно обобщение классического мультициклического генератора с г регистрами, выходная последовательность которого состоит из элементов, образованных произведениями двоичных знаков в регистрах при их циклическом сдвиге друг относительно друга. Знаки, заполняющие каждый регистр, циклически да-зависимы, а регистры независимы между собой. Для случайной величины, равной числу единиц в описанной мультицикличе-ской последовательности, найдены математическое ожидание и дисперсия при помощи формулы, связывающей ее значение с количествами единиц в каждом из регистров. Доказана центральная предельная теорема для числа единиц, когда длины регистров стремятся к бесконечности, а параметры распределений знаков, заполняющих регистры, и число регистров остаются фиксированными. Рассмотрено несколько частных случаев применения предельной теоремы к последовательностям случайных величин специального вида, заполняющих регистры. Для случая независимых и неравновероятных заполнений регистров приведены численные значения скорости сходимости к предельному распределению в равномерной метрике.
Ключевые слова: мультициклическая последовательность; генератор Пола; число единиц; центральная предельная теорема; да-зависимые случайные величины.
Введение
Мультициклический генератор из г регистров над кольцом вычетов по модулю 2 был предложен в работе [1] и определен следующим образом. Пусть (Xq;),...,x2)_1), 7 = 1,...,г —векторы заполнений ячеек регистров взаимно простых длин т1,...,тг, х[;)е{0,1}, £ = 0,...,/я.-1. Знаки мультициклической последовательности образуются по правилу:
z =х(1) © ©х(г) где t{M) = i mod Л/. © — операция сложения по модулю 2. Выходная последовательность такого генератора является чисто периодической с (возможно не минимальной) длиной периода Т = т1...тг. Поэтому при исследовании ее свойств достаточно рассматривать отрезок длины Т. В работе [2] число единиц на цикле выходной последовательности ге-
нератора выражено через числа единиц в его регистрах. В работах [3-5] получены аналогичные формулы для частот знаков в выходной последовательности комбинирующего и фильтрующего генераторов [6]. С помощью этого подхода в [2] получен широкий спектр предельных теорем нормального типа для числа единиц в мультициклической последовательности, когда регистры заполнены независимыми в совокупности двоичными случайными величинами с равномерными распределениями. В [7] результаты работы [2] были обобщены на случай циклически т -зависимых [8, с. 468] заполнений внутри регистров. В [2] и [7] также получены оценки скорости сходимости в предельных теоремах в равномерной метрике [9, с. 475]. Случай неравновероятных независимых заполнений регистров рассмотрен в [10].
1. Постановка задачи
Пусть у = 1,...,г — независимые одинаково распреде-
ленные двоичные случайные величины. Построим по ним последовательности циклически т -зависимых случайных величин [8, с. 468-469]: к = 0,...,т]-1, 7 = 1,...,г. (1)
Здесь /х,...,/г — булевы функции, существенно зависящие хотя бы от одной и не более чем от т переменных, т<т/ = 1,...,г. Случай т = 1
и /х,...,/г Фсот1 соответствует последовательности независимых случайных величин.
Случайные величины имеют одинаковые одномерные
распределения в силу независимости и одинаковой распределенности случайных величин Гд7-1,...,!^^. Пусть
Р]= Р{/Д^),...,7«Й!_1)(^)) = 1}. (2)
Рассмотрим обобщение классического мультициклического генератора [1]. Построим мультициклический генератор с г регистрами, заполненными наборами Х[Г!. к = ().....т, - 1. 7 = 1,...,г, который вырабатывает выходную последовательность по правилу:
В настоящей работе мы изучим вопрос об асимптотическом распределении числа единиц на цикле длины Т = щ ... тг выходной последовательности генератора (3).
2. Число единиц на цикле выходной последовательности мультициклического генератора
т^-1
Обозначим через = — число единиц в / -м регистре,
к=О
Г-1
7 = 1,...,г, 2, = ^7.к —число единиц в цикле последовательности (3).
к=О
Лемма 1. Пусть ^,...,^>1. Тогда имеет место равенство:
^.•А- (4)
Доказательство леммы 1. Докажем формулу (4) по индукции. Пусть сначала г = 2 . Заметим, что ^ = 1} = {^¿¡^ = 1 = 1}. Всего можно
получить ^^ пар единиц при циклическом сдвиге первого и второго регистров друг относительно друга. Таким образом, при г = 2 формула (4) доказана. Остается заметить, что мультициклический генератор с г регистрами может быть представлен как мультициклический генератор с двумя регистрами, первый из которых заполнен выходной последовательностью генератора из первых г-1 регистров исходного генератора, а второй регистр совпадает с г -м регистром исходного генератора. Лемма 1 доказана.
Лемма 2. Пусть закон распределения случайных величин Х^, к = 0,...,т -1, 7 = 1,..., г, задается формулами (1) и (2), т> 1, т1,...,тг >2т — 1. Тогда математическое ожидание и дисперсия числа единиц равны
Щ = ТР1...Рг,
Г ( \
Щ = тЦ\р] + 2 £ = 1} + (|И;. -2т +1 )р)
}=\ V \<1<т-\ У
-(ТР1...РгУ. (5)
Доказательство леммы 2. Начнем с вычисления Е2,. Согласно определению случайной величины SJ :
Ш-, 1 Ш-, 1
= =1}= ¿Р{/Д^),...,7«й!_1)(й!у)) = 1} =
¿=0 ¿=0
= т]П/] (70(Л„. .,7«) = 1} = . (6)
Случайные величины ,..., независимы в совокупности, поэтому Е2, = Е^ .. .Е5,. = т1р1. ..тгрг =Тр1.. ,рг (мы воспользовались определением длины цикла Т и формулой (6)). Теперь вычислим Щ . Заметим, что
В^ = Е(^...^)2-(Е^...^)2=Е^12...Е^2-(7>1...Л)2. (7) Снова воспользуемся определением случайной величины 51.:
( т}-\
ЕЛ'; = Е
Л т}-1
2=1} = Е^ = 1} + Е 2 = 1} =
к=0 ) к=О 0<к,1<тг\-,Ш
0<к,1<тг\,кф1 \<1<тг\
Так как случайная величина Х[]') зависима со случайными величинами ху\...,х^л, иХ(Л Хи\, ТО
1<г<т-1 (
= »7,.
т<1<т;-т
]
р}+ 2 ^ Р^^Х,00 = 1} + (/иу -21И + 1)/?'
V \<1<т-\
Подставив (8) в (7), получим (5). Лемма 2 доказана. Замечание 1. Из формулы (8) нетрудно получить, что
( \
Л
(8)
ЭЛ'; = |иу
(9)
Р}- +2 £ = 1} " (2« " 1)/>.
V \<1<т-\
Замечание 2. При т = 1 и р . =1/2 (для независимых равномерно распределенных случайных величин) получим, что
Е2, = 2"Т, Щ = + (т} -1) Л - (2-Т)2 = +1)
У=1
4
( г-1
)-4<Т2 =
= 4~Т
У=1 Л
1+1 I
^ к=11<Л<.. <у4<г
/я, .../и,
Л
Обозначим
»7
= Р]+2И пхуху=1}-(2т-1)р% 7 = 1.....л
1<г<т-1
Теорема 1. Пусть закон распределения случайных величин Х^-1, к = 0,..., т. -1, у = 1,..., г задается формулами (1) и (2), длины регистров щ,...,тг—>сс, а т> 1 и ^,...,^£(0,1) остаются фиксированными,
г
^ с^ > 0 . Тогда закон распределения случайной величины
.2 У1'2/
£
Тр\ ■■■Рг
--1
(10)
сходится к стандартному нормальному закону распределения. Доказательство теоремы 1. Обозначим 51. = —
центрированное и нормированное число единиц в j -м регистре,
6
ванное и нормированное число единиц в у-м регистре, 7 = 1,...,г. Тогда (4) можно записать в виде:
„ ( \
=П(т,Р, + = ■• ■■ ^П
У=1 7=1
1+5;
РЫт1
Отсюда получим, что
п
Тр1...рг ;=1
, а,
1+5;
РЫт1
у=1
I I гк
] к=2 1<]\<...<]1<г 1=1
р>Атп
Пусть щ,=1шп{ш1,...,шг}. Так как случайные величины 51* /
стремятся по вероятности к нулю при т1,...,тг
( I Л г
[Тр,...рг ) ^
• со, то
+ 0(т-у2)
У=1
(П)
РЫт1
(запись С, = (){!) означает, что случайная величина С/ 1 ограничена по вероятности при ^ —> со). Значит, из (11) предельные при т]... ..т. —»со
ст,-
законы распределения для
Л
— 1
и =
У=1
ЛлГ,
\Тр1...рг
совпадают. Остается заметить, что для каждой из независимых в совокупности случайных величин 51*,...,выполнена центральная предельная теорема (см. теорему 19.2.1 кн. [8, с. 469]). Поэтому предельный закон распределения случайной величины у является нормальным с нулевым
' с2
средним и положительной дисперсией /и» —■ Значит, случайная ве-
/'
личина (10) также имеет в пределе при т1,...,тг —»со стандартное нормальное распределение. Теорема 1 доказана.
3. Частные случаи теоремы 1
Рассмотрим несколько частных случаев.
Пусть т = 1, Х[]>, к = 0,..,,т -1, 7 = 1,...,г, —независимые в совокупности двоичные случайные величины с одинаковым законом распределения: Р{Х^л =\} = \-Р{Х^) =0} = Р]. Тогда 1)Л\ т Р (\ р.) и
о2,=р,(1-р,).
Следствие 1. Пусть Х[]>, к = 0,...,т -1, 7 = 1,...,г, —
независимые в
совокупности двоичные случайные величины с законом распределения:
7
Р{Х[;) = 1} = 1 - Р{Х[Л = 0} = р Пусть
■со.
А»
(ОД)
остаются фиксированными. Тогда закон распределения случайной вели-
(г 1 _ V1'2/ . л
ЧИНЫ Г) =
£
1 -р,
,=1 '"•/'• у
£
Тр\ ■■■Рг
-1
сходится к стандартному нормаль-
ному закону.
В таблице 1 приведены значения расстояния в равномерной метрике с! между функцией распределения случайной величины г| и функцией распределения стандартного нормального закона Ф для мультициклической последовательности, образованной двумя регистрами.
Таблица 1
Значения с!(1\, Ф) для мультициклической последовательности
с двумя регистрами
щ =10; т2= 11 щ = 30; т2= 31 щ = 50; т2= 51 щ =100; т2 =101
рх = 0.2; р2 = 0.2 0.186 0.106 0.076 0.055
Р1=0.ЪР2=0А 0.134 0.071 0.053 0.037
рх = 0.2; р2 = 0.5 0.111 0.061 0.046 0.032
рх = 0.2; р2 = 0.6 0.095 0.054 0.040 0.027
рг =0.2;р2 =0.8 0.103 0.041 0.032 0.022
р,=0А-р2=0А 0.109 0.069 0.052 0.036
Р1=0А;р2=0.5 0.093 0.056 0.039 0.026
Р1=0А;р2=0.6 0.085 0.048 0.031 0.021
Р1=0А;р2=0Я 0.069 0.040 0.031 0.022
рх = 0.5; р2 = 0.5 0.097 0.058 0.044 0.031
рх = 0.5; р2 = 0.6 0.084 0.046 0.033 0.021
рг =0.5;р2 =0.8 0.068 0.030 0.021 0.013
рх = 0.6; р2 = 0.6 0.091 0.051 0.041 0.029
рг =0.6;р2 =0.8 0.060 0.027 0.019 0.011
рг =0.%;р2 =0.8 0.088 0.052 0.040 0.028
Замечание 3. Утверждение следствия 1 для равновероятных распределений случайных величин Х[Г) ( р. =1/2. / = 1...../') получено в [11].
Пусть теперь т = 2, 70(; -1,..., , / = 1,..., г — независимые одинако-
во распределенные двоичные случайные величины:
Р{7«=1} = 1-Р{7« = 0 } = р,
к = 0,...,т}-\, 7 = 1,..„г. (13)
В этом случае
Р{Х« = 1} = = 1} = Р{У0О) е Y™ = 1} = 2р(\ - р),
v{x(kf> = xi» = i} = p{x¿;) = х(хг> = 1} = р{70(;) е y¡s> = y¡s> е г2(;) = 1} = = (1 -р)2р+р2 (1 -р) = р( 1 - р), при этом элементы последовательности X[J'. к = (),.. ..т . - 1. не являющиеся соседними, независимы между собой (считаем элементы Х[]') и соседними). Тогда формула (9) перепишется в виде:
OS2 = т} (Р] + 2F{XPXP = 1} - (4 -1 )р2) =
= т}. (2р(1 -р) + 2р(\ -р)- 3(2/7(1 - р))2) = Ар(\ - р)т} (1 - Зр(1 - р)).
Значит, а2 = 4р{ 1 - р) (1 - 3р{ 1 — р)) и имеет место следующее утверждение.
Следствие 2. Пусть случайные величины Хк]\ к = 7 = 1,...,г определены равенствами (12) и (13). Пусть т1,...,тг^ со, р е (0,1) фиксировано. Тогда закон распределения случайной величины:
Í Л-1/2 . s
(\-Ър(\-р)^\} ( •= >
р(\~р) Ttm,
ТТрг(\-рУ
-1
1=1"'] /
сходится к стандартному нормальному закону.
Замечание 4. При р = 1/2 случайные величины Х[]>, к = 0,... -1, 7 = 1,...,г, определенные равенствами (12) и (13), независимы в совокупности. В этом случае утверждения следствий 1 и 2 совпадают.
Рассмотрим еще один частный случай, когда распределения случайных величин вместо (13) задаются равенством:
X™ к = 0,...,т}.-1, 7 = 1,..,г, (14)
а 7 = 1,...,г — независимые одинаково распределенные
двоичные случайные величины, удовлетворяющие (12). Тогда = 1} = Р{Х« = 1} = Р^Г/» = 1} = р\ = = 1} = = Х[]) = 1} = Р{70°')71°') = У^У^1 = 1 } = р3.
Аналогично следствию 2 имеет место следующее утверждение.
Следствие 3. Пусть случайные величины Хк]\ к = 0,..,,т -1, 7 = 1,..., г определены равенствами (12) и (14). Пусть т1,...,тг^- со, р е (0,1) фиксировано. Тогда закон распределения случайной величины
\ + 2р-Ър2^ 1 W s >
v Р2 [Tp2r
сходится к стандартному нормальному закону.
-1
Заключение
В работе рассмотрено обобщение классического мультициклического генератора с г регистрами над кольцом вычетов по модулю 2. Знаки внутри каждого регистра циклически га-зависимы, а регистры независимы между собой. Доказана центральная предельная теорема для числа единиц в цикле мультициклической последовательности, когда длины регистров стремятся к бесконечности, а параметры распределений знаков, заполняющих регистры, и их число остаются фиксированными. Рассмотрено несколько частных случаев применения предельной теоремы к последовательностям случайных величин специального вида. Получены численные значения скорости сходимости к нормальному распределению в равномерной метрике для независимых и неравновероятных заполнений регистров.
Литература
1. Pohl P. Description of MCV, a pseudo-random number generator // Scand. Actuar. J. 1976. Vol. 1. P. 1-14. DOI: 10.1080/03461238.1976.10405931.
2. Меженная H. M., Михайлов В. Г. О распределении числа единиц в выходной последовательности генератора Пола над полем GF(2) // Математические вопросы криптографии. 2013. Т. 4, № 4. С. 95-107. DOI: 10.4213/mvkl01.
3. Биляк И. Б., Камловский О. В. Частотные характеристики циклов выходных последовательностей комбинирующих генераторов над полем из двух элементов // Прикладная дискретная математика. 2015. Т. 3, № 29(3). С. 17-31. DOI: 10.17223/20710410/29/2.
4. Камловский О. В. Количество появлений векторов на циклах выходных последовательностей двоичных комбинирующих генераторов // Проблемы передачи информации. 2017. Т. 53, № 1. С. 84-91. DOI: 10.1134/S0032946017010070.
5. Камловский О. В. Количество появлений элементов в выходных последовательностях фильтрующих генераторов // Прикладная дискретная математика. 2013. Т. 3, № 21. С. 11-25.
6. Агибалов Г. П. Конечные автоматы в криптографии // Прикладная дискретная математика. Приложение. 2009. Т. 2. С. 43-73.
7. Mezhennaya N. M. Convergence rate estimators for the number of ones in outcome sequence of MCV generator with m-dependent registers items // Sib. Electron. Math. Reports. 2014. Vol. 11. P. 18-25.
8. Ибрагимов И. A., Линник Ю. В. Независимые и стационарно связанные величины. М.: Наука, 1965. 524 с.
9. Ширяев А. Н. Вероятность-1. 4-е изд. М.: Изд-во МЦНМО, 2011. 552 с.
10. Меженная H. М. О распределении числа единиц в двоичной мультициклической последовательности // Прикладная дискретная математика. 2015. Т. 1(27). С. 69-77.
11. Mezhennaya N. M., Mikhailov V. G. Limit theorem for number of ones in the extended Pohl generator outcome sequence // OP&PM Surveys on Applied and Industrial Mathematics. 2018. Vol. 25, № 1. P. 48-50.
ON THE NUMBER OF ONES IN A MULTI-CYCLIC SEQUENCE WITH DEPENDENT SIGNS
Natalya M. Mezhennaya Cand. Sci. (Phys. and Math.), A/Prof., Bauman Moscow State Technical University 5 2nd Bauman St., Moscow 105005, Russia E-mail: [email protected]
The article considers one extension of a classical multi-cyclic generator with r registers, which output sequence consists of elements formed by the products of binits in the registers under their cyclic shift relative to each other. The signs that fill each register are cyclically ///-dependent. and the registers are independent of each other. We have found the mathematical expectation and variance for a random variable equal to the number of ones in the presented multi-cyclic sequence using the formula connecting its value with the number of ones for each registers. The central limit theorem for the number of ones is proved under the conditions when the lengths of registers tend to infinity, and the parameters of signs distributions filling the registers and the number of registers are fixed. We consider several particular cases of the limit theorem application to the sequences of random variables of a special type filling the registers. The numerical values of the convergence rate to the limiting distribution in the uniform metric for the case of independent and nonuniform fillings of registers are given.
Keywords: multi-cyclic sequence; Pohl generator; number of ones; central limit theorem; /»-dependent random variables.
References
1. Pohl P. Description of MCV, a Pseudo-Random Number Generator. Scand. Actuar. J. 1976. V. 1. Pp. 1-14. DOI: 10.1080/03461238.1976.10405931.
2. Mezhennaya N. M., Mikhailov V. G. O raspredelenii chisla edinits v vykhodnoi posledovatel'nosti generatora Pola nad polem GF (2) [On the Distribution of the Number of Ones in the Output Sequence of MCV-Generator over GF(2)]. Matematicheskie voprosy kriptografii —Mathematical Aspects of Cryptography. 2013. V. 4. No. 4. Pp. 95-107. DOI: 10.4213/mvkl01.
3. Bilyak I. B., Kamlovskii O. V. Chastotnye kharakteristiki tsiklov vykhodnykh posledovatel'nostei kombiniruyushchikh generatorov nad polem iz dvukh elementov [Frequency Characteristics of Cycles in Output Sequences Generated by Combining Generators over a Field of Two Elements]. Prikladnaya diskretnaya matematika — Applied Discrete Mathematics. 2015. V. 3, No. 29 (3). Pp. 17-31. DOI: 10.17223/20710410/29/2.
4. Kamlovskii O. V. Kamlovskii O. V. Kolichestvo poyavlenii vektorov na tsik-lakh vykhodnykh posledovatel'nostei dvoichnykh kombiniruyushchikh generatorov [Occurrence Numbers for Vectors in Cycles of Output Sequences of Binary Combining Generators]. Problemy peredachi informatsii — Problems of Information Transmission. 2017. V. 53, No. 1. Pp. 84-91. DOI: 10.1134/S0032946017010070.
5. Kamlovskii O. V. Kolichestvo poyavlenii elementov v vyhodnyh posle-dovatel'nostyah fil'trayuschih generatorov [Distribution Properties of Sequences Produced by Filtering Generators]. Prikladnaya diskretnaya matematika —Applied Discrete Mathematics. 2013. V. 3, No. 21. Pp. 11-25.
6. Agibalov G. P. Konechnye avtomaty v kriptografii [Finite Automata in Cryptography], Prikladnaya diskretnaya matematika. Prilozhenie — Applied Discrete Mathematics. Supplement. 2009. V. 2. Pp. 43-73.
7. Mezhennaya N. M. Convergence Rate Estimators for the Number of Ones in Outcome Sequence of MCV-generator with ///-dependent Registers Items. Sib. Electron. Math. Reports. 2014. V. 11. Pp. 18-25.
8. Ibragimov I. A., Linnik Yu. V. Nezavisimye i statsionarno svyazannye velichiny [Independent and Stationary Related Variables]. Moscow: Nauka Publ., 1965, 524 p.
9. Shiryaev A. N. Veroyatnost'-l [Probability-1]. 4th ed. Moscow: MTsNMO, 2011. 552 p.
10. Mezhennaya N. M. O raspredelenii chisla edinits v dvoichnoi mul'titsik-licheskoi posledovatel'nosti [On Distribution of Number of Ones in Binary Multicycle Sequence]. Prikladnaya diskretnaya matematika — Applied Discrete Mathematics. 2015. V. 1(27). Pp. 69-77.
11. Mezhennaya N. M., Mikhailov V. G. Limit Theorem for Number of Ones in the Extended Pohl Generator Outcome Sequence. OP&PM Surveys on Applied and Industrial Mathematics. 2018. V. 25. No. 1. Pp. 48-50.