傅里叶变换分析信号的缺点

精品文档 1欢迎下载 傅里叶变换分析信号的缺点傅里叶变换分析信号的缺点 基于傅里叶 Fourier 变换的信号频域表示 揭示了时间函数 和频谱函数之间的内在联系 在传统的平稳信号分析和处理中发挥了 极其重要的作用 很多理论研究和应用研究都把傅里叶变换当作最基 本的经典工具来使用 但是傅里叶变换存在着严重的缺点 用傅里叶 变换的方法提取信号频谱时 需要利用信号的全部时域信息 这是一 种整体变换 缺少时域定位功能 因此必须对其加以改进 傅里叶变换的特点及其局限性 设函数 f t 在 内有定义 且使广义积分 1 t 1 2 2 都收敛 则称 1 式定义的广义积分为函数 f t 的傅里叶变换 记为 F f t 2 式定义的广义积分为逆傅里叶变换 记为 F F 1 傅里叶变换可以完成从时域到频域的转换 正变换 也可以完成从频 域到时域的转换 逆变换 但不能同时具有时域和频域信息 其核函 数是 由于三角函数具有填满整个空间的特性 其在物理空间中是 双向无限延伸的正弦波 在积分变换中体现为积分范围从 到 因此 傅里叶变换是先天的非局限性 它对信号 f t 中体现任何局部 信息处理都是相同的 而事实上 工程技术中的许多信号 如 语音信 精品文档 2欢迎下载 号 地震信号 心电图和各种电脉冲 他们的信号值只出现在一个短 暂的时间间隔 t 内 以后快速减为零 t 以外是未知的 可能为零 也可能是背景噪音 如果用 1 式从信号中提取谱信号 F 就要取 无限的时间量 使用过去的及将来的信号只为计算单个频谱 不能反 映出随时间变化的频率 实际上我们需要的是确定的某个时间间隔内 的频谱 这就使人们想到改进傅里叶变换使其能用来处理某个确定 时间范围内的信号 Gabor 提出的窗口傅里叶变换就是一个有效的 方法 另外 傅里叶变换之所得到广泛应用与透镜能实现傅里叶变换是分不 开的 由公式 1 0 2 2 0 0 0 2 0 0 0 0 其中物平面为 焦平面为 d0 为物距 d1 为象平面 x0 0 x 要使 F 即准确实现傅里叶光学变换 只有在 0 x 0 0 f 时才能实现 否则将出现位相弯曲 并且 只有正透镜才能 1 0 实现傅里叶变换 这些限制给工程技术中无疑增加了困难 这使得人 们不得不寻求新得的方法 分数傅立叶变换不要求严频谱面 可根据 需要在既包含空域信息也包括空频域信息的平面上进行处理 这使光 学信息处理更具灵活性 精品文档 3欢迎下载 1 傅里叶变换缺乏时间和频率的定位功能 傅里叶变换及其逆变换表示如下 t 1 2 由以上两式可知 傅里叶变换是一种整体变换 对信号的表征要么完 全在时域内 要么完全在频域内 和 t 是互相排斥的两个变量 用 傅里叶变换的方法得到某一个频率 0 的频谱分量 S 0 必须从 的整个时间轴上进行积分 如果要从频谱得到信号在某一时刻 t0 的值 s t0 则需要对 S X 在整个频率轴上进行积分 因此 傅里 叶变换得到的是信号 s t 在整个时间范围内的频率特性 它不能告 诉人们在某段时间里信号发生了什么变化 也无法获得某一频率出现 的时刻信息 因此 它不具有时间和频率的定位功能 2 傅里叶变换对于非平稳信号的局限性 信号的瞬时频率 表示了信号的谱峰在时间 频率平面上的位置 及其随时间的变化情况 一般平稳信号的瞬时频率为常数 而非平稳 信号的瞬时频率是时间 t 的函数 从傅里叶变换变换的表达式可以看 出 S X 是单变量 X 的函数 信号的傅里叶变换不随时间的变化而变 化 因此 傅里叶变换仅仅适用于平稳信号 但是 在实际工作中 我们 精品文档 4欢迎下载 分析和处理的往往是时变的或非平稳的信号 它们的频率随时间变化 而变化 其相关函数 功率谱等也是时变信号 用傅里叶变换进行分 析 得到的信号频谱反映的是整体信号中包含的某一频率分量的平均 值 所以傅里叶变换不能反映信号瞬时频率随时间的变化情况 仅仅 适用于分析平稳信号 对频率随时间变化的非平稳信号 傅里叶变换 只能给出其总体效果 不能完整地把握信号在某一时刻的本质特征 3 傅里叶变换在时间和频率分辨率上的局限性 分辨率是信号处理的基本概念之一 包括频率分辨率和时间 分辨率 在时域分析中 信号处理的目标是尽可能地同时获得高的时 间分辨率和频率分辨率 然而 可以证明时域窗和频域窗乘积恒定且 大于等于 1 2 也即不可能同时获得高的时频分辨率 这就是著名的 不确定性原理 傅里叶变换在这方面的表现尤其不尽如人意 傅里叶 变换可以改写成内积的形式 即 由于傅里叶变换等效于 s t 和基函数做内积 而对不同的 构成一族正交基 因此 S 精确地反映了 s t 在该频率点的分量 大小 基函数在频域是位于 处的 函数 因此 当用傅里叶变换来 分析信号的频域特性时 具有最好的频率分辨率 但是 在时域对 应的是正弦函数 其在时域的持续时间是 因此 其时域分辨率 最差 对于傅里叶逆变换 分辨率的情况正好相反 这一结果既体现了 精品文档 5欢迎下载 信号的时频不确定性原理 也反映了傅里叶变换在时域和频域分辨率 方面所固有的矛盾 显然 傅里叶变换本身不可能根据信号的特性来 自动调节时域和频域的分辨率 时频分析时频分析 时频分析 JTFA 即时频联合域分析 Joint Time Frequency Analysis 的简称 作为分析时变非平稳信号的有力工具 成为现 代信号处理研究的一个热点 它作为一种新兴的信号处理方法 近 年来受到越来越多的重视 时频分析方法提供了时间域与频率域的 联合分布信息 清楚地描述了信号频率随时间变化的关系 时频分析的基本思想是 设计时间和频率的联合函数 用它同时 描述信号在不同时间和频率的能量密度或强度 时间和频率的这种 联合函数简称为时频分布 利用时频分布来分析信号 能给出各个 时刻的瞬时频率及其幅值 并且能够进行时频滤波和时变信号研究 信号时频分析具有重要的意义 我们很有必要对信号的时频进行研 究分析 精品文档 6欢迎下载 常用的时频分析方法常用的时频分析方法 时间和频率是描述信号的两个最重要的物理量 信号的时域 和频域之间具有紧密的联系 根据时间和频率之间的关系 信号的 时频分析的主要方法有 窗口傅立叶变换 Gabor 变换 小波变换 希尔伯特黄变 Hilbert Huang Transform HHT 窗口傅里叶变换 窗口傅里叶变换亦称短时傅里叶变换 它是由 Gabor 首先系统地使 用的 其基本想法为 傅里叶变换是频域分析的基本工具 为了达 到时间域上局部化 在傅里叶分析中的基本变换函数之前乘上一个 时间上有限的时限函数 即窗口函数 tg 然后再用它们来作傅里 叶分析 这样 tje 起频限作用 tg 起到时限作用 合起来 就 可起到时频双限制作用 其中 tg 是有紧支集 即窗口外数据为零 的函数 tx 为被分析的信号 随着 的位置变动 tg 所 确定的 时间窗 在 t 轴上移动 使 tx 逐步进入被分析的状态 窗口函数 tg 一般为实的偶函数 窗口外数据为零 紧支集 或 精品文档 7欢迎下载 很快趋于零 这时傅里叶变换结果不再为 X 而是 GX 这里 xG 大致反映了 tx 在时刻 时频率为 的 信号成分 的相对含量 时频局部化就是希望找一种信号的表示方法 它能同 时提供时域和频域的局部化信息 而这种变换确实能反映函数在窗 口内部 附近 的频谱特征 窗口傅里叶变换可使信号达到局部 平稳 更好地研究局部范围的特性 窗口函数 tg 的傅里叶变换 它在有限区间之外数据恒等于零 用 tg 乘 tx 即在 附近开 窗口 为窗口傅里叶变换 Gabor 只做了高斯窗 的傅里叶变换 它是窗口傅里变换的一种 尽 管窗口傅里叶变换是一种时频分析 是信号处理的重要工具 并得 到广泛的应用 但是窗口傅里叶变换的一个主要缺点是时域和频域 的采样间隔都是常数 即这种窗口大小和形态与频率无关 是固定 不变的 不能使变换窗口大小随频率而变化 但在处理实际问题 我们希望时域的采样间隔随着频率的增高而减小 同时窗口傅里叶 变换不管如何离散化均不能使它成为一组正交基 为此 J Morlet 等人对窗口傅里叶变换进行了改造 引入了小波变换 连续小波变换 小波变换时今年来在图像处理中受到十分重视的新技术 面向图像 压缩 特征测以及纹理分析等许多方法在时频分析中有重要的应用 线性系统理论中的傅立叶变换是以在两个方向上都无限伸展的正弦 曲线波作为正交基函数的 对于瞬态信号或高度局部化的信号 例 精品文档 8欢迎下载 如边缘 由于这些成分并不类似于任何一个傅立叶基函数 它们的 变换系数 频谱 不是紧凑的 频谱上呈现出一幅相当混乱的构成 这种情况下 傅立叶变换是通过复杂的安排 以抵消一些正弦波的 方式构造出在大部分区间都为零的函数而实现的 为了克服上述缺 陷 使用有限宽度基函数的变换方法逐步发展起来了 这些基函数 不仅在频率上而且在位置上是变化的 它们是有限宽度的波并被称 为小波 wavelet 基于它们的变换就是小波变换 所有小波是通 过对基本小波进行尺度伸缩和位移得到的 基本小波是一具有特殊 性质的实值函数 它是震荡衰减的 而且通常衰减得很快 在数学 上满足积分为零的条件 0 而且其频谱满足条件 2 即基本小波在频域也具有好的衰减性质 有些基本小波实际上在某 个区间外是零 这是一类衰减最快的小波 一组小波基函数是通过 尺度因子和位移因子由本小波来产生 连续小波变换定义为 1 精品文档 9欢迎下载 小波分析的应用是与小波分析的理论研究紧密地结合在一起地 现在 它已经在科技资讯产业领域取得了令人瞩目的成就 电子资 讯技术是六大高新技术中重要的一个领域 它的重要方面是影像和 信号处理 现今 信号处理已经成为当代科学技术工作的重要部分 信号处理的目的就是 准确的分析 诊断 编码压缩和量化 快速 传递或存储 精确地重构 或恢复 从数学的角度来看 信号与影 像处理可以统一看作是信号处理 影像可以看作是二维信号 在小 波分析地许多分析的许多应用中 都可以归结为信号处理问题 现 在 对于其性质随实践是稳定不变的信号 处理的理想工具仍然是 傅立叶分析 但是在实际应用中的绝大多数信号是非稳定的 而特 别适用于非稳定信号的工具就是小波分析 希尔伯特黄变换 希尔伯特特换变换的方法主要由 2 个部分组成 经验模态分解 empirical mode decomposition 简称 EMD 和 Hilbert 谱分析 经验模态分解方法是一种自适应的 高效的数据分解方法 由于这 种分解是以局部时间尺度为基础 因此 它适应于非线性 非平稳 过程 通过经验模型分解 任何复杂的数据集都可以被分解为个数 有限的 而且常常是为数不多的几个固有模函数 intrinsic mode functions 简称 IMF 的线性叠加 一个固有模态函数是满足以下 两个条件的函数 1 1 在整个数据区间内 极值点的数目与过 精品文档 10欢迎下载 零点的数目相等或至多相差 1 个 2 在任意一点处 由局部极大 值点定义的包络以及由局部极小值点定义的包络的均值为零 EMD 方法通过不断的剔出极大值和极小值连接上下包络的均值 将原信号分解为 1 1 其中 为一个 IMF 分量 为残余分量 一般为信号 的平均趋势 为常数序列或单调序列 从基函数理论的角度来看 EMD 对不同信号分解出的基函数 是不同的 它不同傅里叶分解 的基 一系列恒定幅度与频率的正余弦函数 也不同于小波分解的 基函数 预先给定基函数的形式 因此 EMD 分解不仅改进了信号 分解的效率 而且使这种分解方法更有利于非平稳数据处理 通过 分解得到 IMF 后 就可以对每一个分量做希尔伯特变换 得到其瞬 时频率和幅度 设 IMF 分量为 则它的复解析信号为 其中 a t 为幅值函数 表达式为 t 为相位函数 表达式为 t tan 1 精品文档 11欢迎下载 其中幅值函数表示信号每个采样点的瞬时幅度能量 相位函数 表示信号每个采样点的瞬时相位 对其求导就得到瞬时频率 对每 个 IMF 分量做 Hilbert 变换并忽略分解余项 数据可以表示为 t 1 根据式 1 可以将幅度和瞬时相位作为时间的函数表示在三维 平面中 幅度的这种时一频分布被称为希尔伯特幅度谱 简称为希 尔伯特谱 习惯上用幅度的平方来表示能量密度 这里如果用幅度 平方代替希尔伯特幅度谱中的幅度 将得到希尔伯特能量谱 对于 希尔伯特能量谱 如果 EMD 分解得到的 IMF 分量彼此完全正交 那 么信号的平方 式中的第二项为 0 这对于时频能量表示是十分有利的 虽然 对于某些特殊的数据 相邻的分量在不同的时间段内可能含有相同 的频率成分 但从局部意义上说 任何两个分量都是正交的 泄漏 的大小通常与数据长度以及分解结果有着直接的关系 对于有限的 精品文档 12欢迎