clear?all
clc
load?leleccum;%裝在采集的信號leleccum.mat
s=leleccum（1:3920）;%將信號中第1到第3920個采樣點賦給s
ls=length（s）;

%畫出原始信號
subplot（221）;
plot（s）;
title（‘原始信號‘）;grid;

%用db1小波對原始信號進行3層分解并提取系數(shù)
[cL]=wavedec（s3‘db1‘）;
ca3=appcoef（cL‘db1‘3）;
cd3=detcoef（cL3）;
cd2=detcoef（cL2）;
cd1=detcoef（cL1）;

%對信號進行強制性去噪處理并圖示結(jié)果
cdd3=zeros（1length（cd3））;
cdd2=zeros（1length（cd2））;
cdd1=zeros（1length（cd1））;
c1=[ca3?cdd3?cdd2?cdd1];
s1=waverec（c1L‘db1‘）;
subplot（222）;