|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
为了展示HHT波内调制的原理,我需要把原始波形和瞬时频率同时画在图上,并且标注单位。当频率成分只有一个时,同时画出两个波形plotyy可以做到并标注,没问题。可是当HHT分解出多个频率分量时,怎么把一个原始波形和多个频率画在一张图上并分别标注呢?我试过Hold on然后循环plotyy,可是由于前后的单位不一致画出来的图像频率标注根本不对,而且后一个会覆盖前一个的单位标注。下面是我的代码,高手请看一看。
[AX,H1,H2]=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;也是不对的,同样会前后单位不一致并发生覆盖。 |
|