Mercurial > octave
view scripts/plot/draw/comet.m @ 27980:bdad8ca48700
maint: Remove "Created:" lines from code base.
* lin2mu.m, mu2lin.m, bincoeff.m, common_size.m, deal.m, nextpow2.m,
postpad.m, prepad.m, repmat.m, shift.m, xor.m, inpolygon.m, rotx.m, roty.m,
rotz.m, getappdata.m, isappdata.m, rmappdata.m, setappdata.m, colormap.m,
gray.m, gray2ind.m, image.m, imagesc.m, ind2gray.m, ind2rgb.m, ocean.m,
rgb2ind.m, javachk.m, javamem.m, commutation_matrix.m, cross.m,
duplication_matrix.m, gls.m, housh.m, isdefinite.m, ishermitian.m,
issymmetric.m, null.m, ols.m, orth.m, qzhess.m, vech.m, __w2mpth__.m,
expand_rel_paths.m, make_rel_paths.m, daspect.m, pbaspect.m, rticks.m,
thetaticks.m, xticklabels.m, xticks.m, yticklabels.m, yticks.m, zticklabels.m,
zticks.m, comet.m, errorbar.m, loglogerr.m, __errplot__.m, semilogxerr.m,
semilogyerr.m, cla.m, copyobj.m, hdl2struct.m, linkaxes.m, __ghostscript__.m,
__gnuplot_get_var__.m, __gnuplot_has_feature__.m, __gnuplot_has_terminal__.m,
__gnuplot_open_stream__.m, struct2hdl.m, compan.m, conv.m, deconv.m, mpoles.m,
poly.m, polyder.m, polyfit.m, polyint.m, polyout.m, polyreduce.m, polyval.m,
polyvalm.m, residue.m, roots.m, __parse_movargs__.m, detrend.m, fftconv.m,
fftfilt.m, fftshift.m, ifftshift.m, movfun.m, movslice.m, beta.m, betaln.m,
lcm.m, pow2.m, corr.m, kurtosis.m, skewness.m, untabify.m, calendar.m,
datestr.m, datevec.m, eomday.m, now.m, weekday.m:
Remove "Created:" lines from code base.
author | Rik <rik@octave.org> |
---|---|
date | Tue, 21 Jan 2020 13:06:02 -0800 |
parents | bd51beb6205e |
children | 9f9ac219896d |
line wrap: on
line source
######################################################################## ## ## Copyright (C) 2008-2020 The Octave Project Developers ## ## See the file COPYRIGHT.md in the top-level directory of this ## distribution or <https://octave.org/copyright/>. ## ## 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 ## <https://www.gnu.org/licenses/>. ## ######################################################################## ## -*- texinfo -*- ## @deftypefn {} {} comet (@var{y}) ## @deftypefnx {} {} comet (@var{x}, @var{y}) ## @deftypefnx {} {} comet (@var{x}, @var{y}, @var{p}) ## @deftypefnx {} {} comet (@var{hax}, @dots{}) ## Produce a simple comet style animation along the trajectory provided by ## the input coordinate vectors (@var{x}, @var{y}). ## ## If @var{x} is not specified it defaults to the indices of @var{y}. ## ## The speed of the comet may be controlled by @var{p}, which represents the ## time each point is displayed before moving to the next one. The default for ## @var{p} is @code{5 / numel (@var{y})}. ## ## If the first argument @var{hax} is an axes handle, then plot into this axes, ## rather than the current axes returned by @code{gca}. ## @seealso{comet3} ## @end deftypefn ## Author: Ben Abbott bpabbott@mac.com function comet (varargin) [hax, varargin, nargin] = __plt_get_axis_arg__ ("comet", varargin{:}); if (nargin == 0) print_usage (); elseif (nargin == 1) y = varargin{1}; x = 1:numel (y); p = 5 / numel (y); elseif (nargin == 2) x = varargin{1}; y = varargin{2}; p = 5 / numel (y); elseif (nargin == 3) x = varargin{1}; y = varargin{2}; p = varargin{3}; endif oldfig = []; if (! isempty (hax)) oldfig = get (0, "currentfigure"); endif unwind_protect hax = newplot (hax); limits = [min(x), max(x), min(y), max(y)]; num = numel (y); dn = round (num/10); hl = plot (x(1), y(1), "color", "r", "marker", "none", x(1), y(1), "color", "g", "marker", "none", x(1), y(1), "color", "b", "marker", "o"); axis (limits); # set manual limits to speed up plotting ## Initialize the timer t = p; timerid = tic (); for n = 2:(num+dn) m = n - dn; m = max ([m, 1]); k = min ([n, num]); set (hl(1), "xdata", x(1:m), "ydata", y(1:m)); set (hl(2), "xdata", x(m:k), "ydata", y(m:k)); set (hl(3), "xdata", x(k), "ydata", y(k)); pause (t - toc (timerid)); t += p; endfor unwind_protect_cleanup if (! isempty (oldfig)) set (0, "currentfigure", oldfig); endif end_unwind_protect endfunction %!demo %! clf; %! title ("comet() animation"); %! hold on; %! t = 0:.1:2*pi; %! x = cos (2*t) .* (cos (t).^2); %! y = sin (2*t) .* (sin (t).^2); %! comet (x, y, .05); %! hold off;