MATLAB畫雙縱軸曲線。
Matlab plotyy畫雙縱座標圖例項 x = 0:0.01:20; set(AX(1),'XColor','k','YColor','b'); HH1=get(AX(1),'Ylabel'); HH2=get(AX(2),'Ylabel'); set(H1,'LineStyle','-'); legend([H1,H2],{'y1 = 200*exp(-0.05*x).*sin(x)';'y2 = 0.8*exp(-0.5*x).*sin(10*x)'}); Q:右邊用藍色圈起來的tick能去掉嗎?由於用plotyy畫圖,為了使圖儘量地顯示出來,用了set(AX(1),'YLimMode','auto'),但這樣可能會導致左邊AX(1)和右邊AX(2)的tick的間距不一樣,影響美觀。或者說能不能使plotyy畫出的圖兩邊的tick間距是一樣的,這樣在圖形右邊的tick就會重合在一起. A:如果只是想讓plotyy的圖美一些,可以使用其如下形式的呼叫方式: A:也可以用line語句來畫,就沒有左邊和上邊的線了。 Q:plotyy(X1,Y1,X2,Y2,FUN1,FUN2),FUN1和FUN2應該怎麼寫? A:這兩個FUN代表plotyy不一定要用兩個plot,比如下面的例子,一條曲線用plot,一條用semilogy |
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
MATLAB畫雙縱座標
具有兩個縱座標標度的圖形
在MATLAB中,如果需要繪製出具有不同縱座標標度的兩個圖形,可以使用plotyy繪圖函式。呼叫格式為:
plotyy(x1,y1,x2,y2)
其中x1,y1對應一條曲線,x2,y2對應另一條曲線。橫座標的標度相同,縱座標有兩個,左縱座標用於x1,y1資料對,右縱座標用於x2,y2資料對。
雙y軸座標可以用plotyy(x,y1,x,y2)來實現
雙x座標可以用
set(gca,'xaxislocation','bottom','xticklabel',{'0','1','2','3','4'}) (假設x軸的標註為1,2,3,4)
set(gca,'xaxislocation','top','xticklabel',{'0','1','2','3','4'})
進行相應的設定
【 * 例 10.7.3 -1 】製作一個雙座標系用來表現高壓和低溫兩個不同量的過渡過程。
tp=(0:100)/100*5;yp=8+4*(1-exp(-0.8*tp).*cos(3*tp)); % 壓力資料
tt=(0:500)/500*40;yt=120+40*(1-exp(-0.05*tt).*cos(tt)); % 溫度資料
% 產生雙座標系圖形
clf reset,h_ap=axes('Position',[0.13,0.13,0.7,0.75]); %<4>
set(h_ap,'Xcolor','b','Ycolor','b','Xlim',[0,5],'Ylim',[0,15]);
nx=10;ny=6; %<6>
pxtick=0:((5-0)/nx):5;pytick=0:((15-0)/ny):15; %<7>
set(h_ap,'Xtick',pxtick,'Ytick',pytick,'Xgrid','on','Ygrid','on')
h_linet=line(tp,yp,'Color','b'); %<9>
set(get(h_ap,'Xlabel'),'String',' 時間 /rightarrow (分) ')
set(get(h_ap,'Ylabel'),'String',' 壓力 /rightarrow(/times10 ^{5} Pa )')
h_at=axes('Position',get(h_ap,'Position')); %<12>
set(h_at,'Color','none','Xcolor','r','Ycolor','r'); %<13>
set(h_at,'Xaxislocation','top') %<14>
set(h_at,'Yaxislocation','right','Ydir','rev') %<15>
set(get(h_at,'Xlabel'),'String','/fontsize{15}/fontname{ 隸書 } 時間 /rightarrow (分) ')
set(get(h_at,'Ylabel'),'String',' ( {/circ}C )/fontsize{15} /leftarrow /fontname{ 隸書 } 零下溫度 ')
set(h_at,'Ylim',[0,210]) %<18>
line(tt,yt,'Color','r','Parent',h_at) %<19>
xpm=get(h_at,'Xlim'); %<20>
txtick=xpm(1):((xpm(2)-xpm(1))/nx):xpm(2); %<21>
tytick=0:((210-0)/ny):210; %<22>
set(h_at,'Xtick',txtick,'Ytick',tytick) %<23>
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
例項(已驗證):
clc;
clear all;
close all;
x=0:0.1:2*pi;
y1=sin(x);
y2=cos(x);
[AX]=plotyy(x,y1,x,y2);
set(get(gca,'xlabel'),'string','X-axis');
set(get(AX(1),'Ylabel'),'string','left Y-axis');
set(get(AX(2),'Ylabel'),'string','right Y-axis');
set(gca,'xTick',[0:0.5:7]);
set(AX(1),'yTick',[-1:0.2:1]);
set(AX(2),'yTick',[-1:0.5:1]);
尚存在問題:這種設定方法,對各個軸的最小刻度單位可以設定,但是刻度範圍(x取(0~7),y1取(-1~1))不能設定。
2010-12-23修改
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
clc
clear all
close all
runoff=[10700 11400 15800 22900 43100 40700 50500 46000 41800 35000];
sed=[0.105 0.094 0.156 1.264 0.363 0.429 0.731 0.682 0.654 0.290];
m=1:10;
[ax,h1,h2]=plotyy(m,runoff,m,sed); %h-- line handle
set(get(ax(1),'Ylabel'),'string','Runoff (m^3/s))','color','r') %y1
set(get(ax(2),'Ylabel'),'string','Sediment concentration (kg/m^3)','color','k') %y2
xlabel('Month')
set(h1,'linestyle','-','color','r');
set(h2,'linestyle','- -','color','k');
legend([h1 h2],'runoff','sediment concentration') %標註兩條線
legend('boxoff')
% box off
set(ax(:),'Ycolor','k') %設定兩個Y軸的顏色為黑色
set(ax(1),'ytick',[0:10000:100000]); %設定y軸間隔
set(ax(2),'ytick',[0:0.1:1.5])
set(ax,'xlim',[1 12]) % 設定x軸範圍
hold on
scatter(ax(1),4,22900,'r*')
axes(ax(2));
hold on
scatter(4,1.264,'ro')