小波系数实部等值线图和小波系数方差图的绘制
作者：第八星系-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.下图为最终成果图
微信搜索“第八星系人造大气理论爱好者”公众号，关注获取文章数据
页:
[1]