求提取小波系數的matlab程序

裝載leleccum信號

load

leleccum;

s

=

leleccum(1:3920);

%

用db1小波函數對信號進行三尺度小波分解

[C,L]=wavedec(s,2,'db1');

figure(1);

plot(s);

title('leleccum原始信號');

%

提取尺度1的低頻系數

cA1

=

appcoef(C,L,'db1',1);

%用小波分解框架[C.L]計算1層低頻系數的近似值，小波基為db1

%

提取尺度2的低頻系數

cA2

=

appcoef(C,L,'db1',2);

figure(2);

subplot(2,1,1);

plot(cA1);

title('尺度1的低頻系數');

subplot(2,1,2);

plot(cA2);

title('尺度2的低頻系數');

%

提取尺度1的高頻系數

cD1

=

detcoef(C,L,1);

%用小波分解框架[C.L]計算1層高頻系數的近似值，小波基為db1

%

提取尺度2的高頻系數

cD2

=

detcoef(C,L,2);

figure(3);

subplot(2,1,1);

plot(cD1);

title('尺度1的高頻系數');

subplot(2,1,2);

plot(cD2);

title('尺度2的高頻系數');