Mercurial > octave-antonio
view scripts/plot/ribbon.m @ 7165:5cae591f81a8
[project @ 2007-11-12 22:11:46 by jwe]
author | jwe |
---|---|
date | Mon, 12 Nov 2007 22:11:46 +0000 |
parents | 2b5e6c0a9df9 |
children | 4e7e5d51db86 |
line wrap: on
line source
## Copyright (C) 2007 Kai Habel ## ## This file is part of Octave. ## ## Octave is free software; you can redistribute it and/or modify it ## under the terms of the GNU General Public License as published by ## the Free Software Foundation; either version 3 of the License, or (at ## your option) any later version. ## ## Octave is distributed in the hope that it will be useful, but ## WITHOUT ANY WARRANTY; without even the implied warranty of ## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU ## General Public License for more details. ## ## You should have received a copy of the GNU General Public License ## along with Octave; see the file COPYING. If not, see ## <http://www.gnu.org/licenses/>. ## -*- texinfo -*- ## @deftypefn {Function File} ribbon (@var{x}, @var{y}, @var{width}) ## @deftypefnx {Function File} ribbon (@var{y}) ## @deftypefnx {Function File} @var{h} = ribbon (@dots{}) ## Plot a ribbon plot for the columns of @var{y} vs. @var{x}. The ## optional parameter @var{width} specifies the width of a single ribbon ## (default is 0.75). If @var{x} is omitted, a vector containing the ## row numbers is assumed (1:rows(Y)). If requested, return a vector ## @var{h} of the handles to the surface objects. ## @seealso{gca, colorbar} ## @end deftypefn ## Author: Kai Habel <kai.habel at gmx.de> function h = ribbon (X, Y, W) newplot (); if (nargin == 1) Y = X; if (isvector (Y)) Y = Y(:); endif [nr, nc] = size (Y); X = repmat ((1:nr)', 1, nc); W = 0.75; elseif (nargin == 2) W = 0.75; elseif (nargin != 3) print_usage (); end if (isvector (X) && isvector (Y)) if (length (X) != length (Y)) error ("ribbon: in case of vectors, X and Y must have same length") else [X, Y] = meshgrid (X, Y); endif else if (! all (size (X) == size (Y))) error ("ribbon: in case of matrices, X and Y must have same size") endif endif [nr, nc] = size (Y); tmp = zeros (1, nc); for c = nc:-1:1 ZZ = [Y(:,c), Y(:,c)]; y = X(:,c); x = [c - W / 2, c + W / 2]; [XX, YY] = meshgrid (x, y); CC = ones (size (ZZ)) * c; tmp(c) = surface (XX, YY, ZZ, CC); endfor ax = get (tmp(c), "parent"); if (! ishold ()) set (ax, "view", [-37.5, 30], "box", "off", "xgrid", "on", "ygrid", "on", "zgrid", "on"); endif if (nargout > 0) h = tmp; endif endfunction