auting.h

求自相关算法的源代码

double auting(double *x,int nt,int k)
{
	double rl,b;
	int j;
	rl=0.0;
	b=(x[0]*x[k]+x[k+nt-1]*x[nt-1])/2;
	for(j=0;j<nt;j++)
		rl=rl+x[j]*x[j+k];
	rl=(rl-b)/(2.0*nt-2.0);
	return rl;
}