實例2：

首先使用函數(shù)wnoisest獲取噪聲方差，然后使用函數(shù)wbmpen獲取小波去噪閾值，最后使用wdencmp實現(xiàn)信號消噪。

load leleccum;

indx = 1:1024;

x = leleccum（indx）;

%產(chǎn)生含噪信號

init = 2055615886;

randn（‘seed’，init）;

nx = x + 18*randn（size（x））;

%使用小波函數(shù)‘db6’對信號進行3層分解

［c，l］ = wavedec（nx，3，‘db6’）;

%估計尺度1的噪聲標準差

sigma = wnoisest（c，l，1）;

alpha = 2;

%獲取消噪過程中的閾值

thr = wbmpen（c，l，sigma，alpha）;

keepapp = 1;

%對信號進行消噪

xd = wdencmp（‘gbl’，c，l，‘db6’，3，thr，‘s’，keepapp）;

subplot（221）;

plot（x）;

title（‘原始信號’）;

subplot（222）;

plot（nx）;

title（‘含噪信號’）;

subplot（223）;

plot（xd）;

title（‘消噪后的信號’）;