% Comparison of 16-MPSK and 16-QASK
clear all;
clc;

M = 16; x = wavread('ABHI.wav',4000);
x1 = x*128+128 ;% 1000 random quarternary data is generated
y=dec2bin(x1,8);
n=1;i=0;
for i=1:size(x)
    for j=1:4:8
             for l=0:1:3
                t(1+l)=y(i,j+l);
             end
        x1(n)=bin2dec(t);
        n=n+1;
    end
end

y1 = pskmod(x1,M); % 16-MPSK modulation
y1n = awgn(y1,5,'measured'); % SNR of 10 dB
%scatterplot(y1n);
y1r=pskdemod(y1n,M); %16-MPSK demodulation
[num_error,ber]=symerr(x1,y1r) % no semicolon, result echoed  

 n=1;
y2=dec2bin(y1r,2);
for i=1:2:2*size(x)
    for l=0:1:1
            t(1+l*2)=y2(i+l,1);
            t(2+l*2)=y2(i+l,2);
            t(3+l*2)=y2(i+l,3);
            t(4+l*2)=y2(i+l,4);
    end
        x1r(n)=bin2dec(t);
        n=n+1;
           end
x1m = (x1r-128)/128;

wavwrite(x1m,'result_16mpsk.wav');