function [x]=mpsk_impulses(M,nsymbols,nsamples)
u=rand(1,nsymbols);
rinteger=round((M*u)+0.5);
phase=pi/M+((rinteger-1)*(2*pi/M));
x=zeros(1,nsymbols*nsamples);
for m=1:nsymbols
   index=(m-1)*nsamples+round(nsamples/2);
   x(1,index)=exp(i*phase(m));
end