CC BY
26
УДК 519.725
DOI: 10.17586/0021-3454-2020-63-5-405-416
ФОРМИРОВАНИЕ ТРОИЧНЫХ ПОСЛЕДОВАТЕЛЬНОСТЕЙ С ВЫСОКОЙ СТРУКТУРНОЙ СКРЫТНОСТЬЮ В СИСТЕМАХ ПЕРЕДАЧИ ЦИФРОВОЙ ИНФОРМАЦИИ
В. Г. Стародубцев, В. В. Ткаченко, Е. А. Боброва
Военно-космическая академия им. А. Ф. Можайского, 197198, Санкт-Петербург, Россия
E-mail: vka@mil.ru
Представлен алгоритм определения начальных состояний регистров сдвига, входящих в устройство формирования троичных последовательностей Гордона — Миллса — Велча (ГМВ) с периодом N=728. Алгоритм основан на сравнении начальных состояний, полученных в результате решения в конечных полях системы линейных уравнений, и состояний, определенных путем децимации символов базисной М-последовательности. Троичные М-последовательности и ГМВ-последовательности обладают одинаковой двухуровневой периодической автокорреляционной функцией, но различной структурной скрытностью, характеризуемой эквивалентной линейной сложностью. ГМВ-последовательность формируется на основе базисной М-последовательности с аналогичным периодом при ее представлении в виде квазиквадратной матрицы. Показано, что для каждого из 48 примитивных полиномов в конечном поле GF(36) может быть сформировано по три ГМВ-последовательности. Для двоичных ГМВ-последовательностей начальные состояния регистров сдвига образуются путем децимации символов базисной М-последовательности, представленной в каноническом виде; для троичных ГМВ-последовательностей отдельные суммируемые составляющие имеют дополнительный сдвиг на полпериода базисной М-последовательности. Полученные результаты могут найти применение при формировании широкополосных недвоичных сигналов в системах передачи цифровой информации.
Ключевые слова: псевдослучайные последовательности, конечные поля, неприводимые и примитивные полиномы, структурная скрытность, децимация, регистры сдвига
Современные системы передачи цифровой информации (СПЦИ) характеризуются сложной иерархической структурой, территориальной распределенностью, наличием как проводных, так и радиоканалов связи, включающих, в том числе, каналы спутниковой связи, которые могут функционировать в условиях радиоэлектронного противодействия [1, 2]. В таких радиоканалах широкое применение получили сигналы с расширенным спектром (СРС) [3, 4], позволяющие минимизировать воздействие узкополосных и широкополосных преднамеренных помех.
Одним из показателей, характеризующих помехозащищенность систем передачи цифровой информации, является структурная скрытность. В настоящее время в качестве псевдослучайных последовательностей (ПСП), используемых для расширения спектра сигналов, в основном применяются двоичные М-последовательности, последовательности Голда, последовательности малого и большого множеств Касами, а также последовательности Гордона — Миллса — Велча (ГМВ) [5, 6].
Перспективным направлением развития СПЦИ является переход к недвоичным сигналам, формируемым на основе недвоичных ПСП. Среди последовательностей, обладающих двухуровневой периодической автокорреляционной функцией (ПАКФ), можно выделить р-ичные М-последовательности (МП) и ГМВ-последовательности (ГМВП). При этом
в помехозащищенных СПЦИ целесообразно применять ГМВП, структурная скрытность которых, характеризуемая эквивалентной линейной сложностью (ЭЛС), существенно превышает данный показатель в МП [7].
Вопросам формирования и анализа корреляционных и структурных свойств двоичных и недвоичных ПСП посвящено большое количество публикаций [8—18]. Так, в работе [8] проанализированы корреляционные и структурные свойства как двоичных, так и троичных ПСП над полями нечетных характеристик, а также составных троичных последовательностей. Формированию и оценке структурных свойств двоичных ГМВП посвящена работа [9], а в [10] проведен анализ синтеза троичных ГМВП с периодом N=80 с учетом определения начальных состояний регистров сдвига. Процедура формирования фазоманипулированных сигналов с высокой структурной скрытностью рассмотрена в работе [11]. В [12] разработан алгоритм формирования пятеричных ГМВП с периодом N=624 и произведена оценка их линейной сложности. Подробный анализ вопросов формирования недвоичных ПСП и семейств ПСП с заданными корреляционными и структурными свойствами представлен в работах [13—15]. В работах [16—18] приведены результаты по синтезу семейств троичных и р-ичных последовательностей с низким уровнем взаимно корреляционных функций.
Алгоритм формирования троичных ГМВП с периодом N=728 в конечном поле GF(36) рассмотрен в работе [19], где получены проверочные полиномы Иг(х) для трех типов ГМВ-последовательностей с различной ЭЛС, которые представлены в виде произведения трех и девяти неприводимых полиномов-сомножителей ^сг(х). Приведена также структурная схема ГМВП с ЭЛС 4=18, состоящая из трех регистров сдвига с линейными обратными связями. Однако в данной работе не определены начальные состояния регистров сдвига для всех трех типов ГМВП, что не позволяет осуществить их непосредственное формирование.
Цель настоящей статьи — разработка алгоритма определения начальных состояний регистров сдвига, входящих в устройство формирования троичных ГМВП с периодом N=728.
Для двоичных ГМВП начальные состояния регистров сдвига, входящих в устройство формирования, определяются путем децимации символов, начиная с символа d0, базисной М-последовательности, представленной в каноническом виде, по индексам децимации, равным показателям степени корней неприводимых полиномов-сомножителей ^сг(х) проверочного полинома Иг(х) [20]. Данные корни являются элементами конечного поля GF(36). Канонический вид МП определяется через функцию следа из расширенного поля GF(36) в простом поле GF(3) di = tr61a', i = 0, 1, 2,..., N-1, т.е. d0 = tr61a0, d1 = tr^a1,..., где a — примитивный элемент поля GF(36) [5, 12].
13
Например, если корнями некоторого полинома Ис(х) являются элемент a и его различные р-сопряженные элементы, то начальное состояние регистра определяется символами базисной МП d0, d13, d26, d39, d52, d65.
Умножители и сумматоры по mod р (р=2, 3) в цепи обратной связи регистров сдвига расставляются в соответствии с коэффициентами полиномов Ис(х) как для двоичных, так и троичных последовательностей. Искомая ГМВП формируется путем суммирования последовательностей с выходов всех регистров сдвига. Если полином Ъс1(х) является примитивным, то на выходе данного регистра формируется МП. Если же полином И^(х) является только неприводимым, корни которого имеют период меньше максимального, то на выходе регистра формируется ПСП с данным периодом.
Особенность формирования троичных ГМВП с периодом N=728 заключается в том, что некоторые из суммируемых последовательностей начинаются не с символа d0, а имеют сдвиг на полпериода и начинаются с символа d364. В общем случае при формировании р-ичных ГМВ-последовательностей сдвиг начального состояния может принимать значения, кратные величине N/(p-1) [12]. Тогда для рассмотренного выше примера начальное состояние регист-
ра при наличии сдвига на полпериода определялось бы шестью символами базисной МП d364, d377, d390, d403, d416 и d429. Отметим, что номера символов вычисляются по mod 728.
Для разработки алгоритма определения начальных состояний регистров сдвига воспользуемся подходом, изложенным в работах [10, 12] применительно к трем типам ГМВП с периодом N=728, полученным в [19]. Данный подход основан на сравнении начальных состояний, определенных путем децимации символов базисной МП, и состояний, полученных в результате решения в конечных полях системы линейных уравнений.
В работе [19] все три типа ГМВП формируются на основе одной базисной МП с проверочным полиномом ЛМП(х)=х6+х+2, представляемой в виде матрицы размером [JxL] = [26x28], столбцы которой являются различными циклическими сдвигами МП с более коротким периодом J=3 -1=26; такой тип называется характеристической последовательностью (ХП).
В поле GF(33) существует четыре примитивных полинома, на основе которых можно сформировать четыре М-последовательности c периодом J=26, являющиеся одновременно характеристическими последовательностями для МП с периодом N=728. Последовательность циклических сдвигов ХП в матрице МП образует правило формирования. Один столбец матрицы состоит из нулей и называется нулевой последовательностью (НП).
Для получения трех различных ГМВП необходимо в матрице вместо ХП1 поставить в соответствии с правилом формирования другие ХП,, i = 2, 3, 4, проверочные полиномы которых являются примитивными полиномами в поле GF(33), построенном по полиному
31
fx)=x +2x+1. Данные полиномы вместе с их корнями представлены в табл. 1 (здесь а = а). В качестве индексов для полиномов h,(x) в поле GF(33), а также в поле GF(36) используются минимальные показатели степени р-сопряженных корней данных полиномов. Значения показателей степени характеризуются также параметром r, используемым при построении ГМВП [12].
Таблица 1
Элементы поля Параметр r ХП, Минимальные полиномы Корни
а' 2 1 0 а а а
а1 010 r4=17 ХП4 h1(x) = x3 +2x + 1 1 3 9 а а а
а5 212 r3=7 ХП3 h5(x) = x3 + 2x2 + x + 1 5,15,19
а7 122 r2=5 ХП2 h7(x) = x3 + x2 + 2x + 1 7,21,11
а17 210 r1=1 ХП1 h17(x) = x3 + 2x2 + 1 17,25,23
Проверочный полином для ХП1 определяется по алгоритму Берлекэмпа — Месси [5] и в соответствии с табл. 1 имеет вид
32
hxm(x) = h17(x) = x +2x +1.
Этому полиному соответствует параметр r1=1 [10, 19].
Для значений параметров r2 = 5, r3 = 7 и r4 = 17 ХП, будут формироваться с помощью примитивных полиномов, подстрочный индекс в обозначении которых определяется как 17r mod 26:
32
hxm(x) = hnx5 mod 26 (x) = hj(x) = x +x +2x+1;
32
hxro(x) = h17x7 mod 26 (x) = h5(x) = x +2x +x+1,
hxm(x) = h17x17 mod 26 (x) = M(x) = x +2x+1.
Алгоритм определения начальных состояний рассмотрим на примере ГМВП1 с параметром r2=5.
Проверочный полином hrf(x) и сегмент из 18 символов ГМВП1, необходимый для составления и решения системы из 18 линейных уравнений, имеют следующий вид [19]:
hrt(x) = hсl(x)•hс2(x)•hсз(x) = h5(x)h19(x)h31(x) = =(x6+2x5+x4+x3+2x2+x+2) (x6+x4+2x3+x2+x+2) (x6+x5+x4+2x2+2), (1)
Frt = C0C1C2_C17 = 100112200011110102. (2)
Эквивалентная линейная сложность ГМВП1 равна степени проверочного полинома: ls1 = 18. Полиномы Иа(х) являются неприводимыми полиномами 6-й степени поля GF(36), построенного по примитивному полиному _Дх)=х6+х+2, а1=а.
Начальные состояния регистров, получаемые путем децимации символов di базисной МП по индексам децимации /¿1=5, /¿2=19, /¿3=31, характеризуются выражениями
И5(х): d0=0, ds=1, d10=2, d15=1, d20=2, d25=1;
Й19(х): d0=0, d19=0, d38=1, d57=0, d76=1, d95=2; (3)
Й31(х): d0=0, d31=2, d62=2, d93=2, d124=0, d155=0.
Для нахождения сдвигов ПСП, суммируемых по mod 2, необходимо определить значения начальных состояний регистров сдвига путем решения системы линейных уравнений при заданном сегменте троичной ГМВП1 вида (2) длиной 18 символов, т.е. c0 = 1, c1 = 0, c2 = 0, c3 = 1 и т.д. Затем определить сдвиги суммируемых ПСП исходя из условия совпадения начальных состояний регистров, полученных путем децимации символов базисной МП и с помощью решения системы линейных уравнений.
Символы на выходах регистров сдвига определяются коэффициентами полиномов И5(х), И19(х) и И31(х) с помощью рекуррентных выражений
Й5(х): С6+/ = С0+/ + 2c1+/ + С2+/ + 2c3+/ + 2c4+/ + C5+„ i = 0, 1,., 721; (4)
Й19(х): Сб+/ = С0+/ + 2c1+/ + 2c2+i + C3+, + 2c4+„ i = 0, 1,., 721; (5)
Й31(х): Сб+i = C0+i + C2+i + 2c4+/ + 2c5+„ i = 0, 1,., 721. (6)
Составим систему из 18 линейных уравнений вида Yxij = cm, где i = 0, 1,., 5 — номер ячейки в регистре; j = 0, 1, 2 — номер регистра сдвига; m = 0, 1, 2,.. .17 — номер символа в сегменте ГМВП1. Каждое из уравнений определяется с помощью выражений (4)—(6). Для начальных состояний трех регистров 000001, равных начальным символам МП в каноническом виде, первые два уравнения имеют вид
х00 +х01 +х02 =1, х10 +хП +х12 =0.
Полностью система из 18 линейных уравнений в виде коэффициентов при формальных переменных представлена в табл. 2 (для r2=5). Результаты решения этой системы, т.е. начальные состояния Cij трех регистров сдвига, приведены в конце таблицы.
_Таблица 2
х00 х10 х20 х30 х40 х50 х01 х11 х21 х31 х41 х51 х02 х12 х22 х32 х42 х52 ст
1 0 0 0 0 0 1 0 0 0 0 0 1 0 0 0 0 0 1
0 1 0 0 0 0 0 1 0 0 0 0 0 1 0 0 0 0 0
0 0 1 0 0 0 0 0 1 0 0 0 0 0 1 0 0 0 0
0 0 0 1 0 0 0 0 0 1 0 0 0 0 0 1 0 0 1
0 0 0 0 1 0 0 0 0 0 1 0 0 0 0 0 1 0 1
0 0 0 0 0 1 0 0 0 0 0 1 0 0 0 0 0 1 2
1 2 1 2 2 1 1 2 2 1 2 0 1 0 1 0 2 2 2
1 0 0 0 1 0 1 2 2 1 2 2 1 2 1 1 0 0
1 0 0 0 1 2 1 2 1 0 1 0 2 1 2 1 1 0
1 2 2 2 2 1 1 1 0 0 0 0 1 0 0 1 1 0 0
1 0 0 1 1 0 1 1 0 0 0 0 1 0 0 1 1 1
1 0 0 1 1 0 0 1 1 0 0 1 0 2 0 2 0 1
1 2 2 2 2 0 0 0 1 1 0 0 1 0 2 0 2 1
2 1 0 0 1 0 0 0 0 1 1 2 0 0 0 0 1 1
1 1 0 0 2 1 1 2 2 1 2 1 1 2 1 0 2 2 0
1 0 2 2 2 0 1 0 1 0 0 2 2 1 1 1 1 0 1
0 1 0 2 2 2 2 2 1 0 1 0 0 2 1 1 1 1 0
2 1 0 1 0 1 0 2 2 1 0 1 1 0 0 1 0 0 2
Решение системы уравнений
с00 с10 с20 с30 с40 с50 с01 си с21 с31 с41 с51 с02 с12 с22 с32 с42 с52
1 1 0 1 2 0 1 1 2 0 0 2 2 1 1 0 2 0
Для удобства сравнения представим полученные начальные состояния аналогично выражениям (3)—(6), опуская номера регистров сдвига:
^(х): со=1, сх=1, с2=0, сз=1, с4=2, с5=0;
М9(х): со=1, сх=1, с2=2, сз=0, с4=0, с5=2; (7)
Азх(х): со=2, сх=1, с2=1, сз=0, с4=2, с5=0.
Можно показать, что при таких начальных состояниях трех регистров сдвига на выходе устройства будет формироваться троичная ГМВП.
Для определения сдвигов последовательностей с начальными состояниями, полученными путем децимации символов базисной МП и определяемыми согласно выражениям (3)—(6), составим таблицу (табл. 3), в которой представлены сегменты базисной МП в каноническом виде и суммируемые последовательности, являющиеся МП с проверочными полиномами И5(х), И19(х) и Азх(х). При этом каждая из суммируемых последовательностей р формируется для начальных состояний как вида (з), так и вида (7).
Таблица 3
Символы базисной МП с И1(х) и МП с И5(х), ¿19(х) и Из1(х)
do ¿1 ¿2 з 4 5 642 64з 644 645 646 647 648 649 650 651
МП: р 0 0 0 0 0 1 1 1 0 1 2 1 0 1 2 2
И5: ¿1 do d.5 ¿10 15 20 25 298 з0з з08 з1з з18 з2з з28 ззз зз8 з4з
Р5: (з) 0 1 2 1 2 1 1 0 2 0 2 1 2 0 1 1
(7) 1 1 0 1 2 0 2 1 0 1 2 1 2 1 2 2
М9: ^ do ¿19 ¿з8 57 76 95 550 569 588 607 626 645 664 68з 702 721
Р19: (з) 0 0 1 0 1 2 1 1 2 0 0 1 1 1 2 2
Р19: (7) 1 1 2 0 0 2 0 1 0 0 1 0 1 2 1 1
do ¿1 ¿2 з 4 5 278 279 280 281 282 28з 284 285 286 287
МП: Р1 0 0 0 0 0 1 2 2 0 2 1 2 0 2 1 1
НЪ1: ^ do ¿з1 ¿62 9з 124 155 610 641 672 70з 6 з7 68 99 1з0 161
Ръъ (з) 0 2 2 2 0 0 1 2 1 1 0 0 2 1 1 2
Рз1: (7) 2 1 1 0 2 0 1 0 0 2 2 2 0 0 2 2
В строке „^5: (7)" находим соответствующее строке „Р5: (з)" начальное состояние 012121, которое начинается с символа 644 базисной МП. Также с данного символа начинается состояние 001012 в строке „Р19: (7)". Иными словами, последовательности и Р19 суммируются без сдвига относительно друг друга. (Искомые начальные состояния выделены полужирным шрифтом с подчеркиванием.) В строке „Рз1: (7)" находим соответствующее строке „Рз1: (з)" начальное состояние 022200, которое начинается с символа 280 базисной МП. Тогда для формирования ГМВП последовательность Рз1 должна быть сдвинута относительно и Р19 на 644-280=з64 символа, т.е. на полпериода.
Таким образом, для формирования ГМВП с периодом N=728 и параметром г2 = 5 начальные состояния регистров сдвига с учетом (з) определяются через символы базисной МП следующими выражениями:
Й5(х): 4=0, ¿5=1, ¿10=2, ¿15=1, ¿20=2, ¿25=1;
Й19(х): ¿0=0, ¿19=0, ¿з8=1, ¿57=0, ¿76=1, ¿95=2; (8)
Лз1(х): ¿з64=0, ¿з95=1, ¿426=1, ¿457=1, ¿488=0, ¿519=0.
Представим формализованную запись алгоритма определения начальных состояний регистров сдвига для формирования троичных ГМВП с периодом N=728.
Шаг 1. Ввод исходных данных:
— неприводимые полиномы в поле GF(36), _Дх)=х6+х+2;
зз
— проверочные полиномы ХПг- в поле GF(3 ), ,/(х)=х +2х+1;
— параметры г для трех типов ГМВП: г2 = 5, гз = 7 и г4 = 17;
— проверочные полиномы ГМВП Ьп(х)=Ис1(хуИс2(х).. гк^х) [19]:
Йг1(х) = ¿с1(х>йс2(х>йсз(х) = А5(2х)^А19(х)^Аз1(х) =
=(х6+2х5+х4+хз+2х2+х+2) (х6+х4+2хз+х2+х+2) (х6+х5+х4+2х2+2); (9)
к^(х) = ЬС1(хуЬс2(хуИсз(х) = И1(хуИп(хуИз1(х) = = (х6+х3+х2+2х+2) (х6+2х5+2х4+2х3+х2+2) (х6+х4+2х2+х+2); (10)
кй(х) = кс1(хукс2(хуьсз(хукс4(хукс5(хуксб(хукс1(хукс8(х)^кс9(х) =
= к11(хук2з3хУЬ25(б^)^к4з^(хук49(хук95(хук1т(хук1оз(х)^к12^(х) = = (х6+2х5+2х3+2)^(х6+2х5+2х4+2х3+2х+2)^(х6+2х5+2х4+2х2+х+2)х х(х6+х5+2х4+х3+х+2)-(х6+2х5+2х4+х3+х2+х+2)-(х6+х5+х3+2)х
х(х6+2х4+х2+х+2)^(х6+2х5+х4+2х3+х2+2х+2)^(х6+х5+2); (11)
— наборы индексов децимации для трех типов ГМВП:
кг1(х): /¿1=5, /¿2=19, /¿3=31; (12)
кг2(х): /¿1=1, /¿2=11, /¿3=31; (13)
кг3(х): /¿1=11, /¿2=23, /¿3=25, /¿4=43, /¿5=49, /¿6=95, /¿1=101, /¿8=103, /¿9=121; (14)
— сегменты ГМВП;- соответствующей длины (определяется числом уравнений):
кг1(х): ^г1 = с0с1_с11 = 100112200011110102; (15)
кг2(х): ^г2 = с0с1_с11 = 200120200011210112; (16)
кг3(х): ^г3 = с0с1_с53 = 10011220002111002101020202002120220211002001222020121201. (11)
Шаг 2. Составление предварительных наборов начальных состояний вида (3) на основе децимации символов базисной МП по индексам ¡¿^ в соответствии с (12) - (14).
Шаг 3. Составление системы линейных уравнений с учетом рекуррентных выражений (4)—(6) и сегментов ГМВЦ- (15)—(11).
Шаг 4. Формирование начальных состояний регистров сдвига вида (1) на основе решения системы линейных уравнений.
Шаг 5. Определение сдвигов отдельных суммируемых последовательностей на основе их сравнения для начальных состояний, полученных как путем децимации, так и на основе решения системы уравнений.
Шаг 6. Формирование результирующего набора начальных состояний регистров сдвига вида (8).
Начальные состояния регистров сдвига в устройствах формирования ГМВП2 с ЭЛС ¡а =18 и ГМВП3 с ЭЛС ¡ц =54 определяются для параметров г3=1 и г4=11 аналогично процедуре для ГМВП1. В результате выполнения алгоритма с учетом того, что проверочный полином
3 2 3
для ХП3 равен к5(х) = х +2х +х+1, а для ХП4 — к1(х) = х +2х+1, а также с учетом выражений (10)—(11) получены решения систем из 18 и 54 линейных уравнений, приведенные в табл. 4.
Таблица 4
Решение системы у равнений для ГМВП2
с00 с10 с20 с30 с40 с50 с01 сц с21 с31 с41 с51 с02 с12 с22 с32 с42 с52
2 0 0 2 1 2 2 2 0 1 2 2 1 1 0 1 2 2
Решение системы у равнений для ГМВП3
с00 с10 с20 с30 с40 с50 с01 сц с21 с31 с41 с51 с02 с12 с22 с32 с42 с52
1 0 2 2 2 1 1 2 2 2 1 2 1 2 1 1 1 0
с03 с13 с20 с33 с43 с53 с04 с14 с24 с34 с44 с54 с05 с15 с25 с35 с45 с55
2 2 1 2 2 2 2 2 0 1 1 0 2 0 1 2 1 1
с06 с16 с26 с36 с46 с56 с01 с11 с21 с31 с41 сц с08 с18 с28 с38 с48 с58
1 1 0 0 2 2 2 2 2 0 1 0 1 1 0 0 2 0
В результате сравнения начальных состояний, вычисленных путем децимации символов базисной МП по индексам (13), (14) и полученных на основе решения систем уравнений (см. табл. 4), определены сдвиги суммируемых последовательностей и сформированы результирующие наборы начальных состояний ячеек регистров сдвига (Яг), входящих в устройства синтеза ГМВП2 и ГМВП3: см. табл. 5, где также приведены наборы для ГМВП1.
С помощью приведенных в табл. 5 данных можно сформировать две ГМВП с ЭЛС ¡а1= ¡а2 =18 и одну ГМВП с ЭЛС ¡а3=54 на основе базисной МП с полиномом кМП(х) = =к1(х) = х6+х+2. Так как в поле GF(36) существует 48 примитивных полиномов
[19, табл. 2], то аналогичные распределения можно получить для любой из 48 МП, рассматриваемых в качестве базисных.
_Таблица 5
ЭЛС ls hci(x) Начальные состояния ячеек регистра сдвига
Яо Я1 Я2 Я3 Я4 Я5
18 h5(x) d0=0 d5=1 d10=2 d15=1 d20=2 d25=1
hw(x) d0=0 d19=0 d38=1 d57=0 d76=1 d95=2
h31 (x) d364=0 d395=1 d426=1 d457=1 d488=0 d519=0
18 hq(x) d0=0 d7=0 d14=0 d21=0 d28=2 d35=2
h„(x) d0=0 d„=1 d22=0 d33=1 dAA—1 d55=1
h37(x) d364=0 d401 =0 d438=2 d457=0 d512=0 d549—2
54 hn(x) d0=0 dn=1 d34=1 d51=1 d68=2 d85=0
h23(x) d0=0 d23=1 d46=0 d69=1 d92=2 d„5=2
h2s(x) d0=0 d25=1 ds0=0 d75=1 d100=2 d125=2
h43(x) d364=0 d407=1 d450=0 d493— 1 d563=1 d579=2
h49(x) d364=0 d413=2 d462=0 d511=2 d560=1 d609=0
h95(x) d364=0 d459 1 d554=2 d649=1 d16=1 dm=0
h101(x) d0=0 d101=0 d202=2 d303=0 d404=1 d505=1
h103(x) d364=0 d467=2 d570=1 d673=2 d48=1 d151=2
h121(x) d0=0 d121=2 d242=1 d363=2 d484=1 d605=2
При произвольной базисной МП с примитивным полиномом hg(x) для получения новой совокупности полиномов-сомножителей hci(x) необходимо индексы полиномов прежней совокупности умножить по mod 728 на параметр g, являющийся показателем степени корня полинома hg(x). Набор индексов совокупности полиномов-сомножителей h^x) в работе [19] определен как вектор сомножителей.
Пример. Определим проверочные полиномы и начальные состояния ячеек регистров сдвига для трех типов ГМВП с ЭЛС lsi = ls2 =18 и ls3 = 54, если в качестве базисной используется МП с примитивным полиномом hMn(x) = h47(x) = x6+x5+2x3+2x2+2.
В работе [ 19] для hMn(x) = h47(x) определены векторы сомножителей для проверочных полиномов ГМВП:
А5 = (107,29,1), А7 = (49,95,121), Ап = (71,67,149,239,119,97,41,395,197).
Например, проверочный полином hr3(x) ГМВП3 при г4=17 имеет вид
hr3(x) = hсl(x)•hс2(x)•hсз(x)•hс4(x)•hс5(x)•hс6(x)•hс7(x)•hс8(x)•hс9(x) =
= h11(xyh61(xyh149(xyh239(xyhn9(xyh91(xyh41(xyh395(x)^h191(x) =
= (x6+x5+x3+2x2+x+2>(x6+2x4+2x3+x2+2x+2>(x6+x4+x3+x2+2x+2)x
x(x6+2x4+x3+x2+x+2>(x6+2x5+2x4+x3+2>(x6+2x5+2x3+2x2+2x+2)x
x(x6+x5+x4+2x3+2x2+2x+2>(x6+2x5+x4+2x2+2>(x6+x4+2x3+x+2). (18)
При определении начальных состояний ячеек регистров сдвига обязательно соблюдение последовательности полиномов-сомножителей, так как начальные сдвиги определены для конкретных h^x).
С учетом того, что начальные состояния всех ячеек Я0 остаются без изменений, значения остальных ячеек определяются путем децимации символов исходной базисной МП по индексам новых векторов состояний. Например, для ГМВП с ЭЛС ls3 =54 на основе исходной базисной МП с полиномом hMn(x) = h1(x) = x6+x+2 начальные состояния четвертого регистра с полиномом h^(x) = h43(x) равны d364=0, d407=1, d450=0, d493=1, d563=1, d579=2. При замене на базисную МП с полиномом h^(x) = h47(x) = x6+x5+2x3+2x2+2 начальные состояния четвертого регистра с новым полиномом h^(x) = h239(x) равны d364=0, d603=0, d114=1, d353=0, d592=2,
d103=1.
Результаты расчета начальных состояний ячеек регистров сдвига при данной базисной МП для трех типов ГМВП приведены в табл. 6.
_Таблица 6
ЭЛС ls ¿ci(x) Начальные состояния ячеек регистра сдвига
Я0 Я1 Я2 Я3 Я4 Я5
18 ¿107(X) $0—0 $107—2 $214—1 $321—2 $428—1 $535—0
¿29(X) $0—0 $29—1 $58—0 $87—1 $116—0 $145—1
¿1(x) $364—0 $365—0 $366—0 $367—0 $368—0 $369—2
18 ¿49(X) $0—0 $49 1 $98—0 $147—1 $196—2 $245—0
¿95(X) $0—0 $95—2 $190—1 $285—2 $380—2 $475—0
¿121(X) $364—0 $485—1 $606—2 $727—1 $120—2 $241—2
54 ¿71 (x) $0—0 $71—2 $142—1 $213—2 $284—0 $355—2
¿67(X) $0—0 $67—0 $134—2 $201—0 $268—1 $335—1
¿149(X) $0—0 $149 0 $298—1 $447—0 $596—1 $17—1
¿239(X) $364—0 $603—0 $114—1 $353—0 $592—2 $103—1
¿119(X) $364—0 $483—2 $602—0 $721—2 $112—0 $231—2
¿97(X) $364—0 $461—2 $558—2 $655—2 $24—0 $121—2
¿41 (x) $0—0 $41—2 $82—2 $123—2 $164—2 $205—2
¿395(X) $364—0 $31—2 $426—1 $93—2 $488—0 $155—0
¿197(X) $0—0 $197—0 $394 1 $591—0 $60—2 $257—2
На рисунке показана структурная схема устройства формирования ГМВП3 на основе базисной МП с полиномом ¿Mn(X)—¿47(x)—x6+x5+2x3+2x2+2, параметром r2 — 17 и ЭЛС ls3 — 54, состоящая из девяти регистров сдвига. Умножители и сумматоры по mod 3 в цепи обратной связи регистров расставляются в соответствии с коэффициентами неприводимых полиномов (18). Начальные состояния ячеек регистров сдвига соответствуют нижней части табл. 6.
Я5 -Я4 -* Я3 J-*- Я2
^ ^ ^ ^
$355—2 fi?284—0 $213—2 $142—1 $71—2 $0—0 -
¿67 (X)
Я5
T
d335—1 d268—1 $201—0 $134—2 $67—0 $0—0
$17— 1 $596— 1 $447—0 $298— 1 $149—0 $0—0
¿239(X)
-GH-GH-CH--CH--
_ О k—2_Q k—2_G k—2
Я4 - И Я3 ^^ Я2 Я1 ^^ Я0
Я5 Я4
т f
$103—1 $592—2 $353—0 $144—1 $603—0 $364—0
Я3
т
Я2
Г
Я1
т
Я0
т
$231—2 $112—0 $721—2 $602—0 $483—2 $364—0
¿97(X)
Я
Я
Я
о*-
Я
Я1
Я0
Г t t Т Т 7
$205—2 $164—2 $123—2 $82—2 $41 —2 $0—0
-о
$155—0 $488—0 $93—2 $426—1 $31—2 $364—0
-х$и-
Выход
—
ГМВП
$257—2 $60—2 $591 —0 $394—1 $197—0 $0—0
О
$121—2 $24—0 $655—2 $558—2 $461—2 $364—0
■в--0-*-Й^-Ghu-
На выходе устройства формируется ГМВП3 с периодом N=728 и двухуровневой ПАКФ (приведено по 112 начальных и конечных символов):
Fr3=
0 1 2 1 2 1 2 2 0 1 2 0 2 1 0 1 2 1 2 2 0 2 2 0 0 1 0 1
1 0 2 2 2 0 2 1 2 0 1 1 2 2 0 0 1 1 0 2 2 0 0 1 2 1 2 0
0 2 0 1 0 0 1 2 0 0 2 2 1 0 0 2 0 2 2 1 1 0 0 0 0 1 2 2
1 0 0 0 0 2 2 0 2 2 0 2 2 1 0 0 2 1 2 2 1 1 1 1 2 2 0 0
2 1 1 0 1 0 0 0 1 0 0 0 0 1 0 1 2 0 1 0 0 0 0 2 1 1 2 1
2 1 0 2 0 1 0 1 1 1 1 2 0 2 0 1 1 0 2 0 1 2 2 2 1 0 1 1
1 1 2 1 2 0 1 2 2 0 2 2 1 2 0 1 1 2 1 1 1 0 0 1 2 0 0 1
2 2 0 1 0 1 2 2 1 1 2 0 2 2 0 2 1 1 0 2 0 2 2 2 1 2 2 2
соответствии с алгоритмом Берлекэмпа - Месси определяется проверо
(19)
ном, который для компактности представлен в виде коэффициентов при формальной переменной х по убыванию степени от х54 до х0:
^=1220011011120120210222201111020000011212000102211112222. (20)
ЭЛС ГМВП3 /^=54, что в девять раз превышает ЭЛС МП /мп=6.
Легко проверить, что неприводимые полиномы (18) являются делителями полинома (20).
Таким образом, разработан алгоритм определения начальных состояний ячеек регистров сдвига, входящих в устройство формирования троичных ГМВ-последовательностей с периодом N=728, основанный на сравнении начальных состояний, полученных в результате решения в конечных полях системы линейных уравнений, и состояний, определенных путем децимации символов базисной М-последовательности.
В соответствии с алгоритмом определены начальные состояния при формировании трех типов ГМВП для двух базисных МП (с полиномами к\(х) и к47(х)). Всего можно сформировать 144 ГМВП с периодом N=728.
Максимальный выигрыш в структурной скрытности ГМВП составляет 9 раз по сравнению с МП.
Достоинство алгоритма заключается в том, что при любой базисной МП для определения начальных состояний ячеек регистров сдвига, входящих в устройство формирования ГМВП, используются символы только одной базисной МП с полиномом ^мп(х) = М(х) = х6+х+2.
Разработанный алгоритм может найти применение при синтезе устройств формирования ГМВП с основанием р=3, а также при формировании широкополосных недвоичных сигналов в системах передачи цифровой информации, к которым предъявляются повышенные требования по конфиденциальности и помехозащищенности.
список литературы
1. Вишневский В. М., Ляхов А. И., Портной С. Л., Шахнович И. В. Широкополосные беспроводные сети передачи информации. М.: Техносфера, 2005. 592 с.
2. Скляр Б. Цифровая связь. Теоретические основы и практическое применение: Пер. с англ. М.: Изд. дом „Вильямс", 2003. 1104 с.
3. Ипатов В. П. Широкополосные системы и кодовое разделение сигналов. Принципы и приложения: Пер. с англ. М.: Техносфера, 2007. 488 с.
4. CDMA: прошлое, настоящее, будущее / Под ред. Л. Е. Варакина и Ю. С. Шинакова. М.: Международная академия связи, 2003. 608 с.
5. Golomb S. W., Gong G. Signal Design for Good Correlation for Wireless Communication, Cryptography and Radar. Cambridge Univ. Press, 2005. 438 p.
6. Chung H. B., No J. S. Linear span of extended sequences and cascaded GMW sequences // IEEE Transact. on Information Theory. 1999. Vol. 45, N 6. P. 2060—2065.
7. Rizomiliotis P., Kalouptsidis N. Results on the nonlinear span of binary sequences // IEEE Transact. on Information Theory. 2005. Vol. IT—51. P. 1555—1563.
8. Ипатов В. П. Периодические дискретные сигналы с оптимальными корреляционными свойствами. М.: Радио и связь, 1992. 152 с.
9. No Jong-Seon. Generalization of GMW sequences and No sequences // IEEE Transact. on Information Theory. 1996. Vol. 42, N 1. P. 260—262.
10. Стародубцев В. Г., Чернявских А. Е. Формирование троичных последовательностей Гордона — Миллса — Велча на основе регистров сдвига // Изв. вузов. Приборостроение. 2016. Т. 59, № 3. С. 202—210.
11. Tsankov T., Trifonov T., Staneva L. An algorithm for synthesis of phase manipulated signals with high structural complexity // J. Scientific & Applied Research. 2013. Vol. 4. P. 80—87.
12. Стародубцев В. Г. Формирование пятеричных последовательностей Гордона — Миллса — Велча для систем передачи дискретной информации // Тр. СПИИРАН. 2019. Т. 18, № 4. С. 912—948.
13. Lee Wijik, Kim Ji-Youp, No J.S. New families of p-ary sequence of period (p|n-1)/2 with low maximum correlation magnitude // IEICE Transact. on Communications. 2014. Vol. E97-B, N 1. P. 2311—2315.
14. Cho Chang-Min, Kim Ji-Youp, No J. S. New p-ary sequence families of period (pn-1)/2 with good correlation property using two decimated m-sequences // IEICE Transact. on Communications. 2015. Vol. E98, N 7. P. 1268—1275.
15. Tasheva Z. A short survey of p-ary pseudo-random sequences // J. Scientific & Applied Research. 2014. Vol. 2. P. 17—26.
16. Xia Y., Chen S. A new family of p-ary sequences with low correlation constructed from decimated sequences // IEEE Transact. on Information Theory. 2012. Vol. 58, N 9. P. 6037—6046.
17. Helleseth T., Kumar P. V., Martinsen H. A new family of ternary sequences with ideal two-level autocorrelation function // Designs, Codes and Cryptography. 2001. Vol. 23, N 2. P. 157—166.
18. TangX. H., Pingzhi Z. F. A class of pseudonoise sequences over GF(p) with low correlation zone // IEEE Transact. on Information Theory. 2001. Vol. 47, N 4. Pp. 1644—1649.
19. Стародубцев В. Г., Ткаченко В. В., Малышева Е. А. Формирование троичных ГМВ-последовательностей с периодом N=728 в системах передачи измерительной информации // Изв. Тульск. гос. ун-та. Технические науки. 2019. Вып. 6. С. 192—203.
20. Стародубцев В. Г., Бородько Д. Н., Мышко В. В. Алгоритм формирования ГМВ-последовательностей с периодом N=4095 в системах передачи телеметрической информации // Авиакосмическое приборостроение. 2018. № 5. С. 3—15.
Сведения об авторах
Виктор Геннадьевич Стародубцев — канд. техн. наук, доцент; ВКА им. А. Ф. Можайского, преподаватель;
E-mail: vgstarod@mail.ru
Владимир Викторович Ткаченко — канд. техн. наук; ВКА им. А. Ф. Можайского, ст. преподаватель;
E-mail: vik_hohol@mail.ru
Елизавета Алексеевна Боброва — слушатель; ВКА им. А. Ф. Можайского, E-mail: m_ea98@mail.ru
Поступила в редакцию 23.01.2020 г.
Ссылка для цитирования: Стародубцев В. П., Ткаченко В. В., Боброва Е. А. Формирование троичных последовательностей с высокой структурной скрытностью в системах передачи цифровой информации // Изв. вузов.
Приборостроение. 2020. Т. 63, № 5. С. 405—416.
FORMATION OF TERNARY SEQUENCES WITH HIGH STRUCTURAL SECRECY IN DIGITAL INFORMATION TRANSFER SYSTEMS
V. G. Starodubtsev, V. V. Tkachenko, E. A. Bobrova
A. F. Mozhaisky Military Space Academy, 197198, St. Petersburg, Russia
E-mail: vka@mil.ru
An algorithm for determining initial states of shift registers included in the Gordon-Mills-Welch (GMW) ternary sequence generation device with the period of N = 728 is presented. The algorithm is based on comparison of the initial states obtained by solving the system of linear equations in the finite fields and the states determined by decimation of characters of the basic M-sequence. Ternary M-sequences and GMW-sequences have the same two-level periodic autocorrelation function, but different structural secrecy, characterized by equivalent linear complexity. The GMW-sequence is formed using a basic M-sequence with a similar period when it is presented in the form of quasi-square matrix. It is shown that for each of the 48 primitive polynomials in the finite field GF (36), three GMW-sequences can be formed. For binary GMW sequences, the initial states of the shift registers are formed by decimating the symbols of the basic M-sequence, presented in canonical form, by decimation indices corresponding to the roots of indivisible polynomial factor factors. For ternary GMW sequences, the individual summable components have an additional half-period shift of the base M-sequence. It is argued that the obtained results can be used for generating broadband non-binary signals in digital information transmission systems.
Keywords: pseudorandom sequences, finite fields, indivisible and primitive polynomials, structural secrecy, decimation, shift registers
1. Vishnevskij V.M., Lyahov A.I., Portnoj S.L., Shahnovich I.V. Shirokopolosnye besprovodnye seti peredachi informacii (Broadband Wireless Data Transmission Network), Moscow, 2005, 592 p. (in Russ.)
2. Sklar B. Digital Communications: Fundamentals and Applications, Prentice Hall, 2001, 1079 p.
3. Ipatov V.P. Spread Spectrum and CDMA. Principles and Applications, NY, John Wiley and Sons Ltd., 2005, 488 p.
4. Varakin L.E. and Shinakov Yu.S., ed., CDMA: proshloe, nastoyashchee, budushchee (CDMA: Past, Present, Future), Moscow, 2003, 608 p. (in Russ.)
5. Golomb S.W., Gong G. Signal Design for Good Correlation for Wireless Communication, Cryptography and Radar, Cambridge University Press, 2005, 438 p.
6. Chung H.B., No J.S. IEEE Transactions on Information Theory, 1999, no. 6(45), pp. 2060-2065.
7. Rizomiliotis P., Kalouptsidis N. IEEE Transactions on Information Theory, 2005, vol. IT-51, pp. 1555-1563.
8. Ipatov V.P. Periodicheskie diskretnye signaly s optimal'nymi korrelyacionnymi svojstvami (Periodic Discrete Signals with Optimum Correlation Properties), Moscow, 1992, 152 p. (In Russ.).
9. No Jong-Seon. IEEE Transactions on Information Theory, 1996, no. 1(42), pp. 260-262.
10. Starodubtsev V.G., Chernjavskih A.E. Journal of Instrument Engineering, 2016, no. 3(59), pp. 201-210. (in Russ.)
11. Tsankov T., Trifonov T., Staneva L. Journal Scientific & Applied Research, 2013, vol. 4, pp. 80-87.
12. Starodubtsev V.G. Trudy SPIIRAN (SPIIRAS Proceedings), 2019, no. 4(18), pp. 912-948. (in Russ.)
13. Lee Wijik, Kim Ji-Youp, No J.S. IEICE Transactions on Communications, 2014, no. 1(E97-B), pp. 2311-2315.
14. Cho Chang-Min, Kim Ji-Youp, No J.S. IEICE Transactions on Communications, 2015, no. 7(E98), pp. 1268-1275.
15. Tasheva Z. Journal Scientific & Applied Research, 2014, vol. 2, pp. 17-26.
16. Xia Y., Chen S. IEEE Transactions on Information Theory, 2012, no. 9(58), pp. 6037-6046.
17. Helleseth T., Kumar P.V., Martinsen H. Designs, Codes and Cryptography, 2001, no. 2(23), pp. 157-166.
18. Tang X.H., Pingzhi Z.F. IEEE Transactions on Information Theory, 2001, no. 4(47), pp. 1644-1649.
19. Starodubtsev V.G., Tkachenko V.V., Malysheva E.A. Izvestiya Tulskogo gosudarstvennogo universiteta. Tekhnicheskie nauki (Proceedings of Tula State University. Technical science), 2019, no. 6, pp. 192-203 (in Russ.)
20. Starodubtsev V.G., Borod'ko D.N., Myshko V.V. Aviakosmicheskoe priborostroenie (Aerospace Instrumentation), 2018, no. 5, pp. 3-15 (in Russ.)
REFERENCES
Vladimir V. Tkachenko
Victor G. Starodubtsev
Data on authors
PhD, Associate Professor; A. F. Mozhaisky Military Space Academy, Lecturer; E-mail: vgstarod@mail.ru
PhD; A. F. Mozhaisky Military Space Academy, Senior Lecturer; E-mail: vik_hohol@mail.ru
Elizaveta A. Bobrova — Student; A. F. Mozhaisky Military Space Academy,
E-mail: m_ea98@mail.ru
For citation: Starodubtsev V. G., Tkachenko V. V., Bobrova E. A. Formation of ternary sequences with high structural secrecy in digital Information transfer systems. Journal of Instrument Engineering. 2020. Vol. 63, N 5. P. 405—416 (in Russian).
DOI: 10.17586/0021-3454-2020-63-5-405-416
