yulewalk.c

现代谱估计用莱文森-德宾(Levinson-Durbin)算法求解尤利-沃克(Yule-Walker)方程。 形参说明： r:双精度实型一维数组

#include"stdlib.h"
#include"levinson.c"
void yulewalk(double x[],int n,int p,double a[],double *v)
{
	int i,k;
	double s,*r;
	r=(double *)malloc((p+1)*sizeof(double));
	for(k=0;k<=p;k++)
	{
		s=0.0;
		for(i=0;i<(n-k);i++)
		{
			s+=x[i]*x[i+k];
		}
		r[k]=s/n;
	}
	levinson(r,p,a,v);
	free(r);
}