小波系数实部等值线图和小波系数方差图的绘制

文摘   其他   2024-07-15 09:30   重庆  

小波系数实部等值线图

小波系数方差图的绘制

作者:第八星系-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

下图为最终成果图

THE

END

如需获取数据,

请后台回复“小波分析

扫码进群加入我们~

本文编辑:周游


第八星系人造大气理论爱好者
记录与交流python、matlab等科研工具。记录与交流大气科学的学科知识
 最新文章