一维小波去噪重构法，模极大值，软硬阈值

load leleccum;
s=leleccum（1:4000）;
lev=5;
[cl]=wavedec（slev‘db4‘）;
sigma=wnoisest（cl1）;
alpha=2;
thr=wbmpen（clsigmaalpha）;
keepapp1=1;
xd=wdencmp（‘gbl‘cl‘db4‘levthr‘s‘keepapp1）;
figure（1）;
subplot（231）plot（s）title（‘Original signal‘）;
subplot（232）plot（xd）title（‘模极大值去噪图‘）;
ca5=appcoef（cl‘db4‘5）
cd5=detcoef（cl5）;
cd4=detcoef（cl4）;
cd3=detcoef（cl3）;
cd2=detcoef（cl2）;
cd1=detcoef（cl1）;
cdd5=zeros（1length（cd5））;
cdd4=zeros（1length（cd4））;
cdd3=zeros（1length（cd3））;
cdd2=zeros（1length（cd2））;
cdd1=zeros（1length（cd1））;
c1=[ca5 cdd5 cdd4 cdd3 cdd2 cdd1];
s1=waverec（c1