wavelet2.m

来自「对信号进行小波变换」· M 代码 · 共 12 行

load sumsin;
s=sumsin;
swc=swt(s,3,'db1');%用db1小波对信号做静态小波分解，得到的系数存到swc中
[swa,swd]=swt(s,3,'db1');
a0=iswt(swc,'db1');%用swc中的系数重建信号，存入a0中
a1=iswt(swa,swd,'db1');
err0=norm(s-a0,'fro')
err1=norm(s-a1,'fro')
subplot(311);plot(s);
subplot(312);plot(a0);
subplot(313);plot(a1);