例程13-13.m

来自「《小波分析理论与Matlab7实现》书籍源码」· M 代码 · 共 16 行

%装载信号
load leleccum;  
%截取信号中的一段[2600，3100]
s=leleccum(2600:3100);  
%用小波db3对s进行三层分解
[c,l]=wavedec(s,3,'db3');  
%选用全局阈值进行信号压缩处理
thr=40;  
[sd,csd,lsd,perf0,perfl2]=wdencmp('gbl',c,l,'db3',3,thr,'h',1);  
subplot(2,1,1);  
plot(s);  
title('原始信号');  
subplot(2,1,2);  
plot(sd);  
title('压缩后的信号');