一维连续小波变换的matlab基础程序实现_第1页
一维连续小波变换的matlab基础程序实现_第2页
一维连续小波变换的matlab基础程序实现_第3页
一维连续小波变换的matlab基础程序实现_第4页
一维连续小波变换的matlab基础程序实现_第5页
全文预览已结束

下载本文档

版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领

文档简介

小波变换实验二连续小波变换1、 实验目的本实验的目的在于充分理解连续小波变换的算法和作用，利用matlab程序实现对一维信号进行连续小波变换，进而在程序的编辑过程理解一位连续小波变换的小波系数矩阵的含义。同时通过对预算的到的小波系数矩阵进行分析解释，得到原始信号的频谱分布以及了解小波系数在尺度和位移两个分量上的意义。2、 实验原理、实验编程思路1、根据书本的理论知识，知道一维连续小波变换的公式为：实际在编程过程当中，对于上式中积分的求解可以采用将积分函数离散化，通过求和来实现求积分，离散的过程如下式：本实验中，根据题目可以知道采样的时间间隔为0.03s，即上式中t，在实际编程当中为了计算方便可以省略掉这个时间常数，所以在编程过程当中使用的公式实际为：2、 小波函数的选取：使用墨西哥草帽（mexhat）小波来进行小波变换，墨西哥草帽的函数为（支撑区间为-55）：对于连续小波函数的采样间隔，根据不同的尺度参量来进行采样，比如尺度为i，实际对应小波的采样间隔取k/i，以保持和原信号在不同尺度上的同步。3、 程序运算简化：在程序设计过程当中，如果对于小波系数的每一个系数都按照公式来计算，算法的时间复杂度应当为o(n3)。但通过对公式的分析，不难看出，对于同意尺度a，相邻的两个小波系数之间的求和项，只有第一项或者最后一项或者二者都不同，所以在下一个系数求解的时候可以减少一次循环，从而将时间复杂度降到o(n2)，运算效率大大提高。4、 在程序设计的过程当中，还分别对原信号进行傅里叶分析和直接的cwt变换，将得到的结果与设计的连续小波变换程序进行比对分析。3、 实验程序和结果墨西哥草帽小波参数获取函数：mexh.m连续小波变换主函数：mexh-cwt.m傅里叶分析和cwt分析：fft cwt result.m1、 利用mexh-cwt.m对源数据进行分析得到的结果：原信号波形图小波系数矩阵及原信号的三个主成分小波系数矩阵的三维视图功率谱图和功率等值线图从上述的小波变换结果，特别是小波系数矩阵的信息中可以看出，原始信号主要有三个平率不同且时域分布也不同的主成份组成，从图中可以定性地看出，频率最高的成分1始终出现在整个信号段，而频率次之的成分2只在信号刚开始的阶段出现，频率最低的成分3基本上在成分2消失之后开始出现。这些信息都是从单纯的傅里叶频谱分析中无法得到的。2、 利用fft cwt result.m对原信号进行频谱分析和算法校验：从上述傅里叶频谱分析图中同样可以看出信号的三个主成份，且对应频率近似分别为：0.651Hz、1.237Hz、2.604Hz。在小波变换中，每个尺度对应0.03s，上述三个频率对应周期分别为：1.536s、0.808s、0.384。对应到尺度域上容易计算的对应的尺度间隔分别为
人人文库> 全部分类> 行业资料 > 管理策划

温馨提示

  • 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

最新文档

评论

0/150

提交评论