Mercurial > hg > octave-nkf

annotate scripts/polynomial/polyout.m @ 12541:dd2c70b30f28

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add tests for ifftshift.m
author Robert T. Short <octave@phaselockedsystems.com.com>
date Sat, 26 Mar 2011 06:50:12 -0700
parents b0084095098e
children a44ec6d33cc9
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
11523
fd0a3ac60b0e update copyright notices
John W. Eaton <jwe@octave.org>
parents: 10821
diff changeset
1 ## Copyright (C) 1995-2011 Auburn University. All rights reserved.
3424
61e40232a4e8 [project @ 2000-01-13 07:19:43 by jwe]
jwe
parents:
diff changeset
2 ##
61e40232a4e8 [project @ 2000-01-13 07:19:43 by jwe]
jwe
parents:
diff changeset
3 ## This file is part of Octave.
61e40232a4e8 [project @ 2000-01-13 07:19:43 by jwe]
jwe
parents:
diff changeset
4 ##
61e40232a4e8 [project @ 2000-01-13 07:19:43 by jwe]
jwe
parents:
diff changeset
5 ## Octave is free software; you can redistribute it and/or modify it
7016
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe]
jwe
parents: 6046
diff changeset
6 ## under the terms of the GNU General Public License as published by
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe]
jwe
parents: 6046
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: 6046
diff changeset
8 ## your option) any later version.
3424
61e40232a4e8 [project @ 2000-01-13 07:19:43 by jwe]
jwe
parents:
diff changeset
9 ##
7016
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe]
jwe
parents: 6046
diff changeset
10 ## Octave is distributed in the hope that it will be useful, but
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe]
jwe
parents: 6046
diff changeset
11 ## WITHOUT ANY WARRANTY; without even the implied warranty of
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe]
jwe
parents: 6046
diff changeset
12 ## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe]
jwe
parents: 6046
diff changeset
13 ## General Public License for more details.
3424
61e40232a4e8 [project @ 2000-01-13 07:19:43 by jwe]
jwe
parents:
diff changeset
14 ##
61e40232a4e8 [project @ 2000-01-13 07:19:43 by jwe]
jwe
parents:
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: 6046
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: 6046
diff changeset
17 ## <http://www.gnu.org/licenses/>.
3424
61e40232a4e8 [project @ 2000-01-13 07:19:43 by jwe]
jwe
parents:
diff changeset
18
61e40232a4e8 [project @ 2000-01-13 07:19:43 by jwe]
jwe
parents:
diff changeset
19 ## -*- texinfo -*-
10687
a8ce6bdecce5 Improve documentation strings.
Rik <octave@nomad.inbox5.com>
parents: 10549
diff changeset
20 ## @deftypefn {Function File} {} polyout (@var{c})
a8ce6bdecce5 Improve documentation strings.
Rik <octave@nomad.inbox5.com>
parents: 10549
diff changeset
21 ## @deftypefnx {Function File} {} polyout (@var{c}, @var{x})
a8ce6bdecce5 Improve documentation strings.
Rik <octave@nomad.inbox5.com>
parents: 10549
diff changeset
22 ## @deftypefnx {Function File} {@var{str} =} polyout (@dots{})
3457
e031284eea27 [project @ 2000-01-19 08:49:56 by jwe]
jwe
parents: 3426
diff changeset
23 ## Write formatted polynomial
5016
bdbee5282954 [project @ 2004-09-22 02:50:35 by jwe]
jwe
parents: 4898
diff changeset
24 ## @tex
bdbee5282954 [project @ 2004-09-22 02:50:35 by jwe]
jwe
parents: 4898
diff changeset
25 ## $$ c(x) = c_1 x^n + \ldots + c_n x + c_{n+1} $$
bdbee5282954 [project @ 2004-09-22 02:50:35 by jwe]
jwe
parents: 4898
diff changeset
26 ## @end tex
8517
81d6ab3ac93c Allow documentation tobe built for other formats than tex and info
sh@sh-laptop
parents: 8286
diff changeset
27 ## @ifnottex
10821
693e22af08ae Grammarcheck documentation of m-files
Rik <octave@nomad.inbox5.com>
parents: 10687
diff changeset
28 ##
3424
61e40232a4e8 [project @ 2000-01-13 07:19:43 by jwe]
jwe
parents:
diff changeset
29 ## @example
9051
1bf0ce0930be Grammar check TexInfo in all .m files
Rik <rdrider0-list@yahoo.com>
parents: 8920
diff changeset
30 ## c(x) = c(1) * x^n + @dots{} + c(n) x + c(n+1)
3424
61e40232a4e8 [project @ 2000-01-13 07:19:43 by jwe]
jwe
parents:
diff changeset
31 ## @end example
10821
693e22af08ae Grammarcheck documentation of m-files
Rik <octave@nomad.inbox5.com>
parents: 10687
diff changeset
32 ##
8517
81d6ab3ac93c Allow documentation tobe built for other formats than tex and info
sh@sh-laptop
parents: 8286
diff changeset
33 ## @end ifnottex
10687
a8ce6bdecce5 Improve documentation strings.
Rik <octave@nomad.inbox5.com>
parents: 10549
diff changeset
34 ## and return it as a string or write it to the screen (if
a8ce6bdecce5 Improve documentation strings.
Rik <octave@nomad.inbox5.com>
parents: 10549
diff changeset
35 ## @var{nargout} is zero).
a8ce6bdecce5 Improve documentation strings.
Rik <octave@nomad.inbox5.com>
parents: 10549
diff changeset
36 ## @var{x} defaults to the string @code{"s"}.
3426
f8dde1807dee [project @ 2000-01-13 08:40:00 by jwe]
jwe
parents: 3424
diff changeset
37 ## @seealso{polyval, polyvalm, poly, roots, conv, deconv, residue,
10224
f6e0404421f4 point to polyint in @seealso, not polyinteg
John W. Eaton <jwe@octave.org>
parents: 9211
diff changeset
38 ## filter, polyderiv, polyint}
5642
2618a0750ae6 [project @ 2006-03-06 21:26:48 by jwe]
jwe
parents: 5443
diff changeset
39 ## @end deftypefn
3424
61e40232a4e8 [project @ 2000-01-13 07:19:43 by jwe]
jwe
parents:
diff changeset
40
61e40232a4e8 [project @ 2000-01-13 07:19:43 by jwe]
jwe
parents:
diff changeset
41 ## Author: A. S. Hodel <a.s.hodel@eng.auburn.edu>
61e40232a4e8 [project @ 2000-01-13 07:19:43 by jwe]
jwe
parents:
diff changeset
42 ## Created: May 1995
61e40232a4e8 [project @ 2000-01-13 07:19:43 by jwe]
jwe
parents:
diff changeset
43 ## Nov 1998: Correctly handles complex coefficients
3426
f8dde1807dee [project @ 2000-01-13 08:40:00 by jwe]
jwe
parents: 3424
diff changeset
44
3424
61e40232a4e8 [project @ 2000-01-13 07:19:43 by jwe]
jwe
parents:
diff changeset
45 function y = polyout (c, x)
61e40232a4e8 [project @ 2000-01-13 07:19:43 by jwe]
jwe
parents:
diff changeset
46
3457
e031284eea27 [project @ 2000-01-19 08:49:56 by jwe]
jwe
parents: 3426
diff changeset
47 if (nargin < 1) || (nargin > 2) || (nargout < 0) || (nargout > 1)
6046
34f96dd5441b [project @ 2006-10-10 16:10:25 by jwe]
jwe
parents: 5642
diff changeset
48 print_usage ();
3424
61e40232a4e8 [project @ 2000-01-13 07:19:43 by jwe]
jwe
parents:
diff changeset
49 endif
61e40232a4e8 [project @ 2000-01-13 07:19:43 by jwe]
jwe
parents:
diff changeset
50
4030
22bd65326ec1 [project @ 2002-08-09 18:58:13 by jwe]
jwe
parents: 3457
diff changeset
51 if (! isvector (c))
8685
983ac67dc3d4 polyout.m: additional style fixes
John W. Eaton <jwe@octave.org>
parents: 8684
diff changeset
52 error ("polyout: first argument must be a vector");
3424
61e40232a4e8 [project @ 2000-01-13 07:19:43 by jwe]
jwe
parents:
diff changeset
53 endif
3426
f8dde1807dee [project @ 2000-01-13 08:40:00 by jwe]
jwe
parents: 3424
diff changeset
54
3424
61e40232a4e8 [project @ 2000-01-13 07:19:43 by jwe]
jwe
parents:
diff changeset
55 if (nargin == 1)
61e40232a4e8 [project @ 2000-01-13 07:19:43 by jwe]
jwe
parents:
diff changeset
56 x = "s";
8684
6d15bc6c4c15 polyout.m: Replace com2str with num2str.
Thomas D. Dean <tomdean@speakeasy.org>
parents: 8517
diff changeset
57 elseif (! ischar (x))
8685
983ac67dc3d4 polyout.m: additional style fixes
John W. Eaton <jwe@octave.org>
parents: 8684
diff changeset
58 error ("polyout: second argument must be a string");
3424
61e40232a4e8 [project @ 2000-01-13 07:19:43 by jwe]
jwe
parents:
diff changeset
59 endif
61e40232a4e8 [project @ 2000-01-13 07:19:43 by jwe]
jwe
parents:
diff changeset
60
8684
6d15bc6c4c15 polyout.m: Replace com2str with num2str.
Thomas D. Dean <tomdean@speakeasy.org>
parents: 8517
diff changeset
61 n = length (c);
3424
61e40232a4e8 [project @ 2000-01-13 07:19:43 by jwe]
jwe
parents:
diff changeset
62 if(n > 0)
61e40232a4e8 [project @ 2000-01-13 07:19:43 by jwe]
jwe
parents:
diff changeset
63 n1 = n+1;
61e40232a4e8 [project @ 2000-01-13 07:19:43 by jwe]
jwe
parents:
diff changeset
64
4898
8fd9495f5054 [project @ 2004-06-04 13:29:33 by jwe]
jwe
parents: 4030
diff changeset
65 tmp = coeff (c(1));
3457
e031284eea27 [project @ 2000-01-19 08:49:56 by jwe]
jwe
parents: 3426
diff changeset
66 for ii = 2:n
e031284eea27 [project @ 2000-01-19 08:49:56 by jwe]
jwe
parents: 3426
diff changeset
67 if (real (c(ii)) < 0)
10549
95c3e38098bf Untabify .m scripts
Rik <code@nomad.inbox5.com>
parents: 10224
diff changeset
68 ns = " - ";
95c3e38098bf Untabify .m scripts
Rik <code@nomad.inbox5.com>
parents: 10224
diff changeset
69 c(ii) = -c(ii);
3457
e031284eea27 [project @ 2000-01-19 08:49:56 by jwe]
jwe
parents: 3426
diff changeset
70 else
e031284eea27 [project @ 2000-01-19 08:49:56 by jwe]
jwe
parents: 3426
diff changeset
71 ns = " + ";
e031284eea27 [project @ 2000-01-19 08:49:56 by jwe]
jwe
parents: 3426
diff changeset
72 endif
3424
61e40232a4e8 [project @ 2000-01-13 07:19:43 by jwe]
jwe
parents:
diff changeset
73
4898
8fd9495f5054 [project @ 2004-06-04 13:29:33 by jwe]
jwe
parents: 4030
diff changeset
74 tmp = sprintf ("%s*%s^%d%s%s", tmp, x, n1-ii, ns, coeff (c(ii)));
3426
f8dde1807dee [project @ 2000-01-13 08:40:00 by jwe]
jwe
parents: 3424
diff changeset
75
3424
61e40232a4e8 [project @ 2000-01-13 07:19:43 by jwe]
jwe
parents:
diff changeset
76 endfor
61e40232a4e8 [project @ 2000-01-13 07:19:43 by jwe]
jwe
parents:
diff changeset
77 else
61e40232a4e8 [project @ 2000-01-13 07:19:43 by jwe]
jwe
parents:
diff changeset
78 tmp = " ";
61e40232a4e8 [project @ 2000-01-13 07:19:43 by jwe]
jwe
parents:
diff changeset
79 endif
61e40232a4e8 [project @ 2000-01-13 07:19:43 by jwe]
jwe
parents:
diff changeset
80
61e40232a4e8 [project @ 2000-01-13 07:19:43 by jwe]
jwe
parents:
diff changeset
81 if(nargout == 0)
11589
b0084095098e missing semicolons in script files
John W. Eaton <jwe@octave.org>
parents: 11523
diff changeset
82 disp (tmp);
3424
61e40232a4e8 [project @ 2000-01-13 07:19:43 by jwe]
jwe
parents:
diff changeset
83 else
61e40232a4e8 [project @ 2000-01-13 07:19:43 by jwe]
jwe
parents:
diff changeset
84 y = tmp;
61e40232a4e8 [project @ 2000-01-13 07:19:43 by jwe]
jwe
parents:
diff changeset
85 endif
61e40232a4e8 [project @ 2000-01-13 07:19:43 by jwe]
jwe
parents:
diff changeset
86
61e40232a4e8 [project @ 2000-01-13 07:19:43 by jwe]
jwe
parents:
diff changeset
87 endfunction
4898
8fd9495f5054 [project @ 2004-06-04 13:29:33 by jwe]
jwe
parents: 4030
diff changeset
88
8684
6d15bc6c4c15 polyout.m: Replace com2str with num2str.
Thomas D. Dean <tomdean@speakeasy.org>
parents: 8517
diff changeset
89 function str = coeff (c)
4898
8fd9495f5054 [project @ 2004-06-04 13:29:33 by jwe]
jwe
parents: 4030
diff changeset
90 if (imag (c))
8fd9495f5054 [project @ 2004-06-04 13:29:33 by jwe]
jwe
parents: 4030
diff changeset
91 if (real (c))
8684
6d15bc6c4c15 polyout.m: Replace com2str with num2str.
Thomas D. Dean <tomdean@speakeasy.org>
parents: 8517
diff changeset
92 str = sprintf ("(%s)", num2str (c, 5));
4898
8fd9495f5054 [project @ 2004-06-04 13:29:33 by jwe]
jwe
parents: 4030
diff changeset
93 else
8684
6d15bc6c4c15 polyout.m: Replace com2str with num2str.
Thomas D. Dean <tomdean@speakeasy.org>
parents: 8517
diff changeset
94 str = num2str (c, 5);
4898
8fd9495f5054 [project @ 2004-06-04 13:29:33 by jwe]
jwe
parents: 4030
diff changeset
95 endif
8fd9495f5054 [project @ 2004-06-04 13:29:33 by jwe]
jwe
parents: 4030
diff changeset
96 else
8684
6d15bc6c4c15 polyout.m: Replace com2str with num2str.
Thomas D. Dean <tomdean@speakeasy.org>
parents: 8517
diff changeset
97 str = num2str (c, 5);
4898
8fd9495f5054 [project @ 2004-06-04 13:29:33 by jwe]
jwe
parents: 4030
diff changeset
98 endif
8684
6d15bc6c4c15 polyout.m: Replace com2str with num2str.
Thomas D. Dean <tomdean@speakeasy.org>
parents: 8517
diff changeset
99 endfunction