алгоритм БПФ

← →
Valek (2002-11-20 13:51) [0]

Помогите пожалуйста. Как мне написать алгоритм БПФ. Надо для любой функции получить ее спектральное представление и отобразить в виде графика. Заранее спасибо.

← →
MBo © (2002-11-20 13:57) [1]

А что, в Инете не нашел??????

← →
Alx2 © (2002-11-20 14:14) [2]

>MBo © (20.11.02 13:57)
А он там и ищет :)

← →
MBo © (2002-11-20 14:25) [3]

>Alx2
Чужими руками ;)

← →
roadster (2002-11-20 14:35) [4]

Есть у меня, да сидюк дома. Завтра вышлю, не падай духом.

за шо боролись на то и напоролись:

template<class TREAL>
void FFT(int n,TREAL *x,TREAL *y)
{
int n2=n>>1,kd,i1,i;
TREAL Tx,Ty,T;

for(int j=i=0;i<n-1;i++)
{
if(i<j){swap(x[i],x[j]);swap(y[i],y[j]);}
for(kd=n2;j>=kd;kd=(kd+1)>>1)j-=kd;
j+=kd;
}
TREAL Wx,Wy;
for(kd=1;kd<=n2;kd<<=1)
for(j=0;j<kd;j++)
{
Wx=cosl(M_PI/kd*j);
Wy=sinl(M_PI/kd*j);
for(i=j;i<n;i+=kd<<1)
{
i1=i+kd;
Tx=x[i1]*Wx-y[i1]*Wy;
Ty=x[i1]*Wy+y[i1]*Wx;
x[i1]=x[i]-Tx;
y[i1]=y[i]-Ty;
x[i]+=Tx;
y[i]+=Ty;
}
}
for(i=0;i<n;i++)
{
x[i]/=n2;
y[i]/=n2;
}
}

нашел где-то в нете

алгоритм БПФ Найти похожие ветки