小波变换理论及应用..docVIP

2011-2012 学年第学期以图示的方式详细说明连续小波变换（CWT）的运算过程，分析小波变换的内涵；并阐述如何从多分辨率（MRA）的角度构造正交小波基。（20分） 综述小波变换理论与工程应用方面，不少于3000字。（25分） 运用MATLAB中的小波函数和小波工具箱，分别对taobao.wav语音信号在加噪之后的taobao_noise.wav信号进行降噪处理，要求列出程序、降噪结果及降噪的理论依据。（25分）一）连续小波变换（CWT）的运算过程 （ 1.1） 其中，aR且a≠0。式（1.19）定义了连续小波变换，a为尺度因子，表示与频率相关的伸缩，b为时间平移因子。其中为窗口函数也是小波母函数。 从式（1.1）可以得出，连续小波变换计算分以下5个步骤进行。 ① 选定一个小波，并与处在分析时段部分的信号相比较。 ② 计算该时刻的连续小波变换系数C。如图1.5所示，C表示了该小波与处在分析时段内的信号波形相似程度。C愈大，表示两者的波形相似程度愈高。小波变换系数依赖于所选择的小波。因此，为了检测某些特定波形的信号，应该选择波形相近的小波进行分析。 图1.5 计算小波变换系数示意图 ③ 如图1.6所示，调整参数b，调整信号的分析时间段，向右平移小波，重复①～②步骤，直到分析时段已经覆盖了信号的整个支撑区间。 ④ 调整参数a，尺度伸缩，重复①～③步骤。 ⑤ 重复①～④步骤，计算完所有的尺度的连续小波变换系数，如图1.7所示。 图1.6 不同分析时段下的信号小波变换系数计算 图1.7 不同尺度下的信号小波变换系数计算 小波变换的实质是用小波（微小的特定波形）与待分析信号波形分段求内积，所得的系数反映了小波与待分析信号的相似度，相似度越高则系数越高。通过改变平移因子b可以实现对信号时频域的分析。通过改变尺度因子可以改变小波与待分析信号的相似度。最后由得到的系数和所选小波的特性可以知道待分析信号的特性或是待分析信号某一时段或频段的特征。 （二）从多分辨率（MRA）的角度构造正交小波基分解为细节部分和大尺度逼近部分,然后将大尺度逼近部分进一步分解。如此重复就可以得到任意分辨率上的逼近部分和细节部分。在MRA理论中同一尺度下小波函数和尺度函数分别满足。 同一尺度下小波函数同尺度函数正交 小波函数和尺度函数在多分辨率分析中满足方程 这两个方程就是二尺度方程。利用二尺度方程可以构造出小波母函数，通过伸缩平移就得到整个平方可积空间的基。正交尺度函数构造正交小波基，还有当尺度函数为Riesz基是构造的正交小波基函数。所以说MRA不仅为正交小波基的构造提供了一种简单的方法，而且为正交小波变换的快速算法提供了理论依据。 （三）小波变换理论与工程应用方面 对于任意的X中的x,y,z p(x,z)=p(x,y)+p(y,z)都成立 L2(μ(X))表示X 上所有在几乎处处（almost everywhere）意义下平方可积（square-integrable）的复值的可测函数的集合。平方可积表示该函数的绝对值的平方的积分是有限的。　　巴拿赫空间 Banach Space　　巴拿赫空间是一个完备的赋范矢量空间Normed Vector Space，它是希尔伯特空间的推广。巴拿赫空间定义为完备的线性赋范矢量空间。即是说，它是一个实数或复数的矢量空间并且有一个完备的范数||·|| ，即其每个柯西Cauchy序列都是收敛列。作为窗口函数，用同待分析函数相乘，然后在傅里叶变换： （2．1） （2．2） 这一加窗变换使得我们可以分析出一个信号在任意局部范围的频率特征，这是比傅里叶变换优越之处。这一类加窗变换Fourier变换统称为短时傅里叶变换（Short Time Fourier Transform，简称为STFT）。但是其时频窗口不随频率和时间的变化而变化，使它的灵活性与普遍性运用受到限制。 2.2小波变换基本理论 为了使得短时傅里叶变换的时，频窗口均随频率的变化而变化，以实现对低频分量采用大时窗，对高频分量采用小时窗的符合自然规律的分析方法。我们设计一组连续变化的伸缩平移基，称为连续小波基函数，来代替STFT中的。 小波函数的确切定义为：设为一平方可积函数，也即，若傅里叶变换满足条件： （2．3） 则称为一个基本小波或小波母函数，并称式（2.3）为小波函数的可容许性条件。 连续小波变换：将任意平方可积空间中的f（t）在小波基下进行展开，称这种展开为函数f（t）的连续小波变换（Continue Wavelet Transform，简记为CWT）其表达式为 （2．4） 由表达式可知小波变换也是类