www.pudn.com > rfft.rar > rfft.c

 
#include "math.h" 
 
void rfft(double x[],int n) 
{int i,j,k,m,i1,i2,i3,i4,n1,n2,n4; 
 double a,e,cc,ss,xt,t1,t2; 
 for(j=1,i=1;i<16;i++) 
    {m=i; 
     j=2*j; 
     if(j==n)break; 
    } 
    n1=n-1; 
    for(j=0,i=0;i