离散傅里叶变换(DFT) C代码

#include"math.h" void dft(x,y,a,b,n,sign) /*x-双精度实型一维数组，长度为n。存放要变换数据的实部。 y-双精度实型一维数组，长度为n。存放要变换数据的虚部。 a-双精度实型一维数组，长度为n。存放变换结果的实部。 b-双精度实型一维数组，长度为n。存放变换结果的虚部。 n-整型变量。数据长度。 sign-整型变量。当sign=1时，子函数dft()计算离散傅里叶正变换；当sign=-1时，dft()计 算离散傅里叶反变换。 */ int n,sign; double x[],y[],a[],b[]; { int i,k; double c,d,q,w,s; q = 6.28318530715/n; for(k=0;k<n;k++) { w = k*q; a[k] = b[k] = 0.0; for(i=0;i<n;i++) { d = i*w; c = cos(d); s = sin(d)*sign; a[k] += c*x[i] + s*y[i]; b[k] += c*y[i] - s*x[i]; } } if(sign == -1) { c=1.0/n; for(k=0;k<n;k++) { a[k] = c*a[k]; b[k] = c*b[k]; } } }