function m=cyclic_mean(x,alpha)%% CYCLIC_MEAN  Calculates the cyclic-mean of a signal%              m(k*alpha)=E{x(t)exp(-j(k*alpha)t)}%              for k = 0 ... 1/alpha%% USAGE%              m=cyclic_mean(x,alpha)%              % File: cyclic_mean.m% Last Revised: 25/11/97% Created: 24/11/97% Author: Andrew C. McCormick% (C) University of Strathclyde% Simple error checksif nargin~=2  error('Incorrect number of arguments for function cyclic_mean');endif alpha>2*pi  error('Cyclic frequency must be less than 2*pi in function cyclic_mean');endT=ceil(2*pi/alpha)-1;t=0:length(x)-1;m=zeros(1,T+1);for k=0:T  m(k+1)=mean(x.*exp(-j*k*alpha*t));end