第3章-离散傅里叶变换幻灯片.ppt

混 叠 现 象 的 结 论 由?f = fs / N ≥2fmax /N 看出： 在 N 给 定 时, 为 避 免混 叠 失 真 而 一 味 提 高 抽 样 频 率 fs ，必 然 导 致?f增 加, 即 频 率 分 辨 力 下 降; 反 之, 若 要 提 高 频 率 分 辨 力 即 减 小 ?f , 则 导 致 减 小fs, 最 终 必 须 减 小 信 号 的 高 频 容 量. 以 上 两 点 结 论 都 是 在记录长度内抽样点数 N 给 定 的 条 件 下 得 到 的. 所 以 在 高 频 容 量 fmax 与 频 率 分 辨 力 F 参 数 中, 保 持 其 中 一 个 不 变 而 使 另 一 个 性 能 得 以 提 高 的 唯 一 办 法, 就 是 增 加 记 录 长 度 内 的 点 数 N, 即 fmax 和 F 都 给 定 时, 则 N 必 须 满 足 N ≥2fmax /F 这是未采用任何特殊数据处理（例如加窗）情况下，为实现基本DFT算法所必须满足条件。 例子 1 有一频谱分析仪用的FFT处理器，其抽样点数必须是2的整数幂。假定没有采用任何特殊的数据处理措施，已给条件为： (1)频率分辨力≤10Hz (2)信号的最高频率≤4kHz 试确定以下参量： (1) 最小记录长度Tp； (2)抽样点的最大时间间隔T； (3)在一个记录中的最少点数N。 解： (1)由分辨力的要求确定最小记录长度Tp. Tp=1/ ?f =1/10=0.1(s) 故最小记录长度为0.1秒。 (2)从信号的最高频率确定最大的抽样时间间隔T。 fs≥2fmax, T=1/fs ≤1/2fmax=0.125*10-3 (s) (3)最小记录点数N，它应满足 N≥2fmax / ?f =800 该处理器所需最少采样点数为N=210=1024点。（因为N=29=512点不够） 频 谱 泄 漏 在实际中，要把观测的信号x(n)限制在一定的时间间隔之内，即采取截断数据的过程。 时域的截断在数学上的意义为原连续时间信号乘上一个窗函数，使原连续时间函数成为两端突然截断，中间为原信号与窗函数相乘的结果。 时域两函数相乘，在频域是其频谱的卷积.由于窗函数不可能取无限宽,即其频谱不可能为一冲激函数,信号的频谱与窗函数的卷积必然产生拖尾现象。造成 频谱泄漏. 所 以 在 截 取 (即 在 窗 函 数 的 选 取) 时, 应 尽 量 选 择 适 当 形 状 的 窗 函 数 对时域信号进行截断, 使频谱泄漏最小. 频 谱 泄 漏 注 意 点 由于我们无法取无数个点，所以在DFT时，时域的截断是必然的，因而泄漏也是必然存在的。 为了减少频率泄漏可采用： 适当加大窗口宽度，增加M值； 采用适当形状的窗函数截断 例子 2 设信号为x(n)=1/2π,经过矩形窗函数截断，求信号经过矩形窗函数前后的频谱函数。 解：设信号经过矩形窗函数后的信号为x1(n),矩形窗函数为W(n),其频谱函数为X1(ejw) x1(n)=x(n)W(n) X1(ejw)=X(ejw)*W (ejw) 很明显： X1(ejw)≠ X(ejw) 相当于X(ejw)失真，这种失真是由于X(ejw)的频谱泄漏引起，其现象为“拖尾（扩展现象），称之频谱泄漏。 因为X(ejw)=δ(w),矩形窗函数 w X(ejw) X1(ejw) w 产生泄漏 栅 栏 效 应 利 用 DFT 逼 近 连 续 时 间 信 号 的 傅 里 叶 变 换, 其 频 谱 将 不 再 是 连 续 函 数 而 是 基 频 F 的 整 数 倍。 用 DFT 计 算 频 谱, 就 如 通 过 一 个 栅栏观 看 一 个 景 色, 只 能 在 离 散 点 的 地 方 看 到 真 实 的 景 象, 从 而 产 生 栅 栏 效 应. 如 果 在 两 离 散 的 谱 线 间 频 谱 有 很 大 变 化, 不 作 特 殊 处 理, 则 无 法 将 其 检 测 出 来. 减 小 栅 栏 效 应方 法 减 小 栅 栏 效 应 的 一 个 方 法 是 在 所 取 数 据 的 末 端 加 一 些 零 值 点, 使 一 个 周 期 内 点 数 增 加, 但 是 不 改 变 原 有 的 记 录 数 据. 这种方法 等 效 于 加 长 了 周 期 Tp . 因 公 式 F = 1/ Tp (F是 抽 样 间 隔). Tp 增 加, 抽 样 间 隔 变 小, 从 而 能 保 持 原 来 频 谱 形 式 不 变 的 情 况 下 使 谱 线 变 密, 也 就