clc;

clear;

%载入信号

load noisdopp;

s=noisdopp;

%对信号s作3层平稳小波分解，长度仍与原信号长度相同

[swa,swd]=swt(s,3,'db8');
mzero=zeros(size(swd));
A=mzero;

%还原第三层的近似系数；此处不能用swa（3，：）表示，因为由于分解算法的原因，会产生冗余
A(3,:)=iswt(swa,mzero,'db8');
%还原各层的细节系数

D=mzero;
for i=1:3
    swcfs=mzero;
    swcfs(i,:)=swd(i,:);
    D(i,:)=iswt(mzero,swcfs,'db8');
end

%还原第2层，第1层的细节系数
A(2,:)=A(3,:)+D(3,:);
A(1,:)=A(2,:)+D(2,:);

%还原原始信号
A0=A(1,:)+D(1,:);
%计算原信号与重构信号的平方差
norm（A0-s）
plot((0:1023),s,(0:1023),A0);