写个小波变换去噪matlab代码

以下是一个简单的小波去噪的MATLAB代码示例：

% 读取原始信号
load noisysignal.mat;
x = noisysignal;

% 小波去噪
wname = 'sym4';     % 选择小波类型
level = 5;          % 设置小波分解层数
thr = 0.2;          % 设置阈值
[C, L] = wavedec(x, level, wname);
cD = detcoef(C, L, level);
thrC = wthresh(cD, 'h', thr);
C([1:L(level)+1]) = wthresh(C([1:L(level)+1]), 'h', thr);
for k = level:-1:1
    idx = (L(k)+1):(L(k+1)-1);
    C(idx) = wthresh(C(idx), 's', thrC);
end
xrec = waverec(C, L, wname);

% 绘制结果
subplot(2,1,1); plot(x);
title('Original Signal');
subplot(2,1,2); plot(xrec);
title('Denoised Signal');

在这个示例中，我们首先从文件 `noisysignal.mat` 中读取了一个被噪声污染的原始信号。然后，我们使用 `wavedec` 函数对原始信号进行小波分解，将其分解成多个小波系数。接着，我们使用 `wthresh` 函数来设置阈值，并将小于阈值的系数设置为0。然后，我们使用 `waverec` 函数将处理后的小波系数合并成一个信号，得到去噪后的信号。最后，我们使用 `subplot` 函数将原始信号和去噪后的信号绘制在同一张图上进行对比。