|
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}) |
|
|
22 ## @deftypefnx {Function File} @var{h} = ribbon (@dots{}) |
|
|
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 |
|
7165
|
33 function h = ribbon (X, Y, W) |
|
7163
|
34 |
|
|
35 newplot (); |
|
|
36 |
|
|
37 if (nargin == 1) |
|
7165
|
38 Y = X; |
|
7164
|
39 if (isvector (Y)) |
|
7165
|
40 Y = Y(:); |
|
7163
|
41 endif |
|
7165
|
42 [nr, nc] = size (Y); |
|
7164
|
43 X = repmat ((1:nr)', 1, nc); |
|
7165
|
44 W = 0.75; |
|
7163
|
45 elseif (nargin == 2) |
|
7165
|
46 W = 0.75; |
|
7164
|
47 elseif (nargin != 3) |
|
|
48 print_usage (); |
|
7165
|
49 end |
|
7163
|
50 |
|
7165
|
51 if (isvector (X) && isvector (Y)) |
|
|
52 if (length (X) != length (Y)) |
|
|
53 error ("ribbon: in case of vectors, X and Y must have same length") |
|
7163
|
54 else |
|
7165
|
55 [X, Y] = meshgrid (X, Y); |
|
7163
|
56 endif |
|
|
57 else |
|
7165
|
58 if (! all (size (X) == size (Y))) |
|
|
59 error ("ribbon: in case of matrices, X and Y must have same size") |
|
7163
|
60 endif |
|
|
61 endif |
|
|
62 |
|
7165
|
63 [nr, nc] = size (Y); |
|
7164
|
64 tmp = zeros (1, nc); |
|
7163
|
65 |
|
|
66 for c = nc:-1:1 |
|
7165
|
67 ZZ = [Y(:,c), Y(:,c)]; |
|
|
68 y = X(:,c); |
|
|
69 x = [c - W / 2, c + W / 2]; |
|
|
70 [XX, YY] = meshgrid (x, y); |
|
|
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 |
