小波包降噪和提升小波降噪的对比

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