小波动分析及构建单变量时间序列流程

网友‘李明杨艳’指出，我之前三个版本的小波信号分解重构程序中，针对一维信号分解与重构的函数mydwt和myidwt均存在一个严重的错误。当时编写代码时，仅依据haar小波的特性进行设计，未充分考虑在使用其他小波函数滤波器组时，卷积运算输出序列长度的变化问题。后续版本主要关注二维图像处理，忽略了对一维信号的修正。目前该问题已得到解决，更新后的程序代码如下：

function [c,l] = mydwt(x,lpd,hpd,dim);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% 函数 [c,l]=MYDWT(X,LPD,HPD,DIM) 对输入序列x执行一维离散小波分解操作，并返回分解后的系数序列c

% 输入参数：x——输入数据序列；

% lpd——低通滤波器参数；

% hpd——高通滤波器参数；

% dim——设定的小波分解层数。

% 输出参数：c——生成的小波分解系数序列；

% l——记录各级小波分解系数的长度信息。

%

% 版权声明：Zou Yuhua ( chenyusiyuan ) ，创建时间：2007-11-10，

% 修改时间：2008-09-01

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

cA=x;

c=[];

l=length(x);

fo