Документация

График поля точек

Установите фильтр RRC, схему модуляции и параметры графического вывода.

span = 10;          % Filter span
rolloff = 0.2;      % Rolloff factor
sps = 8;            % Samples per symbol
M = 4;              % Modulation alphabet size
k = log2(M);        % Bits/symbol
phOffset = pi/4;    % Phase offset (radians)
n = 1;              % Plot every nth value of the signal
offset = 0;         % Plot every nth value of the signal, starting from offset+1

Создайте коэффициенты фильтра с помощью функции rcosdesign.

filtCoeff = rcosdesign(rolloff,span,sps);

Сгенерируйте случайные символы для размера алфавита M.

rng default
data = randi([0 M-1],5000,1);

Примените модуляцию QPSK.

dataMod = pskmod(data,M,phOffset);

Отфильтруйте модулируемые данные.

txSig = upfirdn(dataMod,filtCoeff,sps);

Вычислите ОСШ для сверхдискретизированного сигнала QPSK.

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

Добавьте AWGN в переданный сигнал.

rxSig = awgn(txSig,snr,'measured');

Применяйтесь RRC получают фильтр.

rxSigFilt = upfirdn(rxSig, filtCoeff,1,sps);

Демодулируйте отфильтрованный сигнал.

dataOut = pskdemod(rxSigFilt,M,phOffset,'gray');

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

h = scatterplot(sqrt(sps)*txSig(sps*span+1:end-sps*span),sps,offset,'g.');
hold on
scatterplot(rxSigFilt(span+1:end-span),n,offset,'kx',h)
scatterplot(dataMod,n,offset,'r*',h)
legend('Transmit Signal','Received Signal','Ideal','location','best')

Индикаторная диаграмма

Отобразите 1 000 точек переданной индикаторной диаграммы сигнала более чем два периода символа.

eyediagram(txSig(sps*span+1:sps*span+1000),2*sps)

Отобразите 1 000 точек полученной индикаторной диаграммы сигнала.

eyediagram(rxSig(sps*span+1:sps*span+1000),2*sps)

Заметьте, что полученная индикаторная диаграмма начинает закрываться из-за присутствия AWGN. Кроме того, фильтр имеет конечную длину, которая также способствует неидеальному поведению.