wden 函数：一维小波去噪

function [xd,cxd,lxd,thrs] = wden(in1,in2,in3,in4,in5,in6,in7)

描述

wden 函数是 MATLAB 中用于执行一维信号小波去噪的函数。它使用不同的阈值规则和阈值类型来对小波系数进行阈值处理，并返回去噪后的信号。

语法

XD = WDEN(X,TPTR,SORH,SCAL,N,WNAME) 返回输入信号 X 的去噪版本 XD，该版本通过对小波系数进行阈值处理获得。

TPTR 是作为字符串指定的阈值选择规则。TPTR 的支持选项如下：

'modwtsqtwolog' 使用最大重叠离散小波变换 (MODWT) 对信号进行去噪，并使用 Donoho 和 Johnstone 的通用阈值和级别相关的阈值处理。
'rigrsure' 使用 Stein 的无偏风险原理。
'heursure' 是 Stein 的无偏风险的启发式变体。
'sqtwolog' 使用 Donoho 和 Johnstone 的通用阈值和 DWT。
'minimaxi' 使用极小极大阈值处理。

SORH 使用 's' 或 'h' 指定软或硬阈值处理。

SCAL 定义阈值重新缩放的类型：

'one' 表示不进行重新缩放。
'sln' 表示使用基于第一级系数的噪声估计进行重新缩放。
'mln' 表示使用级别相关的噪声估计进行重新缩放。

'mln' 是 MODWT 去噪支持的唯一选项。

N 是小波变换的级别，WNAME 是作为字符串指定的小波。对于 MODWT 去噪，WNAME 必须对应于正交小波。

XD = wden(C,L,TPTR,SORH,SCAL,N,WNAME) 返回通过对 DWT 系数向量 C 和每个级别 DWT 系数的数量 L 进行操作获得的去噪信号。C 和 L 是 WAVEDEC 的输出。您必须在 WAVEDEC 和 WDEN 中使用相同的小波。

XD = wden(W,'modwtsqtwolog',SORH,'mln',N,WNAME) 返回通过对 MODWT 变换矩阵 W 进行操作获得的去噪信号。W 是 MODWT 的输出。您必须在 MODWT 和 WDEN 中使用相同的小波。

[XD,CXD] = WDEN(...) 返回去噪的小波系数。对于 DWT 去噪，CXD 是 C 的去噪版本。对于 MODWT 去噪，CXD 是 W 的去噪版本。

关于 SORH 参数的代码

在 wden 函数中，SORH 参数指定了软或硬阈值处理，其中 's' 代表软阈值，'h' 代表硬阈值。在函数中没有专门有关 SORH 的代码段，但在函数的内部实现中，根据 SORH 参数的值分别采用软或硬阈值方法进行阈值处理。具体实现细节可以查看函数的源代码。

示例

% 生成一个带噪声的信号
signal = sin(2*pi*0.1*[1:100]) + randn(1,100);

% 使用 'modwtsqtwolog' 阈值规则和软阈值处理对信号进行去噪
denoisedSignal = wden(signal,'modwtsqtwolog','s','mln',4,'db4');

% 绘制原始信号和去噪后的信号
figure;
subplot(2,1,1);
plot(signal);
title('原始信号');
subplot(2,1,2);
plot(denoisedSignal);
title('去噪后的信号');