一维信号小波降噪matlab程序

%装载采集的信号
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）