循环自相关函数matlab,[求助]循环自相关函数

function R=cyclic_autocorrelation(x,alpha,max_tau)

R=cyclic_cross_correlation(x,x,alpha,max_tau);

function R=cyclic_cross_correlation(x,y,alpha,max_tau)

%

% CYCLIC_CROSS_CORRELATION

%

%              calculates the cyclic cross correlation between

%              two signals x,y at frequency alpha

%

%              R(k*alpha,tau)=E{x(t-tau/2)y(t+tau/2)exp(-jk(alpha)t)}

%              for k=0 ... 2*pi/alpha-1

%

% USAGE

%              R=cyclic_cross_correlation(x,y,alpha,max_tau)

%

%              calculate cross correlation up to max_tau time lags

% File: cyclic_cross_correlation.m

% Simple error checks

if nargin~=4

error('Incorrect number of arguments for function cyclic_cross_correlation');

end

if alpha>2*pi

error(