数字信号处理第二版第三章.pptVIP

3.4 按频率抽取（DIF）的基 -2 FFT算法 k=0, 1, …, N-1 设序列点数N=2M，M为正整数。在把输出X(k)按k的奇偶分 组之前，先把输入序列按前一半、后一半分开（不是按偶数、 奇数分开）， 把N点DFT写成两部分： 3.4.1 算法原理 k=0, 1, …, N-1 （3-23） 当k为偶数时，(-1)k=1；k为奇数时，(-1)k=-1。因此，按k的奇偶可将X(k)分为两部分: （3-24） 又： 故： （3-25） 令: （3-27） 图 3-12 频率抽取法蝶形运算单元 图 3-13 按频率抽取的第一次分解 图 3-14 按频率抽取的第二次分解 图 3-15 按频率抽取的FFT（N=8）信号流图 3.4.2 按频率抽取法的运算特点 式中，m表示m列迭代，k,j为整数所在行数，此式的蝶形运算如图3-16所示，也需要一次复乘和两次复加。 通过(N/2)M个蝶形运算完成，每一个蝶形结构为： DIF法与DIT法的区别： 1、DIF先减后复乘，DIT则是先复乘后再作加减法； 2、频率抽取FFT算法的输入是自然顺序，输出是倒位序。 如果将DIT的基本蝶形加以转置，就得到DIF的基本蝶形; 反过来, 将DIF的基本蝶形加以转置， 就得到DIT的基本蝶形， 因而DIT法与DIF法的基本蝶形是互为转置的。 联系： 3.6 线性调频Z变换（Chirp-Z变换）算法 图 3-17 单位圆与非单位圆采样 （a） 沿单位圆采样; (b) 沿AB弧采样 3.6.1 算法基本原理 （3-28） 沿Z平面上的一段螺线作等分角的采样，z的这些采样点zk为： zk=AW-k k=0, 1, …, M-1 （3-29） M为所要分析的复频率的点数，不一定等于N； A和W都是任意复数，可表示为: 已知x(n)（0≤n≤N-1）是有限长序列，其Z变换为： （3-30） （3-31） 将式（3-30）与式（3-31）代入式（3-29）， 可得 （3-32） 因此有： 图 3-18 螺线采样 （1）A0表示起始采样点z0的矢量半径长度，通常A0≤1; 否则z0将处于单位圆|z|=1 的外部。 结论： （3）φ0表示两相邻采样点之间的角度差。φ0为正时，表示zk的路径是逆时针旋转的；φ0为负时，表示zk的路径是顺时针旋转的。 （2）θ0表示起始采样点z0的相角，它可以是正值或负值。 （4）W0的大小表示螺线的伸展率。W01 时，随着k的增加螺线内缩; W01 时，则随k的增加螺线外伸; W0=1 时，表示是半径为A0的一段圆弧。 若A0=1, 则这段圆弧是单位圆的一部分 当M=N，A=A0ejθ0=1，W=W0·e-jφ0= （W0=1，φ0=2π/N）时，各zk就均匀等间隔地分布在单位圆上，这就是求序列的DFT。 0≤k≤M-1 (3 - 33) 将式（3-29）的zk代入变换表达式（3-28），可得: nk可以用以下表达式来替换: 将式（3-34）代入式（3-33），可得: 如果定义: n=0, 1, …, N-1 则它们的卷积为: 式中，k=0, 1, …, M-1。 因此又可以用卷积的形式表示为: k=0, 1, …, M-1 图 3-1９ Chirp-Z变换的线性系统表示 3.6.2 Chirp-Z变换（CZT）的实现步骤 CZT运算的实现步骤： 0≤n≤N-1 0≤n≤L-1 （3-40） 利用FFT法求此序列的L点DFT: 0≤r≤L-1 （3-41） （2）将 补上零值点， 变为L点序列，有： （1）选择一个最小的整数L，使其满足L≥N+M-1，同时满足 L=2m，以便采用基-2 FFT算法。 0≤n≤M-1 L-N+1≤n≤L-1 0≤n≤M-1 （3-42） 是图3-20（a）的序列 以L为周期的周期延拓序列的主值序列。 （3）形成L点序列h(n) : 用FFT法求其L点DFT: 0≤r≤L-1 （3-43） （4） 将H(r)和G(r)相乘，得Q(r)=H(r)G(r)。 （5）用FFT法求Q(r)的L点IDFT，得h(n)和g(n)的圆周卷积 L （3-44） （6） 求X(zk)： （3-45） （1）形成L点序列 ，需要N次复乘，而系数 可事先准备好，不必在实时分析时计算。 3.6.3 运算量的估计 CZT所需的乘法如下：