附录 2 – C DFT 源程序
　　#include　<math.h>
　　#define　PI　3.1415926535897932384
　　void　dft(double　x[],double　a[],double　phi[],int　*m)
　　{
　　double　y_re[NMAX],　y_im[NMAX],　t,　s,　term_re,　term_im;
　　int　i,j,k,n=*m;
　　for(i=0;i<n;++i)　{
　　y_re[i]=y_im[i]=0;
　　}
　　for(k=0;k<n;++k)
　　{
　　//　compute　y(k),　k-th　DFT　output
　　for(i=0;i<n;++i)
　　{
　　//　compute　i-th　term　of　y(k):
　　//　x(k)*exp(-2*pi*I*(k-1)*(i-1)/n)
　　//　compute　real　and　imaginary　parts　of　i-th　term
　　//　using　exp(I*t)=exp(t)*(cos(t)+I*sin(t))
　　t=-2.*PI*k*i/(double)n;
　　term_re=x[i]*exp(t)*cos(t);
　　term_im=x[i]*exp(t)*sin(t);
　　//　add　term　to　sum
　　y_re[k]+=term_re;
　　y_im[k]+=term_im;
　　}
　　}
　　//　transform　y　to　polar　coordinates
　　for(k=0;k<n;++k)
　　{
　　//　compute　amplitude　of　y(k)
　　a[k]=sqrt(y_re[k]*y_re[k]+y_im[k]*y_im[k]);
　　//　compute　phase　of　y(k)
　　phi[k]=atan2(y_im[k],y_re[k]);
　　}
　　}
　　//　initialize　input　data
　　void　init(double　a[],int　*m)
　　{
　　int　j,n=*m;
　　for(j=0;j<n;++j)
　　{
　　a[j]=sin(1./sqrt((double)j+1.0));
　　}
　　}
　　//　Dummy　function　to　use　result,　preventing　compiler　from
　　//　optimizing　away　the　computation.
　　void　consume(double　a[],double　b[],double　c[])
　　{
　　}

上一页1234下一页