小波分析matlab实例.doc
到小波分析1 背景 传统的信号理论,是建立在Fourier分析基础上的,而Fourier变换作为一种全局性 的变化,其有一定的局限性。在实际应用中人们开始对Fourier变换进行各种改进,小波 分析由此产生了。小波分析是一种新兴的数学分支,它是泛函数、Fourier分析、调和分 析、数值分析的最完美的结晶;在应用领域,特别是在信号处理、图像处理、语音处理以 及众多非线性科学领域,它被认为是继Fourier分析之后的又一有效的时频分析方法。 小波变换是近年发展起来的一种基于时频域的信号分析工具,它具有良好的时频局部 性、选基灵活性和去相关性等优点,可用于光谱信号的噪声滤波和基线校正等。此后,多 位物理、数学家的合作共同奠定了小波变换的理论和应用基础。由于小波变换能够更精确 地分析信号的局部特征,在很多领域得到了越来越多地应用。小波分析的应用领域十分广 泛,它包括:数学领域的许多学科;信号分析、图象处理;量子力学、理论物理;军事电 子对抗与武器的智能化;计算机分类与识别;音乐与语言的人工合成;医学成像与诊断; 地震勘探数据处理;大型机械的故障诊断等方面;例如,在数学方面,它已用于数值分析、 构造快速数值方法、曲线曲面构造、微分方程求解、控制论等。在信号分析方面的滤波、 去噪声、压缩、传递等。在图象处理方面的图象压缩、分类、识别与诊断,去污等。以及 在医学方面的应用,如核磁共振成像时间、提高 CT 、B 超等分辨率。 2 小波变换的产生及去噪的必要性 我们在一维信号分析中,可知傅里叶变换将信号分解成一系列不同频率的正弦或余弦 波的叠加,与之类似,小波变换也可将信号分解成一系列小波函数的叠加,这一系列小波 函数都由某个母小波函数经过平移和尺度变换得来。以不规则的小波信号来逼近局部信号 显然比用光滑的正弦信号逼近程度要好,而用不同尺度小波对同一信号进行逼近又有利于 对信号进行逐步细致的分析,这正是小波分析的基本思想。小波变换采用变化的时频窗, 窗口面积固定,但形状可变。分析低频信号时,采用拉伸的小波和长的时间窗以获取足够 信息,分析高频信号时,采用压缩小波和短时间窗以获取足够精度。常见的小波函数有 Meyer 波、Morlet 波、8阶高斯波等。 传统的去噪方法常使用 Fourier 变换去噪,将含噪信号变换到频域,然后采用低通滤波 器进行滤波,但是基于 Fourier 变换的去噪方法存在着保护信号局部性和抑制噪声之间的矛 盾。Fourier 变换去噪不能有效的将噪声与有用信号的高频部分和有噪声引起的高频干扰加以有效的区分开来。这就使得我们在研究信号去噪课题上注意到小波的好处,小波去噪可 以很好的保护有用信号的尖峰和突变部分的信号。小波变换具有良好的时频局部化性质, 具有以下优点:(1)小波分解可以覆盖整个频域(提供了一个数学上完备的描述);(2)小波变 换通过选取合适的滤波器,可以极大的减小或去除所提取得不同特征之间的相关性;(3)小 波变换具有“变焦”特性,在低频段可用高频率分辨率和低时间分辨率(宽分析窗口),在高 频段,可用低频率分辨率和高时间分辨率(窄分析窗口);(4)小波变换实现上有快速算法 (Mallat 小波分解算法)。因此采用小波去噪是具有必要性的。3 小波变换理论 3.1 小波定义 满足以下条件(1) d 0 2 ) ( 或其等价条件(2) 0 ) ( dt t 的函数 称为基本小波,或母小波。其中 为 的傅里叶变换。式(2)说明母 ) (t ) ( ) (t 小波函数具有一定的振荡性,即包含某种频率特性。(3) ) ( ) ( 2 1 , a b t a t b a 式中 均为常数,且 。显然, 是基本函数 先作移位再作伸缩以后得到 b a, 0 a ) ( , t b a ) (t 的。若 不断地变化,我们可得到一族函数 。a 为伸缩因子,反映函数的尺度, b a, ) ( , t b a a1 波形被拉伸,越大拉伸越多。b为平移因子,表示 沿 t 轴的平移位置。 是母小波经移位和伸缩所产生的一族函数,我们称之为小波基 ) ( , t b a 函数,或简称小波基。 3.2 小波的特性 连续小波 的时频窗口中心和宽度可以精确定位,且都随尺度 a 的变化而伸缩。 ) ( , t b a 若将时、频域窗口综合考虑,根据公式推导可得时频窗口的面积与尺度 a 无关,即时间分辨 率和频率分辨率是相互制约的。小波尺度 a 与频率 相对应。当 变小时,对 的时域 a ) (t x 观察范围变窄,但对 在频率观察的范围变宽,且观察的中心频率向高频处移动;反 ) ( X 之,当 变大时,对 的时域观察范围变宽,频域的观察范围变窄,且分析的中心频率 a ) (t x 向低频处移动。小波变换恒 Q 性质。带宽/中心频率= ,不论 为何 Q a a 0 0 / / / a 值 , 始终保持了和 具有性同的品质因数。恒 Q 性质是小波变换的一个 ) 0 ( a ) ( a t ) (t 重要性质,也是区别于其它类型的变换且被广泛应用的一个重要原因。3.3 连续小波变换和反变换定义 函数 以小波 为基的连续小波变换定义为函数 和 的内积, ) (x f ) (x ) (x f ) ( , x b a dx a b x a x f f b a W b a f ) ( 1 ) ( , ) , ( , 在1984年,A.Grossman和J.Morlet指出,连续小波的逆变换为, dadb a x f C x f b a b a 2 , , ) ( , 2 ) ( 其中, 为母小波y(x)的允许条件(admissible condition), C d C ) ( ˆ 其中, 为 的傅立叶变换,而 是在平方可积的实数空间 。 ) ( ˆ ) (x ) (x ) ( 2 R L 3.4 离散小波变换 在计算连续小波变换时,实际上也是用离散的数据进行计算的,只是所用的缩放因子 和平移参数比较小而已。不难想象,连续小波变换的计算量是惊人的。为了解决计算量的 问题,缩放因子和平移参数都选择 (j>0 的整数)的倍数。使用这样的缩放因子和平移参 j 2 数的小波变换叫做双尺度小波变换,它是离散小波变换的一种形式。执行离散小波变换的 有效方法是使用滤波器。该方法是 Mallat 在 1988年开发的,叫做 Mallat 算法,这种方法 实际上是一种信号的分解方法,在数字信号处理中称为双通道子带编码。 3.5 傅里叶分析与小波包分析的比较 从以上分析中可以看出通过傅立叶分析进行滤波得到的结果与小波分析得到的结果有 些差异,主要是由于信号集中在低频部分,噪声分布在高频部分,所以通过低通滤波器进行滤 波,不能将有用信号的高频部分和由噪声引起的高频干扰加以有效地区分。若低通滤波器太 窄,