郭秀梅《数字信号处理》ch4 FFT.pptVIP

* 常系数线性差分方程 x1(0) x1(1) -1 x1(2) x1(3) -1 x2(0) x2(1) -1 x2(2) x2(3) -1 N/2点 DFT N/2点 DFT x(0) x(7) x(1) x(2) x(3) x(4) x(5) x(6) X1(0)=X(0) X2(0)=X(1) X1(1)=X(2) X1(2)=X(4) X1(3)=X(6) X2(1)=X(3) X2(2)=X(5) X2(3)=X(7) N /2仍为偶数，进一步分解：N /2 N /4 x3(0) x3(1) -1 -1 x4(0) x4(1) N/4点 DFT N/4点 DFT x1(0) x1(1) x1(2) x1(3) X3(0)=X1(0)=X(0) X4(0)=X1(1)=X(2) X3(1)=X1(2)=X(4) X4(1)=X1(3)=X(6) 同理： 其中： 逐级分解，直到2点DFT 当N=8时，即分解到x3(n)，x4(n)，x5(n)，x6(n)，n=0,1 2、算法特点 1）原位计算 -1 L级蝶形运算，每级N/2个蝶形，每个蝶形结构： m表示第m级迭代，k，j表示数据所在的行数 2）蝶形运算 对N=2L点FFT，输入自然序，输出倒位序， 两节点距离：2L-m=N / 2m 第m级运算： 蝶形运算两节点的第一个节点为k值，表示成L位二进制数，左移m-1位，把右边空出的位置补零，结果为r的二进制数。 DIT DIF 基本蝶形不同 DIT: 先复乘后加减 DIF: 先减后复乘 运算量相同 都可原位运算 DIT和DIF的基本蝶形互为转置 § 4.4 DIT与DIF的异同 §4.5 离散傅里叶变换(IDFT)的快速算法 1．改变FFT程序来计算IFFT方法 DFT变换对 比较以上两式知：只要把 再乘以常数 ，则FFT程序可用来计算IDFT 2．不改变FFT程序来计算IFFT方法 ∵ ∴ 共轭 FFT 共轭 乘1/ N §4.8 线性调频Z变换算法 解决的问题：对信号的某一频段感兴趣，也就是只需要计算单位圆上某一段的频谱值，例如对窄带信号，希望在窄带频带内频率的抽样能够非常密集，提高分辨率，带外则不予考虑。如果用DFT方法，则需增加频域抽样点数，增加了窄带之外不需要的计算量。在这种情况下，我们一般采用线性调频 变换(CZT)(chirp-Z变换)。 算法的基本原理： 已知： 是有限长序列，其ZT为： ， ，…， ，… ， 沿 平面上的一段螺线作等分角抽样，抽样点为 ， 、 为任意复数， ， ∴ 抽样点在 平面上所在的周线见右图 ： 由我们所给定的抽样方式及右图，我们有以下结论： ⑴ ：起始抽样点 的矢量半径长度；通常 ，否则 将处于单位圆 的外部。 ⑵ ：表示起始抽样点的相角，可正可负。 ⑶ ：两相邻抽样点之间的角度差， 时， 的路径逆时针旋转， 时， 的路径顺时针旋转。 ⑷ ：其大小表示螺旋线的伸展率 ，随着K的增加螺线内缩 ，表示半径为 的一段圆弧，若 ，则表示单位圆的一部分圆弧。 ，随着 的增加螺线外展（外伸） ⑸最特殊的情况： ， ， （即 ， ），各 均匀等间隔地分布在单位圆上，这就是求序列的DFT。 将 代入ZT表达式： 如果直接计算这一公式，与直接计算DFT相似，总共算 个抽样点，需 次复乘与 次复加。当 、 很大时，运算量是相当可观的。 如何提高运算速度呢？采用布鲁斯坦（Bluestein）提出的等式，可将以上运算转换为卷积和形式。布鲁斯坦提出的等式为： 把它代入 中： 令 则： 上式是 与 的线性卷积，然后再乘上系数 得到，即: 序列 可以想象为频率随时间 成线性增长的复指数序列。在雷达系统中这种信号称为线性调频信号（Chirp Signa