changelog shortlog tags changeset files revisions annotate raw

scripts/plot/feather.m

changeset 10289: 4b124317dc38
parent:16f53d29049f
author: John W. Eaton <jwe@octave.org>
date: Tue Feb 09 20:58:55 2010 -0500 (30 minutes ago)
permissions: -rw-r--r--
description: base_properties::set_children: account for hidden children 
     1## Copyright (C) 2007, 2008, 2009 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} {} feather (@var{u}, @var{v})

    21## @deftypefnx {Function File} {} feather (@var{z})

    22## @deftypefnx {Function File} {} feather (@dots{}, @var{style})

    23## @deftypefnx {Function File} {} feather (@var{h}, @dots{})

    24## @deftypefnx {Function File} {@var{h} =} feather (@dots{})

    25##

    26## Plot the @code{(@var{u}, @var{v})} components of a vector field emanating

    27## from equidistant points on the x-axis.  If a single complex argument

    28## @var{z} is given, then @code{@var{u} = real (@var{z})} and

    29## @code{@var{v} = imag (@var{z})}.

    30##

    31## The style to use for the plot can be defined with a line style @var{style}

    32## in a similar manner to the line styles used with the @code{plot} command.

    33##

    34## The optional return value @var{h} provides a list of handles to the 

    35## the parts of the vector field (body, arrow and marker).

    36##

    37## @example

    38## @group

    39## phi = [0 : 15 : 360] * pi / 180;

    40## feather (sin (phi), cos (phi))

    41## @end group

    42## @end example

    43##

    44## @seealso{plot, quiver, compass}

    45## @end deftypefn

    46

    47function retval = feather (varargin)

    48

    49  [h, varargin, nargin] = __plt_get_axis_arg__ ("feather", varargin{:});

    50

    51  arrowsize = 0.25;

    52

    53  if (nargin == 0)

    54    print_usage ();

    55  elseif (nargin == 1 || (nargin == 2 && ! isnumeric (varargin{2})))

    56    ioff = 2;

    57    z = varargin{1}(:).';

    58    u = real (z);

    59    v = imag (z);

    60  elseif (nargin > 1 && isnumeric (varargin{2}))

    61    ioff = 3;

    62    u = varargin{1}(:).';

    63    v = varargin{2}(:).';

    64  endif

    65

    66  line_spec = "b-";

    67  have_line_spec = false;

    68  while (ioff <= nargin)

    69    arg = varargin{ioff++};

    70    if ((ischar (arg) || iscell (arg)) && ! have_line_spec)

    71      [linespec, valid] = __pltopt__ ("feather", arg, false);

    72      if (valid)

    73	line_spec = arg;

    74	have_line_spec = false;

    75	break;

    76      else

    77	error ("feather: invalid linespec");

    78      endif

    79    else

    80      error ("feather: unrecognized argument");

    81    endif

    82  endwhile

    83

    84  ## Matlab draws feather plots, with the arrow head as one continous 

    85  ## line, and each arrow separately. This is completely different than 

    86  ## quiver and quite ugly.

    87  n = length (u);

    88  xend = [1 : n] + u;

    89  xtmp = [1 : n] + u .* (1 - arrowsize);

    90  yend = v;

    91  ytmp = v .* (1 - arrowsize);

    92  x = [[1 : n]; xend; xtmp  - v * arrowsize; xend; ...

    93       xtmp + v * arrowsize];

    94  y = [zeros(1, n); yend; ytmp  + u * arrowsize / 3; yend; ...

    95       ytmp - u * arrowsize / 3];

    96

    97  oldh = gca ();

    98  unwind_protect

    99    axes (h);

   100    newplot ();

   101    hlist = plot (h, x, y, line_spec, [1, n], [0, 0], line_spec);

   102  unwind_protect_cleanup

   103    axes (oldh);

   104  end_unwind_protect

   105

   106  if (nargout > 0)

   107    retval = hlist;

   108  endif

   109

   110endfunction

   111

   112%!demo

   113%! phi = [0 : 15 : 360] * pi / 180;

   114%! feather (sin (phi), cos (phi))