求助：小波包能量谱图的修改

yuleiwudi

大家注意看这两个图的横坐标，我希望我用matlab画出的图能跟第一个一样，但是我在matlab中直接使用bar命令画出的图就只是第二张图，请大家帮忙，需要怎么修改才能画成第一张图的样子.

[ 本帖最后由 yuleiwudi 于 2008-6-10 13:46 编辑 ]

sjdwjt

不明白你干嘛要画成第一个那样，第二个图很好啊，:@L

ch_j1985

原帖由 yuleiwudi 于 2008-6-10 13:08 发表
大家注意看这两个图的横坐标，我希望我用matlab画出的图能跟第一个一样，但是我在matlab中直接使用bar命令画出的图就只是第二张图，请大家帮忙，需要怎么修改才能画成第一张图的样子.

可以参考此帖的方法：
http://www.chinavib.com/forum/viewthread.php?tid=23755&extra=page%3D2%26amp%3Bfilter%3Dtype%26amp%3Btypeid%3D23

yuleiwudi

哦，我第一个图画的时候没注意，我不是想要具有箭头的坐标轴，而只是想要跟第二张图形式一样但坐标跟第一个一致 的图

ch_j1985

原帖由 yuleiwudi 于 2008-6-11 09:36 发表
哦，我第一个图画的时候没注意，我不是想要具有箭头的坐标轴，而只是想要跟第二张图形式一样但坐标跟第一个一致 的图

这样行不行：先求出百分比，再用plot命令画

yuleiwudi

哦，我已经解决了，程序如下：
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层分解图

wuyunfei001

:handshake