xiaobobao 自己的小波包程序（含消噪、能量特征提取、信号重构等）

自己的小波包程序（含消噪、能量特征提取、信号重构等），希望对大家有帮助

%%小波包程序：能量-特征向量算法一
tic
clear
clc
a=15; b=0.5; f=200; m=pi/3; %设定初始值
N=400; flag=1;
%下面代码是分段绘制信号，时间间隔为0.001，数据点数是1001个
t1=linspace（0.3010.5200）;
xt1=9*sin（100*pi*t1）+5*sin（180*pi*t1）+sin（240*pi*t1）;
t2=linspace（0.5010.7200）;
xt2=9*sin（100*pi*t2）+5*sin（180*pi*t2）+sin（240*pi*t2）+sin（2*pi*f*t2+m.*exp（-a*（t2-b）））;
t=[t1t2];
xt=[xt1xt2];
figure（1）;
plot（txt）;
title（‘原始信号图‘）; xlabel（‘时间（t）‘）; ylabel（‘幅值‘）; 

T=wpdec（xt3‘db2‘）; %进行小波包分解
plot（T）;  figure（2）;
%计算小波包分解系数
x1=wpcoef（T[10]）; 
x2=wpcoef（T[11]）; 
z1=wpcoef（T[20]）;  
z2=wpcoef（T[21]）; 
z3=wpcoef（T[22]）; 
z4=wpcoef（T[23]）; 
y1=wpcoef（T[30]）;  
y2=wpcoef（T[31]）; 
y3=wpcoef（T[32]）; 
y4=wpcoef（T[33]）; 
y5=wpcoef（T[34]）; 
y6=wpcoef（T[35]）; 
y7=wpcoef（T[36]）; 
y8=wpcoef（T[37]）

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        4422  2020-07-27 12:30  ╨┬╜¿╬─╝■╝╨\xiaobobao.m
     文件        1249  2020-07-27 12:30  ╨┬╜¿╬─╝■╝╨\xiaobobao2.m
     文件        4719  2020-07-27 12:30  ╨┬╜¿╬─╝■╝╨\xiaobobao3.m
     文件         547  2020-07-27 12:30  ╨┬╜¿╬─╝■╝╨\xiaobobao4.m