rs = 20e6;
fd = 3;

chan = stdchan('cdmaTUx',rs,fd);
chan.Visualization = 'Frequency response';

data = randi([0 3],10000,1);
txSig = pskmod(data,4,pi/4);

y = chan(txSig);

M = 8; % Modulation order, 8-PSK
Rbit = 9600; % Input bit rate
Rs = Rbit / log2(M); % Symbol rate
Nsamples = 5e2; % Number of samples per frame
Nframes = 10; % Number of frames

v = 10 * 1e3/3600; % Mobile speed (m/s)
fc = 1800e6; % Carrier frequency
c = physconst('LightSpeed'); % Speed of light in free space
fd = v*fc/c; % Maximum Doppler shift of diffuse component

modulator  = comm.PSKModulator(M,'PhaseOffset',0);

channel = stdchan('gsmeqx6',Rs,fd);
channel.RandomStream = 'mt19937ar with seed'; % set for reproducibility
channel.Visualization = 'Impulse and frequency responses';
channel.SamplesToDisplay = '100%';

refC = constellation(modulator);
constDiagram = comm.ConstellationDiagram( ...
    'ReferenceConstellation',refC, ...
    'XLimits',[-3 3],'YLimits',[-3 3]);

for iFrames = 1:Nframes
    msg = randi([0 M-1], Nsamples, 1);
    modSignal = modulator(msg);
    chanOut = channel(modSignal);
    constDiagram(chanOut);
end

release(constDiagram);
release(channel);

v = 120 * 1e3 / 3600; % Mobile speed (m/s)
fd = v*fc/c; % Maximum Doppler shift of diffuse component

channel.MaximumDopplerShift = fd; % Adjust maximum doppler shift

for iFrames = 1:Nframes
    msg = randi([0 M-1], Nsamples, 1);
    modSignal = modulator(msg);
    chanOut = channel(modSignal);
    constDiagram(chanOut);
end