小波分解matlab代码

以下是一个使用 Matlab 实现小波分解的示例代码：

% 定义信号
x = linspace(0, 1, 2^8);
y = sin(2*pi*10*x) + cos(2*pi*20*x) + randn(size(x));

% 小波分解
w = 'db4'; % 选用 Daubechies 4 小波
level = 5; % 分解层数
[C, L] = wavedec(y, level, w); % 小波分解

% 重构信号
A = wrcoef('a', C, L, w, level); % 从近似系数重构信号
D = cell(1, level); % 存储细节系数
for i = 1:level
    D{i} = wrcoef('d', C, L, w, i); % 从各层细节系数重构信号
end

% 显示结果
subplot(level+1, 1, 1);
plot(x, y);
title('Original Signal');
for i = 1:level
    subplot(level+1, 1, i+1);
    plot(x, D{i});
    title(['Level ' num2str(i) ' Detail Coefficients']);
end

代码中使用了 Matlab 自带的 `wavedec` 函数进行小波分解，其中需要指定小波类型和分解层数。然后使用 `wrcoef` 函数从分解后的系数中重构信号。最后将原始信号和各层细节系数绘制在不同的子图中进行展示。