小波去噪matlab代码

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

function [denoised_signal] = wavelet_denoising(signal, wavelet_type, level, noise_std)

    % 小波去噪函数
    % signal: 待去噪信号
    % wavelet_type: 使用的小波类型
    % level: 分解层数
    % noise_std: 噪声标准差

    % 小波分解
    [C, L] = wavedec(signal, level, wavelet_type);

    % 计算噪声标准差的阈值
    sigma = noise_std * sqrt(2 * log(length(signal)));
    % 每层小波系数的阈值
    threshold = sigma * sqrt(2 * log(length(signal) ./ (2.^(0:level-1))));

    % 阈值处理小波系数
    C_den = wthresh(C, 's', threshold);

    % 小波重构
    denoised_signal = waverec(C_den, L, wavelet_type);

end

使用方法：

% 生成带噪声的信号
signal = randn(1, 1000) + 1;
noise_std = 0.5;
noisy_signal = signal + noise_std * randn(1, 1000);

% 小波去噪
wavelet_type = 'db4'; % 小波类型
level = 5; % 分解层数
denoised_signal = wavelet_denoising(noisy_signal, wavelet_type, level, noise_std);

% 画图比较
figure()
subplot(3,1,1)
plot(signal);
title('原始信号');
subplot(3,1,2)
plot(noisy_signal);
title('带噪声信号');
subplot(3,1,3)
plot(denoised_signal);
title('去噪后信号');

这里使用的是db4小波，可以根据需要进行更改。阈值的计算使用的是固定阈值，也可以尝试其他的计算方法。