Количество битовых ошибок и частоты ошибок по битам (BER)
[ сравнивает бинарное представление без знака элементов в number,ratio] = biterr(x,y)x тем в y. Функция возвращает number, количество битов, которые отличаются по сравнению и ratio, отношение number к общему количеству битов. Функция определяет порядок, в котором она сравнивает x и y на основе их размеров. Для получения дополнительной информации смотрите раздел Algorithms.
[ возвращает бинарный результат сравнения number,ratio,individual] = biterr(___)x и y как матричный individual. Можно задать любую комбинацию входных аргументов от предыдущих синтаксисов.
Создайте две бинарных матрицы.
x = [0 0; 0 0; 0 0; 0 0]
x = 4×2
0 0
0 0
0 0
0 0
y = [0 0; 0 0; 0 0; 1 1]
y = 4×2
0 0
0 0
0 0
1 1
Определите количество битовых ошибок.
numerrs = biterr(x,y)
numerrs = 2
Вычислите количество постолбцовых ошибок.
numerrs = biterr(x,y,[],'column-wise')numerrs = 1×2
1 1
Вычислите количество построчных ошибок.
numerrs = biterr(x,y,[],'row-wise')numerrs = 4×1
0
0
0
2
Вычислите количество полных ошибок. Поведение совпадает с поведением по умолчанию.
numerrs = biterr(x,y,[],'overall')numerrs = 2
Демодулируйте шумный 64-QAM сигнал и оцените частоту ошибок по битам (BER) для области значений значений Eb/No. Сравните оценку BER с теоретическими значениями.
Установите параметры симуляции.
M = 64; % Modulation order k = log2(M); % Bits per symbol EbNoVec = (5:15)'; % Eb/No values (dB) numSymPerFrame = 100; % Number of QAM symbols per frame
Инициализируйте вектор результатов.
berEst = zeros(size(EbNoVec));
Основной цикл обработки выполняет эти шаги.
Сгенерируйте двоичные данные и преобразуйте в 64-ary символы.
QAM-модулируйте символы данных.
Передайте модулируемый сигнал через канал AWGN.
Демодулируйте полученный сигнал.
Преобразуйте демодулируемые символы в двоичные данные.
Вычислите количество битовых ошибок.
while цикл продолжает обрабатывать данные, пока с или 200 ошибками не сталкиваются или 1e7, биты передаются.
for n = 1:length(EbNoVec) % Convert Eb/No to SNR snrdB = EbNoVec(n) + 10*log10(k); % Reset the error and bit counters numErrs = 0; numBits = 0; while numErrs < 200 && numBits < 1e7 % Generate binary data and convert to symbols dataIn = randi([0 1],numSymPerFrame,k); dataSym = bi2de(dataIn); % QAM modulate using 'Gray' symbol mapping txSig = qammod(dataSym,M); % Pass through AWGN channel rxSig = awgn(txSig,snrdB,'measured'); % Demodulate the noisy signal rxSym = qamdemod(rxSig,M); % Convert received symbols to bits dataOut = de2bi(rxSym,k); % Calculate the number of bit errors nErrors = biterr(dataIn,dataOut); % Increment the error and bit counters numErrs = numErrs + nErrors; numBits = numBits + numSymPerFrame*k; end % Estimate the BER berEst(n) = numErrs/numBits; end
Определите теоретическую кривую BER при помощи berawgn функция.
berTheory = berawgn(EbNoVec,'qam',M);Отобразите предполагаемые и теоретические данные о BER на графике. Предполагаемые точки данных BER хорошо выравниваются с теоретической кривой.
semilogy(EbNoVec,berEst,'*') hold on semilogy(EbNoVec,berTheory) grid legend('Estimated BER','Theoretical BER') xlabel('Eb/No (dB)') ylabel('Bit Error Rate')
Функция использует размеры x и y определить порядок, в котором это сравнивает их элементы.
Если входные параметры являются матрицами тех же размерностей, то функция сравнивает входные параметры поэлементно. number неотрицательное целое число в этом случае. Например, смотрите случай (a) на рисунке.
Если один вход является матрицей, и другой вход является вектор-столбцом, то функция сравнивает каждый столбец элемента матрицы элементом с вектор-столбцом. Количество строк в матрице должно быть равно длине вектор-столбца. Другими словами, если матрица имеет размерности m-by-n, то вектор-столбец должен иметь размерности m-by-1. Например, смотрите случай (b) на рисунке.
Если один вход является матрицей, и другой вход является вектором-строкой, то функция сравнивает каждую строку элемента матрицы элементом с вектором-строкой. Количество столбцов в матрице должно быть равно длине вектора-строки. Другими словами, если матрица имеет размерности m-by-n, то вектор-строка должен иметь размерности 1 n. Например, смотрите случай (c) на рисунке.
