% 裝載bumps的噪聲信號
load noisbump; x = noisbump;
% 使用sym6執行5層小波分解
wname = 'sym6'; lev = 5;
[c,l] = wavedec(x,lev,wname);
% 由第1層的高頻系數估計噪聲標準差
sigma = wnoisest(c,l,1);
% 信號降噪時，采用調整參數選擇全局閾值
alpha = 2;
thr = wbmpen(c,l,sigma,alpha)
% 使用軟閾值和保存的低頻信號，進行信號降噪
keepapp = 1;
xd = wdencmp('gbl',c,l,wname,lev,thr,'s',keepapp);
% 畫出原始信號和降噪后的信號
figure(1)
subplot(211), plot(x), title('原始信號')
subplot(212), plot(xd), title('降噪后的信號')