原创 用于NIOS软核处理器的FFT程序

 cen=11;                       //如果n=8,则cen=3;
 zu=n;
 for(i=0;i<cen;i++)       
 {
  zu="zu/2";         //组数(4,2,1)
  for(j=0;j<zu;j++)        //第i 层中的第 j组
  { 
   ge=(n/zu)/2;//一组中蝶形的个数(1,2,4)     
   for(k=0;k<ge;k++)       //第 j 组中的第k个
   {                       //第 a 和第 b 个元素作蝶形运算,WNC
    a="j"*ge*2+k;            //(0,2,4,6)(0,1,4,5)(0,1,2,3)
    b="a"+ge;                //(ge=1,2,4)
    c="k"*zu;                //(0,0,0,0)(0,2,0,2)(0,1,2,3)
    qn=(2048/n)*c;
    //注意sign
    tr=(x*costap[qn]-y*sintap[qn])>>15;//由于正余弦表放大32768倍,此处缩小32768倍
    ti=(y*costap[qn]+x*sintap[qn])>>15;
    p="x"[a];
    q="y"[a];
    //蝶形运算
    x[a]=(p+tr)>>1;//防止数据溢出
    y[a]=(q+ti)>>1;
    x=(p-tr)>>1;
    y=(q-ti)>>1;
   }
  }
 }
}