高效FFT算法深度解析与应用探究.pptVIP

§7-3按频率抽取的FFT算法;一、离散傅立叶反变换（IDFT）的快速变换（IFFT）;*;*;二、N为组合数的FFT算法;*;*;*;*;例：讨论一种长N＝18点序列x（n）的DFT快速计算问题。

解：因为N＝18＝3×3×2

令p1＝3，q1＝6

首先将原始序列分成p1＝3个子序列，每个子序列长度q1=6个取样

序列1（L=0）;x(0),x(3),x(6),x(9),x(12),x(15)

序列2（L=1）;x(1),x(4),x(7),x(10),x(13),x(16)

序列3（L=p1-1=2）;x(2),x(5),x(8),x(11),x(14),x(17);*;*;三、实数序列的FFT;一种N点FFT算法可同步计算两个N点实数序列的DFT;*;*;;一种N点FFT一次计算一种2N点实数序列的DFT;*;*;*;;;;k=8; % 所以k=8

n1=[0:1:19]; %向量n1为0,1,2,…,19，截取长度N为20

xa1=sin(2*pi*n1/k); %生成正弦周期序列xa1

subplot(2,2,1) %指定画图左上角

plot(n1,xa1) %画图xa1

xlabel(‘t/T’);ylabel(‘x(n)’); %标x,y轴说明

xk1=fft(xa1);xk