例程10-10.m

来自「《小波分析理论与Matlab7实现》书中例程」· M 代码 · 共 15 行

% 采用补零的扩展模式（参见dwtmode函数）
% 装载一维尺度信号
load leleccum; 
s = leleccum(1:3920); 
% 使用db1在第3层进行分解
[c,l] = wavedec(s,3,'db1');
subplot(4,1,1);plot(s);
title('原始信号');
% 从小波分解结构[c,l]中提取1、2及3层的细节系数
[cd1,cd2,cd3] = detcoef(c,l,[1 2 3]);
% 绘图命令
subplot(4,1,2);plot(cd3);Ylabel('cd3');
subplot(4,1,3);plot(cd2);Ylabel('cd2');
subplot(4,1,4);plot(cd1);Ylabel('cd1');