wavelet_extract.m

来自「小波变换的条件限制的实现的编程」· M 代码 · 共 24 行

%
%小波对信号特征的提取
%
%
load nearbrk;
s=nearbrk;
%一维原始信号
subplot(4,2,1);
plot(s);
ylabel('s');
title('原始信号');
[c,l]=wavedec(s,4,'haar');%利用harr小波信号进行3层分解
for i=1:4
    decmpa=wrcoef('a',c,l,'haar',5-i);
    subplot(5,2,2*i+1);
    plot(decmpa);
    ylabel(['a',num2str(5-i)]);
end
for i=1:4
    decmpd=wrcoef('d',c,l,'haar',5-i);
    subplot(5,2,2*(i+1));
    plot(decmpd);
    ylabel(['d',num2str(5-i)]);
end