% 當前延拓模式是補零
% 產生高斯白噪聲
init = 2055415866; randn('seed',init); 
x = randn(1,1000);
% 使用db3小波對x進行2層分解
[c,l] = wavedec(x,2,'db3');
% 在第1層和第2層估計系數的標準差
% 既然x是方差為1的高斯白噪聲，那么估計值應接近1
sd1=wnoisest(c,l,1:2)
% 現假設x包含10個非正常值
ind = 50:50:500;
x(ind) = 100 * ones(size(ind));
% 使用db3小波對x進行1層分解
[ca,cd] = dwt(x,'db3');
% 普通的cd標準差估計方法
sd2=std(cd)
% cd標準差的魯棒估計方法
sd3=median(abs(cd))/0.6745