转载: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