MATLAB: How to apply the raised cosine filter to an OQPSK signal generated by the OQPSKMOD function using the RCOSFLT function from the Communications Toolbox

Communications Toolboxoqpskpulseshapeshaping

I want to apply a raised cosine filter to pulse-shape an OQPSK signal generated from the OQPSKMOD function. OQPSK signals generated by the OQPSKMOD function inherently up-samples the input samples by 2 to model the half symbol offset. The up-sampling is done by repeating the sample twice.
However, there is a problem when pulse-shaping the OQPSK modulated signal. For example, when I use the following code:
x = randint(100,1,4);
y = real(oqpskmod(x));
rcosflt(y,2,4,'fir/sqrt',0.3,3);
The command pads two zeroes after each y sample and generates an incorrect result.

Best Answer

The OQPSKMOD function produces two samples for each symbol while the RCOSFLT function takes one sample for each symbol. To work around this issue, manually apply a half-symbol offset to the OQPSK signal. Here are the steps:
x = randint(100,1,4);
y = oqpskmod(x);
Fs = 8; Fd = 1;
halfsym = (Fs/Fd)/2;
rolloff = 0.3; delay = 3;
zI = [rcosflt(real(y(1:2:end-1)),Fd,Fs,'normal',rolloff,delay);zeros(halfsym,1)]; % append half a symbol to real part
zQ = [zeros(halfsym,1);rcosflt(imag(y(2:2:end)),Fd,Fs,'normal',0.3,3)]; % offset imaginary part by half a symbol
z = zI + 1j * zQ; % z is the pulse-shaped OQPSK modulated signal
% Plot eye diagram
filtLoading = (delay+1)*Fs/Fd; % Don't display filter loading time
eyediagram(z(filtLoading : end-filtLoading), Fs/Fd);
Related Question