(数字信号处理)第3章快速傅里叶变换.ppt

(3) 最小记录点数N应满足 取 如果我们事先不知道信号的最高频率，可以根据信号的时域波形图来估计它。例如， 某信号的波形如图 3-23 所示。 先找出相邻的波峰与波谷之间的距离，如图中t1，t2，t3，t4。 然后，选出其中最小的一个如t4。这里, t4可能就是由信号的最高频率分量形成的。 峰与谷之间的距离就是周期的一半。 因此，最高频率为 知道fh后就能确定采样频率 图 3-23 估算信号最高频率fh 3.7.2 可能出现的误差 利用FFT对连续信号进行傅里叶分析时可能造成的误差如下。 1. 频谱混叠失真 在图3-21 画出的基本步骤中，A/D变换前利用前置低通滤波器进行预滤波，使xc(t)频谱中最高频率分量不超过fh。假设A/D变换器的采样频率为fs，按照奈奎斯特采样定理，为了不产生混叠， 必须满足 也就是采样间隔T满足 一般应取 fs=(2.5～3.0)fh 如果不满足fs≥2fh，就会产生频谱混叠失真。 对于FFT来说，频率函数也要采样，变成离散的序列，其采样间隔为F(即频率分辨率)。 由式（3-55）可得 （3-58） （3-59） 从以上Ｔ和tp两个公式来看，信号的最高频率分量fh与频率分辨率F存在矛盾关系，要想fh增加，则时域采样间隔T就一定减小， 而fs就增加，由式（3-57）可知，此时若是固定N，必然要增加F， 即分辨率下降。 3.6 线性调频Z变换（Chirp-Z变换）算法 图 3-17 单位圆与非单位圆采样 （a） 沿单位圆采样; (b) 沿AB弧采样 3.6.1 算法基本原理 已知x(n)（0≤n≤N-1）是有限长序列，其Z变换为 （3-28） 为适应z可以沿Z平面更一般的路径取值，故沿Z平面上的一段螺线作等分角的采样，z的这些采样点zk为 zk=AW-k k=0, 1, …, M-1 （3-29） M为所要分析的复频率的点数，即采样点的总数，不一定等于N； A和W都是任意复数，可表示为: 将式（3-30）与式（3-31）代入式（3-29）， 可得 （3-30） （3-31） （3-32） 因此有： 图 3-18 螺线采样 采样点在Z平面上所沿的周线如图3-18所示。由以上讨论和图3-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。 将式（3-29）的zk代入变换表达式（3-28），可得 0≤k≤M-1 (3 - 33) 直接计算这一公式，与直接计算DFT相似，总共算出M个采样点， 需要NM次复数乘法与(N-1)M次复数加法。当N，M很大时， 这个量很大， 这就限制了运算速度。 但是，下面我们将看到，通过一定的变换，以上运算可以转换为卷积形式， 从而可以采用FFT算法，这样就可以大大提高运算速度。 nk可以用以下表达式来替换 将式（3-34）代入式（3-33），可得 如果定义: n=0, 1, …, N-1 则它们的卷积为 式中，k=0, 1, …, M-1。 式（3-38）正好是式（3-35）的一部分， 因此式（3-35）又可以用卷积的形式表示为 k=0, 1, …, M-1 由式（3-39）可看出，如果我们对信号按式（3-36）先进行一次加权处理，加权系数为 ; 然后，通过一个单位脉冲响应为h(n)的线性系统即求g(n)与h(n)的线性卷积; 最后，对该系统的前M点输出再做一次加权，这样就得到了全部M点螺线采样值X(zn)(n=0, 1, …, M-1)。这个过程可以用图3-19 表示。从图中我们看到， 运算的主要部分是由线性系统来完成的。由于系统的单位脉冲响应 可以想象为频率