wavelet9.m

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

data=textread('data.txt');
figure(1);

f=fft(data,2000);
f=abs(f);
subplot(2,1,1);
plot(data);
subplot(2,1,2);
n=10:2000;
plot(n,f(n));
[c,l]=wavedec(data,6,'db6');
figure(2);
subplot(8,1,1);
plot(data);
ylabel('data');
a=wrcoef('a',c,l,'db6',6);
subplot(8,1,2);
plot(a);
ylabel('a6');
for i=1:6,
    d=wrcoef('d',c,l,'db6',7-i);
    subplot(8,1,i+2);
    plot(d);
    ylabel(['d',num2str(7-i)]);
end