%%小波包程序：能量-特征向量算法一
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