使用小波变换重构某个节点的信号

转载：http://blog.163.com/lxfan_0406@126/blog/static/6096992620125713240582/

1.[d,a]=wavedec(y,3,'db5'); %信号分解

1）wrcoef 由多层小波分解重构某一层的分解信号；

2）waverec 直接重构原始信号

a3=wrcoef('a',d,a,'db5',3);

d3=wrcoef('d',d,a,'db5',3);

d2=wrcoef('d',d,a,'db5',2);

d1=wrcoef('d',d,a,'db5',1); %X = wrcoef('type',C,L,'wname',N)

sx=waverec(d,a,'db5');

如果分解三层，则

信号：C=CA3+CD3+CD2+CD1，

长度标示：L=[size(CA3),size(CD3),size(CD2),size(CD1),size(C)]

如果要求解节点CA3，则使用a3=wrcoef('a',d,a,'db5',3);

如果要求解节点CD3，则使用d3=wrcoef('d',d,a,'db5',3);

如果要求解节点CD2，则使用d2=wrcoef('d',d,a,'db5',2);

如果要求解节点CD1，则使用d1=wrcoef('d',d,a,'db5',1);

%////////////////////////////////////////////////////////////////////////

wpt = wpdec(y,3,'db5');   

cfs_a3= wprcoef(wpt,[3 0]);     % 得到CA3

cfs_d1= wprcoef(wpt,[1 1]);     % 得到CD1

cfs_d2= wprcoef(wpt,[2 1]);     % 得到CD2

cfs_d3= wprcoef(wpt,[3 1]);     % 得到CD3