| Новости | Поиск | FAQ | Статьи | Книги | Новости VCL | | Форумы ЧАТ | Орешник | KOL и MCK |
Форумы -» Основная («Начинающим»)/ Базы / WinAPI / Компоненты / Сети / Media / Игры / Corba и COM / KOL / FreePascal / .Net / Прочее / RSDN.ORG

	« предыдущая ветвь | форум | следующая ветвь »

алгоритм БПФ [D5, Win95/98]

Valek ©   (2002-11-20 13:51:00) 

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

---

Ru   (2002-11-22 11:00:00) [1]

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

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

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

---

roadster   (2002-11-20 14:35:00) [2]

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

---

MBo   (2002-11-20 14:25:00) [3]

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

---

Alx2   (2002-11-20 14:14:00) [4]

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

---

MBo   (2002-11-20 13:57:00) [5]

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

---

Наверх