Mercurial > hg > octave-lyh

annotate scripts/polynomial/conv.m @ 13141:e81ddf9cacd5

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
maint: untabify and remove trailing whitespace from source files * bicg.m, gmres.m, pkg.m: Untabify and remove trailing whitespace. * libcruft/Makefile.am, libcruft/blas-xtra/cdotc3.f, libcruft/blas-xtra/cmatm3.f, libcruft/blas-xtra/ddot3.f, libcruft/blas-xtra/dmatm3.f, libcruft/blas-xtra/sdot3.f, libcruft/blas-xtra/smatm3.f, libcruft/blas-xtra/zdotc3.f, libcruft/blas-xtra/zmatm3.f, libcruft/lapack-xtra/crsf2csf.f, libcruft/lapack-xtra/zrsf2csf.f, liboctave/Array.cc, liboctave/DASPK-opts.in, liboctave/DASRT-opts.in, liboctave/DASSL-opts.in, liboctave/LSODE-opts.in, liboctave/Makefile.a,mliboctave/Quad-opts.in, liboctave/Sparse-perm-op-defs.h, scripts/Makefile.a,mscripts/deprecated/glpkmex.m, scripts/general/blkdiag.m, scripts/general/interp1.m, scripts/general/profshow.m, scripts/general/quadl.m, scripts/general/triplequad.m, scripts/help/__makeinfo__.m, scripts/io/strread.m, scripts/io/textread.m, scripts/io/textscan.m, scripts/linear-algebra/rank.m, scripts/miscellaneous/gzip.m, scripts/miscellaneous/private/__xzip__.m, scripts/miscellaneous/tempdir.m, scripts/miscellaneous/unpack.m, scripts/pkg/pkg.m, scripts/plot/allchild.m, scripts/plot/ancestor.m, scripts/plot/cla.m, scripts/plot/clf.m, scripts/plot/findall.m, scripts/plot/findobj.m, scripts/plot/gca.m, scripts/plot/gcf.m, scripts/plot/hggroup.m, scripts/plot/isfigure.m, scripts/plot/ishghandle.m, scripts/plot/legend.m, scripts/plot/line.m, scripts/plot/loglog.m, scripts/plot/patch.m, scripts/plot/print.m, scripts/plot/private/__quiver__.m, scripts/plot/private/__scatter__.m, scripts/plot/rectangle.m, scripts/plot/semilogx.m, scripts/plot/semilogy.m, scripts/plot/surface.m, scripts/plot/text.m, scripts/plot/title.m, scripts/plot/trisurf.m, scripts/plot/view.m, scripts/plot/whitebg.m, scripts/plot/xlabel.m, scripts/plot/xlim.m, scripts/plot/ylabel.m, scripts/plot/ylim.m, scripts/plot/zlabel.m, scripts/plot/zlim.m, scripts/polynomial/mkpp.m, scripts/polynomial/polygcd.m, scripts/polynomial/ppint.m, scripts/polynomial/ppjumps.m, scripts/polynomial/ppval.m, scripts/set/setxor.m, scripts/sparse/bicgstab.m, scripts/sparse/cgs.m, scripts/sparse/spconvert.m, scripts/specfun/nthroot.m, scripts/strings/strmatch.m, scripts/strings/untabify.m, scripts/testfun/demo.m, scripts/testfun/example.m, src/DLD-FUNCTIONS/filter.cc, src/DLD-FUNCTIONS/mgorth.cc, src/DLD-FUNCTIONS/quadcc.cc, src/DLD-FUNCTIONS/str2double.cc, src/Makefile.a,msrc/gl-render.cc, src/gl2ps-renderer.cc, src/graphics.cc, src/octave-config.cc.in, src/octave-config.in, src/ov-class.h, src/ov-fcn.h, src/profiler.cc, src/profiler.h, src/pt-binop.cc, src/pt-unop.cc, src/symtab.cc, src/txt-eng-ft.cc: Remove trailing whitespace.
author John W. Eaton <jwe@octave.org>
date Thu, 15 Sep 2011 12:51:10 -0400
parents fd0a3ac60b0e
children 96c7143304d9
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: 11124
diff changeset
1 ## Copyright (C) 1994-2011 John W. Eaton
2313
5ca126254d15 [project @ 1996-07-11 21:25:22 by jwe]
jwe
parents: 2312
diff changeset
2 ##
5ca126254d15 [project @ 1996-07-11 21:25:22 by jwe]
jwe
parents: 2312
diff changeset
3 ## This file is part of Octave.
5ca126254d15 [project @ 1996-07-11 21:25:22 by jwe]
jwe
parents: 2312
diff changeset
4 ##
5ca126254d15 [project @ 1996-07-11 21:25:22 by jwe]
jwe
parents: 2312
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: 2312
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: 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.
2313
5ca126254d15 [project @ 1996-07-11 21:25:22 by jwe]
jwe
parents: 2312
diff changeset
9 ##
5ca126254d15 [project @ 1996-07-11 21:25:22 by jwe]
jwe
parents: 2312
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: 2312
diff changeset
11 ## WITHOUT ANY WARRANTY; without even the implied warranty of
5ca126254d15 [project @ 1996-07-11 21:25:22 by jwe]
jwe
parents: 2312
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: 2312
diff changeset
13 ## General Public License for more details.
5ca126254d15 [project @ 1996-07-11 21:25:22 by jwe]
jwe
parents: 2312
diff changeset
14 ##
5ca126254d15 [project @ 1996-07-11 21:25:22 by jwe]
jwe
parents: 2312
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/>.
1025
f558749713f1 [project @ 1995-01-11 20:52:10 by jwe]
jwe
parents: 904
diff changeset
18
3368
a4cd1e9d9962 [project @ 1999-11-20 17:22:48 by jwe]
jwe
parents: 3202
diff changeset
19 ## -*- texinfo -*-
11085
2beacd515e09 Update docstrings for convolution family of functions (conv, conv2, fftconv)
Rik <octave@nomad.inbox5.com>
parents: 11084
diff changeset
20 ## @deftypefn {Function File} {} conv (@var{a}, @var{b})
2beacd515e09 Update docstrings for convolution family of functions (conv, conv2, fftconv)
Rik <octave@nomad.inbox5.com>
parents: 11084
diff changeset
21 ## @deftypefnx {Function File} {} conv (@var{a}, @var{b}, @var{shape})
2311
2b5788792cad [project @ 1996-07-11 20:18:38 by jwe]
jwe
parents: 2303
diff changeset
22 ## Convolve two vectors.
3426
f8dde1807dee [project @ 2000-01-13 08:40:00 by jwe]
jwe
parents: 3368
diff changeset
23 ##
11085
2beacd515e09 Update docstrings for convolution family of functions (conv, conv2, fftconv)
Rik <octave@nomad.inbox5.com>
parents: 11084
diff changeset
24 ## @code{c = conv (@var{a}, @var{b})} returns a vector of length equal to
2beacd515e09 Update docstrings for convolution family of functions (conv, conv2, fftconv)
Rik <octave@nomad.inbox5.com>
parents: 11084
diff changeset
25 ## @code{length (@var{a}) + length (@var{b}) - 1}.
2beacd515e09 Update docstrings for convolution family of functions (conv, conv2, fftconv)
Rik <octave@nomad.inbox5.com>
parents: 11084
diff changeset
26 ## If @var{a} and @var{b} are the coefficient vectors of two polynomials, the
2beacd515e09 Update docstrings for convolution family of functions (conv, conv2, fftconv)
Rik <octave@nomad.inbox5.com>
parents: 11084
diff changeset
27 ## returned value is the coefficient vector of the product polynomial.
11082
4558aad4c41d conv.m: handle "same" shape argument
John W. Eaton <jwe@octave.org>
parents: 10224
diff changeset
28 ##
11085
2beacd515e09 Update docstrings for convolution family of functions (conv, conv2, fftconv)
Rik <octave@nomad.inbox5.com>
parents: 11084
diff changeset
29 ## The optional @var{shape} argument may be
11082
4558aad4c41d conv.m: handle "same" shape argument
John W. Eaton <jwe@octave.org>
parents: 10224
diff changeset
30 ##
11084
0f6c5efce96e conv.m: docstring fix
John W. Eaton <jwe@octave.org>
parents: 11082
diff changeset
31 ## @table @asis
11082
4558aad4c41d conv.m: handle "same" shape argument
John W. Eaton <jwe@octave.org>
parents: 10224
diff changeset
32 ## @item @var{shape} = "full"
11085
2beacd515e09 Update docstrings for convolution family of functions (conv, conv2, fftconv)
Rik <octave@nomad.inbox5.com>
parents: 11084
diff changeset
33 ## Return the full convolution. (default)
11084
0f6c5efce96e conv.m: docstring fix
John W. Eaton <jwe@octave.org>
parents: 11082
diff changeset
34 ##
11082
4558aad4c41d conv.m: handle "same" shape argument
John W. Eaton <jwe@octave.org>
parents: 10224
diff changeset
35 ## @item @var{shape} = "same"
11085
2beacd515e09 Update docstrings for convolution family of functions (conv, conv2, fftconv)
Rik <octave@nomad.inbox5.com>
parents: 11084
diff changeset
36 ## Return the central part of the convolution with the same size as @var{a}.
11082
4558aad4c41d conv.m: handle "same" shape argument
John W. Eaton <jwe@octave.org>
parents: 10224
diff changeset
37 ## @end table
4558aad4c41d conv.m: handle "same" shape argument
John W. Eaton <jwe@octave.org>
parents: 10224
diff changeset
38 ##
11085
2beacd515e09 Update docstrings for convolution family of functions (conv, conv2, fftconv)
Rik <octave@nomad.inbox5.com>
parents: 11084
diff changeset
39 ## @seealso{deconv, fftconv, conv2, poly}
3368
a4cd1e9d9962 [project @ 1999-11-20 17:22:48 by jwe]
jwe
parents: 3202
diff changeset
40 ## @end deftypefn
2311
2b5788792cad [project @ 1996-07-11 20:18:38 by jwe]
jwe
parents: 2303
diff changeset
41
3202
44d82b369c78 [project @ 1998-10-29 20:27:57 by jwe]
jwe
parents: 2847
diff changeset
42 ## Author: Tony Richardson <arichard@stark.cc.oh.us>
2312
204cc7db6f4a [project @ 1996-07-11 21:20:36 by jwe]
jwe
parents: 2311
diff changeset
43 ## Created: June 1994
204cc7db6f4a [project @ 1996-07-11 21:20:36 by jwe]
jwe
parents: 2311
diff changeset
44 ## Adapted-By: jwe
204cc7db6f4a [project @ 1996-07-11 21:20:36 by jwe]
jwe
parents: 2311
diff changeset
45
11082
4558aad4c41d conv.m: handle "same" shape argument
John W. Eaton <jwe@octave.org>
parents: 10224
diff changeset
46 function y = conv (a, b, shape = "full")
2325
b5568c31ee2c [project @ 1996-07-15 22:20:21 by jwe]
jwe
parents: 2313
diff changeset
47
11082
4558aad4c41d conv.m: handle "same" shape argument
John W. Eaton <jwe@octave.org>
parents: 10224
diff changeset
48 if (nargin < 2 || nargin > 3)
6046
34f96dd5441b [project @ 2006-10-10 16:10:25 by jwe]
jwe
parents: 5642
diff changeset
49 print_usage ();
787
c5d35bb139b6 [project @ 1994-10-11 00:34:13 by jwe]
jwe
parents:
diff changeset
50 endif
c5d35bb139b6 [project @ 1994-10-11 00:34:13 by jwe]
jwe
parents:
diff changeset
51
4030
22bd65326ec1 [project @ 2002-08-09 18:58:13 by jwe]
jwe
parents: 3457
diff changeset
52 if (! (isvector (a) && isvector (b)))
9433
38a0f9dc0ab4 conv.m: fix Matlab incompatibility; new tests
Robert T. Short <octave@phaselockedsystems.com>
parents: 8920
diff changeset
53 error ("conv: both arguments must be vectors");
787
c5d35bb139b6 [project @ 1994-10-11 00:34:13 by jwe]
jwe
parents:
diff changeset
54 endif
c5d35bb139b6 [project @ 1994-10-11 00:34:13 by jwe]
jwe
parents:
diff changeset
55
c5d35bb139b6 [project @ 1994-10-11 00:34:13 by jwe]
jwe
parents:
diff changeset
56 la = length (a);
c5d35bb139b6 [project @ 1994-10-11 00:34:13 by jwe]
jwe
parents:
diff changeset
57 lb = length (b);
c5d35bb139b6 [project @ 1994-10-11 00:34:13 by jwe]
jwe
parents:
diff changeset
58
c5d35bb139b6 [project @ 1994-10-11 00:34:13 by jwe]
jwe
parents:
diff changeset
59 ly = la + lb - 1;
c5d35bb139b6 [project @ 1994-10-11 00:34:13 by jwe]
jwe
parents:
diff changeset
60
11082
4558aad4c41d conv.m: handle "same" shape argument
John W. Eaton <jwe@octave.org>
parents: 10224
diff changeset
61 if (ly == 0)
4558aad4c41d conv.m: handle "same" shape argument
John W. Eaton <jwe@octave.org>
parents: 10224
diff changeset
62 y = zeros (1, 0);
4558aad4c41d conv.m: handle "same" shape argument
John W. Eaton <jwe@octave.org>
parents: 10224
diff changeset
63 else
4558aad4c41d conv.m: handle "same" shape argument
John W. Eaton <jwe@octave.org>
parents: 10224
diff changeset
64 ## Use the shortest vector as the coefficent vector to filter.
4558aad4c41d conv.m: handle "same" shape argument
John W. Eaton <jwe@octave.org>
parents: 10224
diff changeset
65 ## Preserve the row/column orientation of the longer input.
4558aad4c41d conv.m: handle "same" shape argument
John W. Eaton <jwe@octave.org>
parents: 10224
diff changeset
66 if (la <= lb)
4558aad4c41d conv.m: handle "same" shape argument
John W. Eaton <jwe@octave.org>
parents: 10224
diff changeset
67 if (ly > lb)
4558aad4c41d conv.m: handle "same" shape argument
John W. Eaton <jwe@octave.org>
parents: 10224
diff changeset
68 if (size (b, 1) <= size (b, 2))
4558aad4c41d conv.m: handle "same" shape argument
John W. Eaton <jwe@octave.org>
parents: 10224
diff changeset
69 x = [b, (zeros (1, ly - lb))];
4558aad4c41d conv.m: handle "same" shape argument
John W. Eaton <jwe@octave.org>
parents: 10224
diff changeset
70 else
4558aad4c41d conv.m: handle "same" shape argument
John W. Eaton <jwe@octave.org>
parents: 10224
diff changeset
71 x = [b; (zeros (ly - lb, 1))];
4558aad4c41d conv.m: handle "same" shape argument
John W. Eaton <jwe@octave.org>
parents: 10224
diff changeset
72 endif
8158
15e4a450bf84 conv.m: Correct row/col orientation of output
Ben Abbott <bpabbott@mac.com>
parents: 7411
diff changeset
73 else
11082
4558aad4c41d conv.m: handle "same" shape argument
John W. Eaton <jwe@octave.org>
parents: 10224
diff changeset
74 x = b;
8158
15e4a450bf84 conv.m: Correct row/col orientation of output
Ben Abbott <bpabbott@mac.com>
parents: 7411
diff changeset
75 endif
11082
4558aad4c41d conv.m: handle "same" shape argument
John W. Eaton <jwe@octave.org>
parents: 10224
diff changeset
76 y = filter (a, 1, x);
787
c5d35bb139b6 [project @ 1994-10-11 00:34:13 by jwe]
jwe
parents:
diff changeset
77 else
11082
4558aad4c41d conv.m: handle "same" shape argument
John W. Eaton <jwe@octave.org>
parents: 10224
diff changeset
78 if (ly > la)
4558aad4c41d conv.m: handle "same" shape argument
John W. Eaton <jwe@octave.org>
parents: 10224
diff changeset
79 if (size (a, 1) <= size (a, 2))
4558aad4c41d conv.m: handle "same" shape argument
John W. Eaton <jwe@octave.org>
parents: 10224
diff changeset
80 x = [a, (zeros (1, ly - la))];
4558aad4c41d conv.m: handle "same" shape argument
John W. Eaton <jwe@octave.org>
parents: 10224
diff changeset
81 else
4558aad4c41d conv.m: handle "same" shape argument
John W. Eaton <jwe@octave.org>
parents: 10224
diff changeset
82 x = [a; (zeros (ly - la, 1))];
4558aad4c41d conv.m: handle "same" shape argument
John W. Eaton <jwe@octave.org>
parents: 10224
diff changeset
83 endif
4558aad4c41d conv.m: handle "same" shape argument
John W. Eaton <jwe@octave.org>
parents: 10224
diff changeset
84 else
4558aad4c41d conv.m: handle "same" shape argument
John W. Eaton <jwe@octave.org>
parents: 10224
diff changeset
85 x = a;
4558aad4c41d conv.m: handle "same" shape argument
John W. Eaton <jwe@octave.org>
parents: 10224
diff changeset
86 endif
4558aad4c41d conv.m: handle "same" shape argument
John W. Eaton <jwe@octave.org>
parents: 10224
diff changeset
87 y = filter (b, 1, x);
787
c5d35bb139b6 [project @ 1994-10-11 00:34:13 by jwe]
jwe
parents:
diff changeset
88 endif
11082
4558aad4c41d conv.m: handle "same" shape argument
John W. Eaton <jwe@octave.org>
parents: 10224
diff changeset
89 if (strcmp (shape, "same"))
4558aad4c41d conv.m: handle "same" shape argument
John W. Eaton <jwe@octave.org>
parents: 10224
diff changeset
90 idx = ceil ((ly - la) / 2);
4558aad4c41d conv.m: handle "same" shape argument
John W. Eaton <jwe@octave.org>
parents: 10224
diff changeset
91 y = y(idx+1:idx+la);
787
c5d35bb139b6 [project @ 1994-10-11 00:34:13 by jwe]
jwe
parents:
diff changeset
92 endif
c5d35bb139b6 [project @ 1994-10-11 00:34:13 by jwe]
jwe
parents:
diff changeset
93 endif
c5d35bb139b6 [project @ 1994-10-11 00:34:13 by jwe]
jwe
parents:
diff changeset
94
c5d35bb139b6 [project @ 1994-10-11 00:34:13 by jwe]
jwe
parents:
diff changeset
95 endfunction
7411
83a8781b529d [project @ 2008-01-22 21:52:25 by jwe]
jwe
parents: 7017
diff changeset
96
9433
38a0f9dc0ab4 conv.m: fix Matlab incompatibility; new tests
Robert T. Short <octave@phaselockedsystems.com>
parents: 8920
diff changeset
97 %!test
38a0f9dc0ab4 conv.m: fix Matlab incompatibility; new tests
Robert T. Short <octave@phaselockedsystems.com>
parents: 8920
diff changeset
98 %! x = ones(3,1);
38a0f9dc0ab4 conv.m: fix Matlab incompatibility; new tests
Robert T. Short <octave@phaselockedsystems.com>
parents: 8920
diff changeset
99 %! y = ones(1,3);
38a0f9dc0ab4 conv.m: fix Matlab incompatibility; new tests
Robert T. Short <octave@phaselockedsystems.com>
parents: 8920
diff changeset
100 %! b = 2;
38a0f9dc0ab4 conv.m: fix Matlab incompatibility; new tests
Robert T. Short <octave@phaselockedsystems.com>
parents: 8920
diff changeset
101 %! c = 3;
38a0f9dc0ab4 conv.m: fix Matlab incompatibility; new tests
Robert T. Short <octave@phaselockedsystems.com>
parents: 8920
diff changeset
102 %! assert (conv (x, x), [1; 2; 3; 2; 1]);
38a0f9dc0ab4 conv.m: fix Matlab incompatibility; new tests
Robert T. Short <octave@phaselockedsystems.com>
parents: 8920
diff changeset
103 %! assert (conv (y, y), [1, 2, 3, 2, 1]);
38a0f9dc0ab4 conv.m: fix Matlab incompatibility; new tests
Robert T. Short <octave@phaselockedsystems.com>
parents: 8920
diff changeset
104 %! assert (conv (x, y), [1, 2, 3, 2, 1]);
38a0f9dc0ab4 conv.m: fix Matlab incompatibility; new tests
Robert T. Short <octave@phaselockedsystems.com>
parents: 8920
diff changeset
105 %! assert (conv (y, x), [1; 2; 3; 2; 1]);
38a0f9dc0ab4 conv.m: fix Matlab incompatibility; new tests
Robert T. Short <octave@phaselockedsystems.com>
parents: 8920
diff changeset
106 %! assert (conv (c, x), [3; 3; 3]);
38a0f9dc0ab4 conv.m: fix Matlab incompatibility; new tests
Robert T. Short <octave@phaselockedsystems.com>
parents: 8920
diff changeset
107 %! assert (conv (c, y), [3, 3, 3]);
38a0f9dc0ab4 conv.m: fix Matlab incompatibility; new tests
Robert T. Short <octave@phaselockedsystems.com>
parents: 8920
diff changeset
108 %! assert (conv (x, c), [3; 3; 3]);
38a0f9dc0ab4 conv.m: fix Matlab incompatibility; new tests
Robert T. Short <octave@phaselockedsystems.com>
parents: 8920
diff changeset
109 %! assert (conv (y, c), [3, 3, 3]);
38a0f9dc0ab4 conv.m: fix Matlab incompatibility; new tests
Robert T. Short <octave@phaselockedsystems.com>
parents: 8920
diff changeset
110 %! assert (conv (b, c), 6);
11085
2beacd515e09 Update docstrings for convolution family of functions (conv, conv2, fftconv)
Rik <octave@nomad.inbox5.com>
parents: 11084
diff changeset
111
8158
15e4a450bf84 conv.m: Correct row/col orientation of output
Ben Abbott <bpabbott@mac.com>
parents: 7411
diff changeset
112 %!test
11082
4558aad4c41d conv.m: handle "same" shape argument
John W. Eaton <jwe@octave.org>
parents: 10224
diff changeset
113 %! a = 1:10;
4558aad4c41d conv.m: handle "same" shape argument
John W. Eaton <jwe@octave.org>
parents: 10224
diff changeset
114 %! b = 1:3;
4558aad4c41d conv.m: handle "same" shape argument
John W. Eaton <jwe@octave.org>
parents: 10224
diff changeset
115 %! assert (size(conv(a,b)), [1, numel(a)+numel(b)-1])
4558aad4c41d conv.m: handle "same" shape argument
John W. Eaton <jwe@octave.org>
parents: 10224
diff changeset
116 %! assert (size(conv(b,a)), [1, numel(a)+numel(b)-1])
4558aad4c41d conv.m: handle "same" shape argument
John W. Eaton <jwe@octave.org>
parents: 10224
diff changeset
117
9433
38a0f9dc0ab4 conv.m: fix Matlab incompatibility; new tests
Robert T. Short <octave@phaselockedsystems.com>
parents: 8920
diff changeset
118 %! a = (1:10).';
38a0f9dc0ab4 conv.m: fix Matlab incompatibility; new tests
Robert T. Short <octave@phaselockedsystems.com>
parents: 8920
diff changeset
119 %! b = 1:3;
38a0f9dc0ab4 conv.m: fix Matlab incompatibility; new tests
Robert T. Short <octave@phaselockedsystems.com>
parents: 8920
diff changeset
120 %! assert (size(conv(a,b)), [numel(a)+numel(b)-1, 1])
38a0f9dc0ab4 conv.m: fix Matlab incompatibility; new tests
Robert T. Short <octave@phaselockedsystems.com>
parents: 8920
diff changeset
121 %! assert (size(conv(b,a)), [numel(a)+numel(b)-1, 1])
38a0f9dc0ab4 conv.m: fix Matlab incompatibility; new tests
Robert T. Short <octave@phaselockedsystems.com>
parents: 8920
diff changeset
122
8158
15e4a450bf84 conv.m: Correct row/col orientation of output
Ben Abbott <bpabbott@mac.com>
parents: 7411
diff changeset
123 %!test
9433
38a0f9dc0ab4 conv.m: fix Matlab incompatibility; new tests
Robert T. Short <octave@phaselockedsystems.com>
parents: 8920
diff changeset
124 %! a = 1:10;
11124
e79f59d31a74 Add tests for fftconv.m
Rik <octave@nomad.inbox5.com>
parents: 11085
diff changeset
125 %! b = (1:3).';
e79f59d31a74 Add tests for fftconv.m
Rik <octave@nomad.inbox5.com>
parents: 11085
diff changeset
126 %! assert (size(conv(a,b)), [1, numel(a)+numel(b)-1])
e79f59d31a74 Add tests for fftconv.m
Rik <octave@nomad.inbox5.com>
parents: 11085
diff changeset
127 %! assert (size(conv(b,a)), [1, numel(a)+numel(b)-1])
e79f59d31a74 Add tests for fftconv.m
Rik <octave@nomad.inbox5.com>
parents: 11085
diff changeset
128
e79f59d31a74 Add tests for fftconv.m
Rik <octave@nomad.inbox5.com>
parents: 11085
diff changeset
129 %!test
e79f59d31a74 Add tests for fftconv.m
Rik <octave@nomad.inbox5.com>
parents: 11085
diff changeset
130 %! a = 1:10;
11082
4558aad4c41d conv.m: handle "same" shape argument
John W. Eaton <jwe@octave.org>
parents: 10224
diff changeset
131 %! b = 1:3;
11085
2beacd515e09 Update docstrings for convolution family of functions (conv, conv2, fftconv)
Rik <octave@nomad.inbox5.com>
parents: 11084
diff changeset
132 %! assert (conv(a,b,"full"), conv(a,b))
2beacd515e09 Update docstrings for convolution family of functions (conv, conv2, fftconv)
Rik <octave@nomad.inbox5.com>
parents: 11084
diff changeset
133 %! assert (conv(b,a,"full"), conv(b,a))
2beacd515e09 Update docstrings for convolution family of functions (conv, conv2, fftconv)
Rik <octave@nomad.inbox5.com>
parents: 11084
diff changeset
134
2beacd515e09 Update docstrings for convolution family of functions (conv, conv2, fftconv)
Rik <octave@nomad.inbox5.com>
parents: 11084
diff changeset
135 %!test
2beacd515e09 Update docstrings for convolution family of functions (conv, conv2, fftconv)
Rik <octave@nomad.inbox5.com>
parents: 11084
diff changeset
136 %! a = 1:10;
2beacd515e09 Update docstrings for convolution family of functions (conv, conv2, fftconv)
Rik <octave@nomad.inbox5.com>
parents: 11084
diff changeset
137 %! b = 1:3;
11082
4558aad4c41d conv.m: handle "same" shape argument
John W. Eaton <jwe@octave.org>
parents: 10224
diff changeset
138 %! assert (conv(a,b,'same'), [4, 10, 16, 22, 28, 34, 40, 46, 52, 47])
4558aad4c41d conv.m: handle "same" shape argument
John W. Eaton <jwe@octave.org>
parents: 10224
diff changeset
139 %! assert (conv(b,a,'same'), [28, 34, 40])
4558aad4c41d conv.m: handle "same" shape argument
John W. Eaton <jwe@octave.org>
parents: 10224
diff changeset
140
11085
2beacd515e09 Update docstrings for convolution family of functions (conv, conv2, fftconv)
Rik <octave@nomad.inbox5.com>
parents: 11084
diff changeset
141 %% Test input validation
2beacd515e09 Update docstrings for convolution family of functions (conv, conv2, fftconv)
Rik <octave@nomad.inbox5.com>
parents: 11084
diff changeset
142 %!error conv (1);
2beacd515e09 Update docstrings for convolution family of functions (conv, conv2, fftconv)
Rik <octave@nomad.inbox5.com>
parents: 11084
diff changeset
143 %!error conv (1,2,3,4);
2beacd515e09 Update docstrings for convolution family of functions (conv, conv2, fftconv)
Rik <octave@nomad.inbox5.com>
parents: 11084
diff changeset
144 %!error conv ([1, 2; 3, 4], 3);
2beacd515e09 Update docstrings for convolution family of functions (conv, conv2, fftconv)
Rik <octave@nomad.inbox5.com>
parents: 11084
diff changeset
145 %!error conv (2, []);
11124
e79f59d31a74 Add tests for fftconv.m
Rik <octave@nomad.inbox5.com>
parents: 11085
diff changeset
146