用matlab写一个小波包去噪的代码

好的，下面是一个基于小波包的去噪Matlab代码，其中使用了db4小波基函数进行分解和重构：

function denoised_signal = wavelet_packet_denoise(noisy_signal, level, threshold)

% 小波包去噪函数
% 输入：
%   noisy_signal - 带噪信号
%   level - 分解的层数
%   threshold - 阈值
% 输出：
%   denoised_signal - 去噪后的信号

% 小波包分解
[tree, ~] = wpdec(noisy_signal, level, 'db4');

% 统计每个小波包系数的能量
energy = wpcoef(tree, 'e');

% 计算每个小波包的能量平均值
avg_energy = mean(energy);

% 计算阈值
t = threshold * avg_energy;

% 执行阈值处理
for i = 1:length(tree)
    energy = wpcoef(tree, i, 'e');
    coef = wpcoef(tree, i);
    mask = energy < t;
    coef(mask) = 0;
    tree = wpcoef(tree, i, coef);
end

% 小波包重构
denoised_signal = wprec(tree, 'db4');

end

使用方法：

% 读取带噪信号
noisy_signal = load('noisy_signal.mat');

% 小波包去噪
denoised_signal = wavelet_packet_denoise(noisy_signal, 5, 0.1);

% 显示去噪后的信号
plot(denoised_signal);