帮我写一个小波阈值降噪的MATLAB程序

下面是一个基本的小波阈值降噪的MATLAB程序：

% 读取原始信号
x = load('signal.mat');

% 使用 Daubechies 4 小波进行分解
[c, l] = wavedec(x, 4, 'db4');

% 设置阈值
thr = wthrmngr('sqtwolog', c);

% 应用阈值
s = wdencmp('gbl', c, l, 'db4', 4, thr, 'h');

% 绘制原始信号和降噪后的信号
subplot(2, 1, 1);
plot(x);
title('Original Signal');

subplot(2, 1, 2);
plot(s);
title('Denoised Signal');

在这个程序中，我们首先读取了一个名为 `signal.mat` 的原始信号。然后使用 Daubechies 4 小波进行了分解，并使用 `wthrmngr` 函数计算了阈值。最后，使用 `wdencmp` 函数对信号进行了降噪。最后，我们绘制出原始信号和降噪后的信号的图形。