Модулируйте случайный сигнал

Схема модуляции использует 16-QAM основную полосу, и сигнал проходит через канал аддитивного белого Гауссова шума (AWGN). Основные операции симуляции используют их функции MATLAB® и Communications Toolbox™.

Сгенерируйте случайный поток двоичных данных

Обычный формат для представления сигнала в MATLAB является вектором или матрицей. randi функция создает вектор-столбец, содержащий значения потока двоичных данных. Длина потока двоичных данных (то есть, количество строк в вектор-столбце) произвольно установлена в 30 000.

Задайте параметры.

M = 16; % Modulation order (alphabet size or number of points in signal constellation)
k = log2(M); % Number of bits per symbol
n = 30000; % Number of bits to process
sps = 1; % Number of samples per symbol (oversampling factor)

Установите rng функционируйте к его состоянию по умолчанию или любому статическому значению seed, так, чтобы пример привел к повторяемым результатам. Затем используйте randi функция, чтобы сгенерировать случайные двоичные данные.

rng default;
dataIn = randi([0 1],n,1); % Generate vector of binary data

Используйте диаграмму стебель-листья, чтобы показать двоичные значения для первых 40 битов случайного потока двоичных данных. Используйте двоеточие (:) оператор в вызове stem функционируйте, чтобы выбрать фрагмент бинарного вектора.

stem(dataIn(1:40),'filled');
title('Random Bits');
xlabel('Bit Index');
ylabel('Binary Value');

Преобразуйте двоичный сигнал в сигнал с целочисленным знаком

Настройка по умолчанию для qammod функция ожидает, что данные с целочисленным знаком как вводимые символы будут модулировать. В этом примере поток двоичных данных предварительно обрабатывается в целочисленные значения перед использованием qammod функция. В частности, bi2de функция преобразует каждого с 4 кортежами в соответствующее целое число в области значений [0, (M–1)]. Порядок модуляции, M, 16 в этом примере.

Выполните отображение бита к символу первым изменением данных в бинарный k- кортежи, где k количество битов на символ, заданный $$k=lo{g}_2(M)$$. Затем используйте bi2de функционируйте, чтобы преобразовать каждого с 4 кортежами в целочисленное значение.

dataInMatrix = reshape(dataIn,length(dataIn)/k,k);
dataSymbolsIn = bi2de(dataInMatrix);

Постройте первые 10 символов в диаграмме стебель-листья.

figure; % Create new figure window.
stem(dataSymbolsIn(1:10));
title('Random Symbols');
xlabel('Symbol Index');
ylabel('Integer Value');

Модулируйте Используя 16-QAM

Используйте qammod функция, чтобы применить 16-QAM модуляцию к dataSymbolsIn вектор-столбец для естественно закодированных и серо-закодированных отображений бита к символу.

dataMod = qammod(dataSymbolsIn,M,'bin'); % Binary coding with phase offset of zero
dataModG = qammod(dataSymbolsIn,M); % Gray coding with phase offset of zero

Операция модуляции выходные параметры объединяет вектор-столбцы, содержащие значения, которые являются элементами 16-QAM сигнального созвездия. Позже в этом созвездии в качестве примера схемы показывают естественное и двоичное отображение символа Грэя.

Для получения дополнительной информации о функциях модуляции смотрите Цифровую Модуляцию. Для примера, который использует Грэя, кодирующего с модуляцией манипулирования сдвига фазы (PSK), смотрите Грэя Закодированный 8-PSK.

Добавьте белый Гауссов шум

Модулируемый сигнал проходит через канал при помощи awgn функция с заданным отношением сигнал-шум (SNR). Преобразуйте отношение энергии на бит к шумовой степени спектральная плотность (${\mathit{E}}_{\mathrm{b}}/{\mathit{N}}_0$) к значению ОСШ для использования awgn функция. sps переменная не является значительной в этом примере, но делает расширение примера, чтобы использовать импульс, формирующий легче. Для получения дополнительной информации смотрите, что Импульс Использования Формирует на 16-QAM примере Сигнала.

Вычислите ОСШ, когда канал будет иметь ${\mathit{E}}_{\mathrm{b}}/{\mathit{N}}_0$ из 10 дБ.

EbNo = 10;
snr = EbNo+10*log10(k)-10*log10(sps);

Передайте сигнал через канал AWGN для двоичного файла, и Грэй закодировал отображения символа.

receivedSignal = awgn(dataMod,snr,'measured');
receivedSignalG = awgn(dataModG,snr,'measured');

Создайте схему созвездия

Используйте scatterplot функционируйте, чтобы отобразить синфазные компоненты и квадратурные компоненты модулируемого сигнала, dataMod, и сигнал с шумом, полученный после канала. Эффекты AWGN присутствуют в схеме созвездия.

sPlotFig = scatterplot(receivedSignal,1,0,'g.');
hold on
scatterplot(dataMod,1,0,'k*',sPlotFig)

Демодулируйте 16-QAM

Используйте qamdemod функция, чтобы демодулировать полученные данные и вывести символы данных с целочисленным знаком.

dataSymbolsOut = qamdemod(receivedSignal,M,'bin');
dataSymbolsOutG = qamdemod(receivedSignalG,M);

Преобразуйте сигнал с целочисленным знаком в двоичный сигнал

Используйте de2bi функционируйте, чтобы преобразовать символы данных от демодулятора QAM, dataSymbolsOut, в бинарную матрицу, dataOutMatrix, с размерностями ${\mathit{N}}_{\mathrm{sym}}$–${\mathit{N}}_{\mathrm{bits}/\mathrm{sym}}$. В матрице, ${\mathit{N}}_{\mathrm{sym}}$ общее количество символов QAM, и ${\mathit{N}}_{\mathrm{bits}/\mathrm{sym}}$ количество битов на символ. Для 16-QAM, ${\mathit{N}}_{\mathrm{bits}/\mathrm{sym}}$ = 4. Затем преобразуйте матрицу в вектор-столбец длины, равной количеству входных битов, 30,000. Повторите процесс для серо-закодированных символов данных, dataSymbolsOutG.

Инвертируйте отображение бита к символу, выполняемое ранее в этом примере.

dataOutMatrix = de2bi(dataSymbolsOut,k);
dataOut = dataOutMatrix(:); % Return data in column vector
dataOutMatrixG = de2bi(dataSymbolsOutG,k);
dataOutG = dataOutMatrixG(:); % Return data in column vector

Вычислите системный BER

biterr функция вычисляет статистику битовой ошибки от исходного потока двоичных данных, dataIn, и полученные потоки данных, dataOut и dataOutG. Серое кодирование значительно уменьшает BER.

Используйте функцию коэффициента ошибок, чтобы вычислить ошибочную статистику. Используйте fprintf функция, чтобы отобразить результаты.

[numErrors,ber] = biterr(dataIn,dataOut);
fprintf('\nThe binary coding bit error rate is %5.2e, based on %d errors.\n', ...
    ber,numErrors)

The binary coding bit error rate is 2.40e-03, based on 72 errors.

[numErrorsG,berG] = biterr(dataIn,dataOutG);
fprintf('\nThe Gray coding bit error rate is %5.2e, based on %d errors.\n', ...
    berG,numErrorsG)

The Gray coding bit error rate is 1.33e-03, based on 40 errors.

Постройте сигнальные созвездия

Схема созвездия, показанная ранее, построила точки в созвездии QAM, но это не указало на отображение между значениями символа и точками созвездия. В этом разделе схема созвездия указывает на отображения для естественного кодирования и серо-кодирования двоичных данных к точкам созвездия.

Покажите естественное и серое закодированное бинарное отображение символа для 16-QAM созвездия

Примените 16-QAM модуляцию к полным наборам точек созвездия при помощи естественно закодированного бинарного отображения символа и серо-закодированного отображения символа.

M = 16; % Modulation order
x = (0:15); % Integer input
symbin = qammod(x,M,'bin'); % 16-QAM output (natural-coded binary)
symgray = qammod(x,M,'gray'); % 16-QAM output (Gray-coded)

Используйте scatterplot функционируйте, чтобы построить схему созвездия и аннотировать ее естественным (красным) и Грэем (черные) бинарные представления точек созвездия.

scatterplot(symgray,1,0,'b*');
for k = 1:M
    text(real(symgray(k)) - 0.0,imag(symgray(k)) + 0.3,...
        dec2base(x(k),2,4));
     text(real(symgray(k)) - 0.5,imag(symgray(k)) + 0.3,...
         num2str(x(k)));
    
    text(real(symbin(k)) - 0.0,imag(symbin(k)) - 0.3,...
        dec2base(x(k),2,4),'Color',[1 0 0]);
    text(real(symbin(k)) - 0.5,imag(symbin(k)) - 0.3,...
        num2str(x(k)),'Color',[1 0 0]);
end
title('16-QAM Symbol Mapping')
axis([-4 4 -4 4])

Исследуйте графики

Используя серо-закодированный символ отображение улучшает производительность BER, потому что серо-закодированные точки сигнального созвездия отличаются только на один бит от каждой смежной соседней точки. Где как с естественным бинарным отображением символа, некоторые смежные точки созвездия. Например, естественно закодированные двоичные значения для 1 (0 0 0 1) и 2 (0 0 1 0) отличайтесь на два бита (третье и дальше биты).