Количество битовых ошибок и частота битовых ошибок (BER)
[ сравнивает беззнаковое двоичное представление элементов в number,ratio] = biterr(x,y)x для тех, кто в y. Функция возвращает number, количество битов, которые различаются в сравнении, и ratio, отношение number к общему количеству битов. Функция определяет порядок, в котором она сравнивает x и y исходя из их размеров. Дополнительные сведения см. в разделе Алгоритмы.
[ возвращает результат двоичного сравнения 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-арные символы.
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 является неотрицательным целым числом в данном случае. Например, см. вариант (а) на рисунке.
Если один вход является матрицей, а другой - вектором столбца, то функция сравнивает каждый столбец элемента матрицы по элементам с вектором столбца. Число строк в матрице должно быть равно длине вектора столбца. Другими словами, если матрица имеет размеры m-by-n, то вектор столбца должен иметь размеры m-by-1. Например, см. вариант (b) на рисунке.
Если один вход является матрицей, а другой - вектором строки, то функция сравнивает каждую строку элемента матрицы по элементам с вектором строки. Число столбцов в матрице должно быть равно длине вектора строки. Другими словами, если матрица имеет размеры m-by-n, то вектор строки должен иметь размеры 1-by-n. Например, см. вариант (c) на рисунке.
