matlab编程代做highspeedlogic★小波分析理论

1、小波分析理论简介 刘玉民（一） 傅立叶变换伟大的历史贡献及其局限性1 Fourier 变换 1807年，由当年随拿破仑远征埃及的法国数学、物理学家傅立叶（JeanBaptistle Joseph Fourier ,1786-1830),提出任意一个周期为T (=)的函数，都可以用三角级数表示： = = + + (1) = = (2) (3)对于离散的时程 ，即 N 个离散的测点值 ，0,1，2，,N-1, 为测量时间：=+= （4） 其中 ，0，1，2，, (5) , 1,2, -1 (6) ，0，1，2，,N-1 (7) ， (8) 当T 时，化为傅立叶积分（即 Fourier 变换）： =

2、 （9） (10)傅立叶变换的理论是人类数学发展史上的一个里程碑，从1807 年开始，直到1966年(1807年傅立叶提出任意一个周期函数都可以表示为傅立叶级数的结论是有误的，直到1966年才证明了可积的周期函数才能表示为傅立叶级数)，整整用了一个半世纪多，才发展成熟。她在各个领域产生了深刻的影响，得到了广泛的应用，推动了人类文明的发展。其原因是，傅立叶理论不仅仅在数学上有很大的理论价值，更重要的是傅立叶变换或傅立叶积分得到的频谱信息具有物理意义。所以说，傅立叶理论是万古流芳的。数学上的插值方法。除傅立叶级数外，还有拉格朗日插值，有限元插值，勒让德多项式插值即高斯积分使用的插值方法。遗憾的是，

3、这种理论具有一定的局限性：（1） 傅立叶变换的三种形式中的傅立叶系数都是常数，不随时间 t 变化，因而只能处理频谱成分不变的平稳信号，相反的，在处理非平稳信号时会带来很大误差，甚至与实际情况大相径庭。（举例：无阻尼与有阻尼的单自由度的自由振动、打秋千、座钟、讨论会与大合唱等）。在实际信号中，若高频与低频差别很大，在相同的时间间隔内，高频信号衰减了而低频信号尚未衰减，所以，在不同时刻，信号的频谱成分是不同的。硬要用傅立叶变换找出所有时刻的频谱成分，硬要把幅值的变化用频率的变化来补偿，不仅高频的傅立叶系数有误差，低频的傅立叶系数也有很大误差，包括求出的频率当然也有误差。 （2） 求傅立叶系数是全时

4、间域上的加权平均，这从上面的（5）、（6）、（7）公式可以清楚看到。局部突变信息被平均掉了，局部突变信息的作用很难反映出来（好比吃大锅饭，平均主义）。差别很大的信号，如方波、三角波、正弦波，都可以得到相同的频率，所以，处理、捕捉突变信号如故障信号，灵敏度很差。处理、捕捉突变信号应使用能反映局部信息的变换。 为了克服以上两点局限性，这就要求：（1） 将变换系数视为随时间变化的，级数求和由一重变为两重。（2） 使用能反映局部信息的变换，则函数组不能使用全域上的函数，只能使用有所谓紧支撑的函数，即“小波函数”或 加窗傅立叶变换的窗函数。2 Garbor 变换窗口 Fourier 变换 在时间频率分析

5、中， Fourier 变换公式的不足已经被 D. Garbor 注意到了，在 1946 年的论文中，为了提取信号的Fourier 变换的局部信息，引入了一个时间局部化的 Gaussian 函数作为“窗函数” g(t-b),其中参数 b 用于平移动窗以便覆盖整个时间域。因为一个 Gaussian 函数的Fourier 变换还是Gaussian 函数，所以Fourier 逆变换即频率也是局部的。窗口 Fourier 变换简介。对于时间局部化的“最优”窗，用任一Gaussian 函数 （11）“Garbor 变换”的定义为 （12）由于 1 （13）所以 = （14）令 = （15）利用 Parse

6、val 恒等式， = = = （16）这个等式说明，除去乘数项 之外，在 具有窗函数的的“窗口 Fourier 变换”，与在 具有窗函数的的“窗口 Fourier 逆变换”一致，根据窗函数 的宽度是 2 的 结论，这两个窗的宽度分别是2 和 这两个窗的笛卡儿积是 加窗傅立叶变换的“时间频率窗”的宽度对于观察所有的频率是不变的。在较长的时间窗内，对于高频信号，可能经过了很多周期，因而求出的Fourier 变换系数是很多周期的平均值，局部化性能不能得到体现。若减小时间窗（减小），高频信号局部化性能得到体现，但对于很低的频率信号来讲，检测不到。总上所述，加窗傅立叶变换对于高频与低频差别很大的信号仍不

7、是很有效的。3 窗口 Fourier 变换的测不准原理对于一个非平凡函数 ，若满足 （A）条件，则可作为短时窗口 Fourier 变换的窗函数，若其Fourier 变换也满足上述条件，那么 （B）而且，等号成立，如且仅如 （C）其中 。小结：（1） 傅立叶级数的正弦与余弦系数为常数，不能反映振幅变化的情况；（2） 求傅立叶系数需要所考虑的时间域上所有信息，不能反映局部信息的特征；（3） 加窗傅立叶变换时间窗是固定不变的，高频与低频的时间局部化不能同时满足。由于上述原因，必须进一步改进，克服上述不足，这就导致了小波分析。（二） 小波分析 将时程函数表示为下面的小波级数： = （17） （18）其

8、中， 是小波函数， 是小波系数，且 = ( 19） 由公式（17）到（19） 可以看到，小波级数是两重求和，小波系数的指标不仅有频率的指标，而且还有时间的指标 。也就是说，小波系数不仅像傅立叶系数那样，是随频率不同而变化的，而且对于同一个频率指标 ，在不同时刻 ，小波系数也是不同的。这样就克服了上面所述的第一个不足。由于小波函数具有紧支撑的性质，即某一区间外为零。这样在求各频率水平不同时刻的小波系数时，只用到该时刻附近的局部信息，从而克服了上面所述的第二个不足。与有限元比较。在这一点，小波插值要比有限元高明。有限元虽然是局部的“单元插值”，但单元之间的公共节点上，只能保证阶连续，而导数不连续。

9、小波插值可保证二阶导数连续，只要选三次样条小波就能做到。第三个不足，小波分析是如何克服的呢？通过与加窗傅立叶变换的“时间频率窗”的相似分析，可得到小波变换的“时间频率窗”的笛卡儿积是 （20）其中，时间窗的宽度为 ，随着频率的增大（即 的增大）而变窄，随着频率的减小（即 的减小）而变宽，之所以有这样的结果，关键在于公式（18）中，时间变量 前面乘了个“膨胀系数” 。小波变换的“时间频率窗”的宽度，检测高频信号时变窄，检测低频信号时变宽，这正是时间频率分析所希望的。根据小波变换的“时间频率窗”的宽度可变的特点，为了克服上面所述的第三个不足，只要不同时检测高频与低频信息，问题就迎刃而解了。如，选择

10、从高频到低频的检测次序，首先选择最窄的时间窗，检测到最高频率信息，并将其分离。然后，适当放宽时间窗，再检测剩余信息中的次高频信息。再分离，再放宽时间窗，再检测次次高频信息，依次类推。为了检测到不同频率水平信息，即求出不同频率水平下不同时刻的小波系数，首先要选好小波函数。选择小波函数的“四项原则”。在求小波系数公式（19）中，如果 是 空间的正交基，则的为的复共轭。小波分析的最重要的应用是滤波，为了保证滤波不失真，小波函数必须具有线性相位，至少具有广义线性相位。小波分析的另一重要应用是捕捉、分析突变信号，这就要使用函数的导数，小波函数至少是连续。由前面分析可知，小波函数必须具有紧支撑的性质。所以

11、，正交、线性相位、连续、紧支撑是选择小波函数的“四项原则”。如果选择某个小波函数，同时满足四项指标，那真是人类的福气。遗憾的是，上帝像是有意考验我们的数学家，没有将“四合一”的小波函数“直接”恩赐给人类。数学家们已经证明，具有正交、线性相位、紧支撑的小波函数只有 Harr函数，而Harr函数是间断函数，对于工程应用来说，是不理想的。目前，一种倾向是坚持正交性。另一种倾向是放弃正交性，另辟途径，进行艰辛的长征，前仆后继，花费了将近半个世纪的探索，才使小波分析理论成熟起来，得以在工程中应用。作为后人，我们要忠心地感谢他们。为了进行小波分解与重构，“四合一”的小波函数不存在，数学家们“一分为四”，选

12、择了四个函数，巧妙地解决了这些问题。这四个函数是：尺度函数 ，小波函数 ，对偶尺度函数 ，对偶小波函数 。为什么要选择四个函数呢？由前面小波变换的“时间频率窗”分析可知，小波变换的“时间频率窗”的宽度，当检测高频信号时变窄，检测低频信号时变宽。为了检测到所有频率信号，“时间频率窗”的宽度必须按一定的次序变化，不失一般性，从窄到宽，检测频率信号从高频到低频的次序进行实际上也正是这样的次序。在最高频率水平 （即根据实测数据的时间测量间隔 ，最高能检测到的频率为 Nyquist 频率 ），选择最窄的“时频窗”宽度，检测到原始信号中的最高频率信号，并将这些信号从原始信号中剥离，存放在空间，而将剥离后的

13、剩余低频信号的总合，存放在另一空间。然后，增大“时频窗”的宽度，再检测空间中的高频信息，将这些信号从空间中剥离，存放在空间，而将剥离后的剩余低频信号的总合，存放在另一空间。依次类推。这就要求有两个互相有联系的空间： = = , 子空间性质简介：（1） （2） clos () = （3） , 即 （21）（4）（5） （6） 这样，对于参考子空间，需要单个函数在意义 = （22） 上生成，其中， (23）对于参考子空间，需要单个函数在意义= clos （24）上生成，其中， （25）首先，这就要求有两个函数： 和 ，前者称为尺度函数，后者称为小波函数。并且，它们肯定是有关系的。由（21）中的（1

14、）式可知，由（4）式可知，而：是 的一个基，所以存在唯一的序列、 （26） （27）并引入记号 = （28） = （29）（26）称为尺度函数的“两尺度关系”，（27）称为小波函数的“两尺度关系”，称为尺度函数的“两尺度符号”，称为小波函数的“两尺度符号”。为了由高频到低频逐次检测到不同频率水平的信息，仅有上述两个函数是不够的。由前面分析可知，公式（17）中的小波函数与求小波系数使用的与作内积的函数不是同一个函数，除非使用正交的小波函数。这就要求寻找尺度函数与小波函数的对偶函数：对偶尺度函数 、对偶小波函数 ，以便分解原始信号时能求出小波系数来。对偶关系简介与（26）（29）相对应，有 （30

15、） （31） = （32） = （33）为了满足对偶关系，必须满足以下对偶条件： = (34) = (35)与（34）、（35）等价的条件是 ， （36） ， （37） = det ， (38) = det ， (39)其中 ，G(z)为的共轭，为的共轭， = (40) = (41)由（36）和 （38）便可以得到 的基 到 和 中的基分解关系： （42）由（42）式便可得到频率水平子空间和中向量坐标的分解算法 （43） （44）其中 = （45） （46）根据两尺度关系，便可得到频率水平子空间 中向量坐标的重构算法： （47）有了上述的系数，就可以使用多分辨分析的金字塔算法，快速求出小波系数

16、。以三次样条函数尺度函数 为例，说明其步骤。1 将投影到上 = （48）2 小波分解算法使用多分辨分析的金字塔算法 / / / （49） （50） = （51） （52） （53）而、是符号多项式、的系数：= （54）= (55）= （56）3 小波重构算法 / / / 在我们研制的小波分析软件 “XIAOBO ”中，使用的是三次样条函数尺度函数及三次样条小波（B-小波）。三次样条小波 = （57）其中 = ，（）（58）阶样条函数的定义是 ， （59）是区间0，1的特征函数。 （60）三次样条小波具有连续、紧支撑、广义线性相位和半正交的性质。三次样条小波是阶连续，用它模拟一个信号，二阶导数都

17、是连续的,精度极高。三次样条小波支撑区间为0，7，三次样条尺度函数支撑区间为0，4。三次样条小波是对称函数，因而具有广义线性相位，滤波不会失真。唯一不足的是三次样条小波不是正交的，只是半正交的，在小波分解中，求小波系数必须使用其对偶小波 ，对偶小波 的傅立叶变换为 （61） （62）小波分类简介1 正交小波正交小波的两个族与满足：（1） <> =（2） <> = 0（3） <> = （4） 用振动的语言形象的比喻：是 子空间的低阶振型，是 子空间的高阶振型， ，。2 半正交小波 半正交小波的四个族、与满足：(1) <> = (2) <>

18、; = (3) <> = 0, <> = 0，；j,l,k,m 。(4) 是 子空间的低阶振型的线性组合，即低阶剩余模态。(5) 是 子空间的高阶振型的线性组合，即高阶剩余模态。(6) ， ，(7) 3 小波 小波的与满足：(1) <> = (2) 正交性都是相对于对偶来说的，本身都不是正交族，只是线性无关族。神奇的小波分析为了检验小波分析的能力，在实测数据中叠加上一个低频、幅值为原始最大幅值的十分之一的正弦拍波，经过小波分解，在和两个分辨率水平下都检测到了这个波，尤其是 检测到的正弦拍波曲线，如荷花“出污泥而不染”！（三） 正交小波包与双正交小波包分析从上

19、面介绍的多分辨分析的金字塔算法进行小波分解可以看到，每次分解，都是对进行的 ：对 进行小波分解，将 分解为 和 ；再 对 进行小波分解，将 分解为 和 ，依次类推。而各子空间 不再分解，也就是说，每次都是对低频进行再分解，而高频不再分解！究其原因，便是 的基向量到 和 基、有分解关系（42），而的基向量 没有相应的分解关系。事实上，与、是线性无关的 ！显然，高频的时间频率局部化不是最优的。为了克复这个缺点，必须使用小波包的分解方法。 无论正交小波包还是双正交小波包，与小波的最基本的区别在于它们具有 （63）的关系。根据这个特殊的正交的两尺度关系，可以定义关于尺度函数的“小波包”： （64）并且

20、，小波包族 具有下属正交性： < > = (65) < > = 0 (66)因此有分解关系： ， （67）小波分解中的关系 （68）现在可以改写作 （69）并且，（69）式可以由 推广到任意一 ： （70） = （71）这样，对于正交小波的小波分解得到的 （72）对于每个1，2，3，,可以用小波包再分解： （73）双正交小波包与正交小波包类似，略去。（四） 向量分解小波包简介 上面介绍了正交小波包与双正交小波包分解，与小波分解的最基本的区别在于它们具有（63）式 因而导出了正交小波包和双正交小波包的分解。正如前面提到的，小波分解与小波包分解的共同点，都是高频空间的基向量

21、能向低频空间的基向量分解，若不能分解，便进行不下去。如小波分解中，的基向量 没有相应的分解关系，所以不能再分解。小波包克服了这种弊病，得到了各子空间都能适用的基向量分解关系式（71），才解决了这个问题。我们能不能跳出基向量能分解才能分解的框框呢？即，基向量不能分解时也能分解？答案是肯定的。这就是向量分解小波包的理论。现在，让我们从小波重构关系式（47），k=1,2,3,M出发，导出向量分解小波包的理论。首先，将（47）改写为矩阵的形式 ： = + （74）其中，矩阵、 由两尺度序列、中的数组成。我们把矩阵、的每一列，看成向量空间的一个向量，把也看成向量空间的一个向量，那么，就是这些矩阵、的列向

22、量的一个线性组合： = + （75）问题是，这些矩阵、的列向量是否线性无关？由分解等式（43）和（44）可知，当为零向量时，k=1,2,3,，这就得出了矩阵、的列向量是线性无关的结论。因此，可以把矩阵中的 M 个列向量看成 M 维向量空间的一组基。找到了M 维向量空间的一组基，就可以对任一M 维向量进行分解，找出在这组基下的坐标。只要认准被分解的对象是M 维向量空间的一个向量，而不是再看做 或 中向量的坐标。这样，在小波分解中，不仅得到各频率水平的坐标 可按原来的步骤和方法进行分解，而且得到的各频率水平的坐标 也可按向量分解的理论进行再分解，并且分解的步骤和方法可以和对坐标 进行再分解的完全一

23、样。可能有人会产生疑问：前面讲到， = ，现在 中的向量怎么又变成 中的向量呢 ？其实，这里使用了“移花接木”的技巧。我们现在是对 中的向量分解，这是千真万确的。但，是通过对其坐标 的分解来实现的。也就是说，要求的是等于什么 ？ 而不是 中的向量 等于什么 ？ 中的向量 等于= （76）既然要求的是等于什么，而不是 等于什么，不妨“移花接木”一下，把拿来，虚构 中一向量 ， = （77）不就可以对进行分解了吗？反正在分解过程中，使用的是、 ，求的是新的和，而根本就不出现，管它是什么样子呢！只要保证分解后低频与高频分离即可。既然分解，低频与高频能分离，那么，分解，同样低频与高频可以分离的。其实从

24、（74）式中两个矩阵、的组成，都是两尺度序列、中的数，仔细分析一下，不难看出，中的列相当于低阶振型，中的列相当于高阶振型。以后的工作和小波包的分解类似，只不过这里是直和分解，不是正交分解罢了。前面曾提到，小波分解的第一步，将投影到上（见公式（48）。由于投影到上求小波系数要解方程，很费时间，所以，好多书上都说，可以近似地把原始数据作为小波系数，误差不大。用向量分解小波包的理论来看，这并非是近似的，而是精确的！事实上，使用上面“移花接木”的技巧，虚构 中一向量 = （78） 是原始数据组成的向量。那末，不就是中的向量的坐标吗？只要对它进行分解与重构就得了，而且，重构后，也用不着通过（78）式求，

25、重构得到的，就是原始数据组成的向量，两头省！ 当然，我并非反对投影，而且在捕捉突变奇异信号时，最灵敏的方法是利用导数确定，投影到上先求出作为的中的向量的坐标，能得到快速的求导方法。但这是另外的问题了。有了向量分解小波包新理论，可以在分解过程中，把各种小波组合到一起，例如，求 时，可以用三次样条小波，这样确保滤波重构不失真、求导精度高的长处，在分解的时候，可使用正交分解的分解序列系数，以提高分解效率。这样，上帝虽然没有将“四合一”的小波函数“直接”恩赐给人类，但还是间接地恩赐给人类了。我们研制的小波包信号处理软件，具有三种小波包分解、分解后各频率水平的分析、64 种滤波重构功能及重构后的傅立叶变换和奇异信号大小分析功能。三种小波包分别是“向量分解小波包”、“双正交小波包”和“Daubechies正交小波包”。通过实例对比，向量分解小波包分析效果远比后者好，而且分解速度相当快，4096个点，分解所花 CPU 为0.027 秒。具体实例见后。图 一 原始实测数据（4096=个点）图 二