%裝載采集的信號leleccum.mat
load leleccum;
%將信號中第1160到第1235個采樣點賦給s
index=1160:1235;
s=leleccum(index);
%畫出原始信號
figure(1);
plot(index,s);
xlabel('樣本序號 n');
ylabel('幅值 A');

%用db3小波進行5層分解
[c,l]=wavedec(s,5,'db3');
%重構第1-5層逼近信號
a5 = wrcoef('a',c,l,'db3',5);
a4 = wrcoef('a',c,l,'db3',4);
a3 = wrcoef('a',c,l,'db3',3);
a2 = wrcoef('a',c,l,'db3',2);
a1 = wrcoef('a',c,l,'db3',1);
%重構第1-5層細節信號
d5 = wrcoef('d',c,l,'db3',5);
d4 = wrcoef('d',c,l,'db3',4);
d3 = wrcoef('d',c,l,'db3',3);
d2 = wrcoef('d',c,l,'db3',2);
d1 = wrcoef('d',c,l,'db3',1);

%顯示細節信號
figure(2)
subplot(3,1,1);
plot(index,d3,'LineWidth',2);
ylabel('d3');
subplot(3,1,2);
plot(index,d2,'LineWidth',2);
ylabel('d2');
subplot(3,1,3);
plot(index,d1,'LineWidth',2);
ylabel('d1');
xlabel('樣本序號 n');
%消除奇異點
%%設置1-3層細節信號為零
s0=a5+d5+d4;
%畫出重構信號
figure(3);
plot(index,s0);
xlabel('樣本序號 n');
ylabel('幅值 A');