Exponenta Banner

Уровень 2/3 Сверточный код в AWGN

Открыть скрипт

Этот пример генерирует немного коэффициента ошибок по сравнению с кривой Eb/No для ссылки, которая использует 16-QAM модуляцию и уровень 2/3 сверточный код в AWGN.

Установите порядок модуляции и вычислите количество битов на символ.

M = 16;
k = log2(M);

Создайте решетку для уровня 2/3 сверточный код. Установите traceback и параметры уровня кода.

trellis = poly2trellis([5 4],[23 35 0; 0 5 13]);
traceBack = 16;
codeRate = 2/3;

Создайте сверточный энкодер и его эквивалентный Декодер Витерби.

convEncoder = comm.ConvolutionalEncoder('TrellisStructure',trellis);
vitDecoder = comm.ViterbiDecoder('TrellisStructure',trellis, ...
    'InputFormat','Hard','TracebackDepth',traceBack);

Создайте 16-QAM модулятор и пару демодулятора, имеющую битные вводы и выводы. Установите совокупность нормировать на средней степени.

qamModulator = comm.RectangularQAMModulator('BitInput',true,'NormalizationMethod','Average power');
qamDemodulator = comm.RectangularQAMDemodulator('BitOutput',true,'NormalizationMethod','Average power');

Создайте объект коэффициента ошибок. Установите задержку получателя с дважды traceback глубиной, которая является задержкой через декодер.

errorRate = comm.ErrorRate('ReceiveDelay',2*traceBack);

Установите область значений значений Eb/No моделироваться. Инициализируйте матрицу статистики частоты ошибок по битам.

ebnoVec = 0:2:10;
errorStats = zeros(length(ebnoVec),3);

Моделируйте ссылку путем выполнения этих шагов:

  • Сгенерируйте двоичные данные.

  • Закодируйте данные с уровнем 2/3 сверточный код.

  • Модулируйте закодированные данные.

  • Передайте сигнал через канал AWGN.

  • Демодулируйте полученный сигнал.

  • Декодируйте демодулируемый сигнал при помощи Декодера Витерби.

  • Соберите ошибочные статистические данные.

for m = 1:length(ebnoVec)

    snr = ebnoVec(m) + 10*log10(k*codeRate);

    while errorStats(m,2) <= 100 && errorStats(m,3) <= 1e7
        dataIn = randi([0 1],10000,1);
        dataEnc = convEncoder(dataIn);
        txSig = qamModulator(dataEnc);
        rxSig = awgn(txSig,snr);
        demodSig = qamDemodulator(rxSig);
        dataOut = vitDecoder(demodSig);
        errorStats(m,:) = errorRate(dataIn,dataOut);
    end
    reset(errorRate)

end

Вычислите теоретический BER по сравнению с кривой Eb/No для случая без прямого кодирования с коррекцией ошибок.

berUncoded = berawgn(ebnoVec','qam',M);

Постройте BER по сравнению с кривой Eb/No для моделируемых закодированных данных и теоретических незакодированных данных.

semilogy(ebnoVec,[errorStats(:,1) berUncoded])
grid
legend('Coded','Uncoded')
xlabel('Eb/No (dB)')
ylabel('Bit Error Rate')

В выше значениях Eb/No код с коррекцией ошибок предоставляет выигрыши в производительности.

Документация Communications Toolbox
Поддержка