小波分析MATLAB实例

到小波分析

1 背景

传统的信号理论，是建立在Fourier分析基础上的，而Fourier变换作为一种全局性的变化，其有一定的局限性。在实际应用中人们开始对Fourier变换进行各种改进，小波分析由此产生了。小波分析是一种新兴的数学分支，它是泛函数、Fourier分析、调和分析、数值分析的最完美的结晶；在应用领域，特别是在信号处理、图像处理、语音处理以及众多非线性科学领域，它被认为是继Fourier分析之后的又一有效的时频分析方法。

小波变换是近年发展起来的一种基于时频域的信号分析工具，它具有良好的时频局部性、选基灵活性和去相关性等优点，可用于光谱信号的噪声滤波和基线校正等。此后，多位物理、数学家的合作共同奠定了小波变换的理论和应用基础。由于小波变换能够更精确地分析信号的局部特征，在很多领域得到了越来越多地应用。小波分析的应用领域十分广泛，它包括：数学领域的许多学科；信号分析、图象处理；量子力学、理论物理；军事电子对抗与武器的智能化；计算机分类与识别；音乐与语言的人工合成；医学成像与诊断；地震勘探数据处理；大型机械的故障诊断等方面；例如，在数学方面，它已用于数值分析、构造快速数值方法、曲线曲面构造、微分方程求解、控制论等。在信号分析方面的滤波、去噪声、压缩、传递等。在图象处理方面的图象压缩、分类、识别与诊断，去污等。以及在医学方面的应用，如核磁共振成像时间、提高CT 、B超等分辨率。

2 小波变换的产生及去噪的必要性

我们在一维信号分析中，可知傅里叶变换将信号分解成一系列不同频率的正弦或余弦波的叠加，与之类似，小波变换也可将信号分解成一系列小波函数的叠加，这一系列小波函数都由某个母小波函数经过平移和尺度变换得来。以不规则的小波信号来逼近局部信号显然比用光滑的正弦信号逼近程度要好，而用不同尺度小波对同一信号进行逼近又有利于对信号进行逐步细致的分析，这正是小波分析的基本思想。小波变换采用变化的时频窗，窗口面积固定，但形状可变。分析低频信号时，采用拉伸的小波和长的时间窗以获取足够信息，分析高频信号时，采用压缩小波和短时间窗以获取足够精度。常见的小波函数有Meyer波、Morlet 波、8阶高斯波等。

传统的去噪方法常使用Fourier变换去噪，将含噪信号变换到频域，然后采用低通滤波器进行滤波，但是基于Fourier变换的去噪方法存在着保护信号局部性和抑制噪声之间的矛盾。Fourier变换去噪不能有效的将噪声与有用信号的高频部分和有噪声引起的高频干扰加以有效的区分开来。这就使得我们在研究信号去噪课题上注意到小波的好处，小波去噪可以很好的保护有用信号的尖峰和突变部分的信号。小波变换具有良好的时频局部化性质，具有以下优点：(1)小波分解可以覆盖整个频域(提供了一个数学上完备的描述)；(2)小波变换通过选取合适的滤波器，可以极大的减小或去除所提取得不同特征之间的相关性；(3)小波变换具有“变焦”特性，在低频段可用高频率分辨率和低时间分辨率(宽分析窗口)，在高频段，可用低频率分辨率和高时间分辨率(窄分析窗口)；(4)小波变换实现上有快速算法(Mallat小波分解算法)。因此采用小波去噪是具有必要性的。

3 小波变换理论

小波定义

满足以下条件

∞<ψ⎰

∞ωωωd 02)( （1） 或其等价条件

⎰∞

∞-=0)(dt t ψ （2）

的函数)(t ψ称为基本小波，或母小波。其中)(ωψ为)(t ψ的傅里叶变换。式（2）说明母小波函数具有一定的振荡性，即包含某种频率特性。

)()(21

,a

b t a t b a -=-ψψ （3） 式中b a ,均为常数，且0>a 。显然，)(,t b a ψ是基本函数)(t ψ先作移位再作伸缩以后得到的。若b a ,不断地变化，我们可得到一族函数)(,t b a ψ。a 为伸缩因子，反映函数的尺度，a<1波形被压缩，越小压缩越厉害，a>1波形被拉伸，越大拉伸越多。b 为平移因子，表示沿t 轴的平移位置。)(,t b a ψ是母小波经移位和伸缩所产生的一族函数，我们称之为小波基函数，或简称小波基。

小波的特性

连续小波)(,t b a ψ的时频窗口中心和宽度可以精确定位，且都随尺度a 的变化而伸缩。若将时、频域窗口综合考虑，根据公式推导可得时频窗口的面积与尺度a 无关,即时间分辨率和频率分辨率是相互制约的。小波尺度a 与频率ω相对应。当a 变小时，对)(t x 的时域观察范围变窄，但对)(ΩX 在频率观察的范围变宽，且观察的中心频率向高频处移动；反之，当a 变大时，对)(t x 的时域观察范围变宽，频域的观察范围变窄，且分析的中心频率向低频处移动。小波变换恒Q 性质。带宽/中心频率=Q a

a 00=Ω∆=Ω∆ΩΩ///，不论a 为何值)0(>a ，)(a

t ψ始终保持了和)(t ψ具有性同的品质因数。恒Q 性质是小波变换的一个重要

性质，也是区别于其它类型的变换且被广泛应用的一个重要原因。

连续小波变换和反变换定义

函数)(x f 以小波)(x ψ为基的连续小波变换定义为函数)(x f 和)(,x b a ψ的内积，

⎰∞

∞--ψ>=ψ=

b x a x f f b a W b a f )(1)(,),(, 在1984年，和指出，连续小波的逆变换为，

⎰⎰∞∞-∞∞--ψψ>ψ<=dadb a x f C x f b a b a 2,,)(,2)(

其中，ψC 为母小波y (x )的允许条件(admissible condition)，

⎰∞∞-ψ∞<ψ

=ωωωd C )(ˆ

其中，)(ˆωψ

为)(x ψ的傅立叶变换，而)(x ψ是在平方可积的实数空间)(2R L 。 离散小波变换

在计算连续小波变换时，实际上也是用离散的数据进行计算的，只是所用的缩放因子和

平移参数比较小而已。不难想象，连续小波变换的计算量是惊人的。为了解决计算量的问题，缩放因子和平移参数都选择j

2(j>0的整数)的倍数。使用这样的缩放因子和平移参数的小波变换叫做双尺度小波变换，它是离散小波变换的一种形式。执行离散小波变换的有效方法是使用滤波器。该方法是Mallat 在1988年开发的，叫做Mallat 算法，这种方法实际上是一种信号的分解方法，在数字信号处理中称为双通道子带编码。

傅里叶分析与小波包分析的比较

从以上分析中可以看出通过傅立叶分析进行滤波得到的结果与小波分析得到的结果有

些差异,主要是由于信号集中在低频部分,噪声分布在高频部分,所以通过低通滤波器进行滤波,不能将有用信号的高频部分和由噪声引起的高频干扰加以有效地区分。若低通滤波器太窄,则在滤波后,信号中仍存在大量的噪声,若低通滤波器太宽,则将一部分有用信号当作噪声被滤掉。因此小波分析对非平稳信号的消噪有着傅立叶分析不可比拟的优点。 4 小波去噪