小波去噪的matlab指令应用1（阀值获取）

南风晓晓

　　MATLAB小波阈值去噪，主要包括阈值获取和阈值去噪两方面。(受帖子字数限制，本部分只包括阀值获取，阀值去噪部分在另一个帖子中)
　　1.阈值获取
　　MATLAB中实现阈值获取的函数有ddencmp、thselect、wbmpen和wwdcbm，下面对它们的用法进行简单的说明。
　　①ddencmp的调用格式有以下三种：
　　(1)[THR,SORH,KEEPAPP,CRIT]=ddencmp(IN1,IN2,X)
　　(2)[THR,SORH,KEEPAPP,CRIT]=ddencmp(IN1,'wp',X)
　　(3)[THR,SORH,KEEPAPP,CRIT]=ddencmp(IN1,'wv',X)
　　函数ddencmp用于获取信号在消噪或压缩过程中的默认阈值。输入参数X为一维或二维信号;IN1取值为'den'或'cmp'，'den'表示进行去噪，'cmp'表示进行压缩;IN2取值为'wv'或'wp'，wv表示选择小波，wp表示选择小波包。返回值THR是返回的阈值;SORH是软阈值或硬阈值选择参数;KEEPAPP表示保存低频信号;CRIT是熵名(只在选择小波包时使用)。
　　②函数thselect的调用格式如下：
　　THR=thselect(X,TPTR);
　　THR=thselect(X,TPTR)根据字符串TPTR定义的阈值选择规则来选择信号X的自适应阈值。
　　自适应阈值的选择规则包括以下四种：
　　*TPTR='rigrsure'，自适应阈值选择使用Stein的无偏风险估计原理。
　　*TPTR='heursure'，使用启发式阈值选择。
　　*TPTR='sqtwolog'，阈值等于sqrt(2*log(length(X))).
　　*TPTR='minimaxi'，用极大极小原理选择阈值。
　　阈值选择规则基于模型 y = f(t) + e，e是高斯白噪声N(0,1)。
　　③函数wbmpen的调用格式如下：
　　THR=wbmpen(C,L,SIGMA,ALPHA);
　　THR=wbmpen(C,L,SIGMA,ALPHA)返回去噪的全局阈值THR。THR通过给定的一种小波系数选择规则计算得到，小波系数选择规则使用Birge-Massart的处罚算法。是进行去噪的信号或图像的小波分解结构;SIGMA是零均值的高斯白噪声的标准偏差;ALPHA是用于处罚的调整参数，它必须是一个大于1的实数，一般取ALPHA=2。
　　设t*使crit(t)=-sum(c(k)^2,k<=t) + 2 * SIGMA^2 * t*(ALPHA+log(n/t))的最小值，其中c(k)是按绝对值从大到小排列的小波包系数，n是系数的个数，则THR=|c(t*)|。
　　wbmpen(C,L,SIGMA,ALPHA,ARG)计算阈值并画出三条曲线。
　　2 * SIGMA^2 * t*(ALPHA+log(n/t))
　　sum(c(k)^2, k<=t)
　　crit(t)
　　④wdcbm的调用格式有以下两种：
　　(1)[THR,NKEEP]=wdcbm(C,L,ALPHA);
　　(2)[THR,NKEEP]=wdcbm(C,L,ALPHA,M);
　　函数wdcbm是使用Birge-Massart算法获取一维小波变换的阈值。返回值THR是与尺度无关的阈值，NKEEP是系数的个数。[C,L]是要进行压缩或消噪的信号在j=length(L)-2层的分解结构;LAPHA和M必须是大于1的实数;THR是关于j的向量，THR(i)是第i层的阈值;NKEEP也是关于j的向量，NKEEP(i)是第i层的系数个数。一般压缩时ALPHA取1.5，去噪时ALPHA取3.

和丹

好文，有用

王虎

太谢谢了~~

bzp小鹏

Good！！！！！

Luis_sun

感谢分享