[c,l]=wavedec(f,10,'dmey');%对f进行10尺度小波分解
ca10=appcoef(c,l,'dmey',10);%提取第10尺度上的低频小波分解系数,即逼近分量
cd10=detcoef(c,l,10); %提取各尺度上的高频小波分解系数,即细节分量
cd9=detcoef(c,l,9);
cd8=detcoef(c,l,8);
cd7=detcoef(c,l,7);
cd6=detcoef(c,l,6);
cd5=detcoef(c,l,5);
cd4=detcoef(c,l,4);
cd3=detcoef(c,l,3);
cd2=detcoef(c,l,2);
cd1=detcoef(c,l,1);
caa10=zeros(1,length(ca10));
cdd10=zeros(1,length(cd10));%把第1尺度上的高频小波系数置零
cdd9=zeros(1,length(cd9));
cdd4=zeros(1,length(cd4));%把第4尺度上的高频小波系数置零
cdd3=zeros(1,length(cd3));
cdd2=zeros(1,length(cd2));
cdd1=zeros(1,length(cd1));
c1=[caa10,cd10,cd9,cd8,cd7,cd6,cd5,cdd4,cdd3,cdd2,cdd1];%新的小波分解结构
f1=waverec(c1,l,'dmey'); |