annotate scripts/plot/contour.m @ 7017:a1dbe9d80eee

[project @ 2007-10-12 21:27:11 by jwe]
author jwe
date Fri, 12 Oct 2007 21:27:37 +0000
parents 93c65f2a5668
children 9f38c6293317
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
7017
a1dbe9d80eee [project @ 2007-10-12 21:27:11 by jwe]
jwe
parents: 7016
diff changeset
1 ## Copyright (C) 1993, 1994, 1995, 1996, 1997, 1998, 1999, 2000, 2002,
a1dbe9d80eee [project @ 2007-10-12 21:27:11 by jwe]
jwe
parents: 7016
diff changeset
2 ## 2003, 2004, 2005, 2006, 2007 Shai Ayal
2313
5ca126254d15 [project @ 1996-07-11 21:25:22 by jwe]
jwe
parents: 2311
diff changeset
3 ##
6440
98ee80702bca [project @ 2007-03-23 15:13:19 by jwe]
jwe
parents: 6434
diff changeset
4 ## This file is part of Octave.
98ee80702bca [project @ 2007-03-23 15:13:19 by jwe]
jwe
parents: 6434
diff changeset
5 ##
98ee80702bca [project @ 2007-03-23 15:13:19 by jwe]
jwe
parents: 6434
diff changeset
6 ## Octave is free software; you can redistribute it and/or modify it
2313
5ca126254d15 [project @ 1996-07-11 21:25:22 by jwe]
jwe
parents: 2311
diff changeset
7 ## 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
8 ## 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
9 ## your option) any later version.
2313
5ca126254d15 [project @ 1996-07-11 21:25:22 by jwe]
jwe
parents: 2311
diff changeset
10 ##
6440
98ee80702bca [project @ 2007-03-23 15:13:19 by jwe]
jwe
parents: 6434
diff changeset
11 ## Octave is distributed in the hope that it will be useful, but
2313
5ca126254d15 [project @ 1996-07-11 21:25:22 by jwe]
jwe
parents: 2311
diff changeset
12 ## WITHOUT ANY WARRANTY; without even the implied warranty of
5ca126254d15 [project @ 1996-07-11 21:25:22 by jwe]
jwe
parents: 2311
diff changeset
13 ## 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
14 ## General Public License for more details.
5ca126254d15 [project @ 1996-07-11 21:25:22 by jwe]
jwe
parents: 2311
diff changeset
15 ##
5ca126254d15 [project @ 1996-07-11 21:25:22 by jwe]
jwe
parents: 2311
diff changeset
16 ## 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
17 ## 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
18 ## <http://www.gnu.org/licenses/>.
245
16a24e76d6e0 [project @ 1993-12-03 02:00:15 by jwe]
jwe
parents: 4
diff changeset
19
3368
a4cd1e9d9962 [project @ 1999-11-20 17:22:48 by jwe]
jwe
parents: 3136
diff changeset
20 ## -*- texinfo -*-
6888
b3d286e8a243 [project @ 2007-09-12 06:40:58 by jwe]
jwe
parents: 6604
diff changeset
21 ## @deftypefn {Function File} {@var{c} =} contour (@var{z})
b3d286e8a243 [project @ 2007-09-12 06:40:58 by jwe]
jwe
parents: 6604
diff changeset
22 ## @deftypefnx {Function File} {@var{c} =} contour (@var{z}, @var{vn})
b3d286e8a243 [project @ 2007-09-12 06:40:58 by jwe]
jwe
parents: 6604
diff changeset
23 ## @deftypefnx {Function File} {@var{c} =} contour (@var{x}, @var{y}, @var{z})
b3d286e8a243 [project @ 2007-09-12 06:40:58 by jwe]
jwe
parents: 6604
diff changeset
24 ## @deftypefnx {Function File} {@var{c} =} contour (@var{x}, @var{y}, @var{z}, @var{vn})
6604
c2b700808e0a [project @ 2007-05-02 05:36:43 by jwe]
jwe
parents: 6547
diff changeset
25 ## Plot level curves (contour lines) of the matrix @var{z}, using the
c2b700808e0a [project @ 2007-05-02 05:36:43 by jwe]
jwe
parents: 6547
diff changeset
26 ## contour matrix @var{c} computed by @code{contourc} from the same
6895
76e3d985ae56 [project @ 2007-09-13 18:22:38 by jwe]
jwe
parents: 6888
diff changeset
27 ## arguments; see the latter for their interpretation. The set of
76e3d985ae56 [project @ 2007-09-13 18:22:38 by jwe]
jwe
parents: 6888
diff changeset
28 ## contour levels, @var{c}, is only returned if requested. For example:
6257
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents: 6172
diff changeset
29 ##
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents: 6172
diff changeset
30 ## @example
6604
c2b700808e0a [project @ 2007-05-02 05:36:43 by jwe]
jwe
parents: 6547
diff changeset
31 ## @group
c2b700808e0a [project @ 2007-05-02 05:36:43 by jwe]
jwe
parents: 6547
diff changeset
32 ## x = 0:2;
c2b700808e0a [project @ 2007-05-02 05:36:43 by jwe]
jwe
parents: 6547
diff changeset
33 ## y = x;
c2b700808e0a [project @ 2007-05-02 05:36:43 by jwe]
jwe
parents: 6547
diff changeset
34 ## z = x' * y;
c2b700808e0a [project @ 2007-05-02 05:36:43 by jwe]
jwe
parents: 6547
diff changeset
35 ## contour (x, y, z, 2:3)
6257
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents: 6172
diff changeset
36 ##
6604
c2b700808e0a [project @ 2007-05-02 05:36:43 by jwe]
jwe
parents: 6547
diff changeset
37 ## @end group
6257
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents: 6172
diff changeset
38 ## @end example
6547
4fb053f24fd6 [project @ 2007-04-19 21:47:40 by jwe]
jwe
parents: 6440
diff changeset
39 ## @seealso{contourc, line, plot}
3368
a4cd1e9d9962 [project @ 1999-11-20 17:22:48 by jwe]
jwe
parents: 3136
diff changeset
40 ## @end deftypefn
4
b4df021f796c [project @ 1993-08-08 01:26:08 by jwe]
jwe
parents:
diff changeset
41
6257
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents: 6172
diff changeset
42 ## Author: shaia
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents: 6172
diff changeset
43
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents: 6172
diff changeset
44 function retval = contour (varargin)
4
b4df021f796c [project @ 1993-08-08 01:26:08 by jwe]
jwe
parents:
diff changeset
45
6257
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents: 6172
diff changeset
46 [c, lev] = contourc (varargin{:});
6172
2de853a110df [project @ 2006-11-15 22:53:47 by jwe]
jwe
parents: 6157
diff changeset
47
6257
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents: 6172
diff changeset
48 cmap = get (gcf(), "colormap");
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents: 6172
diff changeset
49
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents: 6172
diff changeset
50 levx = linspace (min (lev), max (lev), size (cmap, 1));
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents: 6172
diff changeset
51
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents: 6172
diff changeset
52 newplot ();
3063
4287b8c06fbf [project @ 1997-06-26 02:36:05 by jwe]
jwe
parents: 2847
diff changeset
53
6434
7caa5942d4da [project @ 2007-03-23 13:59:29 by jwe]
jwe
parents: 6279
diff changeset
54 ## Decode contourc output format.
6257
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents: 6172
diff changeset
55 i1 = 1;
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents: 6172
diff changeset
56 while (i1 < length (c))
6434
7caa5942d4da [project @ 2007-03-23 13:59:29 by jwe]
jwe
parents: 6279
diff changeset
57
6257
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents: 6172
diff changeset
58 clev = c(1,i1);
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents: 6172
diff changeset
59 clen = c(2,i1);
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents: 6172
diff changeset
60
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents: 6172
diff changeset
61 ccr = interp1 (levx, cmap(:,1), clev);
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents: 6172
diff changeset
62 ccg = interp1 (levx, cmap(:,2), clev);
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents: 6172
diff changeset
63 ccb = interp1 (levx, cmap(:,3), clev);
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents: 6172
diff changeset
64
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents: 6172
diff changeset
65 ii = i1+1:i1+clen;
6279
195382e92d51 [project @ 2007-02-08 03:47:19 by jwe]
jwe
parents: 6257
diff changeset
66 line (c(1,ii), c(2,ii), "color", [ccr, ccg, ccb]);
6257
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents: 6172
diff changeset
67
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents: 6172
diff changeset
68 i1 += c(2,i1)+1;
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents: 6172
diff changeset
69 endwhile
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents: 6172
diff changeset
70
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents: 6172
diff changeset
71 if (nargout > 0)
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents: 6172
diff changeset
72 retval = c;
4
b4df021f796c [project @ 1993-08-08 01:26:08 by jwe]
jwe
parents:
diff changeset
73 endif
b4df021f796c [project @ 1993-08-08 01:26:08 by jwe]
jwe
parents:
diff changeset
74
b4df021f796c [project @ 1993-08-08 01:26:08 by jwe]
jwe
parents:
diff changeset
75 endfunction
6257
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents: 6172
diff changeset
76
44c91c5dfe1d [project @ 2007-01-30 19:16:52 by jwe]
jwe
parents: 6172
diff changeset
77