没有合适的资源?快使用搜索试试~ 我知道了~
基于MATLAB小波分析及应用,包括小波工具箱的使用,连续小波的介绍,离散小波的介绍等
资源推荐
资源详情
资源评论
第八章 小波分析理论及应用
第八章 小波分析及应用
8.1 引言
把函数分解成一系列简单基函数的表示,无论是在理论上,还是实际应用中都有重
要意义。
1822 年法国数学家傅里叶(J. Fourier 1768-1830)发表的研究热传导理论的“热的力学
分析”,提出并证明了将周期函数展开为正弦级数的原理,奠定了傅里叶级数理论的基
础
[1]
。傅里叶级数理论研究的是把函数在三角函数系下的展开,使得对信号和系统的研
究归结为对简单的三角函数的研究。傅里叶级数与傅里叶变换共同组成了平常所说的
傅里叶分析
[2]
。傅里叶级数用于分析周期性的函数或分布,理论分析时经常假定周期是
,定义如式(8.1-1)、(8.1-2)
, (8.1-1)
其中 (8.1-2)
然而,被分析函数的性质并不能完整地由傅里叶系数来刻划,这里有一个例子来说
明
[3]
:从任一个平方可和的函数 出发,为了得到一个连续函数 ,只需或者增
大 f(x)的傅里叶系数的模,或者保持它不变并适当地改变系数的位相。因此,不可能仅
根据傅里叶系数大小的阶就预知函数的性质(如大小、正则性)。
傅里叶变换的定义如式(8.1-3)、(8.1-4)
(8.1-3)
(8.1-4)
通过引入广义函数或分布的概念,可获得奇异函数(如冲击函数)的傅里叶变换的
存在。对于时域的常量函数,在频域将表现为冲击函数,表明具有很好的频域局部化
性质。由式(8.1-3)可知,为了得到 ,必须有关于 f(x)的过去和未来的所有知识,
而且 f(x)在时域局部值的变化会扩散到整个频域,也就是 的任意有限区域的信息
都不足以确定任意小区域的 f(x)。在时域,哈尔(Haar)基是一组具有最好的时域分辨能
力的正交基,它在时域上是完全局部化的,但在频域的局部化却很不好,这是由于哈
尔系的两个缺点:缺乏正则性与缺乏振动性。研究者们希望寻找关于空间变量(或时
间变量)与频域变量都同时好的希尔伯特(Hilbert)基,R. Balian 认为:“在通讯理论中,
人们对于在完全给定的时间内,把一个振动信号表示成由其中每一个都拥有足够确定
的位置与有一个频率的小波的叠加这件事感兴趣。事实上,有用的信息常常同时被发
射信号的频率与信号的时间结构(如音乐)所传递。当把一个信号表达成时间的函数
时,其中的频谱表现并不好;相反地,信号的傅里分析却显示不了信号每一分量发射
16
第八章 小波分析理论及应用
信号的瞬时与持续时间。一个适当的表示应结合这两者互补描述的优点,并用一个离
散的刻划来表示,以适应通讯理论
[3]
。”
为此,人们提出了短时傅里叶变换(STFT)的概念:
定义 8.1-1 若 选择得使 W 与它的傅里叶变换 满足:
那么使用 W 作为窗函数,在式(8.1-5)中引入的窗口傅里叶变换称为“短时傅里叶变
换”(STFT):
(8.1-5)
当窗函数选择为高斯(Gaussian)函数时,则为 Gabor 变换
[2]
。
STFT 的缺点是分析窗的大小和形状是恒定的。因为频率与周期成反比,所以反映
信号的高频成份需要窄的时间窗,而反映信号的低频成份需要宽的时间窗, STFT 无法
满足要求,此外,STFT 的冗余很大,增加了不必要的计算量。
小波变换作为能随频率的变化自动调整分析窗大小的分析工具,自八十处代中期以
来得到了迅猛的发展,并在信号处理、计算机视觉、图像处理、语音分析与合成等众
多的领域得到应用。
小波分析方法的出现可以追溯到 1910 年 Haar 提出 Haar 规范正交基,以及 1938 年
Littlewood-Paley 对傅里叶级数建立的 L-P 理论。为克服传统傅里叶分析的不足,在八
十年代初,便有科学家使用“小波”的概念来进行数据处理,比较著名的是 1984 年法国
地球物理学家 Morlet 引入小波的概念对石油勘探中的地震信号进行存贮和表示。在数
学方面所做的探索主要是 R. Coifman 和 G. Weiss 创立的“原子”和“分子”学说,这些“原
子”和“分子”构成了不同函数空间的基的组成部分。L. Carleron 使用了非常象“小波”的函
数构造了 Stein 和 Weiss 的空间 的无条件基。直到 1986 年,法国数学家 Meyer 成功
地 构 造 出 了 具 有 一 定 衰 减 性 的 光 滑 函 数 , 它 的 二 进 伸 缩 与 平 移
构成 的规范正交基。此前,人们普遍认为这是
不可能的,如 Daubechies,Grossman 和 Meyer 都退而研究函数系 构
成 的框架的条件去了。
Lemarie 和 Battle 继 Meyer 之后也分别独立地给出了具有指数衰减的小波函数。
1987 年,Mallat 利用多分辨分析的概念,统一了这之前的各种具体小波的构造,并提
出了现今广泛应用的 Mallat 快速小波分解和重构算法。1988 年 Daubechies 构造了具有
紧支集的正交小波基。Coifman, Meyer 等人在 1989 年引入了小波包的概念。基于样条
函数的单正交小波基由崔锦泰和王建忠在 1990 年构造出来。 1992 年 A. Cohen, I.
Daubechhies 等人构造出了紧支撑双正交小波基。同一时期,有关小波变换与滤波器组
之间的关系也得到了深入研究。小波分析的理论基础基本建立起来。
17
第八章 小波分析理论及应用
近年来,一种简明有效的构造小波基的方法--提升方案(Lifting Scheme)得到很大的
发展和重视
[4,5]
。利用提升方案可把现存的所有紧支撑小波分解成更为基本的步骤
[6]
,另
外,它还为构造非线性小波提供了一种有力的手段,所以,利用提升方案构造的小波
被认为是第二代小波
[5]
。小波理论及其应用仍然处在发展中,其未来将在非线性多尺度
方法、非规则集上的小波构造以及非平稳、非均匀、时变信号处理等方面等到更深入
的研究。
8.2 小波变换及其基本性质
8.2.1 连续小波变换
, 的连续小波变换(有时也称为积分小波变换)定义为:
(8.2-1)
或用内积形式:
(8.2-2)
式中
要使逆变换存在, 要满足允许性条件:
(8.2-3)
式中 是 的傅里叶变换。
这时,逆变换为
(8.2-4)
这个常数限制了能作为“基小波(或母小波)”的属于 的函数 的类,尤其是
若还要求 是一个窗函数,那么 还必须属于 ,即
故 是 R 中的一个连续函数。由式(8.2-3)可得 在原点必定为零,即
18
第八章 小波分析理论及应用
(8.2-5)
从式(8.2-5)可以发现小波函数必然具有振荡性。
连续小波变换具有如下性质:
性质 1(线性):设 ,则
性质 2(平移不变性):若 ,则 。平移不变
性是一个很能好的性质,在实际应用中,尽管离散小波变换要用得广泛一些,但在需
要有平移不变性的情况下,离散小波变换是不能直接使用的。
性 质 3 ( 伸 缩 共 变 性 ) : 若 , 则 , 其 中
c>0。
性质 4(冗余性):连续小波变换中存在信息表述的冗余度。其表现是由连续小波
变换恢复原信号的重构公式不是唯一的,小波变换的核函数 存在许多可能的选
择。尽管冗余的存在可以提高信号重建时计算的稳定性,但增加了分析和解释小波变
换的结果的困难。
8.2.2 连续小波变换的离散化
由于连续小波变换存在冗余,因而有必要搞清楚,为了重构信号,需针对变换域的
变 量 a , b 进 行 何 种 离 散 化 , 以 消 除 变 换 中 的 冗 余 , 在 实 际 中 , 常 取
,这时
常简写为: 。
变换形式为:
为了能重构信号 ,要求 是 的 Riesz 基。
19
第八章 小波分析理论及应用
定义 8.2-1 一个函数 称为一个 R 函数,如果 在下述意义上是一个
Risez 基: 的线性张成在 中是稠密的,并且存在正常数 A 与 B,
,使
对所有二重双无限平方可和序列 成立,即对于 的
成立。
假定 是一个 R 函数,那么存在 的一个唯一的 Riesz 基 ,它在意
义
上与 对偶。这时,每个 有如式(8.2-6)的唯一级数表示:
(8.2-6)
特别地,若 构成 的规范正交基时,有
重构公式为:
(8.2-7)
8.3 多分辨分析与 Mallat 算法
8.3.1 多分辨分析
Mallat 使用多分辨分析的概念统一了各种具体小波基的构造方法,并由此提出了现
今广泛使用的 Mallat 快速小波分解和重构算法,它在小波分析中的地位与快速傅里叶
变换在傅里叶分析中的地位相当
[7]
。
定义 8.3-1 空间 的多分辨分析是指构造该空间内一个子空间列 ,使其具
20
剩余35页未读,继续阅读
资源评论
zoufeng214
- 粉丝: 0
- 资源: 1
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功