|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
小波包的分解和重构 与 小波包的消噪 有什么 联系与区别?
clc;
clear;
fs=12000;
y=load('E:\滚动轴承\风端轴承\0fd normal.txt');
Y=y(:,1);
x=Y(1:1200);
x=x';
t=(0.1/fs:1/fs:0.1);
figure(1);
subplot(411);
plot(t,x);
title('正常信号');
n=length(x);
%3层小波包分解;
T=wpdec(x,3,'db4');
plot(T);%显示小波包3层分解树结构
bt = besttree(T);
plot(bt);%最优小波包树
w我想对最优小波包树的[2,0],[3,2],[3,3]这六个系数进行 消噪 该怎么写?以下是我自己编 的 看看哪错了?大侠们
cfs1=wpcoef(T,[2,0]);
cfs2=wpcoef(T,[3,2]);
cfs3=wpcoef(T,[3,3]);
cfs4=wpcoef(T,[3,4]);
cfs5=wpcoef(T,[3,5]);
cfs6=wpcoef(T,[2,3]);
thr=1.169;
xd1=wpdencmp(cfs1,'s',3,'db4','sure',thr,1);
xd2=wpdencmp(cfs2,'s',3,'db4','sure',thr,1);
xd3=wpdencmp(cfs3,'s',3,'db4','sure',thr,1);
xd4=wpdencmp(cfs4,'s',3,'db4','sure',thr,1);
xd5=wpdencmp(cfs5,'s',3,'db4','sure',thr,1);
xd6=wpdencmp(cfs6,'s',3,'db4','sure',thr,1);
xd=wprec(bt);%为新的分解结构
subplot(412);
plot(t,xd);
xt 信号和x信号怎么一样啊,没起消噪作用?帮看看
|
-
x 和xt
-
|