clc;
clear;
t=0:0.001:2;
%給定一余弦信號
x=sin（2*pi*t）;
subplot（321）;
plot（x）;
axis（[0?3000?-2?2]）;
title（‘原始信號?x=sin（2\bullet\pi\bullett）‘）;
y=sin（2*pi*t）+0.3*randn（1length（t））;
subplot（322）;
plot（y）;
title（‘染噪信號?y=sin（2\bullet\pi\bullett）+0.3\bulletrandn（1length（t））‘）;
%?使用wpdencmp函數進行信號降噪
n?=?length（y）;?
thr?=?sqrt（2*log（n*log（n）/log（2）））;??%計算閾值?
xwpd?=?wpdencmp（y‘s‘4‘sym5‘‘sure‘thr1）;
%?畫出