Модулируйте 8-PSK сигнал, добавляйте белый гауссов шум и постройте график сигнала для визуализации эффектов шума.

Создайте object™ системы модулятора M-PSK. Порядок модуляции по умолчанию для объекта - 8.

pskModulator = comm.PSKModulator;

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

modData = pskModulator(randi([0 7],2000,1));

Добавление белого гауссова шума к модулированному сигналу путем пропускания сигнала через канал аддитивного белого гауссова шума (AWGN).

channel = comm.AWGNChannel('EbNo',20,'BitsPerSymbol',3);

Передача сигнала по каналу AWGN.

channelOutput = channel(modData);

Постройте график бесшумных и шумных данных, используя графики рассеяния для визуализации эффектов шума.

scatterplot(modData)

scatterplot(channelOutput)

Изменить EbNo свойство до 10 дБ для увеличения шума.

channel.EbNo = 10;

Передача модулированных данных через канал AWGN.

channelOutput = channel(modData);

Постройте график выхода канала. Вы можете видеть эффекты повышенного шума.

scatterplot(channelOutput)

Создайте 16-PSK модулятор и демодулятор System objects™ которые используют пользовательское отображение символов. Оцените BER в канале AWGN и сравните производительность с теоретической системой PSK с серым кодированием.

Создайте пользовательское сопоставление символов для схемы модуляции 16-PSK. 16 целых символов должны иметь значения от 0 до 15.

custMap = [0 2 4 6 8 10 12 14 15 13 11 9 7 5 3 1];

Создайте пару модулятора 16-PSK и демодулятора с пользовательским отображением символов, определенным массивом. custMap.

pskModulator = comm.PSKModulator(16,'BitInput',true, ...
    'SymbolMapping','Custom','CustomSymbolMapping',custMap);
pskDemodulator = comm.PSKDemodulator(16,'BitOutput',true, ...
    'SymbolMapping','Custom','CustomSymbolMapping',custMap);

Отображение совокупности модуляторов.

constellation(pskModulator)

Создайте системный объект канала AWGN для использования с 16-кратными данными.

awgnChannel = comm.AWGNChannel('BitsPerSymbol',log2(16));

Создайте объект частоты ошибок для отслеживания статистики BER.

errorRate = comm.ErrorRate;

Инициализируйте векторы моделирования. Величина ${\mathit{}}_{\mathrm{}}{\mathit{}}_{\mathrm{Eb/N0}}$ варьируется от 6 до 18 дБ в 1 дБ.

ebnoVec = 6:18;
ber = zeros(size(ebnoVec));

Оцените BER, модулируя двоичные данные, передавая их через канал AWGN, демодулируя принятый сигнал и собирая статистику ошибок.

for k = 1:length(ebnoVec)
    
    % Reset the error counter for each Eb/No value
    reset(errorRate)
    % Reset the array used to collect the error statistics
    errVec = [0 0 0];
    % Set the channel Eb/No
    awgnChannel.EbNo = ebnoVec(k);
    
    while errVec(2) < 200 && errVec(3) < 1e7
        % Generate a 1000-symbol frame
        data = randi([0 1],4000,1);
        % Modulate the binary data
        modData = pskModulator(data);
        % Pass the modulated data through the AWGN channel
        rxSig = awgnChannel(modData);
        % Demodulate the received signal
        rxData = pskDemodulator(rxSig);
        % Collect the error statistics
        errVec = errorRate(data,rxData);
    end
    
    % Save the BER data
    ber(k) = errVec(1);
end

Генерирование теоретических данных BER для канала AWGN с использованием berawgn функция.

berTheory = berawgn(ebnoVec,'psk',16,'nondiff');

Постройте график смоделированных и теоретических результатов. Характеристики BER 16-PSK модуляции моделируемого пользовательского отображения символов не так хороши, как теоретическая кривая прогнозирования для кодов Грея.

figure
semilogy(ebnoVec,[ber; berTheory])
xlabel('Eb/No (dB)')
ylabel('BER')
grid
legend('Simulation','Theory','location','ne')