Matlab小波工具箱的使用.docx

Matlab 小波工具箱的使用 1 最近想尝试一下小波的用法，就这matlab 的帮助尝试了一下它的例子，顺便翻译了一下帮助的内容，发现matlab 帮助做的确实不错，浅显易懂！现把翻译的文档写出来吧，想学习的共同学习吧！ 小波工具箱简介 小波工具箱包含了图像化的工具和命令行函数，它可以实现如下功能： l 测试、探索小波和小波包的特性 l 测试信号的统计特性和信号的组分 l 对一维信号执行连续小波变换 l 对一维、二维信号执行离散小波分析和综合 l 对一维、二维信号执行小波包分解（参见帮助Using Wavelet Packets） l 对信号或图像进行压缩、去噪 另外，工具箱使用户更方便的展示数据。用户可以做如下选择： l 显示哪个信号 l 放大感兴趣的区域 l 配色设计来显示小波系数细节 工具箱可以方便的导入、导出信息到磁盘或matlab 工作空间。具体详见File Menu Options 一维连续小波分析 这一部分来测试连续小波分析的特性。连续小波分析只需要一个小波函数cwt。在这一部分将学到如下内容： l 加载信号 l 对信号执行连续小波变换 l 绘制小波系数 l 绘制指定尺度的小波系数 l 绘制整个尺度小波系数中的最大值 l 选择显示方式 l 在尺度和伪频率之间切换 l 细节放大 l 在普通或绝对模式下显示系数 l 选择执行小波分析的尺度 使用命令行执行连续小波分析 这个例子是一个包含噪声的正弦波 加载信号 load noissin 可以使用whos 显示信号信息whos Name noissin 执行连续小波变换 Size 1x1000 Bytes 8000 Class double c = cwt(noissin,1:48,db4); 函数 cwt 的参数分别为分析的信号、分析的尺度和使用的小波。返回值 c 包含了在各尺度下的小波系数。对于这里，c 是一个 48x1000 的矩阵，每一行与一个尺度相关。 绘制小波系数 cwt 函数可以接受第四个参数，来指定函数在执行结束后是否绘制连续小波变换系数的绝对值。另外还可以接受更多的参数来定义显示的不同特性，详见 cwt 函数。如下面的语句绘制系数结果 c = cwt(noissin,1:48,db4,plot); 选择分析的尺度 cwt 函数的第二个参数可以设定任意小波分析的尺度，只要这些尺度满足如下要求 l 所有尺幅必须为正实数 l 尺度的增量必须为正 l 最高的尺度不能超过由信号决定的一个最大值如下面的代码可以执行从 2 开始的偶数尺度计算 c = cwt(noissin,2:2:128,db4,plot); 显示结果如下 这幅图像很明确的表示出了信号的周期性。 使用图形接口做连续小波分析 开启一维连续小波工具，只需输入如下命令 wavemenu 出现如下小波工具箱主菜单 选择Continuous Wavelet 1-D 菜单项，出现如下一维信号分析连续小波分析工具 加载信号 选择菜单 File-Load Signal，在 Load Signal 对话框里选择 noissin.mat 文件，它在 matlab 安装目录的 toolbox/wavelet/wavedemo 文件夹下，点击 OK 加载信号。 一维连续小波工具开始加载信号，加载后默认采样频率为1s。 执行连续小波变换 下面来测试使用 db4 小波对尺度 1 到 48 做小波分析，设置如下 点击 Analyze 按钮 在短暂的计算后，工具将绘制小波系数，并在Coefficients line 坐标系中绘制尺度为 24 的小波系数，在 local maxima 坐标系中绘制各尺度的小波系数最大值。 查看小波 Coefficients Line 在小波系数图中右键点击可以选择展示其他尺度的小波系数，选择后点击New Coefficients Line 按钮，Coefficients Line 会相应更新。 查看 Maxima Line 点击 Refresh Maxima Line 按钮，可以显示从尺度 1 到所选尺度的小波系数的最大值。 注意当在系数图中按下鼠标右键并移动时，会在最下面的Info 框中显示当前鼠标位于的 X 位置和尺度。 在尺度和伪频率之间切换 在右边选择 Frequencies，当再在系数图中选择时，在 Info 中显示的将是 Hz。 而关于尺度和频率的转换关系，可以看How to Connect Scale to Frequency? 选择要显示的坐标系 放大细节 在系数框中按鼠标左键可以选择放大的范围。 选择好放大范围后点击最下面的按钮可以实现指定的放大 显示普通系数或系数绝对值 两种显示方式的区别在于，普通模式下，颜色映射是在系数的最大最小之间；而绝对模式， 颜色映射是在