Matlab plotyy画双纵坐标图实例

haidaowang

Matlab plotyy画双纵坐标图实例
) F5 b; t( |( R" s2 t
x = 0:0.01:20;
y1 = 200*exp(-0.05*x).*sin(x);
' o- a5 z/ f2 q: R" f1 Iy2 = 0.8*exp(-0.5*x).*sin(10*x);
[AX,H1,H2] = plotyy(x,y1,x,y2,'plot');

* c) K. z! K& xset(AX(1),'XColor','k','YColor','b');
. r, B: {8 f! Z2 }) e  x( tset(AX(2),'XColor','k','YColor','r');

HH1=get(AX(1),'Ylabel');
set(HH1,'String','Left Y-axis');
0 E1 x  I7 j/ u1 eset(HH1,'color','b');

HH2=get(AX(2),'Ylabel');
set(HH2,'String','Right Y-axis');
set(HH2,'color','r');
  g7 G, |# \& U( A- a( v
set(H1,'LineStyle','-');
set(H1,'color','b');
4 k/ j6 I! S' }8 }8 E/ j/ V. @set(H2,'LineStyle',':');
, S% f9 l) T2 i% b  ]set(H2,'color','r');

; }, \6 E! |) }! ]legend([H1,H2],{'y1 = 200*exp(-0.05*x).*sin(x)';'y2 = 0.8*exp(-0.5*x).*sin(10*x)'});
xlabel('Zero to 20 musec.');
title('Labeling plotyy');


! d+ n0 X6 k' k+ G( ?
Q：右边用蓝色圈起来的tick能去掉吗？由于用plotyy画图，为了使图尽量地显示出来，用了set(AX(1),'YLimMode','auto')，但这样可能会导致左边AX(1)和右边AX(2)的tick的间距不一样，影响美观。或者说能不能使plotyy画出的图两边的tick间距是一样的，这样在图形右边的tick就会重合在一起.
8 e$ w8 R" G* v$ W4 I
A：如果只是想让plotyy的图美一些，可以使用其如下形式的调用方式：
; ^# m: ^% Q0 ~4 r) y4 D[AX,H1,H2] = plotyy(...)
其中AX(2)就是右边Axes对象的句柄，拿到它以后就可以set或者get来处理了，也可以把其ytick关掉。

A：也可以用line语句来画，就没有左边和上边的线了。
- ]. z5 `7 Y: S' P
( i$ ~6 T# [3 n+ wQ：plotyy（X1，Y1，X2，Y2，FUN1，FUN2），FUN1和FUN2应该怎么写？
+ y; _/ d4 n- q' [1 i
# n1 p  {9 N+ A$ {A：这两个FUN代表plotyy不一定要用两个plot，比如下面的例子，一条曲线用plot，一条用semilogy
* Z, S8 R9 A# G
0 L" _" w7 u- a4 Jx1=1:0.1:100;
x2=x1;
y1=x1;
0 n! @5 U7 d4 ~3 k% ky2=x2.^3;
" X  P9 c! k" T8 Fplotyy(x1,y1,x2,y2,@plot,@semilogy)


~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
, L" f3 Z1 F. _/ K) I. j/ R
MATLAB画双纵坐标

具有两个纵坐标标度的图形
在MATLAB中，如果需要绘制出具有不同纵坐标标度的两个图形，可以使用plotyy绘图函数。调用格式为：
plotyy(x1,y1,x2,y2)
其中x1,y1对应一条曲线，x2,y2对应另一条曲线。横坐标的标度相同，纵坐标有两个，左纵坐标用于x1,y1数据对，右纵坐标用于x2,y2数据对。

双y轴坐标可以用plotyy(x,y1,x,y2)来实现
* a7 t* B3 T) U4 s双x坐标可以用
3 B8 f* ~0 E0 j( N# s0 Uset(gca,'xaxislocation','bottom','xticklabel',{'0','1','2','3','4'}) （假设x轴的标注为1，2，3，4）
! i" c6 G5 Z8 j$ D" rset(gca,'xaxislocation','top','xticklabel',{'0','1','2','3','4'})
进行相应的设置

. G, B7 t, _5 s( {: U【 * 例 10.7.3 -1 】制作一个双坐标系用来表现高压和低温两个不同量的过渡过程。
! Y2 o6 I$ M5 l- S* {" }
# v* u7 n6 y: G5 i9 Rtp=(0:100)/100*5;yp=8+4*(1-exp(-0.8*tp).*cos(3*tp)); % 压力数据
& _$ T- V7 J0 ~: M: k+ d- x
tt=(0:500)/500*40;yt=120+40*(1-exp(-0.05*tt).*cos(tt)); % 温度数据
  P4 t2 n. g" R
* k# X' T* Y5 K* B8 E% 产生双坐标系图形
$ n% L$ M, A6 k4 D( M$ v4 H
3 _5 I$ [' q3 Y( @; Q' Y' [clf reset,h_ap=axes('Position',[0.13,0.13,0.7,0.75]); %<4>
, s9 |6 p' p- F) x. }
set(h_ap,'Xcolor','b','Ycolor','b','Xlim',[0,5],'Ylim',[0,15]);
' W* y2 c9 `# c" @; D
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')
& S3 d9 T+ a7 p$ Z
h_linet=line(tp,yp,'Color','b'); %<9>
: g1 y0 h( q& c9 X" _4 m) J, a& m
0 e4 u0 c2 Y. c2 W, }) Eset(get(h_ap,'Xlabel'),'String',' 时间 /rightarrow （分） ')

set(get(h_ap,'Ylabel'),'String',' 压力 /rightarrow(/times10 ^{5} Pa )')
  S  i( S- C- r4 G. }, u5 w
, k+ Q# G& {  V$ Lh_at=axes('Position',get(h_ap,'Position')); %<12>

0 G' Z; {/ d+ F# R# S  V+ e+ Fset(h_at,'Color','none','Xcolor','r','Ycolor','r'); %<13>

set(h_at,'Xaxislocation','top') %<14>

+ v) C5 j7 ~- N% o2 Q8 b* ^set(h_at,'Yaxislocation','right','Ydir','rev') %<15>

set(get(h_at,'Xlabel'),'String','/fontsize{15}/fontname{ 隶书 } 时间 /rightarrow （分） ')
& {1 H0 n) x- D4 |  ^
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>
7 e/ M, x% _) z. a0 z4 R4 r
xpm=get(h_at,'Xlim'); %<20>

txtick=xpm(1)：（(xpm(2)-xpm(1))/nx):xpm(2); %<21>
( ^& Y$ ?7 P- n( ~% m8 b" B% ~+ s
tytick=0：（(210-0)/ny):210; %<22>

. Y; w) R7 E- O4 A6 i) |+ tset(h_at,'Xtick',txtick,'Ytick',tytick) %<23>

: q+ O1 x+ b2 N4 X& u1 ^
" Q0 ]& q! x9 p


( l- N& \# {- f: Y  C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
2 d, B' q, i) k& B/ B
' \0 G4 D2 E; T! r' {3 j& L( e实例（已验证）：
9 `! j0 q3 C. i9 h# G
) n9 d. P& q4 A4 Z2 S: Aclc;
clear all;
0 d8 h! l0 E) @. z' \7 N( {. O/ q7 yclose all;
x=0:0.1:2*pi;
* u$ A; U1 C7 {6 d& p/ ty1=sin(x);
8 r1 @/ ~. p7 Vy2=cos(x);
[AX]=plotyy(x,y1,x,y2);
- F' S- X6 x2 h' Mset(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]);
# Y" k: Q; A7 i5 h# tset(AX(2),'yTick',[-1:0.5:1]);
. m, z4 \4 H  q* J4 B- z7 Y

& j) F8 V5 `3 g% t
3 Q$ s' J: u+ I0 h. ~
尚存在问题：这种设置方法，对各个轴的最小刻度单位可以设置，但是刻度范围（x取(0～7),y1取(-1～1)）不能设置。
! j/ x+ k1 o* h' B3 T) F' P, _

~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

clc
& J  |# L8 O' D6 @6 c6 zclear 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
5 S* d2 k2 R8 g3 v& \set(get(ax(1),'Ylabel'),'string','Runoff (m^3/s))','color','r') %y1
! K# x0 d( B3 U" n7 B. \set(get(ax(2),'Ylabel'),'string','Sediment concentration (kg/m^3)','color','k') %y2
xlabel('Month')
" @5 l, V- v. }0 p$ x6 D4 Cset(h1,'linestyle','-','color','r');   
' g) }: i0 c7 d6 ^! U( Kset(h2,'linestyle','- -','color','k');
legend([h1 h2],'runoff','sediment concentration') %标注两条线
7 l  f: u" l2 `  C% \7 u0 alegend('boxoff')
! @, z3 S+ C! e: O( \% box off
- N$ `% z8 k0 b: f0 T* ], ]7 wset(ax(：）,'Ycolor','k') %设定两个Y轴的颜色为黑色
set(ax(1),'ytick',[0:10000:100000]); %设置y轴间隔
- n2 D8 f. C1 A0 vset(ax(2),'ytick',[0:0.1:1.5])
4 u" B0 H* s% T$ `4 ^6 r- r: Eset(ax,'xlim',[1 12]) % 设置x轴范围
hold on
  G! n- X& b, A% cscatter(ax(1),4,22900,'r*')
axes(ax(2));
hold on
3 Y9 k+ y% v7 O" Uscatter(4,1.264,'ro')

) l, Z8 C  r9 F* L7 v


7 \. i6 f) M( d9 }+ d6 m& W& ^# a
, T8 V  _( v: z! u

ExxNEN

Matlab plotyy画双纵坐标图实例