小波分析1

%含突变点信号freqbrk的连续小波变化mesh图
>> load freqbrk
>> s=freqbrk;
>> coefs=cwt(s,[1;30],'morl');   %coefs为小波系数，s为待分析信号，[1;30]为小波变换尺度，'morl'为选取的小波

>> colormap([0 0 0]) >> mesh(abs(coefs))

突变点信号freqbrk的离散小波变换

>> load freqbrk
>> s=freqbrk;
>> ls=length(s);
>> [C,L]=wavedec(s,4,'db5');
>> A4=wrcoef('a',C,L,'db5',4);
>> D4=wrcoef('d',C,L,'db5',1);
>> subplot(3,1,1);plot(abs(s),'k');
>> subplot(3,1,2);plot(abs(A4),'k');
>> subplot(3,1,3);plot(abs(D4),'k');
>> 

在对信号进行离散变换时，常常先将其利用函数wavedec进行多分辨率分析，进而对结果进行重构，其实质是提取感兴趣的特征。

[C,L]=wavedec(s,4,'db5')，其中C为多分辨率分析结果，L为多分辨率分析结果中矩阵向量对应长度，s为被分析的信号，4为多分辨率分析的chidu ,'db5'为选取的小波。