在Matlab 6.5界面下的Workspace中将cerunoff.mat文件导入，然后双击打开，全部复制到Excel后去掉延伸数据的小波变换系数（本例中去掉前12列和后13列），或只复制原时间序列的小波变换系数到Excel，最后使用Excel中的IMREAL函数计算原时间序列的小波系数实部（图6）。
图6 复小波系数及实部计算示意图
Excel 中IMREAL函数的调用方法为：单击“插入”菜单下的“函数（F）”按钮，弹出图7所示的“插入函数”窗口，在“搜索函数（S）：”框中输入：“IMREAL”后单击“转到”，再单击“确定”，出现函数参数窗口（图8）。在“Inumber”一栏的空白处输入所要计算的数据（图6），单击“确定”即可得到小波系数实部值。
图7 IMEAL 函数调用
图8 IMEAL 函数参数
需要说明的是，从cerunoff.mat文件中转到Excel里的复小波系数，在其实部和虚部中间包含许多“空格”，在计算之前需要先将其去掉。
5. 借助Suffer 8.0，绘制小波系数实部等值线图 5.1 小波系数实部等值线图的绘制
首先，将小波系数实部数据按照图9格式排列，其中列A为时间，列B为尺度，列C为不同时间和 尺度下所对应的小波系数实部值。
图9 小波系数实部数据格式
其次，将图9数据转化成Suffer 8.0识别的数据格式。具体操作为：在Suffer 8.0界面下，单击“网格”菜单下的“数据”按钮，在“打开”窗口选择要打开的文件（小波系数实部.xls），单击“打开”后弹出“网格化数据”对话框（图10）。它给出了多种不同的网格化方法、文件输出路径及网格线索几何学等信息。这里我们选择“克里格“网格方法”，单击“确定”，完成数据格式的转化。
图10 小波系数实部数据格式转化
图11 Suffer8.0中的小系数实部等值线图
最后，绘制小波系数实部等值线图。在Suffer 8.0界面下，单击“地图”菜单下的“等值线图-新建等值线图”按钮，弹出“打开网格”窗口后，选择“小波系数实部.grd”文件，单击“打开”，完成等值线
图的绘制并存盘（图11）。
5.2 小波系数实部等值线图在多时间尺度分析中的作用
小波系数实部等值线图能反映径流序列不同时间尺度的周期变化及其在时间域中的分布，进而能判断在不同时间尺度上，径流的未来变化趋势。为能比较清楚的说明小波系数实部等值线图在径流多时间尺度分析中的作用，我们利用Suffer 8.0对其进一步处理和修
30饰，得到图12显示的小波系数实部等值线图。其中，横坐LHHL标为时间（年份），纵坐标为时间尺度，图中的等值曲线为25L小波系数实部值。当小波系数实部值为正时，代表径流丰水期，在图中我们用实线绘出，“H”表示正值中心；为负时，20表示径流枯水期，用虚线绘出，“L”表示负值中心。
15LHHL由图12可以清楚的看到径流演化过程中存在的多时LHL间尺度特征。总的来说，在流域径流演变过程中存在着10HLHH18~32年，8~17年以及3~7年的3类尺度的周期变化规律。其中，在18~32年尺度上出现了枯-丰交替的准两次震荡；5(a)在8~17年时间尺度上存在准5次震荡。同时，还可以看出年份 1970197519801985199019952000以上两个尺度的周期变化在整个分析时段表现的非常稳
图12 小系数实部等值线图
定，具有全域性；而3~10年尺度的周期变化，在1980s以后表现的较为稳定。
6. 绘制小波系数模和模方等值线图 6.1 小波系数模和模方等值线图的绘制
参考4、5两步，绘制小波系数模和模方等值线图（图13、14）。 说明：在Excel中，复数模的计算函数为“IMABS”。
3030252520201515101055(c)(b)1970197519801985199019952000年份 1970197519801985199019952000年份 图13 小波系数模等值线图 图14 小波系数模方等值线图
6.2 小波系数模等值线图在多时间尺度分析中的作用
Morlet小波系数的模值是不同时间尺度变化周期所对应的能量密度在时间域中分布的反映，系数模值愈大，表明其所对应时段或尺度的周期性就愈强。从图13可以看出，在流域径流演化过程中，18~32年时间尺度模值最大，说明该时间尺度周期变化最明显，18~22年时间尺度的周期变化次之，其他时间尺度的周期性变化较小；
6.2 小波系数模方等值线图在多时间尺度分析中的作用
小波系数的模方相当于小波能量谱，它可以分析出不同周期的震荡能量。由图14知，25~32年时间尺度的能量最强、周期最显著，但它的周期变化具有局部性（1980s前）；10~15年时间尺度能量虽然较弱，但周期分布比较明显，几乎占据整个研究时域(1974~2004年)。
7. 绘制小波方差图 7.1小波方差图的绘制
将不同时间尺度下的小波系数代入式(5)可得径流变化的小波方差，以小波方差为纵坐标，时间尺度a为横坐标，可绘制小波方差图（图15）。
1407.2小波方差图在多时间尺度分析中的作用
小波方差图能反映径流时间序列的波动能量随尺120度a的分布情况。可用来确定径流演化过程中存在的100主周期。
80流域径流的小波方差图中（图15）存在4个较为
60明显的峰值，它们依次对应着28年、14年、8年和4
40年的时间尺度。其中，最大峰值对应着28年的时间尺
度，说明28年左右的周期震荡最强，为流域年径流变20化的第一主周期；14年时间尺度对应着第二峰值，为
005101520253035径流变化的第二主周期，第三、第三峰值分别对应着
(d)时间尺度/a 8年和4年的时间尺度，它们依次为流域径流的第三
图15 小波方差图 和第四主周期。这说明上述4个周期的波动控制着流域径流在整个时间域内的变化特征。
小波方差
气象家园提问截图，小波方差图：
http://bbs.06climate.com/forum.php?mod=viewthread&tid=7906&extra=page=1
8. 主周期趋势图的绘制及其在多时间尺度分析中的作用 根据小波方差检验的结果，我们绘制出了控制流域径流演变的第一和第二主周期小波系数图（图16）。从主周期趋势图中我们可以分析出在不同的时间尺度下，流域径流存在的平均周期及丰-枯变化特征。图
相关推荐：
loading