clear
clc
fs=1000;%采樣頻率
N=1000;%數據點數
n=0:N-1;
t=n/fs;%時間序列
%構造仿真信號
x=[sin（40*pi*t）+cos（10*pi*t）]+0.8*cos（15*pi*t+pi/3）;
s=x+10*randn（1length（x））
figure;
?subplot（221）;
?plot（tx）;
?xlabel（‘樣本序號?n‘）;
?ylabel（‘幅值?A‘）;
?title（‘z原始信號‘）;
?????????????%采用默認閾值，用wdencmp函數進行消噪處理
??????[thrsorhkeepappcrit]=ddencmp（‘den‘‘wp‘s）;
??????
????%用全局閾值選項進行消噪處理
?????[ctreedperf0perfl2]=wpdencmp（ssorh3‘db2‘critthrkeepapp）;
?????subplot（222）;
?????plot（c）;
?????xlabel（‘樣本序號?n‘）;
?????ylabel（‘幅值?A‘）;
?????title（‘默認閾值消噪信號‘）;

????%根據前面的消噪效果，調節閾值大小進行消噪
?????thr=15;
?????[c1treedperf0perfl2]=wpdencmp（ssorh3‘db2‘critthrkeepapp）;
?????subplot（223）;
?????plot（c1）;
?????xlabel（‘樣本序號?n‘）;
?????ylabel（‘幅值?A‘）;
?????title（‘調節閾值后的消噪信號‘）;
%?得到Haar小波的提升方案
lshaar