dft.c

c语言实现的离散傅立叶（DFT）变换。非快速傅里叶变换

void      dft()
{
          extern int     inv;
		  extern long    npt;
		  long           k,n;
		  double         WN,wk,c,s,XR[size],XI[size];
          extern complex x[size];

		  WN=2*pi/npt;

		  if(inv==1)
			  WN=-WN;

		  for(k=0;k<npt;++k)
		  {
		      XR[k]=0.0;XI[k]=0.0;
			  wk=k*WN;

			  for(n=0;n<npt;++n)
			  {
                    c=cos(n*wk);s=sin(n*wk);
					XR[k]=XR[k]+x[n+1].real*c+x[n+1].imag*s;
					XI[k]=XI[k]-x[n+1].real*s+x[n+1].imag*c;
			  }
			  if(inv==1)
			  {
			    XR[k]=XR[k]/npt; 
				XI[k]=XI[k]/npt;
			  }
		  }

		  for(k=1;k<=npt;++k)
		  {
		      x[k].real=XR[k-1];
			  x[k].imag=XI[k-1];
		  }

}