octave-nkf: scripts/plot/polar.m annotate

annotate scripts/plot/polar.m @ 12789:482cc053281b

codesprint: Add 2 demos for polar.m * polar.m : Add 2 rose demos.

author	Rik <octave@nomad.inbox5.com>
date	Sat, 16 Jul 2011 08:16:40 -0700
parents	c792872f8942
children	5f0bb45e615c

rev	line source
11523 fd0a3ac60b0e update copyright notices John W. Eaton <jwe@octave.org> parents: 10730 diff changeset	1 ## Copyright (C) 1993-2011 John W. Eaton
2313 5ca126254d15 [project @ 1996-07-11 21:25:22 by jwe] jwe parents: 2311 diff changeset	2 ##
5ca126254d15 [project @ 1996-07-11 21:25:22 by jwe] jwe parents: 2311 diff changeset	3 ## This file is part of Octave.
5ca126254d15 [project @ 1996-07-11 21:25:22 by jwe] jwe parents: 2311 diff changeset	4 ##
5ca126254d15 [project @ 1996-07-11 21:25:22 by jwe] jwe parents: 2311 diff changeset	5 ## Octave is free software; you can redistribute it and/or modify it
5ca126254d15 [project @ 1996-07-11 21:25:22 by jwe] jwe parents: 2311 diff changeset	6 ## under the terms of the GNU General Public License as published by
7016 93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 6895 diff changeset	7 ## the Free Software Foundation; either version 3 of the License, or (at
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 6895 diff changeset	8 ## your option) any later version.
2313 5ca126254d15 [project @ 1996-07-11 21:25:22 by jwe] jwe parents: 2311 diff changeset	9 ##
5ca126254d15 [project @ 1996-07-11 21:25:22 by jwe] jwe parents: 2311 diff changeset	10 ## Octave is distributed in the hope that it will be useful, but
5ca126254d15 [project @ 1996-07-11 21:25:22 by jwe] jwe parents: 2311 diff changeset	11 ## WITHOUT ANY WARRANTY; without even the implied warranty of
5ca126254d15 [project @ 1996-07-11 21:25:22 by jwe] jwe parents: 2311 diff changeset	12 ## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
5ca126254d15 [project @ 1996-07-11 21:25:22 by jwe] jwe parents: 2311 diff changeset	13 ## General Public License for more details.
5ca126254d15 [project @ 1996-07-11 21:25:22 by jwe] jwe parents: 2311 diff changeset	14 ##
5ca126254d15 [project @ 1996-07-11 21:25:22 by jwe] jwe parents: 2311 diff changeset	15 ## You should have received a copy of the GNU General Public License
7016 93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 6895 diff changeset	16 ## along with Octave; see the file COPYING. If not, see
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 6895 diff changeset	17 ## <http://www.gnu.org/licenses/>.
245 16a24e76d6e0 [project @ 1993-12-03 02:00:15 by jwe] jwe parents: 4 diff changeset	18
3368 a4cd1e9d9962 [project @ 1999-11-20 17:22:48 by jwe] jwe parents: 3063 diff changeset	19 ## -- texinfo --
10730 390d93e20531 Plot commands now print usage information without errors Rik <octave@nomad.inbox5.com> parents: 10549 diff changeset	20 ## @deftypefn {Function File} {} polar (@var{theta}, @var{rho})
390d93e20531 Plot commands now print usage information without errors Rik <octave@nomad.inbox5.com> parents: 10549 diff changeset	21 ## @deftypefnx {Function File} {} polar (@var{theta}, @var{rho}, @var{fmt})
390d93e20531 Plot commands now print usage information without errors Rik <octave@nomad.inbox5.com> parents: 10549 diff changeset	22 ## @deftypefnx {Function File} {} polar (@var{h}, @dots{})
390d93e20531 Plot commands now print usage information without errors Rik <octave@nomad.inbox5.com> parents: 10549 diff changeset	23 ## @deftypefnx {Function File} {@var{h} =} polar (@dots{})
7180 55ff2d9a27c7 [project @ 2007-11-25 04:32:24 by dbateman] dbateman parents: 7017 diff changeset	24 ## Make a two-dimensional plot given the polar coordinates @var{theta} and
3368 a4cd1e9d9962 [project @ 1999-11-20 17:22:48 by jwe] jwe parents: 3063 diff changeset	25 ## @var{rho}.
2311 2b5788792cad [project @ 1996-07-11 20:18:38 by jwe] jwe parents: 2303 diff changeset	26 ##
2b5788792cad [project @ 1996-07-11 20:18:38 by jwe] jwe parents: 2303 diff changeset	27 ## The optional third argument specifies the line type.
6895 76e3d985ae56 [project @ 2007-09-13 18:22:38 by jwe] jwe parents: 6448 diff changeset	28 ## @seealso{plot}
3368 a4cd1e9d9962 [project @ 1999-11-20 17:22:48 by jwe] jwe parents: 3063 diff changeset	29 ## @end deftypefn
4 b4df021f796c [project @ 1993-08-08 01:26:08 by jwe] jwe parents: diff changeset	30
2314 949ab8eba8bc [project @ 1996-07-12 03:58:02 by jwe] jwe parents: 2313 diff changeset	31 ## Author: jwe
949ab8eba8bc [project @ 1996-07-12 03:58:02 by jwe] jwe parents: 2313 diff changeset	32
6302 a5cd8b77e892 [project @ 2007-02-13 08:08:33 by jwe] jwe parents: 6257 diff changeset	33 function retval = polar (varargin)
4 b4df021f796c [project @ 1993-08-08 01:26:08 by jwe] jwe parents: diff changeset	34
7321 359f464342b3 [project @ 2007-12-18 03:48:03 by jwe] jwe parents: 7216 diff changeset	35 [h, varargin, nargs] = __plt_get_axis_arg__ ("polar", varargin{:});
7216 5389a52df87b [project @ 2007-11-29 19:07:29 by jwe] jwe parents: 7215 diff changeset	36
10730 390d93e20531 Plot commands now print usage information without errors Rik <octave@nomad.inbox5.com> parents: 10549 diff changeset	37 if (nargs < 1)
11587 c792872f8942 all script files: untabify and strip trailing whitespace John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	38 print_usage();
10730 390d93e20531 Plot commands now print usage information without errors Rik <octave@nomad.inbox5.com> parents: 10549 diff changeset	39 endif
390d93e20531 Plot commands now print usage information without errors Rik <octave@nomad.inbox5.com> parents: 10549 diff changeset	40
7207 71c03c7239fb [project @ 2007-11-28 02:19:25 by jwe] jwe parents: 7180 diff changeset	41 oldh = gca ();
71c03c7239fb [project @ 2007-11-28 02:19:25 by jwe] jwe parents: 7180 diff changeset	42 unwind_protect
71c03c7239fb [project @ 2007-11-28 02:19:25 by jwe] jwe parents: 7180 diff changeset	43 axes (h);
71c03c7239fb [project @ 2007-11-28 02:19:25 by jwe] jwe parents: 7180 diff changeset	44 newplot ();
3063 4287b8c06fbf [project @ 1997-06-26 02:36:05 by jwe] jwe parents: 2847 diff changeset	45
7207 71c03c7239fb [project @ 2007-11-28 02:19:25 by jwe] jwe parents: 7180 diff changeset	46 if (nargs == 3)
71c03c7239fb [project @ 2007-11-28 02:19:25 by jwe] jwe parents: 7180 diff changeset	47 if (! ischar (varargin{3}))
10549 95c3e38098bf Untabify .m scripts Rik <code@nomad.inbox5.com> parents: 10226 diff changeset	48 error ("polar: third argument must be a string");
7207 71c03c7239fb [project @ 2007-11-28 02:19:25 by jwe] jwe parents: 7180 diff changeset	49 endif
71c03c7239fb [project @ 2007-11-28 02:19:25 by jwe] jwe parents: 7180 diff changeset	50 tmp = __plr2__ (h, varargin{:});
7321 359f464342b3 [project @ 2007-12-18 03:48:03 by jwe] jwe parents: 7216 diff changeset	51 maxr = max (varargin {2} (:));
359f464342b3 [project @ 2007-12-18 03:48:03 by jwe] jwe parents: 7216 diff changeset	52 elseif (nargs == 2)
7207 71c03c7239fb [project @ 2007-11-28 02:19:25 by jwe] jwe parents: 7180 diff changeset	53 if (ischar (varargin{2}))
10549 95c3e38098bf Untabify .m scripts Rik <code@nomad.inbox5.com> parents: 10226 diff changeset	54 tmp = __plr1__ (h, varargin{:});
95c3e38098bf Untabify .m scripts Rik <code@nomad.inbox5.com> parents: 10226 diff changeset	55 if (iscomplex(varargin{1}))
95c3e38098bf Untabify .m scripts Rik <code@nomad.inbox5.com> parents: 10226 diff changeset	56 maxr = max (imag(varargin{1})(:));
95c3e38098bf Untabify .m scripts Rik <code@nomad.inbox5.com> parents: 10226 diff changeset	57 else
95c3e38098bf Untabify .m scripts Rik <code@nomad.inbox5.com> parents: 10226 diff changeset	58 maxr = max (varargin{1}(:));
95c3e38098bf Untabify .m scripts Rik <code@nomad.inbox5.com> parents: 10226 diff changeset	59 endif
7207 71c03c7239fb [project @ 2007-11-28 02:19:25 by jwe] jwe parents: 7180 diff changeset	60 else
10549 95c3e38098bf Untabify .m scripts Rik <code@nomad.inbox5.com> parents: 10226 diff changeset	61 fmt = "";
95c3e38098bf Untabify .m scripts Rik <code@nomad.inbox5.com> parents: 10226 diff changeset	62 tmp = __plr2__ (h, varargin{:}, fmt);
95c3e38098bf Untabify .m scripts Rik <code@nomad.inbox5.com> parents: 10226 diff changeset	63 maxr = max (varargin {2} (:));
7207 71c03c7239fb [project @ 2007-11-28 02:19:25 by jwe] jwe parents: 7180 diff changeset	64 endif
7321 359f464342b3 [project @ 2007-12-18 03:48:03 by jwe] jwe parents: 7216 diff changeset	65 elseif (nargs == 1)
7207 71c03c7239fb [project @ 2007-11-28 02:19:25 by jwe] jwe parents: 7180 diff changeset	66 fmt = "";
71c03c7239fb [project @ 2007-11-28 02:19:25 by jwe] jwe parents: 7180 diff changeset	67 tmp = __plr1__ (h, varargin{:}, fmt);
7321 359f464342b3 [project @ 2007-12-18 03:48:03 by jwe] jwe parents: 7216 diff changeset	68 if (iscomplex(varargin{1}))
10549 95c3e38098bf Untabify .m scripts Rik <code@nomad.inbox5.com> parents: 10226 diff changeset	69 maxr = max (imag(varargin{1})(:));
7321 359f464342b3 [project @ 2007-12-18 03:48:03 by jwe] jwe parents: 7216 diff changeset	70 else
10549 95c3e38098bf Untabify .m scripts Rik <code@nomad.inbox5.com> parents: 10226 diff changeset	71 maxr = max (varargin{1}(:));
7321 359f464342b3 [project @ 2007-12-18 03:48:03 by jwe] jwe parents: 7216 diff changeset	72 endif
7207 71c03c7239fb [project @ 2007-11-28 02:19:25 by jwe] jwe parents: 7180 diff changeset	73 else
71c03c7239fb [project @ 2007-11-28 02:19:25 by jwe] jwe parents: 7180 diff changeset	74 print_usage ();
934 be867ddbe4aa [project @ 1994-11-17 04:23:41 by jwe] jwe parents: 904 diff changeset	75 endif
7207 71c03c7239fb [project @ 2007-11-28 02:19:25 by jwe] jwe parents: 7180 diff changeset	76
7321 359f464342b3 [project @ 2007-12-18 03:48:03 by jwe] jwe parents: 7216 diff changeset	77 set (h, "xlim", [-maxr, maxr], "ylim", [-maxr, maxr],
10549 95c3e38098bf Untabify .m scripts Rik <code@nomad.inbox5.com> parents: 10226 diff changeset	78 "xaxislocation", "zero", "yaxislocation", "zero",
11587 c792872f8942 all script files: untabify and strip trailing whitespace John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	79 "plotboxaspectratio", [1, 1, 1]);
7321 359f464342b3 [project @ 2007-12-18 03:48:03 by jwe] jwe parents: 7216 diff changeset	80
7207 71c03c7239fb [project @ 2007-11-28 02:19:25 by jwe] jwe parents: 7180 diff changeset	81 if (nargout > 0)
71c03c7239fb [project @ 2007-11-28 02:19:25 by jwe] jwe parents: 7180 diff changeset	82 retval = tmp;
934 be867ddbe4aa [project @ 1994-11-17 04:23:41 by jwe] jwe parents: 904 diff changeset	83 endif
7207 71c03c7239fb [project @ 2007-11-28 02:19:25 by jwe] jwe parents: 7180 diff changeset	84 unwind_protect_cleanup
71c03c7239fb [project @ 2007-11-28 02:19:25 by jwe] jwe parents: 7180 diff changeset	85 axes (oldh);
71c03c7239fb [project @ 2007-11-28 02:19:25 by jwe] jwe parents: 7180 diff changeset	86 end_unwind_protect
6302 a5cd8b77e892 [project @ 2007-02-13 08:08:33 by jwe] jwe parents: 6257 diff changeset	87
4 b4df021f796c [project @ 1993-08-08 01:26:08 by jwe] jwe parents: diff changeset	88 endfunction
9896 1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	89
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	90 function retval = __plr1__ (h, theta, fmt)
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	91
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	92 if (nargin != 3)
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	93 print_usage ();
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	94 endif
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	95
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	96 [nr, nc] = size (theta);
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	97 if (nr == 1)
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	98 theta = theta';
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	99 tmp = nr;
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	100 nr = nc;
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	101 nc = tmp;
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	102 endif
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	103 theta_i = imag (theta);
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	104 if (any (theta_i))
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	105 rho = theta_i;
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	106 theta = real (theta);
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	107 else
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	108 rho = theta;
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	109 theta = (1:nr)';
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	110 endif
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	111
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	112 retval = __plr2__ (h, theta, rho, fmt);
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	113
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	114 endfunction
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	115
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	116 function retval = __plr2__ (h, theta, rho, fmt)
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	117
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	118 if (nargin != 4)
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	119 print_usage ();
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	120 endif
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	121
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	122 if (any (imag (theta)))
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	123 theta = real (theta);
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	124 endif
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	125
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	126 if (any (imag (rho)))
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	127 rho = real (rho);
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	128 endif
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	129
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	130 if (isscalar (theta))
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	131 if (isscalar (rho))
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	132 x = rho * cos (theta);
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	133 y = rho * sin (theta);
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	134 retval = __plt__ ("polar", h, x, y, fmt);
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	135 else
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	136 error ("__plr2__: invalid data for plotting");
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	137 endif
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	138 elseif (isvector (theta))
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	139 if (isvector (rho))
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	140 if (length (theta) != length (rho))
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	141 error ("__plr2__: vector lengths must match");
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	142 endif
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	143 if (rows (rho) == 1)
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	144 rho = rho';
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	145 endif
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	146 if (rows (theta) == 1)
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	147 theta = theta';
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	148 endif
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	149 x = rho .* cos (theta);
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	150 y = rho .* sin (theta);
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	151 retval = __plt__ ("polar", h, x, y, fmt);
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	152 elseif (ismatrix (rho))
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	153 [t_nr, t_nc] = size (theta);
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	154 if (t_nr == 1)
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	155 theta = theta';
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	156 tmp = t_nr;
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	157 t_nr = t_nc;
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	158 t_nc = tmp;
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	159 endif
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	160 [r_nr, r_nc] = size (rho);
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	161 if (t_nr != r_nr)
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	162 rho = rho';
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	163 tmp = r_nr;
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	164 r_nr = r_nc;
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	165 r_nc = tmp;
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	166 endif
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	167 if (t_nr != r_nr)
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	168 error ("__plr2__: vector and matrix sizes must match");
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	169 endif
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	170 x = diag (cos (theta)) * rho;
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	171 y = diag (sin (theta)) * rho;
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	172 retval = __plt__ ("polar", h, x, y, fmt);
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	173 else
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	174 error ("__plr2__: invalid data for plotting");
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	175 endif
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	176 elseif (ismatrix (theta))
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	177 if (isvector (rho))
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	178 [r_nr, r_nc] = size (rho);
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	179 if (r_nr == 1)
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	180 rho = rho';
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	181 tmp = r_nr;
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	182 r_nr = r_nc;
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	183 r_nc = tmp;
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	184 endif
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	185 [t_nr, t_nc] = size (theta);
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	186 if (r_nr != t_nr)
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	187 theta = theta';
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	188 tmp = t_nr;
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	189 t_nr = t_nc;
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	190 t_nc = tmp;
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	191 endif
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	192 if (r_nr != t_nr)
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	193 error ("__plr2__: vector and matrix sizes must match");
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	194 endif
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	195 diag_r = diag (rho);
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	196 x = diag_r * cos (theta);
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	197 y = diag_r * sin (theta);
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	198 retval = __plt__ ("polar", h, x, y, fmt);
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	199 elseif (ismatrix (rho))
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	200 if (! size_equal (rho, theta))
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	201 error ("__plr2__: matrix dimensions must match");
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	202 endif
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	203 x = rho .* cos (theta);
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	204 y = rho .* sin (theta);
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	205 retval = __plt__ ("polar", h, x, y, fmt);
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	206 else
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	207 error ("__plr2__: invalid data for plotting");
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	208 endif
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	209 else
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	210 error ("__plr2__: invalid data for plotting");
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	211 endif
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	212
1aeb39118764 convert some plot functions to subfunctions or make some them private John W. Eaton <jwe@octave.org> parents: 7321 diff changeset	213 endfunction
12789 482cc053281b codesprint: Add 2 demos for polar.m Rik <octave@nomad.inbox5.com> parents: 11587 diff changeset	214
482cc053281b codesprint: Add 2 demos for polar.m Rik <octave@nomad.inbox5.com> parents: 11587 diff changeset	215
482cc053281b codesprint: Add 2 demos for polar.m Rik <octave@nomad.inbox5.com> parents: 11587 diff changeset	216 %!demo
482cc053281b codesprint: Add 2 demos for polar.m Rik <octave@nomad.inbox5.com> parents: 11587 diff changeset	217 %! theta = linspace (0, 2*pi, 1000);
482cc053281b codesprint: Add 2 demos for polar.m Rik <octave@nomad.inbox5.com> parents: 11587 diff changeset	218 %! rho = sin (7*theta);
482cc053281b codesprint: Add 2 demos for polar.m Rik <octave@nomad.inbox5.com> parents: 11587 diff changeset	219 %! polar (theta, rho);
482cc053281b codesprint: Add 2 demos for polar.m Rik <octave@nomad.inbox5.com> parents: 11587 diff changeset	220
482cc053281b codesprint: Add 2 demos for polar.m Rik <octave@nomad.inbox5.com> parents: 11587 diff changeset	221 %!demo
482cc053281b codesprint: Add 2 demos for polar.m Rik <octave@nomad.inbox5.com> parents: 11587 diff changeset	222 %! theta = linspace (0, 10*pi, 1000);
482cc053281b codesprint: Add 2 demos for polar.m Rik <octave@nomad.inbox5.com> parents: 11587 diff changeset	223 %! rho = sin (5/4*theta);
482cc053281b codesprint: Add 2 demos for polar.m Rik <octave@nomad.inbox5.com> parents: 11587 diff changeset	224 %! polar (theta, rho);
482cc053281b codesprint: Add 2 demos for polar.m Rik <octave@nomad.inbox5.com> parents: 11587 diff changeset	225

Mercurial > octave-nkf

annotate scripts/plot/polar.m @ 12789:482cc053281b