Количество битовых ошибок и вероятность битовой ошибки (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.
Демодулируйте принятый сигнал.
Преобразуйте демодулированные символы в двоичные данные.
Вычислим количество битовых ошибок.
The 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. Например, см. случай (с) на рисунке.
