2015 Псевдослучайные генераторы № 1(27)
УДК 519.4
О НЕКОТОРЫХ КЛАССАХ РАЗЛОЖИМЫХ ЦЕПЕЙ МАРКОВА НА КОНЕЧНОЙ АБЕЛЕВОЙ ГРУППЕ
М. И. Рожков
Московский институт электроники и математики Национального исследовательского университета Высшая Школа Экономики, г. Москва, Россия
Рассматривается задача разложения простой однородной цепи Маркова Г в сумму s ^ 2 взаимно независимых составляющих цепей Маркова r(i), заданных на конечной абелевой группе G. Данная задача связана с известной процедурой укрупнения состояний цепи Маркова. Описывается широкий класс цепей, допускающих такое разложение в сумму двух цепей бесконечным числом способов. Результаты данной работы могут быть полезны при оценке возможностей восстановления параметров исходных цепей Маркова, поступающих в узел суммирования, по результирующей последовательности. Такого сорта задачи возникают при построении и исследовании свойств генераторов случайных последовательностей.
Ключевые слова: разложение цепи Маркова, сумма цепей Маркова, укрупнение состояний цепи Маркова.
ON SOME CLASSES OF DECOMPOSABLE MARKOV CHAINS ON FINITE ABELIAN GROUP
M. I. Rozhkov
National Research University Higher School of Economics, Moscow, Russia E-mail: [email protected]
The problem for the decomposition of a given simple homogeneous Markov chain into the sum of the components of the chain is considered. The class of chains allowing this decomposition in the infinite number of ways is described.
Keywords: decomposition of a Markov chain, the summation of Markov chains.
1. Постановка задачи
Пусть r(i) = 71 (i), 72 (i), • • ., i = 1,... , s — независимые реализации простых однородных цепей Маркова на конечной абелевой группе G = (G, +) с соответствующими матрицами переходных вероятностей
n(i) = || p(i) (g,h) ||, g,h е G, i = 1,2, ...,s.
Положим Г = Yi, 72,..где
s
7j = Е 7j(i) j = 1 2,---i= 1
В работах [1,2] исследовались условия, при которых элементы результирующей последовательности Г также связаны простой однородной цепью Маркова с некоторой матрицей переходных вероятностей
п = ||p(g,h)||, g,h е G,
которая не зависит от начального распределения исходных цепей (т. е. от распределения величин 7!(г), г = 1, 2,... , в). В этом случае будем говорить, что цепь Г разложима в сумму цепей Маркова. Отметим, что рассматриваемая задача разложения цепи Маркова в сумму цепей тесным образом связана с известной процедурой укрупнения состояний цепи Маркова [3].
С каждой строкой матрицы п(г) = ||р(г)(д, К)||, д, К € С, свяжем элемент группового кольца ОС (определение группового кольца см. в [4])
р%)= Е Р(%,К) ■ К,
ьес
где О — поле действительных (рациональных) чисел. Элемент р(г) (д) будем называть характеристической функцией соответствующей строки матрицы п(г). В [1, теорема 1 и лемма 3] найдено следующее строение матрицы переходных вероятностей п = ||р(д, К) || цепи Г, являющейся суммой исходных цепей Г(г).
Теорема 1. Пусть сумма в ^ 2 цепей Маркова с матрицами переходных вероятностей п(г) = ||р(г)(д, К)||, д, К € С, г = 1, 2,... , в, является цепью Маркова с матрицей переходных вероятностей п = ||р(д, К) ||, д, К € С. Тогда для характеристических функций р(д) строк матрицы п в групповом кольце ОС выполняются соотношения
р(0) = р(1)(0) ■ р(2)(0) ■... ■ рм(0), р(д) = р(0)/(д), д € С, (1)
где / — гомоморфизм группы С в группу мультипликативно обратимых элементов кольца ОС; 0 — нулевой элемент группы С.
Вопрос о том, всегда ли цепь с матрицей п = ||р(д, К)||, характеристические функции строк которой удовлетворяют соотношениям (1), разлагается в сумму в цепей, является открытым. При этом известно следующее достаточное условие марковости для суммы цепей.
Теорема 2 [2, теорема3]. Пусть на конечной абелевой группе С заданы в ^ 2 цепей Маркова Г(г) с матрицами переходных вероятностей
п(г) = ||р(г) (д, К) ||, д, К € С, г = 1, 2,..., в.
Сумма цепей Маркова Г(г) также является цепью Маркова, если характеристические функции строк матриц п(г) удовлетворяют соотношениям
р%) = Р(г)(0)/(д), д € С, г = 1,2,...,в,
где / — гомоморфизм группы С в группу мультипликативно обратимых элементов кольца ОС.
В дальнейшем потребуются два вспомогательных утверждения. Элемент
р = Е р(К) ■К
ьес
группового кольца ОС будем называть вероятностным, если У] р(К) = 1, р(К) ^ 0.
ьес
Лемма 1. Пусть С — конечная абелева группа, О — поле действительных (рациональных) чисел, /(д) = Е /(д, К) ■ К — произвольный гомоморфизм группы С в группе
пу мультипликативно обратимых элементов кольца ОС, такой, что Е /(д, К) = 1 для
ьес
любого элемента g Е G, M — максимальное значение абсолютных величин коэффициентов f (g, h), g, h Е G.
Тогда для любого вероятностного элемента вида
Р = Е h + Е e(h)h, |e(h)| <-!-, Е £(h) = 0
|G| keG keG |G| M keG
произведение p • f (g) также является вероятностным элементом при любом g Е G. Доказательство. Заметим, что
p •f (g) = А Е h + Е Е (e(r)f (g,h - r))h_
|G| keG keGreG Кроме того, справедливо неравенство
| Е e(r)f (g, h - r)| < ^^M • |G| = 1/|G|. reG |G|2M
Отсюда следует, что все коэффициенты элемента pf (g) положительны, и утверждение леммы 1 доказано. ■
Лемма 2. Пусть p = -— Е h + Е e(h)h, Е e(h) = 0,— произвольный вероят-
|G| heG keG keG
ностный элемент. Тогда справедливо разложение
p = L • R,
L = ^E h + E e • h + 5 • 0, e = 0, 5 = 0, e|G| + 5 = 0,
|G| keG keG
R =^E h +E r(h)h, r(h) = e(h)/5.
|G| keG keG
Справедливость леммы 2 проверяется непосредственными расчётами.
2. Основные результаты
Основным результатом настоящей работы является следующее утверждение. Теорема 3. Пусть цепь Маркова Г имеет матрицу переходных вероятностей п = = ||p(g,h)||, для характеристических функций строк которой выполнены следующие условия:
1) p(0) = — Е h + Е e(h)h, |e(h)| < (1/|G|)4(1/M)2, £ e(h) = 0;
|G| keG keG keG
2) p(g)= p(0)f(g), g Е G,
где f (g) = E f (g, h)h — заданный гомоморфизм группы G в группу мультипликатив-keG
но обратимых элементов кольца DG; M — максимальное значение абсолютных величин коэффициентов f (g, h), g, h Е G.
Тогда существует разложение p(0) = p(1)(0)p(2)(0), при котором цепь Г представляется суммой цепей Г = Г(1) + Г(2) c матрицами переходных вероятностей n(i) = = ||p(i)(g, h)||, g, h Е G, i = 1, 2, характеристические функции строк которых имеют вид
p(i)(g)= p(i)(0)f (g), g Е G, i =1, 2.
Доказательство. В соответствии с леммой 2 для характеристической функции p(0) справедливо разложение
p(0) = p(1) (0)p(2) (0) = L • R,
L = ^E h + E e • h + 5 • 0, e = 0, 5 = 0, e|G| + 5 = 0, |G| нес hec
R = ^E h + E r(h)h, r(h) = e(h)/5. |G| hec hec
Положим в данном разложении
|e| = (1/|G|)3(1/M). (2)
Тогда
|5| = |e| • |G| = (1/|G|)2(1/M), |e + 5| = |e|(|G| - 1) < (1/|G|)2(1/M); (3) |r(h)| = |e(h)/5| < (1/|G|)4(1/M)25-1 = (1/|G|)2(1/M). (4)
Из соотношений (2)-(4) следует, что при указанном выборе параметра e каждая из функций p(j)(0), j = 1, 2, удовлетворяет условиям леммы 1. Тем самым при любом g £ G элемент p(j)(0)f (g) является вероятностным. Теперь справедливость теоремы вытекает из теоремы 2. ■
Следующее утверждение показывает, что разложение Г = Г(1) + Г(2) в условиях теоремы 3 является неоднозначным.
Теорема 4. Пусть выполнены условия теоремы 3. Тогда существует бесконечное множество разложений p(0) = p(1)(0)p(2)(0), при каждом из которых цепь Г представляется суммой цепей Г = Г(1) + Г(2) c матрицами переходных вероятностей
n(i) = ||p(i)(g, h)||, g,h £ G, i = 1, 2,
характеристические функции строк которых имеют вид
p(i)(g)= P(i)(0)f (g), g £ G, i =1, 2. Доказательство. Пусть в представлении
p(0) = G E h + E e(h)h, |e(h)| < (1/|G|)4(1/M)2, E e(h) = 0 |G| hec hec hec
имеет место 0 < в < max{|e(h)|}/((1/|G|)4(1/M)2) < 1. Пусть а — произвольное фик-
hec
сированное число, в < а < 1. В соответствии с леммой 2 для характеристической функции p(0) справедливо разложение
p(0)= p(1) (0)p(2) (0) = L • R, (5)
L = ^E h + E e • h + 5 • 0, e = 0, 5 = 0, e|G| + 5 = 0,
|G| hec hec
R =^E h + E r(h)h, r(h) = e(h)/5.
| G| hec hec
Положим в разложении (5)
и = (1/|G|)3(1/M )а. (6)
Тогда
н = И ■ |G| = (1/|G|)2(1/M)a, |е + i| = H(|G| - 1) < (1/|G|)2(1/M)a; (7) |r(h)| = |е(Ь)/*| < e(1/|G|)4(1/M)2^-1 = (e/a)(1/|G|)2(1/M). (8)
Из соотношений (6)-(8) следует, что при указанном выборе параметра е каждая из функций p(j)(0), j = 1, 2, удовлетворяет условиям леммы 1. Тем самым при любом g £ G элемент p(j)(0)f (g) является вероятностным. Следовательно, учитывая теорему 2, исходная цепь Г является суммой цепей Г = Г(1) + Г(2) c матрицами переходных вероятностей
n(i) = ||p(i)(g,h)||, g,h £ G, i = 1, 2, характеристические функции строк которых имеют вид
p(i)(g) = P(i)(0)f (g), g £ G, i =1, 2.
При этом различные параметры а (в < а < 1) отвечают разным матрицам n(i) переходных вероятностей. ■
3. Заключительные замечания
Отметим, что из разложения p(0) = p(1)(0)p(2)(0) и вероятностности многочленов p(0)f (g), g £ G, не всегда «автоматически» следует вероятностность всех многочленов p(j)(0)f (g), j = 1, 2. Об этом свидетельствует следующий пример. Рассмотрим группу G = Zm = {0,1, 2, 3} вычетов по модулю m = 4. В этом случае кольцо DG есть кольцо многочленов D[x]/(x4 — 1) по модулю многочлена x4 — 1. И соответственно вероятностные элементы кольца DG принимают вид многочленов
p = p(h)h = p(j)xj = p(0) + p(1)x + p(2)x2 + p(3)x3. hea jezm
Рассмотрим вероятностные многочлены
p = p(x) = p(1) (x)p(2) (x),
p(1)(x) = p(2)(x) = (1/2 + е) + (1/2 — e)x, 0 < е < 1/2(>/2 — 1),
а также гомоморфизм f группы Zm в группу мультипликативно обратимых элементов кольца D[x]/(x4 — 1) :
f (m) = вт, в =1/2(1+ x + x2 + x3) — 1, m £ {0,1, 2, 3}.
Так как p(1)(x)e = —е + ex + 1/2 ■ x2 + 1/2 ■ x3, многочлен p(1)(x)f(1) не является вероятностным. В то же время любой многочлен p(x)f (m) является вероятностным в силу того, что 1) в2 = 1;
2) все коэффициенты многочлена p(x) = (1/2 + е)2 + (1/2 — 2e2)x + (1/2 — e)2x2 меньше 1/2;
3) p(x)e = 1/2(1 + x + x2 + x3) — p(x).
Результаты данной работы являются развитием аналогичного сорта результатов относительно разложимых законов распределения на конечной абелевой группе (см., например, [5]) и могут быть полезны при оценке возможностей восстановления параметров исходных цепей Маркова, поступающих в узел суммирования, по результирующей последовательности. Такого сорта задачи возникают при построении и исследовании свойств генераторов случайных последовательностей.
ЛИТЕРАТУРА
1. Рожков М. И. О суммировании цепей Маркова на конечной группе // Труды по дискретной математике. 2000. Т.3. С. 195-214.
2. Рожков М. И. Суммирование марковских последовательностей на конечной абелевой группе // Дискретная математика. 2010. Т. 22. №3. С. 44-62.
3. КемениДж., Снелл Дж. Конечные цепи Маркова. М.: Наука, 1970.
4. Кэртис Ч., Райнер И. Теория представлений конечных групп и ассоциативных алгебр. М.: Наука, 1969. 688 с.
5. Воробьев Н. Н. Сложение независимых случайных величин на конечных абелевых группах // Математический сб. 1954. Т. 34(76). №1. С. 89-126.
REFERENCES
1. Rozhkov M. I. O summirovanii cepej Markova na konechnoj gruppe. Trudy po Diskretnoi Matematike, 2000, vol.3, pp. 195-214. (in Russian)
2. Rozhkov M. I. Summirovanie markovskih posledovatel'nostej na konechnoj abelevoj gruppe. Diskretnaya Matematika, 2010, vol. 22, no. 3, pp. 44-62. (in Russian)
3. Kemeni Dzh., Snell Dzh. Konechnye cepi Markova. Moscow, Nauka Publ., 1970. (in Russian)
4. Kjertis Ch., Rajner I. Teorija predstavlenij konechnyh grupp i associativnyh algebr. Moscow, Nauka Publ., 1969. 688 p. (in Russian)
5. Vorob'ev N. N. Slozhenie nezavisimyh sluchajnyh velichin na konechnyh abelevyh gruppah. Matematicheskii Sb., 1954, vol. 34(76), no. 1, pp. 89-126. (in Russian)