|
楼主 |
发表于 2008-6-12 12:02
|
显示全部楼层
哦,我已经解决了,程序如下:
figure
subplot(211);
plot(t,x);
xlabel('t/s')
ylabel('A/mm')
title('原始信号','FontSize',16);
wpt=wpdec(x,2,'db2','shannon');
plot(wpt);
c20=wprcoef(wpt,[2,0]);
subplot(212);
plot(t,c20);
xlabel('t/s')
ylabel('A/mm')
title('重构后小波包节点(2.0)','FontSize',16);
%小波包能量监测
n=2;
L=wpdec(x,n,'db2','shannon');
E=wenergy(L);
xi=1:4;
E=[E(xi)];
SE=E(1)+E(2)+E(3)+E(4);
figure
subplot(211)
w=0.5:1:3.5
bar(w,E/SE,1)
set(gca,'xtick',[1 2 3 4])
set(gca,'xticklabel',{'64','128','192','256'})
xlabel('频带(HZ)')
ylabel('占总能量百分比(%)');
title('原始信号小波包频带能量监测','FontSize',16);
%小波包2层分解图 |
评分
-
1
查看全部评分
-
|