C语言 1024点快速傅里叶变换（FFT）程序，最好经过优化，执行速度快

void fft()

{

int nn,n1,n2,i,j,k,l,m,s,l1

float ar[1024],ai[1024]// 实部 虚部

float a[2050]

float t1,t2,x,y

float w1,w2,u1,u2,z

float fsin[10]={0.000000,1.000000,0.707107,0.3826834,0.1950903,0.09801713,0.04906767,0.02454123,0.01227154,0.00613588,}// 优化

float fcos[10]={-1.000000,0.000000,0.7071068,0.9238796,0.9807853,0.99518472,0.99879545,0.9996988,0.9999247,0.9999812,}

nn=1024

s=10

n1=nn/2 n2=nn-1

j=1

for(i=1i<=nni++)

{

a[2*i]=ar[i-1]

a[2*i+1]=ai[i-1]

}

for(l=1l<n2l++)

{

if(l<j)

{

t1=a[2*j]

t2=a[2*j+1]

a[2*j]=a[2*l]

a[2*j+1]=a[2*l+1]

a[2*l]=t1

a[2*l+1]=t2

}

k=n1

while (k<j)

{

j=j-k

k=k/2

}

j=j+k

}

for(i=1i<=si++)

{

u1=1

u2=0

m=(1<<i)

k=m>>1

w1=fcos[i-1]

w2=-fsin[i-1]

for(j=1j<=kj++)

{

for(l=jl<nnl=l+m)

{

l1=l+k

t1=a[2*l1]*u1-a[2*l1+1]*u2

t2=a[2*l1]*u2+a[2*l1+1]*u1

a[2*l1]=a[2*l]-t1

a[2*l1+1]=a[2*l+1]-t2

a[2*l]=a[2*l]+t1

a[2*l+1]=a[2*l+1]+t2

}

z=u1*w1-u2*w2

u2=u1*w2+u2*w1

u1=z

}

}

for(i=1i<=nn/2i++)

{

ar[i]=a[2*i+2]/nn

ai[i]=-a[2*i+3]/nn

a[i]=4*sqrt(ar[i]*ar[i]+ai[i]*ai[i]) // 幅值

}

}

FFT 只能对 2 的整数次方的点子 做 运算。 一般要重新采样，或在尾部添数值，凑成 1024，2048，4096，....

FFT，通常要用移动窗对输入数字修匀（时域修匀）或对频谱修匀。

1024点FFT的源程序可以修改源程序，变成2的整数次方点。

另一种办法是把 10240 点 分成 10 组。

第0，10，20，..第一组

第1，11，21，..第二组

第2，12，22，..第三组

....

每组1024点

分别 用 1024点FFT

dt 为 原数据 时间步长 的10 倍。

FFT 后，按常规，用 谱窗 截去 高频 和 特低 频。

然后按频段，对10个 FFT 做 算术平均，作为结果。

用1024点FFT的缺点是分辨率降低了。

还有办法，你把时序数据给有现成程序的人，请贺做人家算一下。（轻而易举闷物的事蚂拍液）。

---