УДК 519.214
ПРЕДЕЛЬНЫЕ ТЕОРЕМЫ ДЛЯ ЧИСЛА ПЛОТНЫХ F-РЕКУРРЕНТНЫХ СЕРИЙ И ЦЕПОЧЕК В ПОСЛЕДОВАТЕЛЬНОСТИ НЕ ЗАВИСИМЫХ СЛУЧАЙНЫХ ВЕЛИЧИН
Н.М. Меженная
МГТУ им. Н.Э. Баумана, Москва, Российская Федерация e-mail: [email protected]
Изучены свойства распределения чисел плотных F-рекуррентных серий и цепочек заданной длины в последовательности независимых одинаково распределенных случайных величин над конечным алфавитом. С помощью функционального варианта метода Чена - Стейна получены оценки скорости сходимости распределения вектора из чисел плотных F-рекуррентных серий заданных длин к сопровождающему многомерному пуассоновскому распределению (в метрике расстояния по вариации). Выведены многомерные предельная теорема Пуассона и центральная предельная теорема для чисел плотных F-рекуррентных серий заданной длины и длины не меньше заданной при подходящем изменении параметров схемы (длины последовательности и длины серии). Оценки расстояния по вариации также позволяют получить условия сходимости распределения числа плотных F-рекуррентных цепочек заданной длины к сложному пуассоновскому распределению.
Ключевые слова: плотные серии, рекуррентные цепочки, предельная теорема Пуассона, центральная предельная теорема, метод Чена-Стейна, оценки скорости сходимости.
LIMIT THEOREMS FOR DENSE F-RECCURENT SERIES AND CHAINS NUMBERS IN SEQUENCE OF INDEPENDENT RANDOM VARIABLES
N.M. Mezhennaya
Bauman Moscow State Technical University, Moscow, Russian Federation e-mail: [email protected]
The work is devoted to studying properties of distribution of dense F-recurrent series and chains with given length in the sequence of independent identically distributed random variables over finite alphabet. Using functional modification of Chen-Stein method we investigate estimators of convergence rate of distribution of numbers of dense F-recurrent series with given lengths to accompanying multivariate Poisson distribution (in metric of distance by variance). On the ground of this results multivariate Poisson limit theorem and central limit theorem for numbers of dense F -recurrent series with given lengths and length no less than given were derived under appropriate variation of scheme parameters. The estimators of distance by variance also allow to derive conditions under which distribution ofdense F-recurrent chains with given length converges to compound Poisson distribution.
Keywords: dense runs, recurrent chains, Poisson limit theorem, central limit theorem, Chan-Stein method, convergence rate estimators.
Введение. Пусть Xb X2,..., XT,... — последовательность независимых случайных величин, принимающих значения из алфавита AN = {0,1,... ,N — 1}. Такие случайные величины принято называть случайными знаками, или просто знаками. Пусть f : AlN ^ AN. Согласно работе [1], случайные знаки Xj,..., Xj+i+s—1 образуют
f-рекуррентную цепочку длиной 8 (з > I + 1), если
Хг+1+к = f (Хг+к, . . . , Xi+l+k-1), к = 0, 1, 2, . . . , 8 — 1.
Если это соотношение при каждом значении к выполняется хотя бы для одной функции из набора ^ = {Д,...,/к}, то знаки X,..., образуют ^-рекуррентную цепочку длиной 8. Оче-
видно при К =1 определения f-рекуррентной и ^-рекуррентной цепочек совпадают.
Так, пусть N = 2, А = {0,1} и ^ = ^}, где f : А2 ^ А2, f (а, Ь) = а ф Ь — сложение по модулю 2, а, Ь € А2. Тогда знаки 1, 0,1,1, 0,1,1, 0,... образуют ^-рекуррентную цепочку.
Далее f- и ^-рекуррентные цепочки будем называть сплошными f- и ^-рекуррентными цепочками.
В настоящей работе рассмотрим более общую постановку задачи. Напомним, что знаки Х1, Х2,..., Х8 образуют плотную 1-цепочку длиной 8, если 1 € {Х^,Хт}, г = 1,...,8 — 1 [2-4]. По аналогии с понятием "плотная цепочка" определим понятие "плотная ^-рекуррентная цепочка".
Будем утверждать, что знаки X»,... ,Х^+в+1-1 образуют плотную ^-рекуррентную цепочку длиной 8, если при всех к = 0,1,..., в — 2 найдутся такие функции fk , дк € F, что имеет место хотя бы одно из равенств
Хг+1+к = fk (Хг+к, . . . , Х»+1+к-1)
или
Хг+1+к+1 = Як (Хг+к+1, . . . , Хг+1+к).
Началом F-рекуррентной цепочки Х»,..., Х+3+1-1 будем считать знак
Хг+1.
Пусть N = 2, А2 = {0,1} и F = {f}, где f — сложение по модулю 2. Тогда знаки 1, 0,1,1,1, 0,1,1, 0 не образуют сплошную F-рекуррентную цепочку (курсивом выделен знак, на котором нарушается свойство F-рекуррентности), но формируют плотную F-рекуррентную цепочку.
Пусть N = 3, Аз = {0,1, 2} и F = {Л^}, где £ : А2 ^ Аз, ^(а,Ь) = г, а,Ь € А3, г = 1, 2.
Пусть 1, 2, 1, 2, 0, 1, 0, 2, 0, 0 — первые 10 знаков реализации случайной последовательности Х1, Х2, . . . Первые девять знаков, выделенные курсивом, образуют плотную F-рекуррентную цепочку. Однако в нее нельзя включить последний знак. Отметим также, что выделенные знаки не образуют сплошную F-рекуррентную цепочку, так как на пятом месте нарушается свойство F-рекуррентности (прообразом знака 0 при любом отображении из F будет 0).
В настоящей работе будет доказана многомерная предельная теорема Пуассона для чисел плотных F-рекуррентных серий в после-
довательности независимых случайных величин с оценками скорости сходимости.
Пусть ,..., £r — неотрицательные целочисленные случайные величины. Многомерная предельная теорема Пуассона указывает условия, при которых для всех k1,... ,kr G {0,1,...}:
Г /\ks \
P{6 = ki,...,ir = kre-Xs) ,
где Аь ..., Ar > 0 — параметры предельного пуассоновского распределения.
В работе [1] с помощью многомерной теоремы Севастьянова [5, 6] была доказана предельная теорема пуассоновского типа для чисел f -рекуррентных s-цепочек, а также для ряда связанных с ними случайных величин. В работе [7] распределения этих случайных величин были исследованы с использованием многомерного варианта метода Чена-Стейна [8-11]. Преимущество этого метода по сравнению с большинством классических методов доказательства предельных теорем пуассоновского типа для сумм зависимых случайных индикаторов заключается в том, что он позволяет получить оценки скорости сходимости в предельных теоремах. Для этого достаточно оценок для смешанных моментов второго порядка величин из набора случайных индикаторов (в отличие от других известных методов). Полученные оценки в некоторых случаях позволяют доказать асимптотическую нормальность (как результат сближения с пуассоновским распределением с возрастающим параметром).
Оценки и предельные теоремы. Пусть
Et = {Xt = f (X-,..., Xt-i), f G F}, t = 1, 2,...
Тогда плотной F-рекуррентной цепочке в последовательности X1,X2,...,XT,... соответствует плотная цепочка из единиц в последовательности I{Ei}, I{E1+1},... [2, 3].
Если в последовательности индикаторов I{Ei},I{Ei+1},... возникает плотная цепочка из единиц (плотная 1-цепочка), то на каждом ее конце будет не более чем по одному нулю. Если к такой плотной 1-цепочке (с нулями на концах) приписать еще по одному нулю с двух сторон, то ее нельзя продлить ни в одну сторону. Эту цепочку принято называть плотной 1-серией [2, 3]. Поэтому будем утверждать, что знаки Xi-2,..., Xi+i+s+1 образуют плотную F-рекуррентную серию длиной s, если знаки Xi-1,..., Xi+i+s образуют плотную F-рекуррентную цепочку наибольшей длины, лежащую внутри отрезка Xi-2,..., Xi+i+s+1. Началом серии будем называть знак Xi+i, а концом — знак Xi+i+s-1. Плотной F-рекуррентной серии длиной s, образованной знаками Xi-2,... ,Xi+i+s+1, соответствует плотная 1-серия длиной s, сформированная случайными индикаторами I{Ei+i-2}, I{Ei+i-1},..., I{Ei+i+s+1}.
Пусть событие означает, что знаки Х4—2, Х2,..., ХТ,... образуют плотную F-рекуррентную серию длиной s; для любого события E при c е {0,1}
Ж c = 1;
EC = \E, c = 0. Тогда событие Bi s можно записать как
ТЭ _ II рО 771Ü Т7>1 тра 1 Zpa2 гра.з_2
Bi,s = U Ei+i-2Ei+1-1Ei+1Ei+1+1Ei+1+2 ...Ei+1+s-2X
а1,...,а3_2
хЕ1+г+8—1Ег0+г+8 Ег0+г+8+1, (1)
где объединение ведется по всем числам a1,..., а8—2 е {0,1}, которые образуют плотную цепочку из единиц, т.е. 1 е {а», а^}, i = 1,..., s—3.
Пусть р8 = Р{В4,8}. Определим вероятность появления плотной F-рекуррентной цепочки длиной s — 2 (объединение ведется по тем же наборам, что и в (1)):
^5-2 = Р < У 1+2 . . . ^+1+8-2 (.
Пусть
8+Г
Q8,r ^ ^ q8/ . 8' = 8
Введем событие В^, состоящее в том, что в момент времени t в последовательности X1, X2,..., ХТ,... началась плотная F-рекуррентная серия длиной не меньше s. Тогда
те те
Ps = Р{£М} = ^ Р{В*,8+к} = ^ Ps+k .
к=0 к=0
Т
Пусть ^8>Т = I{В48} — число плотных F-рекуррентных серий
«2 rpas_2
t=1 T
длиной s и £8,т = I{В48} — число плотных F-рекуррентных серий
4=1
длиной не меньше s.
В соответствии с определениями Л8 = Е^8,т = Тр8, Л8 = Е<^8 Т = Тр8 Пусть
" = (бо,Т, ?80+1,Т, . . . , Сзо +г-1,Т, <Л8о+г,^ ;
П = (П80 , П80 + 1, . . . , П80+Г— 1, П80 +Г ),
где п80 ,п80+1,..., п80+г—1, П80+г — независимые случайные величины, каждая из которых распределена по закону Пуассона с параметром
Ав0, Ав0+1,..., Ав0 +г-1, Ав0+Г соответственно, т.е. при к = 0,1, 2,.. .:
дп
Р{п«,+к = п} = е-А^0+к, п = 0,1, 2,... п!
Распределение вектора П принято называть сопровождающим пуассо-новским распределением для вектора
Через р(Х, У) обозначим расстояние по вариации между распределениями случайных величин X и У [12]. Для случайных величин X и У, принимающих значения из счетного множества С, оно определяется по формуле
Р(Х, У ) = 2 Е |Р{Х = с}- Р{У = с}|.
сев
Теорема 1. При 1 < I < в0, г > 1, N > 2
р(2, П) = 2Трзд ((/ + во + г + 2, 5)А0 + (1 + 2^-^). (2)
Замечание 1. Теорема 1 содержит как частный случай теорему, приведенную в работе [7], в которой была получена оценка для вектора из чисел сплошных /-рекуррентных серий. Оценка, данная в работе [7], незначительно отличается от оценки (2) вследствие некоторых различий в доказательстве.
Замечание 2. Частный случай теоремы 1 — оценка, приведенная в работе [2] (формула (8)) для плотных а-серий. Оценка (2) несколько лучше. Это объясняется тем, что в работе [2] метод Чена- Стейна применялся для вывода оценки для вектора (£в0,т,^в0+1,т,..., £в0+г-1,т), а затем из нее получали оценку для вектора 2 = (£в0,т , £в0+1,т,... ..., &ю+г-1,т, 60+г,т). В настоящей работе метод Чена-Стейна использован для вектора 2.
Обозначим через Пв,т число плотных ^-рекуррентных цепочек длиной в, которые начались до момента Т.
Теорема 2. Пусть 1 < I < в0, г > 1, N > 2, и Т, в0 ^ то так, что
А30^ Ак е (0, то), к = 0, 1,. . . ,г - 1; (3)
А^+г ^ Аг е (0, то), ^ 0. (4)
Тогда
1) случайные величины £в0,т,^в0+1,т,... , £в0+г-1,т,<Ав0 +г,т асимптотически независимы в совокупности и распределены в пределе по закону Пуассона с параметрами А0, А1,..., Аг-1, Аг соответственно;
2) распределение случайной величины Пв0,т совпадает в пределе с
распределением выражения + 1)п, где п0,..., пт,... — неза-
^=0
висимые случайные величины, распределенные по закону Пуассона с параметрами А-2, А-1, А0,..., Ат,... соответственно.
Доказательства. Доказательство теоремы 1. Пусть
= ^^ 8 = Зсь . . . , 80 + г — l, ^Мо+г = В4,в0+г;
и = {(г, 8) : г = 1, 2,..., 8 = 80,..., 80 + г};
^ = {/М}(М)еи; п = {пМ}(м)еи, и случайные величины п^ независимы и распределены по закону Пуассона с параметрами Еп^ = Р{/г,в}.
Определим окрестности при 80 < 8 < 80 + г — 1:
0(М) = ШУ):
80 < 8' < 80 + г — 1, шах{1,г — I — 8' — 3} < 3 < г + I + 8 + 3}и и {(3, 80 + г) : шах{1, г — I — 80 — г — 1} < 3 < г + I + 8 + 3};
0(г, 80 + г) = {(3, 8') : 80 < 8' < 80 + г — 1, шах{1, г — I — 8' — 3} <
< 3 < г + I + 80 + г + 1} и {(3, 80 + г) :
шах{1, г — I — 80 — г — 1} < 3 < г + I + 80 + г + 1}.
Построенные окрестности обладают тем свойством, что случайный индикатор /г,5 не зависит от набора /¿^, (3,8') € 0(г,8). С учетом этого многомерный вариант метода Чена-Стейна (теорема 10.А, см. работу [8]) дает оценку расстояния по вариации между наборами W и п следующего вида:
Р(^П) < 5 + Я
где
5 = £ £ Р{/м}Р{4в,}; (5)
(г,в)еи 0>')еО(г,в)
5 = £ £ }. (6)
(г,5)еи (3,8')еО(»,8)\{(.,5)|
Теперь необходимо оценить сверху суммы 51 и 52. Начнем с оценки суммы 51, найденной по (5). С учетом определения окрестностей 0(г, 8) имеем
Т 50+Г
51 = ££ £ р{/,Лр{/^} =
¿=1 5=50 0т>')еО(г,в) Т / 50+Г-1 \
= £ ( £ £ Р{/,5}Р{/^} + £ Р{/^,50+Г }Р{/*-'} I <
¿=1 \ в=в0 (;',в')ёО(»,в) 0>')ёО(г,в0+г) у
Т /в0+г-1 /в0+г-1 г+1+в+Э
< £ £ £ £ р{в,лр{ }+
г=1 \ в=в0 \ в'=в0 ,?=г—1—в'—Э
г+г+в+3 \
+ Е р{^М }р{В,зд+г }) +
,7=г—1—в0 —г—1 /
в0+г-1 г+1+в0+г+1 г+1+в0+г+1 \
+ Е Е Р{Вг,50+г}Р{В>'} + Е Р{В^0+г}Р{В>0+г}1
в'=в0 ,?=г—1—в' —3 ^'=г—1—в0—г—1 /
Т /в0+г—1 /в0+г—1 г+1+в+3 г+1+в+З \
= Е Е Е Е рвРв'+р50+г Е р*) +
г=1 \ в=в0 \ в'=в0 ^'=г—1—в'—3 ^'=г—1—в0—г—1 /
в0+г—1 г+1+в0+г+1 г+1+в0 +г+1
+р5в0+г Е Е рв' + Е (:рв0+г)2
в'=в0 ^'=г—1—в' —3 ^'=г—1—в0—г—1
Во всех внутренних суммах слагаемые не зависят от индекса тогда
''so+r-1/ so+r-1
\ s=so \ s'=so so+r-1
S1<W Е Е (21+s+s'+7)(psPs') + (2/+s+S0+r+5)psPso+r +
+ppso +r E Ps' (2/ + so + r + s' + 5) + (2(1 + So + r) + 3)(p>so+r)2 I <
''so+r-1 so+r-1 so+r-1
< T(2(/ + so+r-1) + 7) E E PsPs' +2j9so+r E Ps + P2o+r <
s=so s'=so s=so
\ 2 N
< 2T(/ + so + r + 2, 5) ((rpso) + 2rpsoPso+r + (Pso+r) ) =
/so+r-1 \ 2
= 2T(/ + so + r + 2, 5)1 Е Pso + Pso+rl =
s=so
= 2T(/ + so + r + 2, 5)p2o. (7)
Перейдем к оценке суммы S2, определяемой по (6). Отметим, что
S2 < 2 Е Е P{Ii,sIj,s'} =
(i,s)eu (j,s')eo(i,s)\{(i,s)}, j>i
T so+r
= 2 ЕЕ Е рил'} =
i=1 s=so (j',s')eO(i,s), j>i
T so+r-1
^E iE E P{Ii,sIj,s'} + E P{Ii
,so+r^s'} I <
i=1 У s=so (j,s')eO(i,s), j>i (j,s')eO(i,so+r),i>i
T /so+r-1 /so+r-1 i+i+s+3
<2 Е Е Е Е m^'}+
i=1 \ s=so \ s'=so j=i+1
г+г+в+3 \
}+
^=¿+1 / в0+г—1 г+1+в0+г+1 г+1+в0+г+1 \
+ £ £ Р{Вг
'} + £ р{В
в'=в0 ^=¿+1 ^=¿+1 /
При (3,8') € 0(г,8) = 0(г, 8), 80 < 8, 8' < 80 + г — 1, 3 > г, события BiíS и В^' совместны, если
3 = г + 8 + 2, г + 8 + 3,..., г + 8 + I + 3
и произведение событий Bг,SBj,S' означает, что знаки Х^2, Х^, ... ..., Х^+5+1 и 2, 1,..., Х^+5'+1 образуют плотные Е-рекуррентные цепочки. Используя обозначения, введенные ранее, запишем
Г- R I I I cA-j+!+s+2 ElÖs'-2
С l U Ei+21+s+2 . . . Ej+Z+s'-2Ej'+i+s'-1
, bi-j+!+3+2>...>b3'-2
С
с Вг,в I У ^«-+21+5+2 . . . -^+5+2+1+5'-1 \«1,...,а3'-г
где 1 € {6г, 6г+1}, г = 1, 2,..., 1 € {аг, аг+1}, г = 1, 2,...
Поскольку событие В^ определяется случайными величинами Хг — 2, . . . , Xг+s+1, а произведение событий Eгa+2г+S+2 ... —¿+5+^+^ 1 — случайными величинами Х^^^,... , Х^+^+г+.з', в силу однородности последовательности Х1,... ,ХТ,... можем написать оценку:
Р{ВМВ.,>'} < Р{В1,5}Р ^ У E¿+21+s+2 . . . Eгa+'s^-2+г+s' —1
а1,...,а3'_г
= раР | и -а1 ...—'—71= р^'—г,
где объединение ведется по всем числам а1,...,а5'—г € {0,1}, 1 € € {а», а^+1}, г = 1, 2,..., 8' — I — 1.
Аналогично можно показать, что при 3 = г + 8 + 2, г + 8 + 3, . . . , г +
+ 8 + I + 3
,50 +г } < р5?50+Г —
P{Bi,so+rBj,s'} < Pso+rqs'-i; P{Bi,so+rBj,so +r} < pso+rqso+r-i.
Тогда
T / so+r-1 / so+r-1 i+i+s+3 i+i+s+3
S2 < 2 E ( ( Ps?s'-i + E Ps^o+r-i ) +
i=1 \ s=so \ s'=so j=i+s+2 j=i+s+2
so+r-1 i+i+so+r+1 i+i+so+r+1
+ .Pso+r qs'-i pso+r qso+r-i |
s'=so j=i+s+2 j=i+s+2
(so+r-1 /so+r-1 \
E ( E Psqs'-H + Psqso+r-i +
s=so s'=so
so+r-1 \
+ ^^ i^so+rqs'-i + Pso+rqso+r-i I =
s'=so
so+r-1
= 2T(/ + 2)Pso Е qs'-i = 2T(/ + 2)PsoQso-i. (8)
s'=so
Из формул (7) и (8) получим
p(W, п) < 2T(/ + so + r + 2, 5)P2o + 2T(/ + 2)PsoQso-i =
= 2Tj5so ((/ + so + r + 2, 5)Pso + (/ + 2)Qso-i).
Следует отметить, что случайные векторы 2 и П получаются из наборов W и п в результате применения к ним одной и той же функции, поэтому
р(2, П) < p(W, п).
Теорема 1 доказана.
Доказательство теоремы 2. Условия (3) и (4) эквивалентны тому, что правая часть оценки (2) стремится к нулю в условиях теоремы. Это доказывает п. 1.
Согласно п. 1, при любом натуральном значении r распределение
r— 1
случайной величины ^ (j + 1)^so+j-2T сходится к распределению
J=o
r-1
случайной величины (j + 1)п.
k=o
При r > 1
P {nso,T = "Т (j + 1)?so+j-2,4 = P( Е ?so+j-2,T > Л = l j=o ) l j=r J
= P{fso +r-2,T > 1} < Efso+r-2,T. (9)
В силу условия (4) ^
Е6о+г,т = £ Е^о+г+к ^ Лг,
к=0
тогда правая часть (9) стремится к нулю при г ^ то. Теорема 2 доказана.
Оценки для равновероятного алфавита. Обратимся к вычислению вероятностей События Е;+1, Е+2, Е+з,... зависимы. Поэтому выписывание вероятности р в явном виде в общем случае невозможно.
Далее ограничимся случаем, когда знаки Х1,Х2,... распределены на алфавите равновероятно.
Пронумеруем функции в наборе Е = {/ь ..., /к}. Пусть при фиксированном значении ] = 1,..., К:
= {X = / (Х-1,...,Х-1)}.
Тогда
K
E = U .
j=1
Пусть для набора индексов i = (z_i, ¿о, ¿i, ¿2,..., is+2), 1 < ij < K, событие
/"Y _ Il 77Ю ipO ipl 77rai
Ci,i = ^J Ei-2,i_i Ei-MoEMiEi+1,i2 . . . ai,...,as_2
Eas-2 E1 E 0 EO
. . . Ei+s-1,is Ei+s-1,is , Ei+s,is + 1 EÎ+S+Ms + 2 ,
где объединение ведется по всем числам
а1,..., as-2 G {0,1}, 1 G {а», ai+1}, i = 1, 2,..., s — 3. Можно записать оценку
Ps = Pj U С1Д < £ P{C1,i},
которая преобразуется в точное равенство, если при 1 < i < j < K
/г(An) П / (AN) = 0,
так как в этом случае события C1i, i G {1,..., K}s+4, попарно несовместны.
В силу однородности последовательности X1, X2,...
,il EÎ+1,Î2 . . . Ei+s-1,is
} = P{E1+1,iiE1+2,i2 . . . E+s,is } =
= P {Xi+1 = fii(хъ... ,X),X+2 = (X2,... ^rn^ ...,
. . . ,X+s = /is (Xs-b . . . ,X1+s-1)} =
= Е - Е = С1> ... Р{Х' = С'}х
хР{Хш = /1(с1,...,сг),Х/+2 =/2(02,..., С',/1(01 ,...,С')),...} =
= N'4- =
N'+5 N5
Это означает, что события ^ ... независимы в сово-
купности [1].
Тогда с учетом свойства, описанного в замечании 1, и результата, полученного в работе [2], имеем
Р{ад = с(д5 - д®); (10)
c = ; q =2 (p + - 3p)) ; qi = 2 (p " " 3p)) где p = N-1. Поскольку
P
c <
- 3р) з'
то
К 4
р < К5+4Р{СМ} < ((Кд)5 - (*?!)-); (11)
< ((Кд)5-2 - (ЗД5-2). (12)
Оценки (11) и (12) для вероятностей содержательны при Кд < 1. При вычислении вероятностей для заданного набора функций можно получить другие оценки. Однако они существенно зависят от этого набора, поэтому их выписывание в общем виде является очень трудоемкой задачей.
Подставляя оценки (11) и (12) в (8), получаем следующую оценку:
8ТК 4( Кд)2в°-'—2
Р(5, П)<87 К -((/+вс+г+2, 5)К4(Кд)'+2 + I + 2) . (13)
Следствие 1. Пусть случайные величины Х1, Х2,..., Хт,... независимы и распределены на алфавите равновероятно. Пусть 1 < I < в0, г > 1, N > 2, дК < 1, и Т, в0 ^ то так, что правая часть (13) стремится к нулю,
Азо+к ^ А^ е (0,то), к = 0,1,...,г - 1, А30+г ^ А е (0,то).
Тогда
1) случайные величины £в0,т, ^^.т,... , £в0+г-1,т, А0 +г,Т асимптотически независимы в совокупности и распределены в пределе по закону Пуассона с параметрами А0, А1,..., Аг-1, А соответственно;
2) распределение случайной величины Пв0,т совпадает в пределе с
те
распределением выражения (3 + 1)п, где п0,..., пт,... — неза-
^=0
висимые случайные величины, распределенные по закону Пуассона с параметрами Л_2, А_1, Л0,..., Лт,... соответственно.
Замечание 3. Можно показать, что условия следствия 1 выполнены, если, например, К =1 и
So =
ln T ln q
+ O(1), T —у то.
Так, пусть семейство Е состоит из одной функции / : А^ ^ А^ (К =1), N = 4, г = 2, Т = 100 000, 8о = 17. Тогда параметры сопровождающего пуассоновского распределения равны Ав0 = 0,735; Ав0+1 = 0,998. Из оценки (13) определим
-к
k!
Лг
P|6° = k, Ло+1 = l}"1^-fе-Лз°
при k, I = 0,1, 2,... В частности,
< 0,016
P{6° = 0, £s°+1 = 0} - 0,177 < 0,016
pie° = 0, £s°+1 = 1} - 0,176 < 0,016
pie° = 0, 65°+1 = 2} - 0,088 < 0,016
pie° = 1, £s°+1 = 0} - 0,130 < 0,016
pie° = 1, 6°+1 = 1} - 0,130 < 0,016
p{e° = 1, 6°+1 = 2} - 0,065 < 0,016
P|6° = 2, 6°+1 = 0} - 0,048 < 0,016
P|6° = 2, 6°+1 = 1} - 0,048 < 0,016
Следствие 2. Пусть случайные величины Х!, Х2,..., Хт,... независимы и распределены на алфавите А^ равновероятно. Пусть 1 < I < в0, г > 1, N > 2, дК < 1, и Т, 80 ^ то так, что правая часть (13) стремится к нулю,
А^+к ^ то, к = 0,1,...,г - 1, Л30+г ^ то.
Тогда случайные величины
бо.Т - Аао <^0+Г-1,Т - А50+Г— 1 <ъ50+г,т - А50+г
Аво Аво+г-1 Ав0+г
асимптотически независимы в совокупности и одинаково распределены в пределе по стандартному нормальному закону.
Замечание 4. Условия следствия 2 выполнены, если, например, К =1 и
1п Т — 1п 1п Т
So = -
+ O(l), T ^то.
1п д
Заключение. Обратимся к схеме плотного вложения [13, 14]. Пусть знаки х,..., х+'+5-1 образуют /-рекуррентную цепочку. Пусть у1 , у2,... — последовательность знаков алфавита . Согласно работе [14], знаки х,..., х+'+5-1 плотно вкладываются в последовательность у1,у2,..., начиная с места > 1, если существуют такие числа
< Л < ^2 < ... < 1, .ь - е {1,2}, и =1,2,..., в +/ - 1,
что = жг+«, и = 0,1, 2,..., в + / - 1.
Это означает, что
= / ,...,Ул+к-1), к = 0,1, 2,в - 1. (14)
Рассмотрим набор Е функций д вида д(уь..., у21-1)=/,..., ), где 1 < ^ < ^2 < ... < ^ < 2/ - 1, ^ - ,7«-1 е {1, 2}, и = ¡2,3,..., 1. Тогда равенство (14) перепишется в виде
уг+к = д(%г+к-2'+1,...,ул+к-l), к = 0,1,2, в - 1 д е е.
Таким образом, знаки у0,..., образуют плотную Е-рекуррентную цепочку.
Задача о свойствах распределений чисел плотных Е-рекуррентных серий и цепочек рассматривается впервые. Она представляет интерес в связи с возможным применением этих статистик при проверке гипотез о свойствах случайных последовательностей. Например, для проверки основной гипотезы о том, что наблюдается последовательность независимых случайных величин с одинаковым равновероятным полиномиальным распределением, против альтернативы о том, что наблюдаемая последовательность получена как результат плотного вложения отрезков двух рекуррентных последовательностей.
В настоящей работе с помощью функционального варианта метода Чена- Стейна получена многомерная предельная теорема Пуассона для вектора из чисел плотных Е-рекуррентных серий заданных длин с оценками скорости сходимости (в метрике расстояния по вариации). Из этих оценок также выведена теорема пуассоновского типа для числа плотных Е-рекуррентных цепочек.
Работа выполнена при поддержке РФФИ (проект 14-01-00318А). Автор выражает глубокую признательность рецензенту за сделанные замечания и внимание к статье.
ЛИТЕРАТУРА
1. Михайлов В.Г. О предельной теореме Б.А. Севастьянова для сумм зависимых случайных индикаторов // Обозрение прикладной и промышленной математики. 2003. Т. 10. Вып. 3. С. 571-578.
2. Меженная Н.М.Предельные теоремы для числа плотных серий в случайной последовательности //Дискретная математика. 2009. Т. 21. Вып. 1. С. 105-116.
3. Меженная НЖПредельная теорема Пуассона для числа плотных серий заданной длины и веса // Вестник МГТУ им. Н.Э. Баумана. Сер. Естественные науки. Спец. вып. 2011. С. 75-82.
4. Меженная Н.М.Предельные теоремы для числа (а, й)-серий заданного веса в последовательности независимых случайных величин // Вестник МГТУ им. Н.Э. Баумана. Сер. Естественные науки. Спец. вып. 2012. С. 20-28.
5. Севастьянов Б.А. Предельный закон Пуассона в схеме сумм зависимых случайных величин // Теория вероятностей и ее применения. 1972. Т. 17. Вып. 4. С. 733-738.
6. Колчин В.Ф., Севастьянов Б.А., Чистяков В.П. Случайные размещения. М.: Наука, 1976. 224 с.
7. Михайлов В.Г. Об асимптотических свойствах числа серий событий // Труды по дискретной математике. 2006. Т. 9. С. 152-163.
8. Barbour A.D., Holst L., Janson S. Poisson approximation. Oxford: Oxford Univ. Press, 1992. 278 p.
9. Lecture note series. Institute of mathematical sciences, National institute of Singapore. In 5 vol. Vol. 4. An introduction to Stein's method / Barbour A.D., Chen L.H.Y. Singapore: Singapore University Press, 2005. 226 p.
10. Lecture note series. Institute of mathematical sciences, National institute of Singapore. In 5 volumes. Vol. 5. Stein's method and applications / Barbour A.D., Chen L.H.Y. Singapore: Singapore University Press, 2005, 297 pp.
11. Михайлов В.Г. Явные оценки в предельных теоремах для сумм случайных индикаторов // Обозрение прикладной и промышленной математики. Т. 1. Вып. 4. С. 580-617.
12. Ширяев А.Н. Вероятность. М.: Наука, 1980. 576 с.
13. Golic J.Dj. Constrained embedding probability for two binary strings // SIAM J. on Discrete Math. Vol. 9. No. 3. 1996. P. 360-364.
14. Михайлов В.Г., Меженная Н.М.Оценки для вероятности плотного вложения одной дискретной последовательности в другую // Дискретная математика. 2005. Т. 17. Вып. 3. C. 19-27.
REFERENCES
[1] Mikhailov V.G. On Sevast'yanov limit theorem for sums of dependent random indicators. Obozrenie prikladnoy i promishllennoy matematiki [OP&PM Surveys in Applied and Industrial Mathematics], 2003, vol. 10, no. 3, pp. 571-578 (in Russ.).
[2] Mezhennaya N.M. Limit theorems for the number of dense series in a random sequence. Diskretnaya matematika [Discrete Mathematics and Applications], 2009, vol. 19, no. 2, pp. 215-228 (in Russ.). DOI: 10.1515/DMA.2009.012
[3] Mezhennaya N.M. Poisson limit theorem for a number of dense series of the given length and weight. Vestn. Mosk. Gos. Tekh. Univ. im. N.E. Baumana, Estestv. Nauki., Spetsvyp. [Herald of the Bauman Moscow State Tech. Univ., Nat. Sci., Spec. Issue], 2011, pp. 75-82 (in Russ.).
[4] Mezhennaya N.M. Limit theorems for the number of (a, d)-series of given weight in a Sequence of Independent Random Quantities. Inzhenernyy zhurnal: nauka i innovatsii. EI. n.-t. izdanie [Engineering Journal: Science and Innovations], 2012, vol. 4, no. 4, p. 20-28 (in Russ.).
[5] Sevast'yanov B.A. Poisson limit law for a scheme of sums of dependent random variables. Teoriya veroyatnostey i ee primeneniya [Theory of Probability and its Applications, 1972, vol. 17, no. 4, pp. 695-699], 1972, vol. 17, no. 4, pp. 733-738 (in Russ.). DOI: 10.1137/1117082
[6] Kolchin V.F., Sevastyanov B.A., Chistyakov V.P. Sluchainye razmescheniya. Teoriya veroyatnostei i matematicheskaya statistika [Random placements. Probability theory and mathematical stattistics]. Moscow, Nauka Publ., 1976. 223 p.
[7] Mikhailov V.G. On asymptotic properties of numbers of event series. Trudy po Diskretnoi Matematike [Proceedings on discrete mathematics], 2006, vol. 9. pp. 152163 (in Russ.).
[8] Barbour A.D., Holst L., Janson S. Poisson Approximation. Oxford: Oxford Univ. Press, 1992. 278 p.
[9] Barbour A.D., Chen L.H.Y., eds. Lecture note series. Institute of mathematical sciences, National institute of Singapore. In 5 volumes. Vol. 4. An introduction to Stein's method. Singapore: Singapore University Press. 2005. 226 p.
[10] Barbour A.D., Chen L.H.Y., eds. Lecture note series. Institute of mathematical sciences, National institute of Singapore. In 5 volumes. Vol. 5. Stein's method and applications. Singapore: Singapore University Press, 2005. 297 p.
[11] Mikhailov V.G. Explicit estimators in limit theorems for sums of random indicators. Obozrenie prikladnoy i promishllennoy matematiki [OP&PM Surveys in Applied and Industrial Mathematics], 1994, vol. 1, no. 4, pp. 580-617 (in Russ.).
[12] Shiryaev A. N. Veroyatnost' [Probability]. Moscow, Nauka Publ., 1980. 576 p.
[13] Golic J.Dj. Constrained embedding probability for two binary strings. SIAM J. Discrete Math., 1996, vol. 9, no. 3, pp. 360-364. DOI: 10.1137/S0895479894246917
[14] Mikhailov V.G., Mezhennaya N.M. Bounds for the probability of a constrained embedding of one discrete sequence into another. Diskretnaya matematika [Discrete Mathematics and Applications], 2005, vol. 15, no. 4, pp. 377-386 (in Russ.). DOI: 10.1515/156939205774464864
Статья поступила в редакцию 25.09.2013
Наталья Михайловна Меженная — канд. физ.-мат. наук, доцент кафедры "Прикладная математика" МГТУ им. Н.Э. Баумана. Автор 10 научных работ в области дискретных задач теории вероятностей, предельных теорем и их применения в математической статистике.
МГТУ им. Н.Э. Баумана, Российская Федерация, 105005, Москва, 2-я Бауманская ул., д. 5.
N.M. Mezhennaya — Cand. Sci. (Phys.-Math.), assoc. professor of "Applied mathematics" department of Bauman Moscow State Technical University. Author of 10 publications in the field of discrete problems of probability theory, limiting theorems and their application in mathematical statistics.
Bauman Moscow State Technical University, Vtoraya Baumanskaya ul. 5, Moscow, 105005 Russian Federation.