Exponenta Banner
exponenta event banner

Передача и получение укороченных кодов Рида-Соломона

Открыть сценарий в реальном времени

Передача и прием стандартных и сокращенных RS-кодированных 64-QAM-modulated данных через канал AWGN. Сравните характеристики стандартного и укороченного кодов.

Задайте параметры кода Рида-Соломона, где N - длина кодового слова, K - номинальная длина сообщения, и S - укороченная длина сообщения. Укажите порядок модуляции, M.

N = 63;  % Codeword length
K = 51;  % Message length
S = 39;  % Shortened message length
M = 64;  % Modulation order

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

numErrors = 200;
numBits = 1e7;
ebnoVec = (8:13)';
[ber0,ber1] = deal(zeros(size(ebnoVec)));

Создание объекта частоты ошибок для сбора статистики об ошибках.

errorRate = comm.ErrorRate;

Создайте пару кодера Рида-Соломона и декодера для кода RS (63,51). Вычислите кодовую скорость.

rsEncoder = comm.RSEncoder(N,K,'BitInput',true);
rsDecoder = comm.RSDecoder(N,K,'BitInput',true);
rate = K/N;

Выполните основной цикл обработки.

for k = 1:length(ebnoVec)
    
    % Convert the coded Eb/No to an SNR. Initialize the error statistics
    % vector.
    snrdB = ebnoVec(k) + 10*log10(rate) + 10*log10(log2(M));
    errorStats = zeros(3,1);
    
    while errorStats(2) < numErrors && errorStats(3) < numBits
        
        % Generate binary data.
        txData = randi([0 1],K*log2(M),1);
        
        % Encode the data.
        encData = rsEncoder(txData);
        
        % Apply 64-QAM modulation.
        txSig = qammod(encData,M, ...
            'UnitAveragePower',true,'InputType','bit');
        
        % Pass the signal through an AWGN channel.
        rxSig = awgn(txSig,snrdB);
        
        % Demodulated the noisy signal.
        demodSig = qamdemod(rxSig,M, ...
            'UnitAveragePower',true,'OutputType','bit');
        
        % Decode the data.
        rxData = rsDecoder(demodSig);
        
        % Compute the error statistics.
        errorStats = errorRate(txData,rxData);
    end
    
    % Save the BER data, and reset the errorRate counter.
    ber0(k) = errorStats(1);
    reset(errorRate)
end

Создайте полином генератора Рида-Соломона для кода RS (63,51).

gp = rsgenpoly(N,K,[],0);

Создание пары кодировщика и декодера Рида-Соломона с использованием укороченной длины сообщения S и генераторный полином gp. Вычислите скорость укороченного кода.

rsEncoder = comm.RSEncoder(N,K,gp,S,'BitInput',true);
rsDecoder = comm.RSDecoder(N,K,gp,S,'BitInput',true);
rate = S/(N-(K-S));

Выполните основной цикл обработки с использованием укороченного кода Рида-Соломона.

for k = 1:length(ebnoVec)
    
    % Convert the coded Eb/No to an SNR. Initialize the error statistics
    % vector.
    snrdB = ebnoVec(k) + 10*log10(rate) + 10*log10(log2(M));
    errorStats = zeros(3,1);
    
    while errorStats(2) < numErrors && errorStats(3) < numBits
        
        % Generate binary data.
        txData = randi([0 1],S*log2(M),1);
        
        % Encode the data.
        encData = rsEncoder(txData);
        
        % Apply 64-QAM modulation.
        txSig = qammod(encData,M, ...
            'UnitAveragePower',true,'InputType','bit');
        
        % Pass the signal through an AWGN channel.
        rxSig = awgn(txSig,snrdB);
        
        % Demodulated the noisy signal.
        demodSig = qamdemod(rxSig,M, ...
            'UnitAveragePower',true,'OutputType','bit');
        
        % Decode the data.
        rxData = rsDecoder(demodSig);
        
        % Compute the error statistics.
        errorStats = errorRate(txData,rxData);
    end
    
    % Save the BER data, and reset the errorRate counter.
    ber1(k) = errorStats(1);
    reset(errorRate)
end

Вычислите приблизительный BER для кода RS (63,51).

berapprox = bercoding(ebnoVec,'RS','hard',N,K,'qam',64);

Сравните кривые BER для кодов RS (63,51) и RS (51,39). Постройте теоретически аппроксимированную кривую BER. Обратите внимание, что сокращение кода не влияет на производительность.

semilogy(ebnoVec,ber0,'o-',ebnoVec,ber1,'c^-',ebnoVec,berapprox,'k--')
legend('RS(63,51)','RS(51,39)','Theory')
xlabel('Eb/No (dB)')
ylabel('Bit Error Rate')
grid

Figure contains an axes. The axes contains 3 objects of type line. These objects represent RS(63,51), RS(51,39), Theory.

Документация по инструментам связи

Поддержка