请问谁做过滤波器，能否提供源代码参考？

[yangbill]

最近正在做数据分析系统，
需要设计数字滤波器，
请问有哪位高手做过，IIR滤波和FIR滤波，
还有窗函数的选择，功率谱分析等数据处理方法
能否提供源代码给我参考一下，
不胜感激。

[NetSpider9804040]

没做过，只能帮楼主顶一下！

[yangbill]

多谢，我查了好久的资料了，
哪位大侠帮帮忙啊。

[constantine]

楼主搞语音的？

[yangbill]

不是。是振动数据采集分析
不知道是否需要数字滤波，
有哪位大侠会的麻烦告诉偶一下啊
摆脱

[GeoPhoenix]

简单的带同滤波就可以了了下面是一段带通滤波的代码，你简单修改一下就可以用了
void __fastcall TMainForm::BandFilter(float HighCut,float LowCut)
{
int i,ch,key;
float *px,*ppx,*py,*ppy;
float tmp;
float f1,f2,w,v,a0,b1,b2,b3,b4;
double Dt;
double e1;
float PI = 3.1415926;
float PP = 1.41421;
int length;
try
{
length = Length;
f2 = HighCut;
f1 = LowCut;
px = new float[length];
if(!px)
{
return ;
}
py = new float[length];
if(!py)
{
delete px;
return;
}
ppx = px;
ppy = py;

if(DataFormat ==8)
Dt = dt*1000000.0;
else
Dt = dt;
e1=-cos(PI*Dt*(f2+f1))/cos(PI*dt*(f2-f1));
w=tan(PI*Dt*(f2-f1));
v=1+PP*w+w*w;
a0=w*w/v;
b1=2*e1*(2+PP*w)/v;
b2=2*(2*e1*e1+1-w*w)/v;
b3=2*e1*(2-PP*w)/v;
b4=(1-PP*w+w*w)/v;
for(ch=0; ch<TotalScan; ch++)
{
px = ppx;
py = ppy;
for(i=0;i<length;i++)
px[i]= dat[ch][i];
px=ppx;
*(py+0)=a0*(*(px+0));
*(py+1)=a0*(*(px+1))-*(py+0)*b1;
*(py+2)=a0*(*(px+2))-2*a0*(*(px+0))-*(py+0)*b2-*(py+1)*b1;
*(py+3)=a0*(*(px+3))-2*a0*(*(px+1))-*(py+0)*b3-*(py+1)*b2-*(py+2)*b1;
py +=4;
px +=4;
for(;py<=ppy+length-1;py++,px++)
{
*py=a0*(*px)-2*a0*(*(px-2))+a0*(*(px-4))-b1*(*(py-1))-b2*\
(*(py-2))-b3*(*(py-3)) -b4*(*(py-4));
}
px=ppx;
py=ppy;
*(px+length-1)=a0*(*(py+length-1));
*(px+length-2)=a0*(*(py+length-2))-b1*(*(px+length-1));
*(px+length-3)=a0*(*(py+length-3))-2*a0*(*(py+length-1))-b1*(*(px+length-2))-b2\
*(*(px+Length-1));
*(px+Length-4)=a0*(*(py+Length-4))-2*a0*(*(py+Length-2))-b1*(*(px+Length-3))-b2\
*(*(px+Length-2))-b3*(*(px+Length-1));
py +=Length-5;
px +=Length-5;
for(;px>=ppx;px--,py--)
{
*px=a0*(*py)-2*a0*(*(py+2))+a0*(*(py+4))-b1*(*(px+1))-b2*
(*(px+2))- b3*(*(px+3))-b4*(*(px+4));
}
px=ppx;
for(i=0;i<Length;i++)
{
tmp = *px++ ;
dat[ch][i] = tmp;
}
}
px = ppx;
py = ppy;
delete px;
delete py;
}
catch(Exception &E)
{
Application->ShowException(&E);
}
}

//这是俺写的一个多道数据滤波的程序，你可根据自己的需要适当修改

[yangbill]

多谢了～～