annotate scripts/special-matrix/toeplitz.m @ 14868:5d3a684236b0

maint: Use Octave coding conventions for cuddling parentheses in scripts directory * lin2mu.m, loadaudio.m, wavread.m, accumarray.m, bicubic.m, celldisp.m, colon.m, cplxpair.m, dblquad.m, divergence.m, genvarname.m, gradient.m, int2str.m, interp1.m, interp1q.m, interp2.m, interpn.m, loadobj.m, nthargout.m, __isequal__.m, __splinen__.m, quadgk.m, quadl.m, quadv.m, rat.m, rot90.m, rotdim.m, saveobj.m, subsindex.m, triplequad.m, delaunay3.m, griddata.m, inpolygon.m, tsearchn.m, voronoi.m, get_first_help_sentence.m, which.m, gray2ind.m, pink.m, dlmwrite.m, strread.m, textread.m, textscan.m, housh.m, ishermitian.m, issymmetric.m, krylov.m, logm.m, null.m, rref.m, compare_versions.m, copyfile.m, dump_prefs.m, edit.m, fileparts.m, getappdata.m, isappdata.m, movefile.m, orderfields.m, parseparams.m, __xzip__.m, rmappdata.m, setappdata.m, swapbytes.m, unpack.m, ver.m, fminbnd.m, fminunc.m, fsolve.m, glpk.m, lsqnonneg.m, qp.m, sqp.m, configure_make.m, copy_files.m, describe.m, get_description.m, get_forge_pkg.m, install.m, installed_packages.m, is_architecture_dependent.m, load_package_dirs.m, print_package_description.m, rebuild.m, repackage.m, save_order.m, shell.m, allchild.m, ancestor.m, area.m, axes.m, axis.m, clabel.m, close.m, colorbar.m, comet.m, comet3.m, contour.m, cylinder.m, ezmesh.m, ezsurf.m, findobj.m, fplot.m, hist.m, isocolors.m, isonormals.m, isosurface.m, isprop.m, legend.m, mesh.m, meshz.m, pareto.m, pcolor.m, peaks.m, plot3.m, plotmatrix.m, plotyy.m, polar.m, print.m, __add_datasource__.m, __add_default_menu__.m, __axes_limits__.m, __bar__.m, __clabel__.m, __contour__.m, __errcomm__.m, __errplot__.m, __ezplot__.m, __file_filter__.m, __fltk_print__.m, __ghostscript__.m, __gnuplot_print__.m, __go_draw_axes__.m, __go_draw_figure__.m, __interp_cube__.m, __marching_cube__.m, __patch__.m, __pie__.m, __plt__.m, __print_parse_opts__.m, __quiver__.m, __scatter__.m, __stem__.m, __tight_eps_bbox__.m, __uigetdir_fltk__.m, __uigetfile_fltk__.m, __uiputfile_fltk__.m, quiver.m, quiver3.m, rectangle.m, refreshdata.m, ribbon.m, scatter.m, semilogy.m, shading.m, slice.m, subplot.m, surface.m, surfl.m, surfnorm.m, text.m, uigetfile.m, uiputfile.m, whitebg.m, deconv.m, mkpp.m, pchip.m, polyaffine.m, polyder.m, polygcd.m, polyout.m, polyval.m, ppint.m, ppjumps.m, ppval.m, residue.m, roots.m, spline.m, splinefit.m, addpref.m, getpref.m, setpref.m, ismember.m, setxor.m, arch_fit.m, arch_rnd.m, arch_test.m, autoreg_matrix.m, diffpara.m, fftconv.m, filter2.m, hanning.m, hurst.m, periodogram.m, triangle_sw.m, sinc.m, spectral_xdf.m, spencer.m, stft.m, synthesis.m, unwrap.m, yulewalker.m, bicgstab.m, gmres.m, pcg.m, pcr.m, __sprand_impl__.m, speye.m, spfun.m, sprandn.m, spstats.m, svds.m, treelayout.m, treeplot.m, bessel.m, factor.m, legendre.m, perms.m, primes.m, magic.m, toeplitz.m, corr.m, cov.m, mean.m, median.m, mode.m, qqplot.m, quantile.m, ranks.m, zscore.m, logistic_regression_likelihood.m, bartlett_test.m, chisquare_test_homogeneity.m, chisquare_test_independence.m, kolmogorov_smirnov_test.m, run_test.m, u_test.m, wilcoxon_test.m, z_test.m, z_test_2.m, bin2dec.m, dec2base.m, mat2str.m, strcat.m, strchr.m, strjust.m, strtok.m, substr.m, untabify.m, assert.m, demo.m, example.m, fail.m, speed.m, test.m, now.m: Use Octave coding conventions for cuddling parentheses in scripts directory.
author Rik <octave@nomad.inbox5.com>
date Tue, 17 Jul 2012 07:08:39 -0700
parents 72c96de7a403
children b81b9d079515
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
14138
72c96de7a403 maint: update copyright notices for 2012
John W. Eaton <jwe@octave.org>
parents: 13889
diff changeset
1 ## Copyright (C) 1993-2012 John W. Eaton
9092
38c3a0f8c6d0 optimize toeplitz
Jaroslav Hajek <highegg@gmail.com>
parents: 9041
diff changeset
2 ## Copyright (C) 2009 VZLU Prague
2313
5ca126254d15 [project @ 1996-07-11 21:25:22 by jwe]
jwe
parents: 2311
diff changeset
3 ##
5ca126254d15 [project @ 1996-07-11 21:25:22 by jwe]
jwe
parents: 2311
diff changeset
4 ## This file is part of Octave.
5ca126254d15 [project @ 1996-07-11 21:25:22 by jwe]
jwe
parents: 2311
diff changeset
5 ##
5ca126254d15 [project @ 1996-07-11 21:25:22 by jwe]
jwe
parents: 2311
diff changeset
6 ## 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
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: 6046
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: 6046
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 ##
5ca126254d15 [project @ 1996-07-11 21:25:22 by jwe]
jwe
parents: 2311
diff changeset
11 ## 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
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: 6046
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: 6046
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
3369
f37ca3017116 [project @ 1999-11-21 16:26:02 by jwe]
jwe
parents: 3140
diff changeset
20 ## -*- texinfo -*-
10791
3140cb7a05a1 Add spellchecker scripts for Octave and run spellcheck of documentation
Rik <octave@nomad.inbox5.com>
parents: 10549
diff changeset
21 ## @deftypefn {Function File} {} toeplitz (@var{c})
3140cb7a05a1 Add spellchecker scripts for Octave and run spellcheck of documentation
Rik <octave@nomad.inbox5.com>
parents: 10549
diff changeset
22 ## @deftypefnx {Function File} {} toeplitz (@var{c}, @var{r})
12639
4d777e05d47c doc: Review and update documentation for "Matrix Manipulation" chapter.
Rik <octave@nomad.inbox5.com>
parents: 11587
diff changeset
23 ## Return the Toeplitz matrix constructed from the first column @var{c},
4d777e05d47c doc: Review and update documentation for "Matrix Manipulation" chapter.
Rik <octave@nomad.inbox5.com>
parents: 11587
diff changeset
24 ## and (optionally) the first row @var{r}. If the first element of @var{r}
4d777e05d47c doc: Review and update documentation for "Matrix Manipulation" chapter.
Rik <octave@nomad.inbox5.com>
parents: 11587
diff changeset
25 ## is not the same as the first element of @var{c}, the first element of
3369
f37ca3017116 [project @ 1999-11-21 16:26:02 by jwe]
jwe
parents: 3140
diff changeset
26 ## @var{c} is used. If the second argument is omitted, the first row is
f37ca3017116 [project @ 1999-11-21 16:26:02 by jwe]
jwe
parents: 3140
diff changeset
27 ## taken to be the same as the first column.
3426
f8dde1807dee [project @ 2000-01-13 08:40:00 by jwe]
jwe
parents: 3408
diff changeset
28 ##
5016
bdbee5282954 [project @ 2004-09-22 02:50:35 by jwe]
jwe
parents: 3426
diff changeset
29 ## A square Toeplitz matrix has the form:
3369
f37ca3017116 [project @ 1999-11-21 16:26:02 by jwe]
jwe
parents: 3140
diff changeset
30 ## @tex
f37ca3017116 [project @ 1999-11-21 16:26:02 by jwe]
jwe
parents: 3140
diff changeset
31 ## $$
5016
bdbee5282954 [project @ 2004-09-22 02:50:35 by jwe]
jwe
parents: 3426
diff changeset
32 ## \left[\matrix{c_0 & r_1 & r_2 & \cdots & r_n\cr
bdbee5282954 [project @ 2004-09-22 02:50:35 by jwe]
jwe
parents: 3426
diff changeset
33 ## c_1 & c_0 & r_1 & \cdots & r_{n-1}\cr
bdbee5282954 [project @ 2004-09-22 02:50:35 by jwe]
jwe
parents: 3426
diff changeset
34 ## c_2 & c_1 & c_0 & \cdots & r_{n-2}\cr
bdbee5282954 [project @ 2004-09-22 02:50:35 by jwe]
jwe
parents: 3426
diff changeset
35 ## \vdots & \vdots & \vdots & \ddots & \vdots\cr
bdbee5282954 [project @ 2004-09-22 02:50:35 by jwe]
jwe
parents: 3426
diff changeset
36 ## c_n & c_{n-1} & c_{n-2} & \ldots & c_0}\right]
3369
f37ca3017116 [project @ 1999-11-21 16:26:02 by jwe]
jwe
parents: 3140
diff changeset
37 ## $$
f37ca3017116 [project @ 1999-11-21 16:26:02 by jwe]
jwe
parents: 3140
diff changeset
38 ## @end tex
8517
81d6ab3ac93c Allow documentation tobe built for other formats than tex and info
sh@sh-laptop
parents: 7411
diff changeset
39 ## @ifnottex
3426
f8dde1807dee [project @ 2000-01-13 08:40:00 by jwe]
jwe
parents: 3408
diff changeset
40 ##
3369
f37ca3017116 [project @ 1999-11-21 16:26:02 by jwe]
jwe
parents: 3140
diff changeset
41 ## @example
f37ca3017116 [project @ 1999-11-21 16:26:02 by jwe]
jwe
parents: 3140
diff changeset
42 ## @group
9041
853f96e8008f Cleanup documentation file matrix.texi
Rik <rdrider0-list@yahoo.com>
parents: 8920
diff changeset
43 ## c(0) r(1) r(2) @dots{} r(n)
853f96e8008f Cleanup documentation file matrix.texi
Rik <rdrider0-list@yahoo.com>
parents: 8920
diff changeset
44 ## c(1) c(0) r(1) @dots{} r(n-1)
853f96e8008f Cleanup documentation file matrix.texi
Rik <rdrider0-list@yahoo.com>
parents: 8920
diff changeset
45 ## c(2) c(1) c(0) @dots{} r(n-2)
10821
693e22af08ae Grammarcheck documentation of m-files
Rik <octave@nomad.inbox5.com>
parents: 10791
diff changeset
46 ## . . . . .
693e22af08ae Grammarcheck documentation of m-files
Rik <octave@nomad.inbox5.com>
parents: 10791
diff changeset
47 ## . . . . .
693e22af08ae Grammarcheck documentation of m-files
Rik <octave@nomad.inbox5.com>
parents: 10791
diff changeset
48 ## . . . . .
9041
853f96e8008f Cleanup documentation file matrix.texi
Rik <rdrider0-list@yahoo.com>
parents: 8920
diff changeset
49 ## c(n) c(n-1) c(n-2) @dots{} c(0)
3369
f37ca3017116 [project @ 1999-11-21 16:26:02 by jwe]
jwe
parents: 3140
diff changeset
50 ## @end group
f37ca3017116 [project @ 1999-11-21 16:26:02 by jwe]
jwe
parents: 3140
diff changeset
51 ## @end example
10821
693e22af08ae Grammarcheck documentation of m-files
Rik <octave@nomad.inbox5.com>
parents: 10791
diff changeset
52 ##
8517
81d6ab3ac93c Allow documentation tobe built for other formats than tex and info
sh@sh-laptop
parents: 7411
diff changeset
53 ## @end ifnottex
12639
4d777e05d47c doc: Review and update documentation for "Matrix Manipulation" chapter.
Rik <octave@nomad.inbox5.com>
parents: 11587
diff changeset
54 ## @seealso{hankel}
3369
f37ca3017116 [project @ 1999-11-21 16:26:02 by jwe]
jwe
parents: 3140
diff changeset
55 ## @end deftypefn
4
b4df021f796c [project @ 1993-08-08 01:26:08 by jwe]
jwe
parents:
diff changeset
56
9092
38c3a0f8c6d0 optimize toeplitz
Jaroslav Hajek <highegg@gmail.com>
parents: 9041
diff changeset
57 ## Author: jwe && jh
2314
949ab8eba8bc [project @ 1996-07-12 03:58:02 by jwe]
jwe
parents: 2313
diff changeset
58
2311
2b5788792cad [project @ 1996-07-11 20:18:38 by jwe]
jwe
parents: 2303
diff changeset
59 function retval = toeplitz (c, r)
4
b4df021f796c [project @ 1993-08-08 01:26:08 by jwe]
jwe
parents:
diff changeset
60
13889
aaefd6b28188 toeplitz.m: Update coding style
Rik <octave@nomad.inbox5.com>
parents: 12639
diff changeset
61 if (nargin < 1 || nargin > 2)
6046
34f96dd5441b [project @ 2006-10-10 16:10:25 by jwe]
jwe
parents: 5781
diff changeset
62 print_usage ();
4
b4df021f796c [project @ 1993-08-08 01:26:08 by jwe]
jwe
parents:
diff changeset
63 endif
b4df021f796c [project @ 1993-08-08 01:26:08 by jwe]
jwe
parents:
diff changeset
64
13889
aaefd6b28188 toeplitz.m: Update coding style
Rik <octave@nomad.inbox5.com>
parents: 12639
diff changeset
65 if (nargin == 1)
aaefd6b28188 toeplitz.m: Update coding style
Rik <octave@nomad.inbox5.com>
parents: 12639
diff changeset
66 if (! isvector (c))
aaefd6b28188 toeplitz.m: Update coding style
Rik <octave@nomad.inbox5.com>
parents: 12639
diff changeset
67 error ("toeplitz: C must be a vector");
aaefd6b28188 toeplitz.m: Update coding style
Rik <octave@nomad.inbox5.com>
parents: 12639
diff changeset
68 endif
aaefd6b28188 toeplitz.m: Update coding style
Rik <octave@nomad.inbox5.com>
parents: 12639
diff changeset
69
aaefd6b28188 toeplitz.m: Update coding style
Rik <octave@nomad.inbox5.com>
parents: 12639
diff changeset
70 r = c;
aaefd6b28188 toeplitz.m: Update coding style
Rik <octave@nomad.inbox5.com>
parents: 12639
diff changeset
71 nr = length (c);
aaefd6b28188 toeplitz.m: Update coding style
Rik <octave@nomad.inbox5.com>
parents: 12639
diff changeset
72 nc = nr;
aaefd6b28188 toeplitz.m: Update coding style
Rik <octave@nomad.inbox5.com>
parents: 12639
diff changeset
73 else
aaefd6b28188 toeplitz.m: Update coding style
Rik <octave@nomad.inbox5.com>
parents: 12639
diff changeset
74 if (! (isvector (c) && isvector (r)))
aaefd6b28188 toeplitz.m: Update coding style
Rik <octave@nomad.inbox5.com>
parents: 12639
diff changeset
75 error ("toeplitz: C and R must be vectors");
aaefd6b28188 toeplitz.m: Update coding style
Rik <octave@nomad.inbox5.com>
parents: 12639
diff changeset
76 elseif (r(1) != c(1))
aaefd6b28188 toeplitz.m: Update coding style
Rik <octave@nomad.inbox5.com>
parents: 12639
diff changeset
77 warning ("toeplitz: column wins anti-diagonal conflict");
aaefd6b28188 toeplitz.m: Update coding style
Rik <octave@nomad.inbox5.com>
parents: 12639
diff changeset
78 endif
aaefd6b28188 toeplitz.m: Update coding style
Rik <octave@nomad.inbox5.com>
parents: 12639
diff changeset
79
aaefd6b28188 toeplitz.m: Update coding style
Rik <octave@nomad.inbox5.com>
parents: 12639
diff changeset
80 nr = length (c);
aaefd6b28188 toeplitz.m: Update coding style
Rik <octave@nomad.inbox5.com>
parents: 12639
diff changeset
81 nc = length (r);
4
b4df021f796c [project @ 1993-08-08 01:26:08 by jwe]
jwe
parents:
diff changeset
82 endif
b4df021f796c [project @ 1993-08-08 01:26:08 by jwe]
jwe
parents:
diff changeset
83
9092
38c3a0f8c6d0 optimize toeplitz
Jaroslav Hajek <highegg@gmail.com>
parents: 9041
diff changeset
84 if (nr == 0 || nc == 0)
38c3a0f8c6d0 optimize toeplitz
Jaroslav Hajek <highegg@gmail.com>
parents: 9041
diff changeset
85 ## Empty matrix.
38c3a0f8c6d0 optimize toeplitz
Jaroslav Hajek <highegg@gmail.com>
parents: 9041
diff changeset
86 retval = zeros (nr, nc, class (c));
38c3a0f8c6d0 optimize toeplitz
Jaroslav Hajek <highegg@gmail.com>
parents: 9041
diff changeset
87 return;
4
b4df021f796c [project @ 1993-08-08 01:26:08 by jwe]
jwe
parents:
diff changeset
88 endif
b4df021f796c [project @ 1993-08-08 01:26:08 by jwe]
jwe
parents:
diff changeset
89
2303
5cffc4b8de57 [project @ 1996-06-24 09:15:24 by jwe]
jwe
parents: 1887
diff changeset
90 ## If we have a single complex argument, we want to return a
5cffc4b8de57 [project @ 1996-06-24 09:15:24 by jwe]
jwe
parents: 1887
diff changeset
91 ## Hermitian-symmetric matrix (actually, this will really only be
5cffc4b8de57 [project @ 1996-06-24 09:15:24 by jwe]
jwe
parents: 1887
diff changeset
92 ## Hermitian-symmetric if the first element of the vector is real).
9092
38c3a0f8c6d0 optimize toeplitz
Jaroslav Hajek <highegg@gmail.com>
parents: 9041
diff changeset
93 if (nargin == 1 && iscomplex (c))
1016
58ce441b5426 [project @ 1995-01-06 19:18:43 by jwe]
jwe
parents: 1015
diff changeset
94 c = conj (c);
58ce441b5426 [project @ 1995-01-06 19:18:43 by jwe]
jwe
parents: 1015
diff changeset
95 c(1) = conj (c(1));
58ce441b5426 [project @ 1995-01-06 19:18:43 by jwe]
jwe
parents: 1015
diff changeset
96 endif
58ce441b5426 [project @ 1995-01-06 19:18:43 by jwe]
jwe
parents: 1015
diff changeset
97
13889
aaefd6b28188 toeplitz.m: Update coding style
Rik <octave@nomad.inbox5.com>
parents: 12639
diff changeset
98 if (issparse (c) && issparse (r))
aaefd6b28188 toeplitz.m: Update coding style
Rik <octave@nomad.inbox5.com>
parents: 12639
diff changeset
99 c = c(:).'; ## enforce row vector
aaefd6b28188 toeplitz.m: Update coding style
Rik <octave@nomad.inbox5.com>
parents: 12639
diff changeset
100 r = r(:).'; ## enforce row vector
aaefd6b28188 toeplitz.m: Update coding style
Rik <octave@nomad.inbox5.com>
parents: 12639
diff changeset
101 cidx = find (c);
aaefd6b28188 toeplitz.m: Update coding style
Rik <octave@nomad.inbox5.com>
parents: 12639
diff changeset
102 ridx = find (r);
9126
5780b3b80425 optimize toeplitz for sparse args
Jaroslav Hajek <highegg@gmail.com>
parents: 9092
diff changeset
103
5780b3b80425 optimize toeplitz for sparse args
Jaroslav Hajek <highegg@gmail.com>
parents: 9092
diff changeset
104 ## Ignore the first element in r.
5780b3b80425 optimize toeplitz for sparse args
Jaroslav Hajek <highegg@gmail.com>
parents: 9092
diff changeset
105 ridx = ridx(ridx > 1);
4
b4df021f796c [project @ 1993-08-08 01:26:08 by jwe]
jwe
parents:
diff changeset
106
9126
5780b3b80425 optimize toeplitz for sparse args
Jaroslav Hajek <highegg@gmail.com>
parents: 9092
diff changeset
107 ## Form matrix.
14868
5d3a684236b0 maint: Use Octave coding conventions for cuddling parentheses in scripts directory
Rik <octave@nomad.inbox5.com>
parents: 14138
diff changeset
108 retval = spdiags (repmat (c(cidx),nr,1),1-cidx,nr,nc) + ...
5d3a684236b0 maint: Use Octave coding conventions for cuddling parentheses in scripts directory
Rik <octave@nomad.inbox5.com>
parents: 14138
diff changeset
109 spdiags (repmat (r(ridx),nr,1),ridx-1,nr,nc);
11587
c792872f8942 all script files: untabify and strip trailing whitespace
John W. Eaton <jwe@octave.org>
parents: 11523
diff changeset
110 else
9126
5780b3b80425 optimize toeplitz for sparse args
Jaroslav Hajek <highegg@gmail.com>
parents: 9092
diff changeset
111 ## Concatenate data into a single column vector.
5780b3b80425 optimize toeplitz for sparse args
Jaroslav Hajek <highegg@gmail.com>
parents: 9092
diff changeset
112 data = [r(end:-1:2)(:); c(:)];
4
b4df021f796c [project @ 1993-08-08 01:26:08 by jwe]
jwe
parents:
diff changeset
113
9126
5780b3b80425 optimize toeplitz for sparse args
Jaroslav Hajek <highegg@gmail.com>
parents: 9092
diff changeset
114 ## Get slices.
5780b3b80425 optimize toeplitz for sparse args
Jaroslav Hajek <highegg@gmail.com>
parents: 9092
diff changeset
115 slices = cellslices (data, nc:-1:1, nc+nr-1:-1:nr);
4
b4df021f796c [project @ 1993-08-08 01:26:08 by jwe]
jwe
parents:
diff changeset
116
9126
5780b3b80425 optimize toeplitz for sparse args
Jaroslav Hajek <highegg@gmail.com>
parents: 9092
diff changeset
117 ## Form matrix.
5780b3b80425 optimize toeplitz for sparse args
Jaroslav Hajek <highegg@gmail.com>
parents: 9092
diff changeset
118 retval = horzcat (slices{:});
5780b3b80425 optimize toeplitz for sparse args
Jaroslav Hajek <highegg@gmail.com>
parents: 9092
diff changeset
119 endif
13889
aaefd6b28188 toeplitz.m: Update coding style
Rik <octave@nomad.inbox5.com>
parents: 12639
diff changeset
120
4
b4df021f796c [project @ 1993-08-08 01:26:08 by jwe]
jwe
parents:
diff changeset
121 endfunction
7411
83a8781b529d [project @ 2008-01-22 21:52:25 by jwe]
jwe
parents: 7017
diff changeset
122
13889
aaefd6b28188 toeplitz.m: Update coding style
Rik <octave@nomad.inbox5.com>
parents: 12639
diff changeset
123
aaefd6b28188 toeplitz.m: Update coding style
Rik <octave@nomad.inbox5.com>
parents: 12639
diff changeset
124 %!assert (toeplitz (1), [1])
aaefd6b28188 toeplitz.m: Update coding style
Rik <octave@nomad.inbox5.com>
parents: 12639
diff changeset
125 %!assert (toeplitz ([1, 2, 3], [1; -3; -5]), [1, -3, -5; 2, 1, -3; 3, 2, 1])
aaefd6b28188 toeplitz.m: Update coding style
Rik <octave@nomad.inbox5.com>
parents: 12639
diff changeset
126 %!assert (toeplitz ([1, 2, 3], [1; -3i; -5i]), [1, -3i, -5i; 2, 1, -3i; 3, 2, 1])
7411
83a8781b529d [project @ 2008-01-22 21:52:25 by jwe]
jwe
parents: 7017
diff changeset
127
13889
aaefd6b28188 toeplitz.m: Update coding style
Rik <octave@nomad.inbox5.com>
parents: 12639
diff changeset
128 %% Test input validation
aaefd6b28188 toeplitz.m: Update coding style
Rik <octave@nomad.inbox5.com>
parents: 12639
diff changeset
129 %!error toeplitz ()
aaefd6b28188 toeplitz.m: Update coding style
Rik <octave@nomad.inbox5.com>
parents: 12639
diff changeset
130 %!error toeplitz (1, 2, 3)
aaefd6b28188 toeplitz.m: Update coding style
Rik <octave@nomad.inbox5.com>
parents: 12639
diff changeset
131 %!error <C must be a vector> toeplitz ([1, 2; 3, 4])
aaefd6b28188 toeplitz.m: Update coding style
Rik <octave@nomad.inbox5.com>
parents: 12639
diff changeset
132 %!error <C and R must be vectors> toeplitz ([1, 2; 3, 4], 1)
aaefd6b28188 toeplitz.m: Update coding style
Rik <octave@nomad.inbox5.com>
parents: 12639
diff changeset
133 %!error <C and R must be vectors> toeplitz (1, [1, 2; 3, 4])
7411
83a8781b529d [project @ 2008-01-22 21:52:25 by jwe]
jwe
parents: 7017
diff changeset
134