morlet小波分析matlabM文件

function [ra2vcoef]=cmor（x）
%Morlet复值小波分析，母小波的选取参考百度文库《宜昌水文站年径流量演变多时间尺度分析》
%输入部分：x为输入原始连续时间序列，例如年降水量序列
%输出部分：coef为输出系数矩阵，r为系数实部矩阵，a为系数模矩阵，a2为模方，v为小波方差
%%
%fi0为确定母小波的常数
%t为周期序列，T0为周期初始值，Tn为周期末值默认取序列的3/5，Tp为步长。时间尺度a的取值已在程序中换算
%时间间隔dt默认取1
%平移因子b默认按时间序列的序号依次取值
l=length（x）;
fi0=6.2;
T0=2;
Tn=2*fix（l*3/10）;
Tp=2;
dt=1;
%%
%时间序列x中心化，使其均值为0，标准差为1
m=mean（x）;
s=std（x1）;
for k=1:l
    x（k）=（x（k）-m）/s
end
%两端对称延长时间序列x。平移原序列至中间位置
for k=（l+1）:2*l
    x（k）=x（k-l）;
end 
%首端对称延长时间序列x
for k=1:l
    x（k）=x（2*l+1-k）;
end
%末端对称延长时间序列x
for k=（2*l+1）:3*l
    x（k）=x（4*l+1-k）;
end
%%
%构建系数空矩阵
coef=zeros（（Tn-T0）/Tp+13*l）;
%小波系数的计算
for b=1:3*l
    for t=T0:Tp:Tn
        a=t/（4*pi/（fi0+（fi0^2+2）^0.5））;
        for k=1:3*l
            fimi=（k*dt-b）/a;
            fimu=c