clear?all;
close?all;
clc;

tic;
N=1000;
t=1:1000;
f=sin（0.03*t）;
load?noissin;
e1=noissin;
init=2055615866;
randn（‘seed‘init）;
e=e1+0.5*randn（size（e1））;
subplot（221）;
plot（tf）;
xlabel（‘1?樣本序列‘）;
ylabel（‘原始信號幅值‘）;
grid
subplot（222）;
plot（e）
xlabel（‘2?測試樣本序列‘）
ylabel（‘含有已加噪聲的信號幅值‘）
grid
s1=wden（e‘minimaxi‘‘s‘‘one‘5‘db12‘）;
subplot（223）;
plot（s1）;
xlabel（‘3?db10降噪后信號‘）
ylabel（‘db10小波降噪后的信號幅值‘）;
grid;
s2=wden（e‘heursure‘‘s‘