расчет коэффициента ошибок по битам 15,11 граф Хемминга - PullRequest
Новая
       98

расчет коэффициента ошибок по битам 15,11 граф Хемминга

0 голосов
/ 18 октября 2018

Ниже в моем коде предполагается найти симулятор BER.Но я получаю ошибку в этом коде, где nErrors = biterr(dataIn,dataDec2) эта строка дает несоответствие размеров матрицы.

Кто-нибудь может мне помочь в этом ?? 

close all;
clear all;

M = 2;                 % Modulation order
k = log2(M);            % Bits per symbol
EbNoVec = -4:2:0;      % Eb/No values (dB)
No = -10;
numSymPerFrame = 100;   % Number of PSK symbols per frame
berEst2 = zeros(size(EbNoVec));
G2=[1 1 0 0 1 0 0 0 0 0 0 0 0 0 0; 
0 1 1 0 0 1 0 0 0 0 0 0 0 0 0; 
0 0 1 1 0 0 1 0 0 0 0 0 0 0 0; 
1 0 1 0 0 0 0 1 0 0 0 0 0 0 0;
1 0 0 1 0 0 0 0 1 0 0 0 0 0 0;
0 1 0 1 0 0 0 0 0 1 0 0 0 0 0;
1 1 1 0 0 0 0 0 0 0 1 0 0 0 0;
0 1 1 1 0 0 0 0 0 0 0 1 0 0 0;
1 0 1 1 0 0 0 0 0 0 0 0 1 0 0;
1 1 0 1 0 0 0 0 0 0 0 0 0 1 0;
1 1 1 1 0 0 0 0 0 0 0 0 0 0 1   ];
H2= gen2par(G2);
decoding2 = syndtable(H2); 
Pt2 = zeros(size(EbNoVec));
for n = 1:length(EbNoVec)
% Convert Eb/No to SNR
snrdB = EbNoVec(n) + 10*log10(15/11);
% Reset the error and bit counters
numErrs = 0;
numBits = 0;
Pt2(n)= 10^((snrdB-10)/10);
while numErrs < 100
    % Generate binary data and convert to symbols
    dataIn = randi([0 1],numSymPerFrame,k)
    dataSym = bi2de(dataIn)
    dataEnc2 = encode(dataIn,15,11,'linear/binary',G2)
    % PSK modulation

    txSig = pskmod(dataEnc2,M);

    % Pass through AWGN channel
    rxSig = awgn(txSig,snrdB,'measured');

    % Demodulate the noisy signal
    rxSym = pskdemod(rxSig,M);
    % Convert received symbols to bits
    dataOut = de2bi(rxSym,k);
    dataDec2 =decode(rxSym,15,11,'linear/binary',G2, decoding2);
    % Calculate the number of bit errors
    nErrors = biterr(dataIn,dataDec2);

    % Increment the error and bit counters
    numErrs = numErrs + nErrors;
    numBits = numBits + numSymPerFrame*k;
end

% Estimate the BER
berEst2(n) = numErrs/numBits;


end
berTheory2 = berawgn(EbNoVec,'psk',M,'nondiff');
hold on 
semilogy(EbNoVec,berEst1,'r','LineWidth',2);

1 Ответ

0 голосов
/ 20 октября 2018

Я отладил ваш код с точкой останова в строке 46 [nErrors = biterr (dataIn, dataDec2)].Кажется, массив dataIn имеет размер 100 x 1, а dataDec2 - массив 110 x 1.Функция biterr вычисляет количество мест, где два вектора одинаковой длины различаются.

Это даст вам лучшее объяснение.https://www.mathworks.com/help/comm/ref/biterr.html

...