小波去噪代码.docVIP

- PAGE PAGE 2 欢迎下载 例1： load leleccum; index = 1:1024; x = leleccum(index); %产生噪声信号 init = 2055615866; randn(seed,init); nx = x + 18*randn(size(x)); %获取消噪的阈值 [thr,sorh,keepapp] = ddencmp(den,wv,nx); %对信号进行消噪 xd = wdencmp(gbl,nx,db4,2,thr,sorh,keepapp); subplot(221); plot(x); title(原始信号); subplot(222); plot(nx); title(含噪信号); subplot(223); plot(xd); title(消噪后的信号); 例2： 本例中，首先使用函数wnoisest获取噪声方差，然后使用函数wbmpen获取小波去噪阈值，最后使用wdencmp实现信号消噪。 load leleccum; indx = 1:1024; x = leleccum(indx); %产生含噪信号 init = 2055615886; randn(seed,init); nx = x + 18*randn(size(x)); %使用小波函数db6对信号进行3层分解 [c,l] = wavedec(nx,3