数字信号处理--第3章__离散傅里叶变换(DFT).pptVIP

* 数字信号处理 例 3.4.1 对实信号进行谱分析， 要求谱分辨率F≤10 Hz，信号最高频率fc=2.5 kHz， 试确定最小记录时间TPmin， 最大的采样间隔Tmax， 最少的采样点数Nmin。 如果fc不变， 要求谱分辨率增加一倍， 最少的采样点九和最小的记录时间是多少? 解： 因此TPmin=0.1 s， 因为要求fs≥2fc， 所以 * 数字信号处理 2. 用DFT对序列进行谱分析 我们已知道单位圆上的Z变换就是序列傅里叶变换， 即 为使频率分辨率提高一倍， F=5 Hz， 要求 * 数字信号处理 对周期为N的周期序列 ， 由(2.3.10)式知道， 其频谱函数为 用DFT的隐含周期性知道， 截取 的主值序列x(n)= (n)RN(n)， 并进行N点DFT得到 其中 * 数字信号处理 如果截取长度M等于 (n)的整数个周期， 即M=mN， m为正整数， 则 令n=n′+rN, r=0, 1, …, m-1, n′=0, 1, …, N-1,则 * 数字信号处理 因为 k/m=整数 k/m≠整数 * 数字信号处理 如果 的周期预先不知道， 可先截取M进行DFT， 即 k/m=整数 k/m≠整数 再将截取长度扩大一倍， 截取 * 数字信号处理 图 3.4.7 单位圆与非单位圆采样 * 数字信号处理 例如， 要求计算序列在半径为r的圆上的频谱， 那么N个等间隔采样点为 , k=0, 1, 2, …, N-1， zk点的频谱分量为 令 则 (3.4.12) * 数字信号处理 3. ChirpZ变换 设序列x(n)长度为N， 要分析z平面上M点频谱采样值， 分析点为zk, k=0, 1, 2, …, M-1。 设 zk=AW-k, 0≤k≤M-1 式中A和W为复数， 用极坐标形式表示为 (3.4.13) 式中A0和W0为实数。 当k=0时有 * 数字信号处理 将zk代入Z变换公式得到 利用下面的关系式： 得到: 令 * 数字信号处理 (3.4.14) 图 3.4.8 Chrip-Z变换分析频率点分布图 * 数字信号处理 图 3.4.9 Chirp z变换计算框图 * 数字信号处理 图 3.4.10 Chirp-Z变换中hL(n)序列的形成 * 数字信号处理 由(3.4.3)式知， y(n)○h(n)是V(n)的周期延拓序列的主值序列， 延拓周期为L， 即 综上所述， 可归纳出具体计算步骤如下： (1) 形成hL(n)序列 (1) (2) (3) * 数字信号处理 (4) (5) 计算 (6) (7) * 数字信号处理 与标准DFT(FFT)算法相比较， Chirp-Z变换有以下特点： (1) 输入序列长度N和输出序列长度不需要相等， 且二者均可以素数。 (2) 分析频率点zk的起始点z0及相邻两点的夹角φ0是任意的(即频率分辨率是任意的)， 因此可从任意频率上开始， 对输入数据进行窄带高分辨率的谱分析。 (3) 谱分析路径可以是螺旋形的。 (4) 当A=1，M=N， 时， zk均匀分布在单位圆上， 此时Chirp-Z变换就是序列的DFT。 * 数字信号处理 4. 用DFT进行谱分析的误差问题 DFT(实际中用FFT计算)可用来对连续信号和数字信号进行谱分析。 (1) 混叠现象。 (2) 栅栏效应。 (3) 截断效应。 根据傅里叶变换的频域卷积定理有 * 数字信号处理 幅度谱RN(ω)~ω曲线如图3.4.11所示(RN（ω）以2π为周期， 只画低频部分)。 图中，|ω|2π/N的部分称为主瓣， 其余部分称为旁