N为复合数的FFT算法课件.ppt

* 第10讲 快速傅里叶变换 4.5 N为复合数的FFT算法 4.6 线性调频Z变换（Chirp-Z变换）算法 FFT算法用于计算有限长序列的z变换X(z)在z平面单位圆上N个等间隔抽样点zk上的采样值 实际应用中在很多情况下并不一定需要计算全部频谱值,而仅需对某一频带内的信号频谱作较密集的分析 另外，采样也不一定局限于单位圆上，而需要计算出某一螺旋线上的等角度间隔的采样值 例如语音信号分析时，往往在靠近语音信号序列z变换的极点的螺旋线上进行采样，可以使语音信号的共振峰变得更尖锐，便于精确确定共振峰频率 4.6 线性调频Z变换（Chirp-Z变换）算法 线性调频z变换就是利用FFT快速计算螺旋线采样的算法 沿单位圆采样 沿AB弧采样 4.6.1 算法基本原理 为适应z可以沿Z平面更一般的路径取值，沿Z平面上的一段螺线作等分角的采样，采样点zk为 设有限长序列x(n)的Z变换为 式中: A和W都是任意复数, 设 zk=AW-k k=0, 1, …, M-1 M为采样点的总数 综合上式得到 zk参数的物理意义 螺线采样 A0: 起始采样点z0的矢量半径 θ0: 起始采样点z0的相角 φ0: 两相邻采样点之间的角度差 W0: 表示螺线的伸展率 序列的DFT是Chirp-Z变换的特例 0≤k≤M-1 将zk=AW-k代入x(n)的Z变换变换表达式中,得到 Chirp-Z变换公式 Chirp-Z变换的FFT 算法 直接计算Chirp-Z变换公式的计算量很大, 但可以将其转换为卷积形式, 从而利用FFT算法, 提高运算速度 将公式中的nk做一变换 Chirp-Z变换公式的卷积形式 0≤k≤M-1 代入公式, 整理后得到 (卷积形式 ) 与卷积公式比较 当 X(zk)可表示成卷积形式 k=0, 1, …, M-1 卷积形式的Chirp-Z公式可采用FFT算法，从而提高运算速度 Chirp-Z变换的计算框图如下 4.6.2 Chirp-Z变换的特点 输入序列和输出序列长度不需要相等, 且二者均可为素数 分析频率点zk的起始点z0及相邻两点的夹角φ0是任意的(即频率分辨率是任意的)， 因此可从任意频率上开始， 对输入数据进行窄带高分辨率的谱分析 谱分析路径可以是螺旋形的 Chirp-Z变换在一定条件下就是序列的DFT 与标准DFT(FFT)算法相比较， Chirp-Z变换有以下特点: 4.7 利用FFT分析时域连续信号频谱 4.7.1 基本步骤 时域连续信号离散傅里叶分析的处理步骤如下图示 频谱分析是指计算信号各个频率的幅值, 相位和功率 处理过程分析 LPF : 避免在模拟信号转换成序列时, 可能出现的频谱混叠现象 A/D: 时域离散化,得到采样序列x(n)，其频谱用X(ejω)表示 DFT(FFT)运算: 加窗后的DFT是 0≤k≤N-1 w(n): 窗函数.为进行FFT，须对x(n)加窗处理，即v(n)=x(n)w(n) 对连续信号进行谱分析时，主要关心两个问题：即谱分析范围与频率分辨率 谱分析范围与频率分辨率 谱分析范围(所能分析的最高频率)受采样频率的限制 频率分辨率用频率采样间隔F描述，表示谱分析中能够分辨的两个频谱分量最小间隔。F较小时，称频率分辨率较高 有限长序列v(n)=x(n)w(n)的DFT相当于v(n)FT的等间隔采样 V(k)是sc(t)的离散频率函数， V(k)的第k点对应的模拟频率为: 谱分析范围与频率分辨率 显然，数字域频率间隔Δω=2π/N 对应的模拟域谱线间距应为 ( fs 为采样频率) 谱线间距: 即频谱分辨率，指可分辨两频率的最小间距 ( fs 为采样频率) T为采样时间间隔(s) fs为采样频率(Hz) F为谱线间距,频谱分辨率(Hz) tp为截取连续时间信号的样本长度, 又称记录长度(s) DFT对连续信号谱分析时，参数间的关系及选择原则 例 时间序列v(n)=(1.1)nR16(n)与16点的DFT 的如下图所示 参数间的关系 参数的选择(T、tp和N) 实际应用中, 根据信号最高频率fh和频谱分辨率F的要求来确定三个参数 DFT对连续信号谱分析时，参数间的关系及选择原则 采样周期T的选择： 为保证采样信号不失真应满足，fs≥2fh，即周期T 说明：信号的谱分析范围 fh与频率分辨率F存在矛盾关系 解释：如果保持采样点数N不变，而提高谱的分辨率(F减小)， 必须降低采样速率fs ，从而引起谱分析范围 fh 减少 解决方法：如维持fs不变, 为提高分辨率，可以考虑增加采样点数N与观察时间Tp N的选择（由频谱分辨率F和T确定） 最小记录长度Tp的选择（由N,