Exponenta Banner

Сконфигурируйте EB/No для каналов AWGN с кодированием

Открыть Live Script

Этот пример показывает, как установить отношение энергии битов к плотности шума (Eb/No) для ссылок связи, использующих кодирование канала.

Укажите кодовое слово и длину сообщения для кода Рида-Соломона. Задайте порядок модуляции.

N = 15;       % R-S codeword length in symbols
K = 9;        % R-S message length in symbols
M = 16;       % Modulation order

Создайте (15,9) энкодер Рида-Соломона и модулятор 16-PSK. Задайте объекты так, чтобы они принимали битовые входы.

rsEncoder = comm.RSEncoder('CodewordLength',N,'MessageLength',K, ...
    'BitInput',true);
pskModulator = comm.PSKModulator('ModulationOrder',M,'BitInput',true);

Создайте соответствующий декодер кода Рида-Соломона и 16-PSK объекты демодулятора.

rsDecoder = comm.RSDecoder('CodewordLength',N,'MessageLength',K, ...
    'BitInput',true);
pskDemodulator = comm.PSKDemodulator('ModulationOrder',M,'BitOutput',true);

Вычислите скорость кода Рида-Соломона на основе отношения символов сообщения к длине кодового слова. Определите биты на символ для PSK модулятора.

codeRate = K/N;
bitsPerSymbol = log2(M);

Задайте незакодированные значения Eb/No в дБ. Преобразуйте незакодированный Eb/No в соответствующий закодированный Eb/No, используя кодовую скорость.

UncodedEbNo = 6;
CodedEbNo = UncodedEbNo + 10*log10(codeRate);

Создайте канал AWGN с учетом количества бит на символ. Установите EbNo свойство channel к кодированному Eb/No.

channel = comm.AWGNChannel('BitsPerSymbol',bitsPerSymbol);
channel.EbNo = CodedEbNo;

Установите общее количество ошибок и бит для симуляции. Для точности симуляция должно выполняться до тех пор, пока не встретится достаточное количество битовых ошибок. Количество общих бит используется, чтобы убедиться, что симуляция не работает слишком долго.

totalErrors = 100;
totalBits = 1e6;

Создайте калькулятор частоты ошибок System object™ и инициализируйте вектор частоты ошибок.

errorRate = comm.ErrorRate;
errorVec = zeros(3,1);

Запустите симуляцию, чтобы определить BER. 

while errorVec(2) < totalErrors && errorVec(3) < totalBits
    % Generate random bits
    dataIn = randi([0,1],360,1);
    % Use the RS (15,9) encoder to add error correction capability
    dataEnc = rsEncoder(dataIn);
    % Apply 16-PSK modulation
    txSig = pskModulator(dataIn);
    % Pass the modulated data through the AWGN channel
    rxSig = channel(txSig);
    % Demodulate the received signal
    demodData = pskDemodulator(rxSig);
    % Decode the demodulated data with the RS (15,9) decoder
    dataOut = rsDecoder(demodData);
    % Collect error statistics
    errorVec = errorRate(dataIn,demodData);
end

Отобразите результирующую частоту битовой ошибки.

ber = errorVec(1)
ber = 0.0935

Документация Communications Toolbox

Поддержка