Основы теории информации и криптографии - -Учебное пособие (В.В. Лидовский)

8. лекция: математическая модель системы связи

Рассматриваются такие классы кодов, как коды с исправлением и обнаружением ошибок. Хорошее математическое обоснование материала лекции. Описываются последовательные коды и их применение на практике. Матричное кодирование позволяет использовать меньший объем памяти при кодировании информации. Практические задания помогут лучше разобраться в сложном материале лекции

Коды делятся на два больших класса. Коды с исправлением ошибок имеют целью восстановить с вероятностью, близкой к единице, посланное сообщение. Коды с обнаружением ошибок имеют целью выявить с вероятностью, близкой к единице, наличие ошибок.

Простой код с обнаружением ошибок основан на схеме проверки четности, применимой к сообщениям a_1ldots a_mлюбой фиксированной длины m. Схема кодирования определяется следующими формулами: E(a_1ldots a_m)=a_1ldots a_ma_{m+1},

a_{m+1}=egin{cases}0,&	ext{если $sum^m_{i=1}a_i$ --- четная;}\
1,&	ext{если $sum^m_{i=1}a_i$ --- нечетная.}
end{cases}

Таким образом, sum^{m+1}_{i=1}a_iдолжна быть четной.

Соответствующая схема декодирования тривиальна:

D(a_1ldots a_ma_{m+1})=
egin{cases}a_1ldots a_m,&	ext{если $sum^{m+1}_{i=1}a_i$ --- четна;}\
         hbox{it \<ошибка\>},&	ext{если $sum^{m+1}_{i=1}a_i$ --- нечетна.}
end{cases}

Разумеется, что четность sum^{m+1}_{i=1}a_iне гарантирует безошибочной передачи.

Пример. Проверка четности при m=2реализуется следующим кодом (функцией E): 00
ightarrow000, 01
ightarrow011, 10
ightarrow101, 11
ightarrow110. В двоичном симметричном канале доля неверно принятых сообщений для этого кода (хотя бы с одной ошибкой) равна q^3+3pq^2+3p^2q(три, две или одна ошибка соответственно). Из них незамеченными окажутся только ошибки точно в двух битах, не изменяющие четности. Вероятность таких ошибок 3pq^2. Вероятность ошибочной передачи сообщения из двух бит равна 2pq+q^2. При малых qверно, что 3pq^2ll2pq+q^2.

Рассмотрим (m,3m)-код с тройным повторением. Коды с повторениями очень неэффективны, но полезны в качестве теоретического примера кодов, исправляющих ошибки. Любое сообщение разбивается на блоки длиной mкаждое и каждый блок передается трижды - это определяет функцию E. Функция Dопределяется следующим образом. Принятая строка разбивается на блоки длиной 3m. Бит с номером i(1le ile m)в декодированном блоке получается из анализа битов с номерами i, i+m, i+2mв полученном блоке: берется тот бит из трех, который встречается не менее двух раз. Вероятность того, что бит в данной позиции будет принят трижды правильно равна p^3. Вероятность одной ошибки в тройке равна 3p^2q. Поэтому вероятность правильного приема одного бита равна p^3+3p^2q. Аналогичным образом получается, что вероятность приема ошибочного бита равна q^3+3pq^2.

Пример. Предположим q=0.1. Тогда вероятность ошибки при передачи одного бита - 0.028, те. этот код снижает вероятность ошибки с 10\% до 2.8\%. Подобным образом организованная передача с пятикратным повторением даст вероятность ошибки на бит q^5+5pq^4+10p^2q^3=0.00856=0.856\%,т.е. менее 1\%. В результате вероятность правильной передачи строки длиной 10 возрастет с 0.9^{10}approx35\%до 0.972^{10}approx75\%при тройных повторениях и до 0.99144^{10}approx92\%при пятикратных повторениях.

Тройное повторение обеспечивает исправление одной ошибки в каждой позиции за счет трехкратного увеличения времени передачи.

Рассмотрим (2048,2313)-код, используемый при записи данных на магнитофонную ленту компьютерами Apple II. К каждому байту исходных данных прибавляется бит четности и, кроме того, после каждых таких расширенных битом четности 256 байт добавляется специальный байт, также расширенный битом четности. Этот специальный байт, который называют контрольной суммой (check sum), есть результат применения поразрядной логической операции "исключающее ИЛИ" (XOR) к 256 предшествующим расширенным байтам. Этот код способен как обнаруживать ошибки нечетной кратности в каждом из отдельных байт, так и исправлять до 8 ошибок в блоке длиной 256 байт. Исправление ошибок основано на том, что если в одном из бит одного из байт 256 байтового блока произойдет сбой, обнаруживаемый проверкой четности, то этот же сбой проявится и в том, что результат операции "исключающее ИЛИ" над всеми соответствующими битами блока не будет соответствовать соответствующему биту контрольной суммы. Сбойный бит однозначно определяется пересечением сбойных колонки байта и строки бита контрольной суммы. На рис. 8.1 изображена схема участка ленты, содержащего ровно 9 ошибок в позициях, обозначенных p_1, p_2, ldots, p_9. Расширенный байт контрольной суммы обозначен CS, а бит паритета (в данном случае четности) - PB (parity bit). Ошибка в позиции p_1может быть исправлена. Ошибки в позициях p_4, p_5, p_6, p_7можно обнаружить, но не исправить. Ошибки в позициях p_2, p_3, p_8, p_9невозможно даже обнаружить.

Рис. 8.1. 

Приведенные ранее примеры простейших кодов принадлежат к классу блочных. По определению, блочный код заменяет каждый блок из mсимволов более длинным блоком из nсимволов. Следовательно, (m,n)-коды являются блочными. Существуют также древовидные или последовательные коды, в которых значение очередного контрольного символа зависит от всего предшествующего фрагмента сообщения. Работа с древовидным шумозащитным кодом имеет сходство с работой с арифметическим кодом для сжатия информации.

Расстоянием (Хэмминга) между двоичными словами длины nназывается количество позиций, в которых эти слова различаются. Это одно из ключевых понятий теории кодирования. Если обозначить двоичные слова как a=a_1ldots a_nи b=b_1ldots b_n, то расстояние между ними обозначается d(a,b).

Весом двоичного слова a=a_1ldots a_nназывается количество единиц в нем. Обозначение w(a). Можно сказать, что w(a)=sum^n_{i=1}a_i.

Пример. Пусть a=1001и b=0011, тогда w(a)=w(b)=2, d(a,b)=2.

Далее операция +при применении к двоичным словам будет означать поразрядное сложение без переноса, т.е. сложение по модулю 2 или "исключающее ИЛИ" (XOR).

Расстояние между двоичными словами aи bравно весу их поразрядной суммы, т.е. d(a,b)=w(a+b).

Если два слова различаются в каком-либо разряде, то это добавит единицу к весу их поразрядной суммы.

Следовательно, если aи b- слова длины n, то вероятность того, что слово aбудет принято как b, равна p^{n-d(a,b)}q^{d(a,b)}.

Наример, вероятность того, что слово 1011 будет принято как 0011, равна p^3q.

Для возможности обнаружения ошибки в одной позиции минимальное расстояние между словами кода должно быть большим 1.

Иначе ошибка в одной позиции сможет превратить одно кодовое слово в другое, что не даст ее обнаружить.

Для того, чтобы код давал возможность обнаруживать все ошибки кратности, не большей k, необходимо и достаточно, чтобы наименьшее расстояние между его словами было k+1.

Достаточность доказывается конструктивно: если условие утверждения выполнено для E, то в качестве декодирующей функции Dследует взять функцию, сообщающую об ошибке, если декодируемое слово отличается от любого из слов из образа E. Необходимость доказывается от противного: если минимальное расстояние k'<k+1, то ошибка в k'позициях сможет превратить одно кодовое слово в другое.

Для такого кода вероятность того, что ошибки в сообщении останутся необнаруженными, равна

sum^n_{i=k+1}C^i_np^{n-i}q^i=C^{k+1}_np^{n-k-1}q^{k+1}+cdots+
C^{n-1}_npq^{n-1}+q^napprox

[при малых qи не слишком маленьких k]approx C^{k+1}_np^{n-k-1}q^{k+1}.

Для того, чтобы код давал возможность исправлять все ошибки кратности, не большей k, необходимо и достаточно, чтобы наименьшее расстояние между его словами было 2k+1.

Достаточность доказывается конструктивно: если условие утверждения выполнено для E, то в качестве декодирующей функции Dследует взять функцию, возвращающую ближайшее к декодируемому слово из образа E. Необходимость доказывается от противного. Пусть расстояние между выбранными словами в коде равно 2k. Тогда если при передаче каждого из этих слов случится kошибок, которые изменят биты, в которых различаются эти слова, то приемник получит два идентичных сообщения, что свидетельствует о том, что в данной ситуации исправление kошибок невозможно. Следовательно, минимальное расстояние между словами кода должно быть большим 2k.

Пример. Рассмотрим (1,3)-код, состоящий из E, задающей отображение 0
ightarrow000и 1
ightarrow111, и D, задающей отображение 000
ightarrow0, 001
ightarrow0, 010
ightarrow0, 011
ightarrow1,
100
ightarrow0, 101
ightarrow1, 110
ightarrow1, 111
ightarrow1. Этот код (с тройным повторением) исправляет ошибки в одной позиции, т.к. минимальное расстояние между словами кода равно 3.

Если код исправляет все ошибки кратности kи меньшей, то вероятность ошибочного приема слова длины nочевидно не превосходит sum^n_{i=k+1}C^i_np^{n-i}q^i. Вероятность правильного приема в этом случае не меньше, чем

sum^k_{i=0}C^i_np^{n-i}q^i=p^n+C^1_np^{n-1}q+cdots+C^k_np^{n-k}q^k.

Передачу данных часто удобно рассматривать следующим образом. Исходное сообщение a=a_1ldots a_mкодируется функцией Eв кодовое слово b=b_1ldots b_n. Канал связи при передаче добавляет к нему функцией Tстроку ошибок e=e_1ldots e_nтак, что приемник получает сообщение r=r_1ldots r_n, где r_i=b_i+e_i,(1le ile n). Система, исправляющая ошибки, переводит rв некоторое (обычно ближайшее) кодовое слово. Система, только обнаруживающая ошибки, лишь проверяет, является ли принятое слово кодовым, и сигнализирует о наличии ошибки, если это не так.

Пример. Пусть передаваемое слово a=01кодируется словом b=0110, а строка ошибок - e=0010. Тогда будет принято слово r=0100. Система, исправляющая ошибки, переведет его в 0110 и затем восстановит переданное слово 01.

Если система только обнаруживает ошибки и расстояние между любыми кодовыми словами kge2, то любая строка ошибок eс единственной единицей приведет к слову r=b+e, которое не является кодовым.

Пример. Рассмотрим (2,3)-код с проверкой четности. Множество кодовых слов - \{000, 011, 101, 110\}. Ни одна из строк ошибок 001, 010, 100, 111 не переводит одно кодовое слово в другое. Поэтому однократная и тройная ошибки могут быть обнаружены.

Пример. Следующий (2,5)-код обнаруживает две ошибки:

a_1=00
ightarrow00000=b_1,qquad
a_2=01
ightarrow01011=b_2,

a_3=10
ightarrow10101=b_3,qquad
a_4=11
ightarrow11110=b_4.

Этот же код способен исправлять однократную ошибку, потому что любые два кодовых слова отличаются по меньшей мере в трех позициях. Из того, что d(b_i,b_j)ge3при i
e j, следует, что однократная ошибка приведет к приему слова, которое находится на расстоянии 1 от кодового слова, которое было передано. Поэтому схема декодирования, состоящая в том, что принятое слово переводится в ближайшее к нему кодовое, будет исправлять однократную ошибку. В двоичном симметричном канале вероятность правильной передачи одного блока будет не меньше чем p^5+5p^4q.

Установлено 1), что в (n-r,n)-коде, минимальное расстояние между кодовыми словами которого 2k+1, числа n, r(число дополнительных разрядов в кодовых словах) и kдолжны соответствовать неравенству

r ge log_2(C_n^k+C_n^{k-1}+cdots+C_n^1+1),

называемому неравенством или нижней границей Хэмминга. Кроме того, если числа n, rи kсоответствуют неравенству

r >
log_2(C_{n-1}^{2k-1}+C_{n-1}^{2k-2}+cdots+C_{n-1}^1+1),

называемому неравенством или верхней границей Варшамова - Гильберта, то существует (n-r,n)-код, исправляющий все ошибки веса kи менее2).

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

Упражнение 37 Имеется (8,9)-код с проверкой четности. Вычислить вероятность того, что в случае ошибки этот код ее не обнаружит, если вероятность ошибки при передаче каждого бита равна 1\%. Вычислить также вероятность ошибочной передачи без использования кода. Сделать аналогичные расчеты для случая, когда вероятность ошибки в десять раз меньше.

Упражнение 38 Вычислить минимальную и максимальную оценки количества дополнительных разрядов rдля кодовых слов длины n, если требуется, чтобы минимальное расстояние между ними было d. Рассмотреть случаи n=32, d=3и n=23, d=7.

Матричное кодирование

Ранее каждая схема кодирования описывалась таблицами, задающими кодовое слово длины nдля каждого исходного слова длины m. Для блоков большой длины этот способ требует большого объема памяти и поэтому непрактичен. Например, для (16,33)-кода потребуется 33*2^{16}=2,162,688бит.

Гораздо меньшего объема памяти требует матричное кодирование. Пусть Eматрица размерности m	imes n, состоящая из элементов e_{ij}, где i- это номер строки, а j- номер столбца. Каждый из элементов матрицы e_{ij}может быть либо 0, либо 1. Кодирование реализуется операцией b=aEили b_j=a_1e_{1j}+a_2e_{2j}+cdots+a_me_{mj}, где кодовые слова рассматриваются как векторы, т.е как матрицы-строки размера 1	imes n.

Пример. Рассмотрим следующую 3	imes6-матрицу:

E=leftlbrackmatrix{1&0&0&1&1&0cr
                        0&1&0&0&1&1cr
                        0&0&1&1&1&1cr}
ight
brack.

Тогда кодирование задается такими отображениями: 000
ightarrow000000, 001
ightarrow001111, 010
ightarrow010011, 011
ightarrow011100, 100
ightarrow100110, 101
ightarrow101001, 110
ightarrow110101, 111
ightarrow111010.

Рассмотренный пример показывает преимущества матричного кодирования: достаточно запомнить mкодовых слов вместо 2^mслов. Это общий факт.

Кодирование не должно приписывать одно и то же кодовое слово разным исходным сообщениям. Простой способ добиться этого состоит в том, чтобы mстолбцов (в предыдущем примере - первых) матрицы Eобразовывали единичную матрицу. При умножении любого вектора на единичную матрицу получается этот же самый вектор, следовательно, разным векторам-сообщениям будут соответствовать разные вектора систематического кода.

Матричные коды называют также линейными кодами. Для линейных (n-r,n)-кодов с минимальным расстоянием Хэмминга dсуществует нижняя граница Плоткина (Plotkin)1) для минимального количества контрольных разрядов rпри n ge 2d-1,

r ge 2d-2-log_2d.

Упражнение 39 Вычислить минимальную оценку по Плоткину количества дополнительных разрядов rдля кодовых слов матричного кода, если требуется, чтобы минимальное расстояние между ними было d. Рассмотреть случаи из предыдущего упражнения.