小波系数实部等值线图和
小波系数方差图的绘制
作者:第八星系-Bella
邮箱:13072806080@163.com
本次以贵阳市近60年年平均气温数据为例,为大家介绍小波系数实部等值线图和小波系数方差图的绘制。
一、将数据转化为MATLAB格式
1
数据准备,如表中所示,排成一列
2
打开MATLAB,将分析数据导入其中。修改数据类型为“列向量”,并导入所选内容,然后关闭导入页面
3
将已经转化为MATLAB形式的数据另存备用
二、边界效应的消除
1
使用waveletAnalyzer语句调用小波分析工具箱,回车,然后单击Signal Extension,进入延伸界面
2
将刚刚保存的MATLAB格式的数据导入
修改“Desired Length”为“Next Power of 2”的2倍
“Extension Mode”选择“Symmetric (Whole-Point)”,
点击Extend,此时,数据分别向前后以对称形式延伸34[(128-60)/2]个数据
保存延长后的数据
关闭数据延长界面
三、计算复小波系数、实部、模和模的平方
1
单击小波分析工具箱中的“Complex Continuous Wavelet 1-D”
进入小波分析界面计算复小波系数,导入延长后的数据
“Wavelet”修改为“comr”
然后单击“Analyze”
保存分析结果(结果为复小波系数)
关闭小波分析工具箱
2
将复小波系数导入MATLAB,coefs中的数据就是复小波系数,需要注意的是:此时的复小波系数包含延长后的数据,需要删除前后34[(128-60)/2]列数据
3
采用real函数计算复小波系数的实部
采用abs函数计算复小波系数的模
计算小波系数模的平方、方差
计算出的数据全部保存在工作区
4
将数据全部复制粘贴到excel表格里以备绘图使用
四、绘制实部等值线图和方差图(其他图同理绘制)
1
实部等值线图
为小波系数实部添加横纵坐标,本例横坐标为年份,纵坐标为周期尺度
复制全部数据,粘贴到origin之中
调整细节即可完成绘图
2
方差图为折线图,可用Excel、Origin、MATLAB等绘制,较为简单,不再赘述
3
下图为最终成果图
如需获取数据,
请后台回复“小波分析”
或扫码进群加入我们~
本文编辑:周游