小波包分解

小波包分解包括重构信号问题的matlab程序，测试可用，采用自己的数据。

y1=awgn（v_ac70）;   %%注意使用measured参
figure（1）
plot（ty1）
y2=awgn（v_bc70）;   %%注意使用measured参
figure（2）
plot（ty2）


depth=4;%分解层数
[T11]=wpdec（y1depth‘db8‘）;%小波包分解
[thr1sorh1keepapp1crit1]=ddencmp（‘den‘‘wp‘y1）;    %%基于小波包分解获取阈值的值
sorh1=‘s‘;
%thr1 = thselect（y1‘minimaxi‘）

T1 = wpthcoef（T11keepapp1sorh1thr1）
%plot（T1）
x1=wprcoef（T1[4 0]）;%重构最后一层小波包系数
x2=wprcoef（T1[4 1]）;
x3=wprcoef（T1[4 2]）;
x4=wprcoef（T1[4 3]）;
x5=wprcoef（T1[4 6]）;
x6=wprcoef（T1[4 7]）;
x7=wprcoef（T1[4 12]）;
x8=wprcoef（T1[4 13]）;
yy1=x1+x2+x3+x4;
%%%%%%%%%%%%%%%%%%%%%%%%%%
[T22]=wpdec（y2depth‘db8‘）;%小波包分解
[thr2sorh2keepap2crit2]=ddencmp（‘den‘‘wp‘y2）;    %%基于小波包分解获取阈值的值
sorh2=‘s‘;
%thr2 = thselect（y2‘minimaxi‘）
T2 = wpthcoef（T22keepapp1sorh2thr1）
g1=wprcoef（T2[4 0]）;%重构最后一层小波包系数
g2=wprcoef（T2[4 1]）;
g3=wprcoef（T2[4 3]）;
g4=wprcoef（T2[4 4]）;
g5=wprcoef（T2[4 6]）;
g6=wprcoef（T2[4 7]）;
g7=wprcoef（T2[4 8]）;
yy2=g1+g2+g3+g4;
%E（i）=wentropy（x‘shannon‘i）;%求小波包能量熵 x1+x2+x3+x4+x5+x6+x7+
%y=x1+x2+x3+x4+x5+x6+x7+x8+x9+x10+x11+x12+x13+x14+x15+x16;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%  yy2=g1+g3+g2+g4+g5+g6;
figure（5）
plot（yy1）
figure（9）
plot（yy2）
%%%%%%%%%%%%%%%%%%%%相关
a=length（yy1）-1;             %%进行直接互相关
%[yya]=xcorr（xwpd2xwpd1）;  %%a为位置向量，y为幅值向量
[Rxya]=xcorr（yy1yy2）;           %%互相关
Rxx0=sum（yy2.*yy2）;           %%自相关
Ryy0=sum（yy1.*yy1）;
rm=Rxy/sqrt（Rxx0*Ryy0）;   %%求出相关系数   两个图形Rxy与rm很相似的

dt=1.0*10^-8;
%%e=find（y==c）;
u=dt*a;
figure
plot（urm）
[YI]=max（rm）; 
text（u（I）Ynum2str（Y））   %%%显示最大幅值的横坐标u（I）最大幅值Y

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1721  2010-03-27 17:42  小波包分解包括重构信号问题程序.m

----------- ---------  ---------- -----  ----

                 1721                    1