myxiaobo.m

来自「进行时频变换的工具箱」· M 代码 · 共 17 行

subplot(311)
plot(x)
ylabel('原始信号')
[c,l]=wavedec(x,3,'db2');%用DB2小波对信号进行3 层分解
ca3=appcoef(c,l,'db2',3);%提取小波分解的低频系数
cd3=detcoef(c,l,3);%提取第三层的高频系数
cd2=detcoef(c,l,2);
cd1=detcoef(c,l,1);
subplot(3,2,3);plot(ca3);
ylabel('ca3');
subplot(3,2,4);plot(cd3);
ylabel('cd3');
subplot(3,2,5);plot(cd2);
ylabel('cd2')
subplot(3,2,6);plot(cd1);
ylabel('cd1');