MATLAB中地阈值获取和阈值去噪(超级有用)
Matlab中的信号降噪与滤波技术详解
Matlab中的信号降噪与滤波技术详解正文部分：在信号处理的领域中，信号的降噪和滤波是非常重要的步骤。
Matlab作为一种常用的工具，提供了丰富的信号处理函数和工具箱，可以帮助我们实现高效的信号降噪和滤波。
本文将详细介绍Matlab中的信号降噪和滤波技术。
一、信号降噪技术1.1 经典的降噪方法在信号降噪的过程中，最常用的方法之一是使用滑动平均法。
该方法通过计算信号在一定窗口内的平均值来消除噪声的影响。
在Matlab中，可以使用smooth函数来实现该方法。
使用该函数时，需要指定窗口的大小。
较大的窗口可以平滑信号，但会导致信号的平均值偏移。
而较小的窗口则可以更有效地去除高频噪声，但可能会保留一些低频噪声。
此外，还可以使用中值滤波法来降噪，该方法能够消除信号中的离群值。
在Matlab中，可以使用medfilt1函数实现中值滤波。
该函数需要指定一个窗口大小，并对信号进行中值滤波处理。
较大的窗口可以更好地降噪，但可能会导致信号的细节信息丢失。
1.2 基于小波变换的降噪方法除了经典的降噪方法外，基于小波变换的降噪方法也是一种常用的技术。
小波变换是一种多分辨率分析方法，可以将信号分解为不同尺度的子信号。
在降噪过程中，可以通过滤除高频子信号中的噪声来实现降噪效果。
在Matlab中，可以使用wdenoise函数来实现基于小波变换的降噪。
该函数需要指定小波族，降噪方法和阈值等参数。
1.3 基于自适应滤波的降噪方法自适应滤波是一种根据信号的统计特性进行滤波的方法，它可以根据信号的自相关矩阵来调整滤波器的参数。
在Matlab中，可以使用wiener2函数来实现自适应滤波。
该函数需要指定一个噪声估计器，通过估计信号和噪声的自相关函数来调整滤波器的参数。
二、信号滤波技术2.1 无限脉冲响应滤波器无限脉冲响应滤波器（IIR滤波器）是一种常用的滤波器，它可以对信号进行低通、高通、带通或带阻滤波。
在Matlab中，可以使用butter函数来设计和应用IIR滤波器。
在Matlab中进行噪声抑制和降噪处理的方法
在Matlab中进行噪声抑制和降噪处理的方法引言：噪声是信号处理中的一个常见问题，它可以由多种因素引起，如传感器本身的噪声、电磁干扰等。
噪声的存在会影响到信号的质量和准确性，因此在许多应用中，我们需要进行噪声抑制和降噪处理。
对于Matlab来说，它提供了多种方法和工具来实现这一目标。
本文将介绍在Matlab中进行噪声抑制和降噪处理的方法。
一、频域滤波方法在Matlab中，频域滤波方法是一种常见且有效的噪声抑制和降噪处理方法。
该方法的基本思想是将信号从时域转换到频域，在频域中对信号进行滤波，并将滤波后的信号再转换回时域。
Matlab提供了丰富的频域滤波函数和工具，如fft、ifft、fftshift等。
通过这些函数，我们可以实现低通滤波、高通滤波、带通滤波等各种滤波操作，从而有效抑制和降噪信号。
二、时域滤波方法时域滤波方法是另一种常用的噪声抑制和降噪处理方法。
该方法的基本思想是在时域中对信号进行滤波，直接对信号进行抽样和滤波处理。
与频域滤波不同的是，时域滤波方法更加直观和易于理解。
在Matlab中，我们可以使用filter函数和fir1函数实现时域滤波。
其中，filter函数可以对信号进行FIR滤波，而fir1函数可以设计并生成FIR滤波器。
三、小波变换方法小波变换是一种非常有用的信号处理方法，它可以将信号在时间和频率上进行局部分析。
在噪声抑制和降噪处理中，小波变换可以帮助我们将信号分解成不同的频率成分，并对噪声进行抑制。
在Matlab中，我们可以使用wavelet函数和wdenoise函数来实现小波变换。
通过这些函数，我们可以选择不同的小波基函数，并设置适当的阈值来实现噪声抑制和降噪处理。
四、自适应滤波方法自适应滤波是一种根据信号特性自动调整滤波器参数的滤波方法。
它可以自动识别和适应信号中的噪声，并对其进行抑制和降噪处理。
在Matlab中，自适应滤波可以通过nlms函数和rls函数来实现。
这些函数基于LMS算法和RLS算法，可以快速、准确地对信号进行自适应滤波。
如何使用MATLAB进行频域滤波与去噪
如何使用MATLAB进行频域滤波与去噪使用MATLAB进行频域滤波与去噪引言：在数字信号处理领域，频域滤波与去噪是一项重要而常见的任务，在实际应用中有很多场景需要对信号进行去除噪声或者滤波处理。
MATLAB作为一种功能强大的科学计算软件，提供了丰富的工具和函数，可以帮助我们完成频域滤波与去噪的任务。
本文将介绍如何使用MATLAB进行频域滤波与去噪，并给出一些实用的例子。
一、频域滤波频域滤波是一种常用的信号处理方法，它通过将信号从时域转换到频域，对频域上的信号进行滤波处理，再将滤波后的信号转换回时域得到最终结果。
MATLAB提供了丰富的函数和工具箱，可以方便地进行频域滤波。
1. FFT（快速傅里叶变换）傅里叶变换是将信号从时域转换到频域的一种数学方法，而快速傅里叶变换（FFT）是对离散信号进行傅里叶变换的一种快速算法。
在MATLAB中，可以使用fft函数完成FFT变换，如下所示：```matlabY = fft(X);```其中，X为输入信号，Y为傅里叶变换后的结果。
通过FFT变换，我们可以将信号转换到频域进行进一步的处理。
2. 频域滤波器设计MATLAB提供了fir1、fir2、butter等函数用于设计常见的滤波器，根据滤波器的需求选择合适的函数进行滤波器设计。
以fir1函数为例，它可以设计出一种FIR （有限脉冲响应）滤波器，实现对频域信号的滤波。
下面是一个示例代码：```matlaborder = 32; % 滤波器阶数cutoff = 0.2; % 截止频率b = fir1(order, cutoff); % FIR滤波器设计```在上述代码中，我们指定了滤波器的阶数和截止频率，通过调用fir1函数进行滤波器设计，并得到滤波器的系数b。
将滤波器系数应用到信号上，可以实现对信号的频域滤波。
3. 频域滤波器应用设计好滤波器后，我们可以将滤波器应用到信号上，实现频域滤波。
MATLAB 提供了fftfilt函数用于对信号进行频域滤波，如下所示：```matlabY = fftfilt(b, X);```其中，b为滤波器系数，X为输入信号，Y为滤波后的结果。
利用Matlab进行图像去噪和图像增强
利用Matlab进行图像去噪和图像增强随着数字图像处理技术的不断发展和成熟，图像去噪和图像增强在各个领域都有广泛的应用。
而在数字图像处理的工具中，Matlab凭借其强大的功能和易于使用的特点，成为了许多研究者和工程师首选的软件之一。
本文将介绍如何利用Matlab进行图像去噪和图像增强的方法和技巧。
一、图像去噪图像去噪是指通过一系列算法和技术，将图像中的噪声信号去除或减弱，提高图像的质量和清晰度。
Matlab提供了多种去噪方法，其中最常用的方法之一是利用小波变换进行去噪。
1. 小波变换去噪小波变换是一种多尺度分析方法，能够对信号进行时频分析，通过将信号分解到不同的尺度上，实现对图像的去噪。
在Matlab中，可以使用"dwt"函数进行小波变换，将图像分解为低频和高频子带，然后通过对高频子带进行阈值处理，将噪声信号滤除。
最后通过逆小波变换将去噪后的图像重构出来。
这种方法能够有效抑制高频噪声，保留图像的细节信息。
2. 均值滤波去噪均值滤波是一种基于平均值的线性滤波方法，通过计算像素周围邻域内像素的平均值，替代原始像素的值来去除噪声。
在Matlab中，可以使用"imfilter"函数进行均值滤波，通过设置适当的滤波模板大小和滤波器系数，实现对图像的去噪。
二、图像增强图像增强是指通过一系列算法和技术，改善图像的质量、增强图像的细节和对比度，使图像更容易被观察和理解。
Matlab提供了多种图像增强方法，以下将介绍其中的两种常用方法。
1. 直方图均衡化直方图均衡化是一种通过对图像像素值的分布进行调整，增强图像对比度的方法。
在Matlab中，可以使用"histeq"函数进行直方图均衡化处理。
该函数能够将图像的像素值分布拉伸到整个灰度级范围内，提高图像的动态范围和对比度。
2. 锐化增强锐化增强是一种通过增强图像边缘和细节来改善图像质量的方法。
在Matlab中，可以使用"imsharpen"函数进行图像的锐化增强处理。
matlab数据选取阈值
matlab数据选取阈值
摘要：
1.引言
2.Matlab 数据选取阈值方法
a.基于最大值和最小值的方法
b.基于均值和标准差的方法
c.基于经验法则的方法
d.基于自定义函数的方法
3.各种方法的优缺点分析
4.结论
正文：
Matlab 作为一款功能强大的数学软件，在数据处理方面有着广泛的应用。
在数据处理过程中，选取合适的阈值对结果具有重要的影响。
本文将介绍几种在Matlab 中选取阈值的方法。
1.基于最大值和最小值的方法
该方法是最简单的阈值选取方法，通过找到数据中的最大值和最小值，将数据分成两部分。
优点是简单易行，适用于数据量较小的情况。
缺点是对于具有多个峰值的数据，该方法可能无法准确划分。
2.基于均值和标准差的方法
该方法通过计算数据的均值和标准差，将数据分成两部分。
当数据分布近似于正态分布时，该方法效果较好。
优点是适用于大部分数据分布情况，缺点
是对异常值较敏感。
3.基于经验法则的方法
经验法则根据数据中的大部分数据点来确定阈值。
常用的经验法则有三种：Q1 和Q3 方法、IQR（四分位距）方法和Median 方法。
这些方法的优点是适用于各种数据分布情况，缺点是对于极端值较多的数据，阈值可能不太准确。
4.基于自定义函数的方法
该方法是通过编写自定义函数来选取阈值。
例如，可以编写一个函数，根据数据的形状、波动程度等因素来确定阈值。
这种方法的优点是可以针对具体数据集进行优化，缺点是需要一定的编程基础。
总之，在Matlab 中选取阈值的方法有多种，各有优缺点。
matlab中的去噪函数
matlab中的去噪函数
MATLAB中有多种用于去噪的函数和工具，其中一些常用的包括：
1. `medfilt1`，这个函数用于对一维信号进行中值滤波，可以
有效地去除椒盐噪声和其他类型的噪声。
2. `medfilt2`，类似于`medfilt1`，这个函数用于对二维图像
进行中值滤波，对于去除图像中的斑点噪声和其他类型的噪声效果
很好。
3. `wiener2`，这个函数实现了维纳滤波器，可以用于图像的
去噪。
维纳滤波器是一种自适应滤波器，可以根据图像的局部特性
进行滤波，适用于各种类型的噪声。
4. `imfilter`，这个函数可以实现各种类型的滤波操作，包括
高斯滤波、均值滤波等，可以根据具体的需求选择合适的滤波器进
行去噪处理。
5. `denoiseWavelet`，MATLAB中还提供了基于小波变换的去
噪函数，可以通过小波阈值处理来去除信号中的噪声成分。
除了以上提到的函数，MATLAB还提供了一些图像处理工具箱，
其中包含了更多高级的去噪算法和工具，比如基于深度学习的去噪
方法、非局部均值去噪（NL-means denoising）等。
这些工具可以
根据具体的应用场景和需求选择合适的去噪方法进行处理。
总的来说，MATLAB提供了丰富的去噪函数和工具，可以根据具
体的信号或图像特性选择合适的方法进行去噪处理。
在实际应用中，需要根据噪声类型、信噪比以及对信号质量的要求来选择合适的去
噪方法。
matlab去噪声方法
matlab去噪声方法MATLAB中去噪声的方法有很多，以下列举了一些常用的方法：1. 均值滤波：均值滤波是一种简单且有效的去噪声方法，它通过计算邻域内像素的平均值来减小噪声。
具体操作是创建一个与输入图像大小相同的零矩阵，然后遍历图像的每个像素，将邻域内的像素值求和，最后除以邻域内像素的数量，得到滤波后的像素值。
2. 中值滤波：中值滤波主要用于去除图像中的脉冲噪声和椒盐噪声。
该方法的核心思想是将图像中相邻像素的灰度值进行排序，然后取中间值作为滤波后的像素值。
3. 高斯滤波：高斯滤波是一种广泛应用于图像去噪的方法，它通过在图像上滑动一个高斯核函数来降低噪声。
高斯核函数的宽度决定了滤波的效果，宽度越大，去噪效果越明显，但同时也会导致图像变得模糊。
4. 双边滤波：双边滤波是一种基于邻域的滤波方法，它同时考虑了像素的空间距离和灰度差异。
通过在图像上滑动一个双边核函数，可以有效地去除噪声并保留图像的细节。
5. 基于小波变换的方法：小波变换可以将图像分解为不同尺度、方向和频率的子带，通过对子带的处理，可以有效地去除图像中的噪声。
常用的有小波分解、小波合成和小波去噪等方法。
6. 基于深度学习的方法：深度学习方法，如卷积神经网络（CNN）和循环神经网络（RNN），在图像去噪领域也取得了很好的效果。
通过训练神经网络，可以学习到图像的复杂特征，从而在去噪的同时保留图像的细节。
在MATLAB中实现这些方法，可以利用内置的图像处理函数或自行编写代码。
例如，使用imfilter函数实现均值滤波，使用im2double函数将图像转换为双精度值等。
同时，可以借助图像处理工具箱中的各种滤波器和图像读取、显示函数，如sobel、roberts、prewitt算子等，来实现特定的去噪效果。
matlab阈值处理
matlab阈值处理MATLAB中的阈值处理是图像处理中常用的一种方法，用于将图像转换为二值图像。
阈值处理通常基于图像的像素强度，并在给定阈值的情况下将像素分为两个类别。
本文将从阈值处理的原理、MATLAB的阈值处理函数、阈值选择方法和实际案例分析等方面详细介绍MATLAB中的阈值处理技术。
一、阈值处理原理阈值处理是一种基于像素强度的图像分割方法。
其基本原理是根据像素灰度值是否超过阈值，将像素分为两类：一类属于背景，一类属于目标。
阈值的选择对于最终的图像分割结果至关重要。
通过选择适当的阈值，我们可以从图像中分割出想要的目标或者特定区域，为后续的图像处理任务提供良好的基础。
二、MATLAB中的阈值处理函数在MATLAB中，有多种用于阈值处理的函数可供使用。
其中，最常用的是`imbinarize`函数。
该函数的基本语法如下：BW = imbinarize(I, level)其中，`I`表示输入的灰度图像，`level`表示阈值。
该函数将返回一个二值图像`BW`，其中超过阈值的像素值设置为1，否则设置为0。
除了`imbinarize`函数，MATLAB还提供了其他常用的阈值处理函数，如`graythresh`、`im2bw`等函数。
这些函数可以根据用户的需求和应用场景选择合适的阈值处理方法。
三、阈值选择方法阈值的选择是阈值处理中一个关键的问题。
合理选择阈值能够有效地分离目标和背景。
下面介绍几种常用的阈值选择方法。
1. 固定阈值法：即将阈值固定为一个预先确定的数值。
该方法简单直观，但需要根据具体问题进行调整。
根据图像的特点和经验，可以选择一个合适的阈值作为全局阈值，并将图像中超过该阈值的像素设置为目标。
2. 基于直方图的阈值法：通过分析图像的直方图，可以得到图像的灰度分布特征。
根据直方图的波峰和波谷，可以选择一个合适的阈值。
常用的方法有Otsu法、Yen法等。
3. 自适应阈值法：该方法根据图像的局部特征，进行局部阈值的选择。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
1.阈值获取MATLAB中实现阈值获取的函数有ddencmp、thselect、wbmpen和wwdcbm，下面对它们的用法进行简单的说明。
（1）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是熵名（只在选择小波包时使用）。
（2）函数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)。
（3）函数wbmpen的调用格式如下：THR=wbmpen(C,L,SIGMA,ALPHA);THR=wbmpen(C,L,SIGMA,ALPHA)返回去噪的全局阈值THR。
THR 通过给定的一种小波系数选择规则计算得到，小波系数选择规则使用Birge-Massart的处罚算法。
{C,L]是进行去噪的信号或图像的小波分解结构；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)（4）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.2.信号的阈值去噪MATLAB中实现信号的阈值去噪的函数有wden、wdencmp、wthresh、wthcoef、wpthcoef以及wpdencmp。
下面对它们的用法作简单的介绍。
（5）函数wden的调用格式有以下两种：[XD,CXD,LXD]=wden(X,TPTR,SORH,SCAL,N,'wname')[XD,CXD,LXD]=wden(C,L,TPTR,SORH,SCAL,N,'wname')函数wden用于一维信号的自动消噪。
X为原始信号，[C,L]为信号的小波分解，N为小波分解的层数。
THR为阈值选择规则：*TPTR='rigrsure'，自适应阈值选择使用Stein的无偏风险估计原理。
*TPTR='heursure'，使用启发式阈值选择。
*TPTR='sqtwolog'，阈值等于sqrt(2*log(length(X))).*TPTR='minimaxi'，用极大极小原理选择阈值。
SORH是软阈值或硬阈值的选择（分别对应's'和'h'）。
SCAL指所使用的阈值是否需要重新调整，包含下面三种：*SCAL='one'不调整；*SCAL='sln'根据第一层的系数进行噪声层的估计来调整阈值。
*SCAL='mln'根据不同的噪声估计来调整阈值。
XD为消噪后的信号，[CXD,LXD]为消噪后信号的小波分解结构。
格式（1）返回对信号X经过N层分解后的小波系数进行阈值处理后的消噪信号XD和信号XD的小波分解结构[CXD,LXD]。
格式（2）返回参数与格式（1）相同，但其结构是由直接对信号的小波分解结构[C,L]进行阈值处理得到的。
（6）函数wdencmp的调用格式有以下三种：(1)[XC,CXC,LXC,PERF0,PERFL2]=wdencmp('gbl',X,'wname',N,TH TR,SORH,KEEPAPP);(2)[XC,CXC,LXC,PERF0,PERFL2]=wdencmp('lvd',X,'wname',N,TH TR,SORH);(3)[XC,CXC,LXC,PERF0,PERFL2]=wdencmp('lvd',C,L,'wname',N, THTR,SORH);函数wdencmp用于一维或二维信号的消噪或压缩。
wname是所用的小波函数，gbl(global的缩写)表示每一层都采用同一个阈值进行处理，lvd表示每层采用不同的阈值进行处理，N表示小波分解的层数，THR为阈值向量，对于格式（2）和（3）每层都要求有一个阈值，因此阈值向量THR的长度为N，SORH表示选择软阈值或硬阈值（分别取值为's'和'h'），参数KEEPAPP取值为1时，则低频系数不进行阈值量化，反之，低频系数要进行阈值量化。
XC是要进行消噪或压缩的信号，[CXC,LXC]是XC的小波分解结构，PERF0和PERFL2是恢复或压缩L^2的范数百分比。
如果[C,L]是X的小波分解结构，则PERFL2=100*(CXC向量的范数/C向量的范数)^2；如果X是一维信号，小波wname是一个正交小波，则PERFL2=100||XC||^2/||X||^2。
（7）函数wthresh的调用格式如下：Y=wthresh(X,SORH,T)Y=wthresh(X,SORH,T) 返回输入向量或矩阵X经过软阈值（如果SORH='s'）或硬阈值（如果SORH='h'）处理后的信号。
T是阈值。
Y=wthresh(X,'s',T)返回的是Y=SIG(X)*(|X|-T)+，即把信号的绝对值与阈值进行比较，小于或等于阈值的点变为零，大于阈值的点为该点值与阈值的差值。
Y=wthresh(X,'h',T)返回的是Y=X*1(|X|>T)，即把信号的绝对值和阈值进行比较，小于或等于阈值的点变为零，大于阈值的点保持不变。
一般来说，用硬阈值处理后的信号比用软阈值处理后的信号更粗糙。
（8）函数wthcoef的调用格式下面四种：（1）NC=wthcoef('d',C,L,N,P)（2）NC=wthcoef('d',C,L,N)（3）NC=wthcoef('a',C,L)（4）NC=wthcoef('t',C,L,N,T,SORH)函数wthcoef用于一维信号小波系数的阈值处理。
格式（1）返回小波分解结构[C,L]经向量N和P定义的压缩率处理后的新的小波分解向量NC，[NC,L]构成一个新的小波分解结构。
N包含被压缩的细节向量，P是把较小系数置0的百分比信息的向量。
N和P的长度必须相同，向量N必须满足1<=N(i)<=length(L)-2。
格式（2）返回小波分解结构[C,L]经过向量N中指定的细节系数置0后的小波分解向量NC。
格式（3）返回小波分解结构[C,L]经过近似系数置0后的小波分解向量NC。
格式（4）返回小波分解结构[C,L]经过将向量N作阈值处理后的小波分解向量NC。
如果SORH=’s‘，则为软阈值；如果SORH='h'则为硬阈值。
N包含细节的尺度向量，T是N相对应的阈值向量。
N 和T的长度必须相等。
（9）函数wpdencmp的调用格式有以下两种：[XD,TREED,PERF0,PERFL2]=wpdencmp(X,SORH,N,'wname',CRIT,PAR, KEEPAPP)[XD,TREED,PERF0,PERFL2]=wpdencmp(TREE,SORH,CRIT,PAR,KEEPAPP )函数wpdencmp用于使用小波包变换进行信号的压缩或去噪。
格式（1）返回输入信号X（一维或二维）的去噪或压缩后的信号XD。
输出参数TREED是XD的最佳小波包分解树；PERFL2和PERF0是恢复和压缩L2的能量百分比。
PERFL2=100*(X的小波包系数范数/X的小波包系数)^2；如果X是一维信号，小波wname是一个正交小波，则PERFL2=100*||XD||^2/||X||^2。
SORH的取值为's'或'h'，表示的是软阈值或硬阈值。
输入参数N是小波包的分解层数，wname是包含小波名的字符串。
函数使用由字符串CRIT定义的熵和阈值参数PAR实现最佳分解。
如果KEEPAPP=1，则近似信号的小波系数不进行阈值量化；否则，进行阈值量化。
格式（2）与格式（1）的输出参数相同，输入选项也相同，只是它从信号的小波包分解树TREE进行去噪或压缩。
（1）Thselect函数：选取用于小波消噪处理的阈值调用方式Thr=thselect(x,tptr);根据信号x和阈值选择标准tptr来确定一个消噪处理过程中所采用的自适应阈值。