一维连续小波变换的matlab基础程序实现

1、小波变换实验二连续小波变换一、实验目的本实验的目的在于充分理解连续小波变换的算法和作用，利用 matlabmatlab 程序实现对一维信号进行连续小波变换，进而在程序的编辑过程理解一位连续小波变换 的小波系数矩阵的含义。同时通过对预算的到的小波系数矩阵进行分析解释，得到原始信号的频谱分布以及了解小波系数在尺度和位移两个分量上的意义。二、实验原理、实验编程思路1 1 根据书本的理论知识，知道一维连续小波变换的公式为：AXWf(a,b) f,a,bf(t)()dtviaa实际在编程过程当中，对于上式中积分的求解可以采用将积分函数离散化，通过求和来实现求积分，离散的过程如下式：(k 1) t1t b

2、Wf(a, b)f(t)a2() tkk ta1上ktba2tf(k t)(-)ka本实验中，根据题目可以知道采样的时间间隔为0.03s0.03s，即上式中厶 t t，在实际编程当中为了计算方便可以省略掉这个时间常数，所以在编程过程当中使用的公式实际为：1k bWf(a,b) a2f(k)()ka2 2、小波函数的选取：使用墨西哥草帽(mexhatmexhat)小波来进行小波变换，墨 西哥草帽的函数为(支撑区间为-5-55 5)：2Xmexh(1 x )e 2 , x 50,else对于连续小波函数的采样间隔，根据不同的尺度参量来进行采样，比如尺度 为 i i，实际对应小波的采样间隔取 k/i

3、k/i，以保持和原信号在不同尺度上的同步。3 3、 程序运算简化：在程序设计过程当中，如果对于小波系数的每一个系数都按照公式来计算，算法的时间复杂度应当为 o(o( n3)n3)。但通过对公式的分析，不难看出，对于同意尺 度 a,a,相邻的两个小波系数之间的求和项，只有第一项或者最后一项或者二者都 不同，所以在下一个系数求解的时候可以减少一次循环，从而将时间复杂度降到o(no(n2)2)，运算效率大大提高。4 4、 在程序设计的过程当中，还分别对原信号进行傅里叶分析和直接的cwtcwt变换，将得到的结果与设计的连续小波变换程序进行比对分析。三、实验程序和结果墨西哥草帽小波参数获取函数：mexh

4、.mmexh.m连续小波变换主函数：mexh-cwt.mmexh-cwt.m傅里叶分析和 cwtcwt 分析：fftfft cwtcwt result.mresult.m1 1 利用 mexh-cwt.mmexh-cwt.m 对源数据进行分析得到的结果：原信号波形图0100150290250300350400小波系数矩阵及原信号的三个主成分小波系数矩阵的三维视图功率谱图和功率等值线图从上述的小波变换结果，特别是小波系数矩阵的信息中可以看出，原始信号 主要有三个平率不同且时域分布也不同的主成份组成，从图中可以定性地看出， 频率最高的成分 1 1 始终出现在整个信号段，而频率次之的成分 2 2 只

5、在信号刚开始 的阶段出现，频率最低的成分 3 3 基本上在成分 2 2 消失之后开始出现。这些信息都 是从单纯的傅里叶频谱分析中无法得到的。2 2、利用 fftfft cwtcwt result.mresult.m 对原信号进行频谱分析和算法校验:从上述傅里叶频谱分析图中同样可以看出信号的三个主成份，且对应频率近似分别为：0.651Hz0.651Hz、1.237Hz1.237Hz、2.604Hz2.604Hz。在小波变换中，每个尺度对应 0.03s,0.03s,上述三个频率对应周期分别为：1.536s1.536s、0.808s0.808s、0.3840.384。对应到尺度域上容易计算的对应的尺度间隔分别为：51.251.2、26.926.9、12.12.8 8。再次观察小波系数，分析得到的三