基于CooleyTukey的FFT算法的递归C实现Linux认证考试

文章作者 100test 发表时间 2010:04:10 18:43:38
来源 100Test.Com百考试题网

　　信号处理课程需要用到FFT变换，于是找资料写了一个C程序，网上的资料不少，但是感觉没什么大用，基本上大部分是MATLAB写的，在 Introduction to Algorithm的第三十章有一个算法伪代码图，小心wn这个是每一次都变的而不是一直是2pi/N,我犯了这个错误，一直不大对，后来在 cprogramming上的以为朋友指出了这个错误，之后犯的错误是结果一直都是不对称，找来找去不知道错在哪里，第二天才发现关键错在

　　re_w = re_w*W_RE - im_w*W_IM. im_w = im_w*W_RE re_w*W_IM.

　　应该替换为

　　save_re = re_w.save_im = im_w.re_w = save_re*W_RE(len_x) - save_im*W_IM(len_x).im_w = save_im*W_RE(len_x) save_re*W_IM(len_x).

　　好了，贴上我的代码，大家有什么问题欢迎讨论

　　#include