【文件属性】:
文件名称:burg 滤波器,用于滤去噪声
文件大小:5KB
文件格式:C
更新时间:2013-06-09 14:09:46
burg
#include
#include
#include
#define N 110000
double x[N];
double xout[N];
double f0[N];
double g0[N];
double f00[N];
double g00[N];
double p0;
double km;
double am[1000];
double amchange[1000];
int m;
double error;
double error1;
double amtry[20]={-1.561735,1.436875,-1.228909,1.114787,-1.286210,1.066792,-0.818686,0.850963,-0.777519,0.624041,-0.474253,0.475683,-0.404298,0.217569,-0.185316,0.154146,-0.082471,-0.043769,0.069651};
//****************************************初始化程序用于计算PO,f0,g0以及确定阶数m**********************************//
double initial ()
{
double z;
int i1,j1;
for(i1=0;i1-0.00005)
{
//controller=0;
printf("success! m = %d\n",m);
}
else
{
/*m=m+1;
kmnew();
update();
update2();*/
}
// }
//printf("m=%d\n",m);
}
/*for(i=20;i