|
7163
|
1 ## Copyright (C) 2007 Kai Habel |
|
|
2 ## |
|
7164
|
3 ## This file is part of Octave. |
|
|
4 ## |
|
|
5 ## Octave is free software; you can redistribute it and/or modify it |
|
7163
|
6 ## under the terms of the GNU General Public License as published by |
|
7164
|
7 ## the Free Software Foundation; either version 3 of the License, or (at |
|
|
8 ## your option) any later version. |
|
7163
|
9 ## |
|
7164
|
10 ## Octave is distributed in the hope that it will be useful, but |
|
7163
|
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 |
|
7164
|
16 ## along with Octave; see the file COPYING. If not, see |
|
|
17 ## <http://www.gnu.org/licenses/>. |
|
7163
|
18 |
|
|
19 ## -*- texinfo -*- |
|
7164
|
20 ## @deftypefn {Function File} ribbon (@var{x}, @var{y}, @var{width}) |
|
|
21 ## @deftypefnx {Function File} ribbon (@var{y}) |
|
7650
|
22 ## @deftypefnx {Function File} {@var{h} =} ribbon (@dots{}) |
|
7164
|
23 ## Plot a ribbon plot for the columns of @var{y} vs. @var{x}. The |
|
|
24 ## optional parameter @var{width} specifies the width of a single ribbon |
|
|
25 ## (default is 0.75). If @var{x} is omitted, a vector containing the |
|
|
26 ## row numbers is assumed (1:rows(Y)). If requested, return a vector |
|
|
27 ## @var{h} of the handles to the surface objects. |
|
|
28 ## @seealso{gca, colorbar} |
|
7163
|
29 ## @end deftypefn |
|
|
30 |
|
|
31 ## Author: Kai Habel <kai.habel at gmx.de> |
|
|
32 |
|
7166
|
33 function h = ribbon (x, y, width) |
|
7163
|
34 |
|
|
35 newplot (); |
|
|
36 |
|
|
37 if (nargin == 1) |
|
7166
|
38 y = x; |
|
|
39 if (isvector (y)) |
|
|
40 y = y(:); |
|
7163
|
41 endif |
|
7166
|
42 [nr, nc] = size (y); |
|
|
43 x = repmat ((1:nr)', 1, nc); |
|
|
44 width = 0.75; |
|
7163
|
45 elseif (nargin == 2) |
|
7166
|
46 width = 0.75; |
|
7164
|
47 elseif (nargin != 3) |
|
|
48 print_usage (); |
|
7165
|
49 end |
|
7163
|
50 |
|
7166
|
51 if (isvector (x) && isvector (y)) |
|
|
52 if (length (x) != length (y)) |
|
7165
|
53 error ("ribbon: in case of vectors, X and Y must have same length") |
|
7163
|
54 else |
|
7166
|
55 [x, y] = meshgrid (x, y); |
|
7163
|
56 endif |
|
|
57 else |
|
7292
|
58 if (! size_equal(x, y)) |
|
7165
|
59 error ("ribbon: in case of matrices, X and Y must have same size") |
|
7163
|
60 endif |
|
|
61 endif |
|
|
62 |
|
7166
|
63 [nr, nc] = size (y); |
|
7164
|
64 tmp = zeros (1, nc); |
|
7163
|
65 |
|
|
66 for c = nc:-1:1 |
|
7166
|
67 zz = [y(:,c), y(:,c)]; |
|
|
68 yy = x(:,c); |
|
|
69 xx = [c - width / 2, c + width / 2]; |
|
|
70 [xx, yy] = meshgrid (xx, yy); |
|
|
71 cc = ones (size (zz)) * c; |
|
|
72 tmp(c) = surface (xx, yy, zz, cc); |
|
7163
|
73 endfor |
|
|
74 |
|
|
75 ax = get (tmp(c), "parent"); |
|
|
76 |
|
7164
|
77 if (! ishold ()) |
|
|
78 set (ax, "view", [-37.5, 30], "box", "off", "xgrid", "on", |
|
|
79 "ygrid", "on", "zgrid", "on"); |
|
7163
|
80 endif |
|
|
81 |
|
|
82 if (nargout > 0) |
|
|
83 h = tmp; |
|
|
84 endif |
|
7164
|
85 |
|
|
86 endfunction |
|
7269
|
87 |
|
|
88 %!demo |
|
|
89 %! [x, y, z] = sombrero (); |
|
|
90 %! [x, y] = meshgrid (x, y); |
|
|
91 %! ribbon (y, z); |
