Matlab程序小波分析时间序列数据

以美国某气象站1894~2010年连续的年降水量为例，试应用小波分析，完成如下任务：①小波变换系数；②绘制小波系数实部等值线图；③绘制小波系数模和模方等值线图；④绘制小波方差图；以及⑤绘制不同时间尺度的小波实部过程线。所谓年降水量时间序列的多时间尺度是指：年降水量在演化过程中，并不存在真正意义上的变化周期，而是其变化周期随着研究尺度的不同而发生相应的变化，这种变化一般表现为小时间尺度的变化周期往往嵌套在大尺度的变化周期之中。也就是说，年降水量变化在时间域中存在多层次的时间尺度结构和局部变化特征。Details: http://blog.sciencenet.cn/blog-1148346-794768.html

% Created by LI Xu
% 15 May 2014
% Version 1.0

% Retrieve the entended numer ahead and tail


function [st ed]=entendvalue（ShortArr LongArr）

% ShortArr=Prec;
% LongArr=ePrec;

num_match=length（ShortArr）;
num_matched=length（LongArr）;

for ii=1:num_matched-num_match+1
    a=（ShortArr==LongArr（ii:ii+num_match-1））;
    count_same=sum（a）;
    if count_same==num_match
        st=ii;
        ed=ii+num_match-1;
        fprintf（‘Original Arrary:%d~%d\n‘ st ed）;
        break;
    end
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件     102984  2015-03-23 09:39  WavletAnalysis\6_var.jpg

     文件     126349  2014-05-15 10:11  WavletAnalysis\cePrec.mat

     文件        532  2014-05-15 09:52  WavletAnalysis\entendvalue.m

     文件        734  2014-05-15 09:24  WavletAnalysis\ePrec.mat

     文件       1133  2014-05-15 16:47  WavletAnalysis\LI_contourf.m

     文件       1047  2014-05-15 20:39  WavletAnalysis\LI_period_plot.asv

     文件        630  2014-05-15 20:40  WavletAnalysis\LI_period_plot.m

     文件        889  2014-05-15 15:19  WavletAnalysis\LI_plot.m

     文件        369  2014-05-15 17:06  WavletAnalysis\LI_Var.m

     文件       2526  2014-05-15 20:39  WavletAnalysis\main.asv

     文件       3262  2015-03-23 09:37  WavletAnalysis\main.m

     文件      41975  2014-05-15 15:50  WavletAnalysis\pictures\0_trend.png

     文件     147343  2014-05-15 09:20  WavletAnalysis\pictures\1.png

     文件     204914  2014-05-15 10:10  WavletAnalysis\pictures\2.png

     文件     513234  2014-05-15 21:03  WavletAnalysis\pictures\3_real_contourf.jpg

     文件     457762  2014-05-15 21:05  WavletAnalysis\pictures\4_mod_contourf.jpg

     文件     456891  2014-05-15 21:05  WavletAnalysis\pictures\5_mod_square_contourf.jpg

     文件     115342  2014-05-15 21:06  WavletAnalysis\pictures\7_55a_coefficient.jpg

     文件     131500  2014-05-15 21:06  WavletAnalysis\pictures\8_30a_coefficient.jpg

     文件        718  2014-05-15 08:47  WavletAnalysis\Prec.mat

     文件        363  2014-05-15 08:48  WavletAnalysis\Year.mat

     目录          0  2015-03-23 09:39  WavletAnalysis\pictures

     目录          0  2015-03-23 09:39  WavletAnalysis

----------- ---------  ---------- -----  ----

              2310497                    23