mqam-demodulation.zip_MQAM

function [ t,x ] = dmqam( y,fs,fd,fc,M,ini_phase ) FsDFd=fs/fd;offset=0; [num,den] = butter(5,fc*2/fs ); t = (0:1/fs:(size(y,1)-1)/fs)'; x = [y .* cos(pi*2 * fc * t + ini_phase) y.*sin(pi*2*fc*t + ini_phase)]; wid_x = size(x, 2); y = x(:,1:2); for i = 1 : wid_x y(:, i) = filter(num, den, x(:, i)) *2; end [I, Q] = qaskenco(M); sizey = size(y); y = integ(y, FsDFd, offset); y = repmat(y(:), 1, FsDFd); y = reshape(y.', sizey(1), sizey(2)); x = demodmap(y, fd, fs, 'qask/arb', I, Q); end