init=2055615866;
[xrefx]=wnoise（5117init）;
xref=xref（1:1000）;
x=x（1:1000）
subplot（221）plot（xref）
title（‘原始信號‘）
subplot（222）plot（x）
title（‘噪聲信號‘）
n=length（x）;
thr=sqrt（2*log（n*