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');