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