自己要解决的问题是,用小波分解对时间序列进行多尺度分解,再分别对各系数进行预测,最后重构得到最后的预测结果。
因为对系数进行了新的处理,所以不能直接用waverec函数重构原始信号,而需要先构造C,再用waverec函数重构原始信号。
所以这里主要解决的问题是如何用新的系数构造C,进而用waverec函数重构原始信号。
如果有哪个地方理解错误,还望指正
%%首先是不对系数进行另外处理的直接重构原始信号
%小波分解与重构
gasload=importdata('load_new.xls');loaddata=gasload.data.Sheet2;
a=loaddata(:,4);
a=a';
%多尺度一维分解
[C,L]=wavedec(a,3,'db1');
%提取系数
cA3=appcoef(C,L,'db1',3);cD3=detcoef(C,L,3);
cD2=detcoef(C,L,2);
cD1=detcoef(C,L,1);
%重构系数
A3=wrcoef('a',C,L,'db1',3);
D1=wrcoef('d',C,L,'db1',1);
D2=wrcoef('d',C,L,'db1',2);
D3=wrcoef('d',C,L,'db1',3);
A0=waverec