合工大 数字信号处理 第四章.ppt

第六节 N为复合数的FFT算法  ——混合基算法 第九节 线性调频 z变换（CZT)算法 FFT不适用于： 1、算法原理 N点有限长序列，其z变换： 抽样点： 求抽样点处的z变换： 2、CZT的实现步骤及运算量的估算 3）形成L点序列h(n): 总运算量: 4）求乘积 3、CZT算法的优点 1) N，M可为任意数，可不等 1、 整数的多基多进制表示形式 （1）二进制： （2）r进制： （3）多基多进制（混合基）： 例： 2、 的快速算法 行 列 行序号 列序号 行 列 行变量 列变量 的DFT 算法 (1) 改写 成 (2) 做 个 点DFT ，得 为参量，输入变量 ，输出变量 的 点 DFT (3) N个 （旋转因子） (4) 做 个 点DFT，得 为参量，输入变量 ，输出变量 的 点DFT (5) 整序 例 当N为高组合素数时： 个 点DFT，乘以旋转因子 个 点DFT 个 点DFT，乘以旋转因子 个 点DFT，乘以旋转因子 L级r点DFT 称基 算法， 基 算法 混合基算法（基 算法） 基 算法 混合基算法的运算量 不计译序、整序工作量 （2）乘N个旋转因子 复乘 N 总计： （1） 个 点DFT 复乘 复加 （3） 个 点DFT 复乘 复加 混合基节省的运算量 次乘N个旋转因子 个 点DFT 个 点DFT 个 点DFT 第七节 基 -4FFT算法 当混合基FFT算法中 时， 即为基-4FFT算法，n、k都为4进制数 个 点DFT 乘N个旋转因子 个 点DFT 乘N个旋转因子 个 点DFT 1) 的4点DFT 的四进制数 按二进制倒位序排列成 3) 的4点 DFT 一个4点FFT不需乘法，只需3次乘旋转因子（ 除外） 而基 -2FFT 基- 4FFT运算量： 每级有N/4个4点FFT，共L级（L-1级要乘旋转因子） 第八节 分裂基FFT算法 对偶序号使用基-2FFT算法，对奇序号使用基-4FFT算法，称分裂基FFT算法 针对 的算法中具有最少乘法次数，且同址运算。 将 分 成三个序列。 偶序号的 点DFT 奇序号的 点DFT 利用周期性 分成四段： 的第一级分解得4个分裂基 同样对 、 、 作进一步分解。 和 的第二级分解分别是基-4的4点DFT。 的第二级分解得2个分裂基。 一个基-4的4点DFT和2个基-2的4点DFT。 基-2，基-4等基本碟形结都没有乘法，只有每个分裂基有两次复乘。 运算量： 分裂基碟形数： 只研究信号的某一频段，要求对该频段抽样密集，提高分辨率； 研究非单位圆上的抽样值； 需要准确计算N点DFT，且N为大的素数；等等。 CZT算法：对z变换采用螺线抽样 chirp-z变换：线性调频 z变换 沿z平面上的一段螺线作等分角抽样，抽样点zk： 其中： M为要分析的复频谱点数 则 ：起始抽样点的矢量半径长度 ：起始抽样点的相角 ：相邻抽样点的角度差 ： 逆时针 ：顺时针 ：螺线的伸展率 W01：螺线内缩 W01： 螺线外伸 当W0=1，则表示半径为A0的一段圆弧 若又有A0=1，则表示单位圆上的一段圆弧 若又有 ，M=N ，即为序列的DFT。 NM次复乘 (N-1) M次复加 FFT算法与直接计算DFT所需乘法次数的比较曲线 4、DIT―FFT的运算规律及编程思想 FFT的每级（列）计算都是由N个复数数据（输入）两两构成一个蝶型（共N/2个蝶形）运算而得到另外N个复数数据（输出）。 当数据输入到存储器以后，每一组运算的结果，仍