%双横坐标绘制

​x1 = [0:.1:40];

y1 = 4.*cos(x1)./(x1+2);

x2 = [1:.2:20];

y2 = x2.^2./x2.^3;

hl1 = line(x1,y1,'Color','r');

ax1 = gca;

set(ax1,'XColor','r','YColor','r')

ax2 =

axes('Position',get(ax1,'Position'),...

'XAxisLocation','top',...

'YAxisLocation','right',...

'Color','none',...

'XColor','k'

,'YColor','k');

hl2 =

line(x2,y2,'Color','k','Parent',ax2);

 双横坐标绘制

​%双纵坐标轴

x = 0:0.01:20;

y1 = 200*exp(-0.05*x).*sin(x);

y2 = 0.8*exp(-0.5*x).*sin(10*x);

[AX,H1,H2] = plotyy(x,y1,x,y2,'plot');

set(get(AX(1),'Ylabel'),'String','SlowDecay')

set(get(AX(2),'Ylabel'),'String','FastDecay')

xlabel('Time (\musec)')

title('Multiple Decay Rates')

双纵坐标轴

​function [ax,hlines] =

plotyyy(x1,y1,x2,y2,x3,y3,ylabels)

%PLOTYYY - Extends plotyyto include a third

y-axis

%

%Syntax:  [ax,hlines]

=plotyyy(x1,y1,x2,y2,x3,y3,ylabels)

%

%Inputs: x1,y1 are thexdata and ydata for the first axes'

line

%

x2,y2 are the xdata and ydata for thesecond axes' line

%

x3,y3 are the

xdata and ydata for the thirdaxes' line

%

ylabels is a 3x1 cell array containingthe ylabel strings

%

%Outputs: ax

-

3x1 double array containing the axes'handles

%

hlines - 3x1 double array containingthe lines' handles

%

%Example:

%x=0:10;

%y1=x;

y2=x.^2;  y3=x.^3;

%ylabels{1}='Firsty-label';

%ylabels{2}='Secondy-label';

%ylabels{3}='Thirdy-label';

%[ax,hlines] =plotyyy(x,y1,x,y2,x,y3,ylabels);

%legend(hlines, 'y =x','y = x^2','y = x^3',2)

%

%m-files required: none

%Author: Denis Gilbert,Ph.D., physical

oceanography

%Maurice LamontagneInstitute

%email:gilbertd@dfo-mpo.gc.ca

%Web:http://www.qc.dfo-mpo.gc.ca/iml/

%April 2000; Lastrevision: 14-Nov-2001

if nargin==6

%Useempty

strings for the ylabels

ylabels{1}=' '; ylabels{2}='

'; ylabels{3}='

';

elseif nargin > 7

error('Too many input

arguments')

elseif nargin < 6

error('Not enough input

arguments')

end

figure('units','normalized',...

'DefaultAxesXMinorTick','on','DefaultAxesYminorTick','on');

% Plot the first twolines with plotyy

[ax,hlines(1),hlines(2)] =plotyy(x1,y1,x2,y2);

cfig =

get(gcf,'color');

pos = [0.1

0.1 0.7  0.8];

offset = pos(3)/5.5;

% Reduce width of the twoaxes generated by

plotyy

pos(3) = pos(3) - offset/2;

set(ax,'position',pos);

% Determine the positionof the third axes

pos3=[pos(1) pos(2)pos(3)+offset pos(4)];

�termine the properx-limits for the

third axes

limx1=get(ax(1),'xlim');

limx3=[limx1(1)

limx1(1) + 1.2*(limx1(2)-limx1(1))];

% Bug fix 14 Nov-2001:the 1.2 scale factor in the line

above

% was contributed byMariano Garcia (BorgWarner Morse TEC

Inc)

ax(3)=axes('Position',pos3,'box','off',...

'Color','none','XColor','k','YColor','r',...

'xtick',[],'xlim',limx3,'yaxislocation','right');

hlines(3) =

line(x3,y3,'Color','r','Parent',ax(3));

limy3=get(ax(3),'YLim');

% Hide unwanted portionof the x-axis line that

lies

% between the end of thesecond and third axes

line([limx1(2)limx3(2)],[limy3(1)

limy3(1)],...

'Color',cfig,'Parent',ax(3),'Clipping','off');

axes(ax(2))

% Label all three y-axes

set(get(ax(1),'ylabel'),'string',ylabels{1})

set(get(ax(2),'ylabel'),'string',ylabels{2})

set(get(ax(3),'ylabel'),'string',ylabels{3})

x=0:10;

y1=x;

y2=x.^2;  y3=x.^3;

ylabels{1}='First

y-label';

ylabels{2}='Second

y-label';

ylabels{3}='Third

y-label';

[ax,hlines] =plotyyy(x,y1,x,y2,x,y3,ylabels);

legend(hlines,'y =

x','y = x^2','y

= x^3',2)

%双纵双横坐标

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]);%设定左边的X、Y轴

nx=10;ny=6;%<6>定义两个变量。给下一语句使用

pxtick=0:((5-0)/nx):5;pytick=0:((15-0)/ny):15;%<7>

X轴被分了10格，Y轴被分了6格，pxtick，pytick分别表示向量

set(h_ap,'Xtick',pxtick,'Ytick',pytick,'Xgrid','on','Ygrid','on')%画出虚线格，'Xgrid','on','Ygrid','on'起到画虚线的作用

h_linet=line(tp,yp,'Color','b');%<9>画一yp波形，对波形进行颜色设置

set(get(h_ap,'Xlabel'),'String','时间

/rightarrow(分)

')%

X横坐标命名

set(get(h_ap,'Ylabel'),'String','压力/rightarrow(/times10

^{5} Pa

)')%Y坐标命名

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>%以上四句语句实现了上X坐标，右Y坐标

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>画yt波形，切记不能缺少Parent

xpm=get(h_at,'Xlim');%<20>定义新变量

txtick=xpm(1):((xpm(2)-xpm(1))/nx):xpm(2);%<21>

xpm(1)为0，xpm(2)为40，在哪条语句定义的？

tytick=0:((210-0)/ny):210;

%<22>%Y坐标分格

set(h_at,'Xtick',txtick,'Ytick',tytick)

%<23>