playfish 发表于 2007-5-7 21:54

多分量plotyy画图问题

为了展示HHT波内调制的原理,我需要把原始波形和瞬时频率同时画在图上,并且标注单位。当频率成分只有一个时,同时画出两个波形plotyy可以做到并标注,没问题。可是当HHT分解出多个频率分量时,怎么把一个原始波形和多个频率画在一张图上并分别标注呢?我试过Hold on然后循环plotyy,可是由于前后的单位不一致画出来的图像频率标注根本不对,而且后一个会覆盖前一个的单位标注。下面是我的代码,高手请看一看。
=plotyy(t,X,ft,fq(1,:));hold on;
   set(H1,'LineStyle','--');
   title('HHT spectrum');
   set(get(AX(1),'Ylabel'),'String','diaphragm displacement(m)');
   set(get(AX(2),'Ylabel'),'String','instantaneous frequency(Hz)');
   xlabel('time(s)');
for i=2:n
    plotyy(ft,fq(i-1,:),ft,fq(i,:));hold on;
end
grid;

其中X是原始波形,fq为瞬时频率。如果把循环里的画图改为plotyy(t,X,ft,fq(i,:));hold on;也是不对的,同样会前后单位不一致并发生覆盖。
页: [1]
查看完整版本: 多分量plotyy画图问题