Главная Помехоустойчивое кодирование



Oi - = -2; 02 = Уг/уо =1; 03 = уз/уо = -2 и

a(z) =z - 22 +z - 2.

2.3. ВЫЧИСЛЕНИЕ ОСОБЫХ ПРОДОЛЖЕНИЙ ГАНКЕЛЕВЫХ МАТРИЦ

Дробно-рациональные функции и особые продолжения ганкелевых матриц. Начнем со следующего утверждения [23]. Бесконечная ганкелева матрица (2.3) имеет конечный ранг г тогда н только тогда, когда сумма формального ряда

2 Oj/zj является дробно-рацнональной функцией а;(2)/ст(2). В этом случае / = 1

ранг г ганкелевой матрицы совпадает с числом корней полинома a{z) с учетом нх кратностей. Из этого утверждения следует, что элементы Д2» • • • * матрицы f. вместе с компонентами Д2г+г = 1, 2, . . . особого продолжения удовлетворяют равенству

2 а{2)

(2.34)

Пусть и > 2 н расширенная ганкелева матрица А задана своими 2п компонентами Д1, 02у -. . , 02п Усекая ряд, стоящий в левой части равенства (2.34), можем представить его в виде

01 02 °2п

+ -г + • • • + -

2 .2П

(2.35)

где полином Q{z} = «2« + i/" + ги + г/ " + • • * определяется продолжениями 2+12и+2 • • • (причем особыми) ганкелевой матрицы >4„„ + j.

Подставим сумму (2.35) в равенство (2.34) и осуществим приведение к общему знаменателю:

Л{2) а;(л) 2«

2 2" а(2)

-(2) или A{z)a{z) =2"[а;(2) -(2)0(2)], (2.36)

где /4(2) = Ji2" +Д2" + . . . + • • • + 2и-1 +2и-

Предположим, что ранг г расширенной ганкелевой матрицы Л. равен и. Составим формальный полином степени п:

СТ(2) = CTqz" + CTj 2""* + . . . + CT„ j2 + ст„. (2.37)

Если же ранг г матрицы Л„ „ + меньше ин«-г = s,tos младших коэффициентов полинома (2.37) приравниваются нулю. Рассмотрим произведение

Л(2)а(2) = (fli22"-*+fl2"~ + -.- + «2«-2+2«-l+2«) 54



X (OqZ + oiz • • • + <-2 n-i n - i = 2П . i = n

V ,2rt-i V ,« - 7 = 2/ az Li a.z • .

1 = 1 7 = 0

(2.38)

Произведение (2.38) можно разбить на три части: A{z)o{z) - P{z) + Q{z) + + R{z)y где P{z), Q{z) н R{z) - соответственно старшая (избыток), средняя и младшая (остаток) части, а именно:

P{z) aiOQZ~ + (j2ao + JiCTi)2"~ + . . . +

+ {aia +a2a 2 + . . . + a oi+aaQ)z ; (2.39)

2h 1

Q{Z) = (fli CT„ + J2<J«-1 + • • • + rtl +rt+l o) +

2w - 2

+ (2 "3 <п- 1 + • • • + fl+I П + 2оУ + . . . +

+ « + 1 1 + • • • + 2«-1 2«<o)"; (2-40)

+ («2«<«-2 + 2«-l«-I 2«-2*«>--+ (fl2«<2 + «2„ ja3 + ...+fl„ + 2*«)"~ + + («2«<l+2«-l*2 + -. + =

где символ - означает "принято по обозначению".

В выражении (2.41) коэффициенты перед степенями г обозначены через г, 1 = 1, 2, . , «. Среди коэффициентов могут быть н нулевые. Ясно, что средняя часть Q{z) в произведении A{z)o{z) равна нулю на основании рекуррентного соотношения (2.29) для вычисления особых продолжений, так как выражение в каждой скобке (2.40) удовлетворяет формуле (2.29).

Остаток R{z) является полу-НОД полиномов z к А{2). Методы его нахождения будут рассмотрены ниже, а здесь покажем как по остатку R{z) и известному полиному A{z) может быть рассчитан полином а(2), а также вычислены без рекуррентных процедур компоненты 2и+/» / = 1, 2, . .. , продолжения расширенной ганкелевой матрицы А n + i-

Из представления (2.41) следует, что полином о {z) может быть вычислен непосредственно по полиномам R{z) hA{z) на основании соотношения

СТ = Rz : А, остаток -

(2.42)

где для простоты полиномы заменены соответствующими векторами.

Отметим также, что обращенный полином P{z) для суммы старших членов представления (2.38) можно рассматривать как "остаток" (избыток), а точнее ~



как полу-НОД для полиномов и /1(2). Тогда вместо соотношения (2.42) можно использовать

о = Pz : А, остаток ~ Л,

(2.43)

где Р определяется из соотношения (2.39).

Поясним сказанное. Ранее рассматривался пример с полиномами a(z) = -2z + z - 2и A{z) = 2z5 + + 3 + + ц2 + 21 при = 3. Умножив эти полиномы, получим

(2г + 3z + 2 + 32 + 1U +21) (2 - 22 + 2 2)

- 22® + 3Z +2 + 3Z + 112+ 212

- 42 - 6z - 22 - 62** - 222 42z +

+ 22 + 3Z +2** + 32 + 112 + 2I2 -

- 42 - 62** 2z - 62 - 222 - 42

= 22® - z - 32 + 02 + 02** + 02 - 37z - 2 - 42

> .--/ .-

zP(z) Qiz) = 0 R{z)

Такимобразом, P(2) = 22 - 2 - 3; Л (2) = -372 - 2 42.

Для упрощения записи вместо полиномов будем использовать их коэффициенты, помня, что разряды должны обрабатываться независимо друг от друга, т. е. без заемов и переносов.

Согласно соотношению (2-42) имеем

42, -1, -37, О, О, О, О, О, О : 21, U, 3, 1, 3. 2 = -2,1,-2, 1 42,-22, -6, -2, -6, -4

21, -31. 2, 6. 4, О 21. П. 3, 1, 3, 2

(2.44)

-42, -1. 5. 1, -2. О -42,-22, -в, -2, -6, -А

21, И, 3, 4, 4, О 21. 11, 3, 1, 3, 2

3. 1, ~2 = -Р. Аналогично, используя соотношение (2.43), получаем

* Определение полу-НОД приведено далее



0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 [17] 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94


0.0253