wavelet3.m

来自「小波在信号处理中的应用」· M 代码 · 共 22 行

data=textread('data.txt');
figure(1);
plot(data);
[C,L]=wavedec(data,2,'db1');
cA1=appcoef(C,L,'db1',1);
cA2=appcoef(C,L,'db1',2);
figure(2);
subplot(2,1,1);
plot(cA1);
title('1');
subplot(2,1,2);
plot(cA2);
title('2');
cD1=detcoef(C,L,1);
cD2=detcoef(C,L,2);
figure(3);
subplot(2,1,1);
plot(cD1);
title('3');
subplot(2,1,2);
plot(cD2);
title('4');