附錄 2 – C DFT 源程序

　　#include　

　　#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

　　y_re[i]=y_im[i]=0;

　　}

　　for(k=0;k

　　{

　　//　compute　y(k),　k-th　DFT　output

　　for(i=0;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