yulewalk.c

来自「现代谱估计用莱文森-德宾(Levinson-Durbin)算法求解尤利-沃克(Y」· C语言 代码 · 共 19 行

#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);
}