matlab 一维小波软阈值降噪程序

以下是一个简单的 MATLAB 一维小波软阈值降噪程序的示例：

function x_denoised = wavelet_soft_thresholding(x, wavelet, level, threshold)
% x: 待降噪信号
% wavelet: 选用的小波基
% level: 分解层数
% threshold: 阈值

% 小波分解
[c, l] = wavedec(x, level, wavelet);

% 软阈值处理
for i = 1:level
    start_index = sum(l(1:i-1)) + 1; % 当前层系数起始下标
    end_index = sum(l(1:i)); % 当前层系数结束下标
    c(start_index:end_index) = sign(c(start_index:end_index)) .* max(abs(c(start_index:end_index)) - threshold, 0);
end

% 小波重构
x_denoised = waverec(c, l, wavelet);
end

使用示例：

% 生成含噪声的信号
x = sin(linspace(0, 10*pi, 1024)) + 0.5*randn(1, 1024);

% 调用小波软阈值降噪函数
x_denoised = wavelet_soft_thresholding(x, 'db4', 4, 0.5);

% 绘制原始信号和降噪后的信号
subplot(211), plot(x), title('Original signal');
subplot(212), plot(x_denoised), title('Denoised signal');

你可以根据自己的需求修改阈值大小和选用的小波基。