program_10_09.m

来自「小波分析理论与MATLABR2007实现一书的原代码」· M 代码 · 共 35 行

% 当前扩展模式是补零（参见dwtmode函数）
% 低频信号由1～6层系数获得
cfs = [1];  
essup = 10; 
figure(1) 
for i=1:6 
    rec = upcoef('a',cfs,'db6',i);
    % essup 是重构信号必须的
    % 当j等于essup时，rec(j) 非常小
    ax = subplot(6,1,i),h = plot(rec(1:essup)); 
    set(ax,'xlim',[1 325]); 
    essup = essup*2; 
end 
subplot(611) 
title(['尺度1到6，由惟一的系数获得的低频信号'])

% 同样可以获得高频信号
% 高频信号可以由惟一系数从尺度1～6获得
cfs = [1]; 
mi = 12; ma = 30;   % db6小波滤波器必须的
rec = upcoef('d',cfs,'db6',1); 
figure(2) 
subplot(611), plot(rec(3:12)) 
for i=2:6 
    rec = upcoef('d',cfs,'db6',i); 
    subplot(6,1,i), plot(rec(mi*2^(i-2):ma*2^(i-2))) 
end 
subplot(611) 
title(['由惟一系数从尺度1到6获得的高频信号'])