第四章 离散傅里叶变换及其快速算法.ppt

?求有限长序列x(n)的DFT的实质是: 将有限长序列x(n)作周期延拓x((n))N, 求其DFS, 取其主值序列,即可得X（K） DFT的引出: 对序列的Z变换在单位圆上取值，可得到该序列的傅立叶变换FT，即 ：若令Z=ejw 则可得??所以有限长序列 x(n)在Z平面单位圆上的Z变换是该序列x(n)的傅立叶变换 2. 时域圆周移位定理  y(n)也是一个长度为N的序列,记为: 2、时域循环卷积定理 当N为偶数时， 将上式中的n换成N/2-n可得到 如图所示。 图中*表示对应点为序列取共轭后的值。  圆周共轭对称与共轭反对称序列示意图 ? 任何有限长序列x(n)都可以表示成其共轭对称分量和共轭反对称分量之和， 即   其对应的DFT：  即 XR(k)= XR(N-k) XI(k)= -XI(N-k)  所以 X1(k)=DFT［x1(n)］=1/2［X(k)+X*(N-k)］ X2(k)=DFT［x2(n)］=-j1/2［X(k)-X*(N-k)］ ? 对X（k）求IDFT，得xN(n)=IDFT［X(k)］ ?说明：xN(n)为原序列x (n)以N为周期 周期延拓后的主值序列。 所以，若x (n)为M长，则只有当频域采样点数N≥M时，才有 即可由频域采样X（K）恢复原序列x (n)。否则产生时域混叠现象。这就是所谓的频域采 样定理? 例：已知序列 对x(n)的Z变换X(z)在单位圆上等间隔采样N点，采样值为： 求有限长序列IDFT[X(k)]解：??????? 例：设序列x(n)={8,7,6,5,4,3,2,1},现对x(n)的DTFT 在一个周期内作N=6点的均匀抽样,得XN(K)。求XN(K)的IDFTxN(n)   ????  ? 问题 ：由X（K）能否完整地表达X（Z） X（Z）的内插公式，它说明已知X（K），可根据内插公式求任意Z点的X（Z）。所以X（Z）的N个取样值X(K），包含了X（Z）的全部信息。令： 进一步化简可得 ? 由此可见， 循环卷积既可在时域直接计算，也可以在频域计算。 由于DFT有快速算法FFT， 当N很大时， 在频域计算的速度快得多， 因而常用DFT(FFT)计算循环卷积。 ? 将x(n)、h(n)的长度都变为N1+N2-1,则： 可看出，每一项 的长度为N+M-1，大于序列分段 的长度M 频谱分析存在的问题 例如， x(n)=cos(ω0n), ω0=π/4其频谱为 ??? 1、 直接计算DFT的特点设 x(n)长度为N，其DFT为：X(K),长度也为N 运算时间：设： x(n)为复数序列，计算一个X(K)值需要N次复数乘及N-1次复数加。则计算一个X(K)序列需要N2次复数乘及N（N-1）次复数加。所以直接计算DFT的计算量与N的平方成正比，当N较大时，计算量太大。  应用了上述性质后，独立的值只有N/2个，为简化DFT的运算提供了有力的依据。?  所以一个N点的DFT可分解为两个N/2点的DFT。由于X1(k)和X2(k)均以N/2为周期，且 ，而X(k)为N点所以：? 用信号流图表示：? 那么， X1(k)又可表示为 ?  用同样的方法将x2(r)分解成两个N/4长的子序列x5(l)和x6(l)， N点DFT的第二次时域抽取分解图(N=8)  ??????????   DIT―FFT算法与直接计算DFT运算量的比较 从图可看出当N=2M时，要经过M级蝶算，每一级蝶算包含N/2个蝶形运算，所以总共需要的蝶形运算为：每个蝶形运算需要一次复数乘和两次复数加法。所以N点的FFT需要： 个复数乘， 个复数加。例如，N=210=1024时，DFT与FFT复数乘法运算之比为：? 其中：? 可看出第 m-1级的两节点xm-1(p), xm-1(q) ，经蝶算后，在 m级中的节点序号是不变的，即xm(p), xm(q) 。且蝶算自成独立单元。即xm(p), xm(q) 只与xm-1(p), xm-1(q)有关，而与其它节点值无关，同时xm-1(p), xm-1(q)也不再参加其它的蝶算。所以可以将计算