%裝載采集的信號
load?leleccum;
%將信號中第2000到第3450個采樣點賦給s
indx=2000:3450;
s=leleccum（indx）;
%畫出原始信號
subplot（221）;
plot（s）;
title（‘yuanxinhao‘）;
ylabel（‘幅值A‘）;
%用db1小波對信號進行3層分解并提取系數
[cl]=wavedec（s3‘db1‘）;
a3=appcoef（cl‘db1‘3）;
d3=detcoef（cl3）;
d2=detcoef（cl2）;
d1=detcoef（cl1）;
%強制消噪處理
dd3=zeros（1length（d3））;
dd2=zeros（1length（d2））;
dd1=zeros（1length（d1））;
c1=[a3?dd3?dd2?dd1];
s1=waverec（c1l‘db1‘）;
subplot（222）