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