wavelet10.m

来自「该例子为大家提供了一个很好的」· M 代码 · 共 16 行

data=textread('data.txt');
figure(1);
[c,l]=wavedec(data,2,'db4');
subplot(4,1,1);
plot(data);
ylabel('data');
a=wrcoef('a',c,l,'db4',2);
subplot(4,1,2);
plot(a);
ylabel('a2');
for i=1:2,
    d=wrcoef('d',c,l,'db6',3-i);
    subplot(4,1,i+2);
    plot(d);
    ylabel(['d',num2str(3-i)]);
end