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