快速Fourier变换教案.pptVIP

* * * IDFT的快速算法 算法 1 在DFT的FFT流程图中: IDFT的快速算法 算法 2 return FFT应用举例 分析过去300年的太阳黑子活动 具有周期性，每11年有一次最大值出现 数据：苏黎世太阳黑子相对数，同时记录了太阳黑子的数量和大小。 * load sunspot.dat year=sunspot(:,1); relNums=sunspot(:,2); plot(year,relNums) title(Sunspot Data) * 前50组数据 plot(year(1:50),relNums(1:50),b.-); * 利用FFT做频域分析 Y = fft(relNums); Y(1)=[]; %First component is the sum of data plot(Y,ro) title(Fourier Coefficients in the Complex Plane); xlabel(Real Axis); ylabel(Imaginary Axis); * 难以理解，考虑其功率(Complex magnitude squared Y)和频率之间的关系，即画出周期图 * plot(freq(1:40),power(1:40)) xlabel(cycles/year)