purefrequencies_detection.m

来自「MATLAB 小波分析高级技术——书中的一些应用程序」· M 代码 · 共 29 行

load sumsin;
x = sumsin;
figure(1)
subplot(611);
plot(x);
ylabel('x');
title('原始信号及各层近似信号');
% 使用DB3小波进行5层分解
[c,l]=wavedec(x,5,'db3');
for i=1:5
    % 对分解的第5层到第1层的低频系数分别进行重构
    a=wrcoef('a',c,l,'db3',6-i);
    subplot(6,1,i+1);
    plot(a);
    ylabel(['a',num2str(6-i)]);
end
figure(2);
subplot(611);
plot(x);
ylabel('x');
title('原始信号及各层细节信号');
for i=1:5
    % 对分解的第5层到第1层的高频系数分别进行重构
    d=wrcoef('d',c,l,'db3',6-i);
    subplot(6,1,i+1);
    plot(d);
    ylabel(['d',num2str(6-i)]);
end