我写的dft程序 不是fft哦

yangzf1986

#include<stdio.h>
#include<math.h>

main()
{int i,j,k,n;
float s;
float ss;
float pi=3.141592653;
float x[128],re[128],im[128],rre[128],rim[128],result[128],real,imag;
for(j=0;j<128;j++)
{x[j]=0.5*sin(2*pi*15*j/100)+2*sin(2*pi*40*j/100);
re[j]=0.0;im[j]=0.0;
rre[j]=0.0;rim[j]=0.0;}
for(k=0;k<128;k++)
{for(n=0;n<128;n++)
{real=cos(2*pi*k*n/128);
imag=-sin(2*pi*k*n/128);
re[n]=x[n]*real;
im[n]=x[n]*imag;
}
s=0.0;ss=0.0;
for(i=0;i<128;i++)
{
s=s+re;
ss=ss+im;}
rre[k]=s;rim[k]=ss;
}
for(i=0;i<128;i++)
  result=sqrt(pow(rre,2)+pow(rim,2));
for(i=0;i<128;i++)
printf("%.4f\n",result);}