小波阈值去噪的matlab代码

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

% 读取图像
img = imread('lena.png');
% 转换为灰度图像
img = rgb2gray(img);
% 加入高斯噪声
img_noise = imnoise(img, 'gaussian', 0, 0.01);

% 小波变换
[c, l] = wavedec2(img_noise, 2, 'db4');
% 计算噪声标准差
sigma = median(abs(c))/0.6745;
% 计算阈值
threshold = sigma * sqrt(2*log(numel(img_noise)));
% 求解软阈值
c_soft = wthresh(c, 's', threshold);
% 重构图像
img_denoised = waverec2(c_soft, l, 'db4');

% 显示结果
subplot(1, 2, 1);
imshow(img_noise);
title('原图像加噪声');
subplot(1, 2, 2);
imshow(img_denoised, []);
title('小波阈值去噪');

其中，`imread` 函数用于读取图像，`rgb2gray` 函数用于将图像转换为灰度图像，`imnoise` 函数用于给图像加入高斯噪声。`wavedec2` 函数用于进行二维小波变换，`median` 函数用于计算中位数，`wthresh` 函数用于求解软阈值，`waverec2` 函数用于进行二维小波逆变换。

你可以根据自己的需要调整参数和阈值，对不同类型的图像进行去噪处理。