1.代码
clear
clc
%Load a signal.
load noisdopp
%Obtain the optimal wavelet packet tree based on an entropy criterion.
wpt=wpdec(noisdopp,3,'db1','shannon');
BstTree=besttree(wpt);
plot(BstTree)
%Reconstruct the signal from the optimal wavelet packet tree
xrec=wprec(BstTree);
plot(xrec)
title('最优树重建信号')
%小波包降噪
[thr1,sorh1,keepapp1,crit]=ddencmp('den','wp',noisdopp);
xwpd=wpdencmp(noisdopp,sorh1,3,'db2',crit,thr1,keepapp1);
%小波降噪
[thr2,sorh2,keepapp2]=ddencmp('den','wv',noisdopp);
xcwv=wdencmp('gbl',noisdopp,'db2',3,thr2,sorh2,keepapp2);
figure(3)
subplot(3,1,1)
plot(noisdopp)
title('原始信号');
subplot(3,1,2)
plot(xwpd)
title('小波包降噪');
subplot(3,1,3)
plot(xcwv)
title('小波降噪');2.运行结果
