annotate scripts/plot/plotyy.m @ 7215:dd88d61d443f

[project @ 2007-11-29 16:44:45 by jwe]
author jwe
date Thu, 29 Nov 2007 16:46:40 +0000
parents a730e47fda4d
children 5389a52df87b
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
7195
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
1 ## Copyright (C) 2007 David Bateman
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
2 ##
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
3 ## This file is part of Octave.
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
4 ##
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
5 ## Octave is free software; you can redistribute it and/or modify it
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
6 ## under the terms of the GNU General Public License as published by
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
7 ## the Free Software Foundation; either version 3 of the License, or (at
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
8 ## your option) any later version.
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
9 ##
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
10 ## Octave is distributed in the hope that it will be useful, but
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
11 ## WITHOUT ANY WARRANTY; without even the implied warranty of
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
12 ## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
13 ## General Public License for more details.
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
14 ##
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
15 ## You should have received a copy of the GNU General Public License
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
16 ## along with Octave; see the file COPYING. If not, see
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
17 ## <http://www.gnu.org/licenses/>.
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
18
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
19 ## -*- texinfo -*-
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
20 ## @deftypefn {Function File} {} plotyy (@var{x1}, @var{y1}, @var{x2}, @var{y2})
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
21 ## @deftypefnx {Function File} {} plotyy (@dots{}, @var{fun})
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
22 ## @deftypefnx {Function File} {} plotyy (@dots{}, @var{fun1}, @var{fun2})
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
23 ## @deftypefnx {Function File} {} plotyy (@var{h}, @dots{})
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
24 ## @deftypefnx {Function File} {[@var{ax}, @var{h1}, @var{h2}] =} plotyy (@dots{})
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
25 ## Plots two sets of data with independent y-axes. The arguments @var{x1} and
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
26 ## @var{y1} define the arguments for the first plot and @var{x1} and @var{y2}
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
27 ## for the second.
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
28 ##
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
29 ## By default the arguments are evaluated with
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
30 ## @code{feval (@@plot, @var{x}, @var{y})}. However the type of plot can be
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
31 ## modified with the @var{fun} argument, in which case the plots are
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
32 ## generated by @code{feval (@var{fun}, @var{x}, @var{y})}. @var{fun} can be
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
33 ## a function handle, an inline function or a string of a function name.
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
34 ##
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
35 ## The function to use for each of the plots can be independently defined
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
36 ## with @var{fun1} and @var{fun2}.
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
37 ##
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
38 ## If given, @var{h} defines the principal axis in which to plot the @var{x1}
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
39 ## and @var{y1} data. The return value @var{ax} is a two element vector with
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
40 ## the axis handles of the two plots. @var{h1} and @var{h2} are handles to
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
41 ## the objects generated by the plot commands.
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
42 ##
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
43 ## @example
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
44 ## @group
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
45 ## x = 0:0.1:2*pi;
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
46 ## y1 = sin (x);
7196
9092375e3936 [project @ 2007-11-26 23:39:26 by jwe]
jwe
parents: 7195
diff changeset
47 ## y2 = exp (x - 1);
9092375e3936 [project @ 2007-11-26 23:39:26 by jwe]
jwe
parents: 7195
diff changeset
48 ## ax = plotyy (x, y1, x - 1, y2, @@plot, @@semilogy);
7195
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
49 ## xlabel ("X");
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
50 ## ylabel (ax(1), "Axis 1");
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
51 ## ylabel (ax(2), "Axis 2");
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
52 ## @end group
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
53 ## @end example
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
54 ## @end deftypefn
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
55
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
56 function [Ax, H1, H2] = plotyy (varargin)
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
57
7215
dd88d61d443f [project @ 2007-11-29 16:44:45 by jwe]
jwe
parents: 7208
diff changeset
58 [ax, varargin] = __plt_get_axis_arg__ ("plotyy", varargin{:});
dd88d61d443f [project @ 2007-11-29 16:44:45 by jwe]
jwe
parents: 7208
diff changeset
59
dd88d61d443f [project @ 2007-11-29 16:44:45 by jwe]
jwe
parents: 7208
diff changeset
60 if (nargin < 4)
dd88d61d443f [project @ 2007-11-29 16:44:45 by jwe]
jwe
parents: 7208
diff changeset
61 print_usage ();
dd88d61d443f [project @ 2007-11-29 16:44:45 by jwe]
jwe
parents: 7208
diff changeset
62 endif
dd88d61d443f [project @ 2007-11-29 16:44:45 by jwe]
jwe
parents: 7208
diff changeset
63 oldh = gca ();
dd88d61d443f [project @ 2007-11-29 16:44:45 by jwe]
jwe
parents: 7208
diff changeset
64 unwind_protect
dd88d61d443f [project @ 2007-11-29 16:44:45 by jwe]
jwe
parents: 7208
diff changeset
65 axes (ax);
7195
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
66 newplot ();
7215
dd88d61d443f [project @ 2007-11-29 16:44:45 by jwe]
jwe
parents: 7208
diff changeset
67 [ax, h1, h2] = __plotyy__ (ax, varargin{:});
dd88d61d443f [project @ 2007-11-29 16:44:45 by jwe]
jwe
parents: 7208
diff changeset
68 unwind_protect_cleanup
dd88d61d443f [project @ 2007-11-29 16:44:45 by jwe]
jwe
parents: 7208
diff changeset
69 axes (oldh);
dd88d61d443f [project @ 2007-11-29 16:44:45 by jwe]
jwe
parents: 7208
diff changeset
70 end_unwind_protect
7195
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
71
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
72 if (nargout > 0)
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
73 Ax = ax;
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
74 H1 = h1;
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
75 H2 = h2;
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
76 endif
7196
9092375e3936 [project @ 2007-11-26 23:39:26 by jwe]
jwe
parents: 7195
diff changeset
77
7195
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
78 endfunction
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
79
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
80 function [ax, h1, h2] = __plotyy__ (ax, x1, y1, x2, y2, varargin)
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
81 if (nargin > 5)
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
82 fun1 = varargin{1};
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
83 else
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
84 fun1 = @plot;
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
85 endif
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
86 if (nargin > 6)
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
87 fun2 = varargin{2};
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
88 else
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
89 fun2 = fun1;
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
90 endif
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
91
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
92 xlim = [min([x1(:); x2(:)]), max([x1(:); x2(:)])];
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
93
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
94 h1 = feval (fun1, x1, y1);
7206
6e4ceeeb1940 [project @ 2007-11-27 22:40:20 by jwe]
jwe
parents: 7196
diff changeset
95 set (ax(1), "ycolor", get (h1(1), "color"));
7195
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
96 set (ax(1), "position", get (ax(1), "outerposition"));
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
97 set (ax(1), "xlim", xlim);
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
98
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
99 cf = gcf ();
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
100 set (cf, "nextplot", "add");
7206
6e4ceeeb1940 [project @ 2007-11-27 22:40:20 by jwe]
jwe
parents: 7196
diff changeset
101 ax(2) = axes ();
7195
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
102 colors = get (ax(1), "colororder");
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
103 set (ax(2), "colororder", [colors(2:end,:); colors(1,:)]);
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
104
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
105 h2 = feval (fun2, x2, y2);
7206
6e4ceeeb1940 [project @ 2007-11-27 22:40:20 by jwe]
jwe
parents: 7196
diff changeset
106 set (ax(2), "yaxislocation", "right");
6e4ceeeb1940 [project @ 2007-11-27 22:40:20 by jwe]
jwe
parents: 7196
diff changeset
107 set (ax(2), "ycolor", get (h2(1), "color"));
6e4ceeeb1940 [project @ 2007-11-27 22:40:20 by jwe]
jwe
parents: 7196
diff changeset
108 set (ax(2), "position", get (ax(1), "outerposition"));
7195
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
109 set (ax(2), "xlim", xlim);
97d2da8ed746 [project @ 2007-11-26 23:13:40 by dbateman]
dbateman
parents:
diff changeset
110 endfunction