**3.1.1.1 Defining parameters**

M = 16; % Size of signal constellation k = log2(M); % Number of bits per symbol n = 3e4; % Number of bits to process nsamp = 1; % Oversampling rate
