/*输入序列对应的值*/
      for(i=0;i<size_x;i++)
         scanf("%lf %lf",&x[i].real,&x[i].img);
      initW();                                     
                                       //计算W(0)~W(size_x-1)的值                                                    
      fft();                          //利用fft快速算法进行DFT变化
     output();                       //顺序输出size_x个fft的结果
 
     return   0;  
  }              /*进行基-2 FFT运算,蝶形算法。这个算法的思路就是先把计算过程分为log(size_x)/log(2)-1级(用i控制级数)然后把每一级蝶形单元分组(用j控制组的第一个元素起始下标)最后算出某一级某一组每一个蝶形单元(用k控制个数,共l个)。
*/
  void fft()  
  {  
     int   i=0,j=0,k=0,l=0;                           
     complex   up,down,product;                      
      {                                        
              //算出j组中第k个蝶形单元
                mul(x[j+k+l],W[(size_x/2/l)*k],&product); /*size/2/l是该级W的相邻上标差,l是该级该组取的W总个数*/
                add(x[j+k],product,&up);  
                sub(x[j+k],product,&down);  
                x[j+k]=up;                          
                                         //up为蝶形单元右上方的值
上一篇:Asp.net可视化程序设计课程网站建设+ER图
下一篇:C#+SQL Server物流公司系统管理和人事管理设计

Android手机考勤平台的设计与实现

基于android的环境信息管理系统设计

java+mysql班级评优系统的设计实现

Python+mysql宠物领养平台的设计与实现

ASP.NET飞翔租贷汽车公司信...

基于激光超声检测金属材...

多频激励下典型非线性系统的振动特性研究

基于Joomla平台的计算机学院网站设计与开发

浅谈高校行政管理人员的...

上海居民的社会参与研究

AES算法GPU协处理下分组加...

酵母菌发酵生产天然香料...

STC89C52单片机NRF24L01的无线病房呼叫系统设计

压疮高危人群的标准化中...

提高教育质量,构建大學生...

从政策角度谈黑龙江對俄...

浅论职工思想政治工作茬...