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