Здравствуйте! Появилась необходимость в реализации сверточного кода со скоростью 2/3. Для образования кода использую генерацию со скоростью 1/2 с помощью двух полиномов 133 и 171. Далее провожу перфорацию данного кода путем отбрасывания кадого четвертого бита. Результирующая скорость кода при этом возрастает до требуемых 2/3. Для проверки работы кода создал модель в Matlab. Проверил - декодирует, даже исправляет ошибки). Но радость была недолгой - обнаружилось, что хотя код в состоянии справить и 5 ошибок, он не может исправить даже одну, если она расположена в "хвосте" закодированной последовательности. Размер исходного кадра - 40 бит, закодированного - 60 бит. Длина "хвоста" примерно 5-8 бит. Привожу код модели.
Код
dataIn = randi([0, 1], 1, 40); % входной вектор
rate = 2/3;
tblen = 20;
ber = 2; % вероятность ошибок в канале, в процентах
trellis = poly2trellis(7,[133,171]); % создание решетки
%% encoder - кодируем на скорости 1/2
dataEncoded = convenc(dataIn, trellis, 0);
%% perforation - выбрасываем каждый 4-ый бит
dataEncoded(4:4:length(dataEncoded))=0.5;
dataEncoded = dataEncoded(find(dataEncoded ~= 0.5));
%% distortion - вносим искажения в случайные позиции исходя из заданной вероятности ошибки
numer = round((length(dataIn) / rate) * ber / 100);
errors = randi([1, length(dataIn) / rate], 1, numer);
%errors = 57; % ФАТАЛИТИ! Можно закомментировать верхнюю строчку и определить позицию ошибки вручную.
dataEncoded(errors) = 1 - dataEncoded(errors);
%% decoder - декодируем перфорированный сигнал
[dataDecoded] = vitdec(dataEncoded, trellis, tblen, 'trunc', 'hard', [1 1 1 0]);
%% result - подсчитываем ошибки
numErrors = sum(dataIn ~= dataDecoded);
res = find(dataIn ~= dataDecoded);
fprintf('ber = %d\n', numErrors);
fprintf('pos = ');
fprintf('%d, ', res);
fprintf('\n');
Может кто-нибудь подсказать, в чем может быть причина такого поведения?
Сверточный код 2/3 не может исправить одну ошибку, Перфорированный сверточный код 2/3 не может исправить одну ошибку
Oct 20 2015, 13:27
