% function pmf_eval(M, mu)

function pmf_eval(M, mu)

% Symbol rate
fs = 6000;
Ts = 1/fs;

% Samples per symbol
N = 4;   

% Number of polyphase taps
Nh1 = 31

% Polyphase upconversion
Nh2 = M*(Nh1+0)

h1 = firrcos(Nh1, 1/Ts, 0.35, N*fs, 'rolloff');
h2 = M*firrcos(Nh2, 1/Ts, 0.35, M*N*fs, 'rolloff');
index = mod(mu,M);

h2a = h2(1+index:M:Nh2);
nh2a = length(h2a)
close all;
sum(h1)
plot(1:Nh1, h1(1:Nh1), '-x', 1:nh2a, h2a(1:nh2a), '-o');
grid;