小波降噪的函数.docVIP

小波降噪的函数 小波变换——压缩与去噪算法程序2009-10-04 14:27 最近忙着做方案设计，还有国庆的杂事比较多，也好久没上空间了。呵呵今天拿出来点东西晒晒，有关小波变换的一些应用，压缩和滤除噪声方面的。 ＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃ ddencmp函数 % 当前延拓模式时补零 % 产生高斯白噪声 init = 2055415866; randn(seed,init); x = randn(1,1000); % 寻找使用小波降噪输出的3个默认值 disp(小波降噪输出的3个默认值) [thr,sorh,keepapp] = ddencmp(den,wv,x) % 寻找使用小波压缩输出的3个默认值 disp(小波压缩输出的3个默认值) [thr,sorh,keepapp] = ddencmp(cmp,wv,x) % 寻找使用小波包降噪输出的3个默认值 disp(小波包降噪输出的4个默认值) [thr,sorh,keepapp,crit] = ddencmp(den,wp,x) % 寻找使用小波包压缩输出的3个默认值 disp(小波包压缩输出的4个默认值) [thr,sorh,keepapp,crit] = ddencmp(cmp,wp,x) ＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃ thselect函数 % 当前的延拓模式是补零 % 产生高斯白噪声 init = 2055415866; randn(seed,init); x = randn(1,1000); % 寻找各个阈值选择规则的阈值 % 使用SURE的自适应阈值 thr1 = thselect(x,rigrsure) % 固定形式阈值 thr2 = thselect(x,sqtwolog) % 启发式阈值 thr3 = thselect(x,heursure) % 极小极大阈值 thr4 = thselect(x,minimaxi) ＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃ wbmpen函数 % 装载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(降噪后的信号) ＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃ wdcbm2函数 % 装载原始图像 load belmont2; nbc = size(map,1); % 使用sym4小波对图像进行3层分解 wname = sym4; lev = 3; [c,s] = wavedec2(X,lev,wname); % 图像压缩，使用推荐参数选择各层不同的阈值 alpha = 1.5; m = 2.7*prod(s(1,:)); [thr,nkeep] = wdcbm2(c,s,alpha,m) % 使用上面的阈值和硬阈值处理方式进行图像压缩 [xd,cxd,sxd,perf0,perfl2] =wdencmp(lvd,c,s,wname,lev,thr,h); % 画出原始图像和压缩后的图像 colormap(pink(nbc)); subplot(121), image(wcodemat(X,nbc)), title(原始图像) subplot(122), image(wcodemat(xd,nbc)), title(压缩后的图像) disp(恢复百分比: ); num2str(perfl2) disp(压缩百分比: ); num2str(perf0) ＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃＃ wdcbm函数 % 装载信号并选择其中一段 load leleccum; indx = 2600:3100; x = leleccum(indx); % 使用db3在第5层执行信号的小波分解 wname = db3; lev =