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

Установите фильтр 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);
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. Кроме того, фильтр имеет конечную длину, что также способствует неидеальному поведению.