Уровень 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;
Создайте сверточный энкодер и его эквивалентный декодер Viterbi.
convEncoder = comm.ConvolutionalEncoder('TrellisStructure',trellis); vitDecoder = comm.ViterbiDecoder('TrellisStructure',trellis, ... 'InputFormat','Hard','TracebackDepth',traceBack);
Создайте объект частоты ошибок. Установите задержку приемника в два раза больше глубины следа, чем задержка через декодер.
errorRate = comm.ErrorRate('ReceiveDelay',2*traceBack);Установите область значений значений Eb/No, которые будут моделироваться. Инициализируйте матрицу статистики вероятности битовой ошибки.
ebnoVec = 0:2:10; errorStats = zeros(length(ebnoVec),3);
Симулируйте ссылку следующими шагами:
Сгенерируйте двоичные данные.
Закодируйте данные сверточным кодом со скоростью 2/3.
16-QAM модулировать закодированные данные, сконфигурируйте битовые входы и среднюю степень модуля.
Передайте сигнал через канал AWGN.
16-QAM демодулируйте принятый сигнал, сконфигурируйте битовые выходы и среднюю степень модуля.
Декодируйте демодулированный сигнал с помощью декодера Viterbi.
Соберите статистику ошибок.
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 = qammod(dataEnc,M, ... 'InputType','bit','UnitAveragePower',true); rxSig = awgn(txSig,snr); demodSig = qamdemod(rxSig,M, ... 'OutputType','bit','UnitAveragePower',true); dataOut = vitDecoder(demodSig); errorStats(m,:) = errorRate(dataIn,dataOut); end reset(errorRate) end
Вычислите теоретическую кривую BER по сравнению с кривой Eb/No для случая без кодирования прямой коррекции ошибок.
berUncoded = berawgn(ebnoVec','qam',M);Постройте график BER по сравнению с кривой Eb/No для моделируемых кодированных данных и теоретических незакодированных данных. При более высоких значениях Eb/No код исправления ошибок обеспечивает преимущества эффективности.
semilogy(ebnoVec,[errorStats(:,1) berUncoded]) grid legend('Coded','Uncoded') xlabel('Eb/No (dB)') ylabel('Bit Error Rate')
