在MATLAB中绘制多纵坐标图

有网友问到这个问题。感谢bainhome 和shunfly 提供答案。SimWe个人空间Di\&o'i(Q:}

function [ax,hlines] = plotyyy(x1,y1,x2,y2,x3,y3,ylabels)
x%Yr!eu_T2t0%PLOTYYY - Extends plotyy to include a third y-axis SimWe个人空间-nxYk`8n]/x
% SimWe个人空间5sER-m7d
QtY
%Syntax:  [ax,hlines] = plotyyy(x1,y1,x2,y2,x3,y3,ylabels) SimWe个人空间%Sk%I'L mL9S
%
fT?`5NeG;fr0%Inputs: x1,y1 are the xdata and ydata for the first axes' line
C PR`$TH0D0%        x2,y2 are the xdata and ydata for the second axes' line SimWe个人空间5au hC"~G}hs`
%        x3,y3 are the xdata and ydata for the third axes' line SimWe个人空间*|j-w/a.Ha2F
%        ylabels is a 3x1 cell array containing the ylabel strings
+QR)q*k'Q,O0% SimWe个人空间\oI%d B
%Outputs: ax -     3x1 double array containing the axes' handles SimWe个人空间;cd SAj;y+}O[u
%         hlines - 3x1 double array containing the lines' handles SimWe个人空间%^}luW
%
7@e0J3o5M'w~:h0%Example: SimWe个人空间.tY |-t-Rx
% x=0:10;  SimWe个人空间&rcU"]-p
gZb*n*z
% y1=x;  y2=x.^2;   y3=x.^3;
ip&TH3CI _;s0% ylabels{1}='First y-label';
'd%Ij6e|"p0% ylabels{2}='Second y-label'; SimWe个人空间?/n^[h$L V Qp
% ylabels{3}='Third y-label';
k
x Y0\3FL5x0% [ax,hlines] = plotyyy(x,y1,x,y2,x,y3,ylabels);
vm-k&RD8t;C M0% legend(hlines, 'y = x','y = x^2','y = x^3',2) SimWe个人空间9{%eC"xX8BrQ Z
% SimWe个人空间
M[4FUl'Ya
%m-files required: none SimWe个人空间(l4UovA1o~N#r@

?DZy,QF-N0%Author: Denis Gilbert, Ph.D., physical oceanography SimWe个人空间:]P{ls+D3f
%Maurice Lamontagne Institute
Ny ~G6qfD0%Dept. of Fisheries and Oceans Canada SimWe个人空间_o9tOy4|_y
%email: [email]gilbertd@dfo-mpo.gc.ca[/email]   
8S7Oc
jg1H`6YU*k0%Web: [url]http://www.qc.dfo-mpo.gc.ca/iml/[/url] SimWe个人空间'X(T+efn
%April 2000; Last revision: 06-Nov-2001 SimWe个人空间9{X/D*` P"n
SimWe个人空间t(OAN};c
if nargin==6 SimWe个人空间!o-X6a!R,s0d zj-k
   %Use empty strings for the ylabels
I^;Zw7j?`+X0   ylabels{1}=' '; ylabels{2}=' '; ylabels{3}=' '; SimWe个人空间5}OP%H M.G4p el
elseif nargin > 7
Y_J/_+B
E0   error('Too many input arguments')
r ~&k n)]
w0elseif nargin < 6
A}:L0aF0   error('Not enough input arguments') SimWe个人空间"L1kq.{7?*^8D5I
end SimWe个人空间1o1J%C!a0Ugqy*KH
SimWe个人空间0L4};x5Ag,c/kK/H7v0}
figure('units','normalized',...
+^v2UvZ0       'DefaultAxesXMinorTick','off','DefaultAxesYminorTick','off');  
m][4fZvr'xr0
TC._9L!z0%Plot the first two lines with plotyy SimWe个人空间-X"iQ_*aA
[ax,hlines(1),hlines(2)] = plotyy(x1,y1,x2,y2); SimWe个人空间
i0Ce0z;] @e-f0[
cfig = get(gcf,'color');
1F$HWhy ^-b1a0pos = [0.1  0.1  0.7  0.8];
Z6j
V$N~%G0offset = pos(3)/5.5;
0v*u!t ^%sj.@}Qj0SimWe个人空间U/n2R1j!n^
D
%Reduce width of the two axes generated by plotyy  
t Y P3[(M1P l-@0pos(3) = pos(3) - offset/2;  SimWe个人空间&FdVZ}C;Td2\:\
set(ax,'position',pos);   SimWe个人空间7n
X@jW
SimWe个人空间~W&J/D1O'S/V \
%Set the position of the third axes
`}*{+MB(U;CFoa/Q0pos3=[pos(1) pos(2) pos(3)+offset pos(4)]; SimWe个人空间,R
pF0Vi
SimWe个人空间;_'iCI Okd,v
limx1=get(ax(1),'xlim');
{Y WH(d1_0xticks=get(ax(1),'xtick'); SimWe个人空间g'hl&G4_\ F[
delta_xtick = xticks(end) - xticks(end-1); SimWe个人空间!r#E,E hR S'y;gJ
limx3=[limx1(1)    limx1(2)+delta_xtick]; SimWe个人空间EqZ&A.Y tb6a6cs
SimWe个人空间wf?ED
ax(3)=axes('Position',pos3,'box','off',... SimWe个人空间V2e\UfM+D3e
   'Color','none','XColor','k','YColor','r',...   
9w%r@5tfn9S0   'xtick',[],'xlim',limx3,'yaxislocation','right'); SimWe个人空间L(j*b3NH"DV

#h Bn$s B {d@3b0hlines(3) = line(x3,y3,'Color','r','Parent',ax(3));
rH,k$Wr(O@ o0limy3=get(ax(3),'YLim'); SimWe个人空间S#p y&K/t[g1YK

4X%N.uxC7g4nj#^Z_0%Hide unwanted portion of the x-axis line that lies
PE#eO[B)t(Ot0%between the end of the second and third axes
5y!^BP"QRwcH0line([limx1(2) limx3(2)],[limy3(1) limy3(1)],...
FA*TcB'Gi \}5J0   'Color',cfig,'Parent',ax(3),'Clipping','off');%'off' SimWe个人空间0Dz)?wq9J2i
line([limx1(2) limx3(2)],[limy3(2) limy3(2)],... SimWe个人空间.ho/E_[H&CHwU
   'Color',cfig,'Parent',ax(3),'Clipping','off');%'off'
p&cQ P}6F\0axes(ax(2)) SimWe个人空间4}I.a}3[
SimWe个人空间lI6Ig^X&pB
%Label all three y-axes SimWe个人空间)hDe2w R2`
set(get(ax(1),'ylabel'),'string',ylabels{1})
6B3on&M'S0set(get(ax(2),'ylabel'),'string',ylabels{2}) SimWe个人空间J"pdqv/T
set(get(ax(3),'ylabel'),'string',ylabels{3})

0N1R Z3L4NR@G0效果图：

效果图