wavlet.m

对信号进行小波变换

load noissin;%读入白噪声
s=noissin(1:1000);%去信号前1000个采样点
[c,l]=wavedec(s,3,'db4');%对信号做层数为3的多尺度分解
[cd1,cd2,cd3]=detcoef(c,l,[1,2,3]);%得到三个尺度的细节数
ca3=appcoef(c,l,'db4',3);%得到尺度3的近似数
subplot(511);plot(1:1000,s);title('s');
subplot(512);plot(1:l(1),ca3);title('ca3');
subplot(513);plot(1:l(2),cd3);title('cd3');%l(i)存放的是相应参数的长度
subplot(514);plot(1:l(3),cd2);title('cd2');
subplot(515);plot(1:l(4),cd1);title('cd1');
%将原始信号和分解后得到的以一组近似系数和3组细节数的波形显示出来