%單尺度一維離散小波變換（dwt函數）

%裝載一維原始信號
load noissin;  
s=noissin(1:1000);
%畫出原始信號的波形
subplot(411);plot(s);  %函數subplot的作用是在標定位置上建立坐標系
title('原始信號');
%下面用haar小波函數進行一維離散小波變換
[ca1,cd1]=dwt(s,'haar');
subplot(4,2,3);plot(ca1);
ylabel('haar(ca1)');
subplot(4,2,4);plot(cd1);
ylabel('haar(cd1)');
%給定一個小波db2,計算與之相關的分解濾波器
[Lo_D,Hi_D]=wfilters('db2','d');
%用分解濾波器Lo_D,Hi_D計算信號s的離散小波分解系數
[ca2,cd2]=dwt(s,Lo_D,Hi_D);
subplot(4,2,5);plot(ca2);
ylabel('db2(ca2)');
subplot(4,2,6);plot(cd2);
ylabel('db2(cd2)');