欢迎来这里逛逛，前些时间我的空间总是跟别人的重叠，搞得我也没心思来这里打理，今后可能好一点儿^_^

【分享】三轴绘图程序

上一篇 / 下一篇 2005-12-16 13:00:29

查看( 478 ) / 评论( 5 )

上次见到有人问怎么画三轴图形，印象中最后好像没结果了，今天碰巧在files exchange看到这个程序，搜索功能现在不好使，只有发新贴了，希望对这位朋友还有用。
function [ax,hlines] = plotyyy(x1,y1,x2,y2,x3,y3,ylabels)
%PLOTYYY - Extends plotyy to include a third y-axis
%
%Syntax:  [ax,hlines] = plotyyy(x1,y1,x2,y2,x3,y3,ylabels)
%
%Inputs: x1,y1 are the xdata and ydata for the first axes' line
%       x2,y2 are the xdata and ydata for the second axes' line
%       x3,y3 are the xdata and ydata for the third axes' line
%       ylabels is a 3x1 cell array containing the ylabel strings
%
%Outputs: ax -    3x1 double array containing the axes' handles
%       hlines - 3x1 double array containing the lines' handles

%Author: Denis Gilbert, Ph.D., physical oceanography
%Maurice Lamontagne Institute
%Dept. of Fisheries and Oceans Canada
%email: gilbertd@dfo-mpo.gc.ca
%Web: http://www.qc.dfo-mpo.gc.ca/iml/
%April 2000; Last revision: 14-Nov-2001

if nargin==6
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');

[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;

pos(3) = pos(3) - offset/2;
set(ax,'position',pos);

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

limx1=get(ax(1),'xlim');
limx3=[limx1(1) limx1(1) + 1.2*(limx1(2)-limx1(1))];

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');

line([limx1(2) limx3(2)],[limy3(1) limy3(1)],...
'Color',cfig,'Parent',ax(3),'Clipping','off');
axes(ax(2))

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

example:
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)

figure shown as below:

fomula.jpg

全部脚印 不留脚印留下脚印:

江风

论坛模式推荐收藏分享给好友推荐到圈子管理

TAG:

pengpai_1978 发布于2006-03-07 14:39:39: Re:【分享】三轴绘图程序
请问：
这句SimWe个人空间H8R w6kx7?0]W
line([limx1(2) limx3(2)],[limy3(1) limy3(1)],...
'Color',cfig,'Parent',ax(3),'Clipping','off');SimWe个人空间8I,T)K9t,l9mk?n

作用是裁剪掉比plotyy横坐标多出来的那段坐标轴吧？
tB$[qudSimWe个人空间那么'Clipping','on'时又有什么作用？
我怎么自己画图时试不出来？（想剪掉一段直线而未果）SimWe个人空间/\rb
_CpxU
究竟怎样才是裁剪？

_9y"L`+px3@YUSimWe个人空间例如：
>> x=1:1:5;
>> y=x;
.D2c3b*p{I>> plot(x,y)e.simwe.com-L+X GV$[
>> cfig=get(gcf,'color')
;jn9i+C,[J't/`d.k7qe.simwe.com>> line([1 3],[1 3],'Color',cfig,'clipping','off')%本意是想剪掉x=1:3这段;结果好像变成虚线了。怎样才能定义成无色？或真正剪掉？

bainhome 发布于2006-03-10 16:58:17: Re:【分享】三轴绘图程序
其实早已经看到这个问题，一直抽不出时间做，今天正好得空儿，写一个程序说明一下，你所说的clipping其实不叫“裁剪”而是自适应于axes或者figure的长度，所谓裁剪是根本就没有了，下面这个程序则可看出您的理解是错误的——线仍然存在，图形裁剪用nan赋值即可：
CODE:
function ClippingProterty clc close all x=-5:15; hFig1=figure; subplot(221) LineHandles=plot(x,x+5,'--r',x,x-3,'g'); TextHandles(1)=text(6.5,5,... 'This String will have clipping off'); TextHandles(2)=text(-1,3.5,... 'This String will have clipping off too'); axis([0 10 0 10]); subplot(222) LineHandles=plot(x,x+5,'--r',x,x-3,'g'); set(LineHandles(1),'clipping','off') TextHandles(1)=text(6.5,5,... 'This String will have clipping on'); TextHandles(2)=text(-1,3.5,... 'This String will have still been clipping off'); axis([0 10 0 10]); subplot(223) LineHandles=plot(x,x+5,'--r',x,x-3,'g'); hSubLine1=findobj(LineHandles,'color','r'); set(hSubLine1,'linewidth',2.5) subplot(224) y1=x+5; x([4:11])=NaN; LineHandles=plot(x,y1,'--r',[-5:15],[-5:15]-3,'g'); hSubLine1=findobj(LineHandles,'color','r'); set(hSubLine1,'linewidth',2.5)

拟和公式抓图1.jpg