频率细化分析方法及其在虚拟仪器中的应用.pdfVIP

书 2005 年 8 月 重庆大学学报（自然科学版） Aug. 2005 第 28 卷第 8 期 JournaI of Chongging University（NOturOI Science Edition） VoI. 28 No. 8 文章编号：1000 - 582X（2005）08 - 0001 - 03 频率细化分析方法及其在虚拟仪器中的应用 ＊ 周 传 德，秦 树 人，汤 宝 平，尹 爱 军 （重庆大学 机械工程学院，重庆 400030） 摘 要：频率细化分析是一种重要的频谱分析工具，是信号处理和虚拟仪器的主要功能之一. 介绍 了几种典型的频率细化分析方法，包括连续 DFT 和复调制方法等，并研究了它们的优缺点，提出了一种 结合子带分解 DFT 和复调制方法的新频率细化方法，它继承了复调制细化的优点，且无需设计特定的 滤波器，减少了运算量，适用于实时细化分析，最后以虚拟仪器中的应用实例展示了几种频率细化的算 法和效果. 关键词：频率细化；虚拟仪器；子带分解 DFT 中图分类号：Th115 文献标识码：A 频率细化分析是 20 世纪 70 年代发展起来的一种 重要的频谱分析手段，也是目前虚拟仪器中 FFT 分析 仪必备的一项功能，进行频率细化分析的主要目的在 于提高频谱的分辨力，得到比一般 FFT 分析更加详尽 频率成分 信 息，或 者 更 为 准 确 的 频 率、幅 值 和 相 位 信息. 1 FFT 误差分析 FFT 是对式（1）的快速算法，在式（1）中 I 标识时 间系列，I 标识频率系列，且 I，I 皆为 I，I（0，N - 1）内 的整数. 这就会产生如下 2 个问题： X（I）= z N -1 I = 0 x（I）WII N . （1） 1）频率分辨力. 当多频信号中的 2 个频率 f 1 、f 2 之差 I f 1 - f 2 I ! f 时，信号 FFT 幅值谱线在 f 1 和 f 2 之 间不可能有谷点，所以分辨不出 2 个频率. 特别对于窄 带信号，这一现象表现得尤为明显. FFT 变换中与频率 相关的变量 K 的步长最小为 1，存在频率分辨率不高 的问题. 2）栅栏效应. 当信号频率对不准某一条谱线时， 通常用邻近的一条谱线来近似. 这会对幅值和相位造 成很大误差，特别是相位误差可达 90 . 栅栏效应主要 受频率分辨率和窗函数的影响. 克服栅栏效应的措施主要包括加窗处理和校正， 文献［1 - 2］研究了频率 - 幅值 - 相位的校正办法；提 高频率分辨力的理论途径主要有 2 种：降低采样频率 f S 和增大采样点数 N，常见的信号处理方法为频率细 化分析，其中复调制细化和相位补偿细化是 2 种目前 普遍应用的频率细化技术. 2 连续 DFT 频率细化 对式（1），X（I）是对连续傅里叶变换频域采样的 结果. 如果对 I 取“连续”变化，则 X（I）也该是“连续” 变化，即实现细化，比如 A I = 0. 1，输出 X（I）系列则是 DFT 直接算法的基础上细化 10 倍. 根据傅里叶级数的 概念也易得到相同的结果［3］. 连续 DFT 频率细化的物理意义明确，对采样数据 长度没有严格要求，可以进行任意倍数的细化. 但细化 倍数太大或者采样点数太少，受加窗的影响严重，效果 较差. 因此适用于数据样本适中，细化倍数不大的离线 或在线测试. 3 复调频 FFT 频率细化 复调制细化是根据傅里叶变换的频移性质，将采 样时 间 序 列 通 过 特 定 的 带 通 滤 波 器 后，与 复 指 数 e - 2 ! If c / f S相乘，将中心频率 f c 的频段信息移到相应频谱 的原点处，再经 M 倍下采样并作 N M 点 FFT，实现 M 倍 频率细化. 算法流程如图 1. 图 1 ZOOM - FFT 频率细化处理流程 ＊ 收稿日期：2005 - 04 - 28 基金项目：国家自然科学基金资助项目 作者简介：周传德（1978 - ），男，湖北宜昌人，重庆大学博士研究生，研究方向信号处理与智能测试. 复调制法 ZFFT 利用带通滤波加频移获得任何感 兴趣频段的细化谱思路，是先进和实用的，但此法要求 事先设计带通滤波器，且一旦设计好以后，其参数不好 随时变更，滤波器也会使相位和幅值发生改变，计算量 也较大. 可以通过特定滤波器（ 零相位滤波器）和幅 值 / 相位校正改善复调制 FFT 细化算法的性能［4］，它 适用于离线信号处理. 若数据较长或感兴趣频段很窄， 即感兴趣频率线数 K — MN，可以通过子带分解 DFT （SD - DFT）实现频率细化. 4 子带分解 SD - FFT 及其改进 SD - DFT 最大特点是可以单独计算出不同带宽 的子带谱，而无需设