第4章 快速傅里叶变换(FFT)--补充.pptVIP

Bb11-14 Xl11-13 说明 推 导 分 裂 基 FFT 时, 是 按 n 的 奇 偶 划 分 的, 所 以 其 输 入 信 号 排 序 是 基--2 算 法 的 倒 位 序. 分 裂 基 的 基 本 蝶 形 运 算: ( 倒 L 形 ) 运算量：一个分裂基蝶形有两个复乘 分裂基蝶形个数 4.9 线性调频Z变换（Chirp-Z变换）算法 利用FFT算法，可以很快地计算出有限长序列的 DFT 值，也即是Z 变换在单位圆上的全部等间隔采样值。然而，在许多场合，并不一定需要计算全部频谱值，而仅需对某一频带内的信号频谱作较密集的分析。另外，采样也不一定局限于单位圆上，而需要计算出某一螺旋线上的等角度间隔的采样值。例如，在语音信号分析时K采用靠近语音信号序列 Z 变换的极点的螺旋线上进行采样，可以使语音信号的共振峰变得更尖锐，便于精确确定共振峰频率。 Chirp-Z变换算法，在螺旋线上进行采样，可 适 用 于 更 一 般 情 况 下，由 x(n) 求 的 快 速 算 法，这 种 变 换 称 为 线 性 调 频 Z 变 换 ( 简 称 CZT ). 图 4-17 单位圆与非单位圆采样 （a） 沿单位圆采样; (b) 沿AB弧采样 4.6.1 算法基本原理 已知x(n)（0≤n≤N-1）是有限长序列，其Z变换为 （4-28） 为适应z可以沿Z平面更一般的路径取值，故沿Z平面上的一段螺线作等分角的采样，z的这些采样点zk为 zk=AW-k k=0, 1, …, M-1 （4-29） M为所要分析的复频率的点数，即采样点的总数，不一定等于N； A和W都是任意复数，可表示为: 将式（4-30）与式（4-31）代入式（4-29）， 可得 （4-30） （4-31） （4-32） 因此有： 图 4-18 螺线采样 采样点在Z平面上所沿的周线如图4-18所示。由以上讨论和图4-18可以看出： （1）A0表示起始采样点z0的矢量半径长度，通常A0≤1; 否则z0将处于单位圆|z|=1 的外部。 （2）θ0表示起始采样点z0的相角，它可以是正值或负值。 （3）φ0表示两相邻采样点之间的角度差。φ0为正时，表示zk的路径是逆时针旋转的；φ0为负时，表示zk的路径是顺时针旋转的。 （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。 将式（4-29）的zk代入变换表达式（4-28），可得 0≤k≤M-1 (3 - 33) 直接计算这一公式，与直接计算DFT相似，总共算出M个采样点， 需要NM次复数乘法与(N-1)M次复数加法。当N，M很大时， 这个量很大， 这就限制了运算速度。 但是，下面我们将看到，通过一定的变换，以上运算可以转换为卷积形式， 从而可以采用FFT算法，这样就可以大大提高运算速度。 nk可以用以下表达式来替换 将式（4-34）代入式（4-33），可得 如果定义: n=0, 1, …, N-1 则它们的卷积为 式中，k=0, 1, …, M-1。 式（4-38）正好是式（4-35）的一部分， 因此式（4-35）又可以用卷积的形式表示为 k=0, 1, …, M-1 由式（4-39）可看出，如果我们对信号按式（4-36）先进行一次加权处理，加权系数为 ; 然后，通过一个单位脉冲响应为h(n)的线性系统即求g(n)与h(n)的线性卷积; 最后，对该系统的前M点输出再做一次加权，这样就得到了全部M点螺线采样值X(zn)(n=0, 1, …, M-1)。这个过程可以用图4-19 表示。从图中我们看到， 运算的主要部分是由线性系统来完成的。由于系统的单位脉冲响应 可以想象为频率随时间(n)呈线性增长的复指数序列。在雷达系统中,这种信号称为线性调频信号（Chirp Signal），因此，这里的变换称为线性调频Z变换。 图 4-1９ Chirp-Z变换的线性系统表示 k=0, 1, …, M-1 n=0, 1, …, N-1 4.6.2 Chirp-Z变换（CZT）的实现步骤 由式（4-37）可看出，线性系统h(