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