pskmc.m

it will make the sinosoidal wave forms after you input the values of phase,freq,delay and sampling i

function y=pskmc(x);
ph=0;
y=[];
for k=1:2:length(x)
    if x(k)==0 && x(k+1)==0
       [a,ph]=sino1(2*pi*10,1000,1,0+ph);
    end
    if x(k)==0 && x(k+1)==1
       [a,ph]=sino1(2*pi*10,1000,1,(pi/2)+ph);
    end
    if x(k)==1 && x(k+1)==0
       [a,ph]=sino1(2*pi*10,1000,1,pi+ph);
    end
    if x(k)==1 && x(k+1)==1
       [a,ph]=sino1(2*pi*10,1000,1,(3*pi/2)+ph);
    end
    y=[y a];
end
subplot(2,1,1),stem(x);
subplot(2,1,2),stem(y);