program_13_02.m

这是基于MATLAB系统上开发的,应用于小波分析方面,很有用的程序模块,书名为"MATLAB在小波分析中的应用"

% 当前延拓模式是补零
% 产生高斯白噪声
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