График рассеяния

Установите фильтр 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);

Вычислите SNR для сигнала 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);
hold on
scatterplot(rxSigFilt(span+1:end-span),n,offset,'bx',h)
scatterplot(dataMod,n,offset,'r+',h)
legend('Transmit Signal','Received Signal','Ideal','location','best')

Диаграмма глаз

Отображение 1000 точек диаграммы переданного сигнального глаза в течение двух символьных периодов.

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

Отображение 1000 точек диаграммы принимаемого сигнального глаза.

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

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