signal_compression1.m

来自「MATLAB 小波分析高级技术——书中的一些应用程序」· M 代码 · 共 15 行

load nelec;
indx = 1:1024; 
x = nelec(indx);
%用小波函数haar对信号进行3层分解
[c,l] = wavedec(x,3,'haar');      
alpha = 1.5; 
%获取信号压缩的阀值
[thr,nkeep] = wdcbm(c,l,alpha)    
%对信号进行压缩
[xd,cxd,lxd,perf0,perfl2] = wdencmp('lvd',c,l,'haar',3,thr,'s');   
subplot(211), plot(indx,x); 
title('原始信号');
subplot(212), plot(indx,xd); 
title('压缩后的信号');