Mercurial > octave
annotate scripts/signal/blackman.m @ 23219:3ac9f9ecfae5 stable
maint: Update copyright dates.
author | John W. Eaton <jwe@octave.org> |
---|---|
date | Wed, 22 Feb 2017 12:39:29 -0500 |
parents | e9a0469dedd9 |
children | 092078913d54 |
rev | line source |
---|---|
23219
3ac9f9ecfae5
maint: Update copyright dates.
John W. Eaton <jwe@octave.org>
parents:
23083
diff
changeset
|
1 ## Copyright (C) 1995-2017 Andreas Weingessel |
3426 | 2 ## |
3922 | 3 ## This file is part of Octave. |
4 ## | |
5 ## Octave is free software; you can redistribute it and/or modify it | |
6 ## under the terms of the GNU General Public License as published by | |
7016 | 7 ## the Free Software Foundation; either version 3 of the License, or (at |
8 ## your option) any later version. | |
3426 | 9 ## |
3922 | 10 ## Octave is distributed in the hope that it will be useful, but |
3191 | 11 ## WITHOUT ANY WARRANTY; without even the implied warranty of |
12 ## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU | |
3426 | 13 ## General Public License for more details. |
14 ## | |
3191 | 15 ## You should have received a copy of the GNU General Public License |
7016 | 16 ## along with Octave; see the file COPYING. If not, see |
17 ## <http://www.gnu.org/licenses/>. | |
3191 | 18 |
3449 | 19 ## -*- texinfo -*- |
20852
516bb87ea72e
2015 Code Sprint: remove class of function from docstring for all m-files.
Rik <rik@octave.org>
parents:
20231
diff
changeset
|
20 ## @deftypefn {} {} blackman (@var{m}) |
516bb87ea72e
2015 Code Sprint: remove class of function from docstring for all m-files.
Rik <rik@octave.org>
parents:
20231
diff
changeset
|
21 ## @deftypefnx {} {} blackman (@var{m}, "periodic") |
516bb87ea72e
2015 Code Sprint: remove class of function from docstring for all m-files.
Rik <rik@octave.org>
parents:
20231
diff
changeset
|
22 ## @deftypefnx {} {} blackman (@var{m}, "symmetric") |
3449 | 23 ## Return the filter coefficients of a Blackman window of length @var{m}. |
3191 | 24 ## |
19826
26fb4bfa4193
blackman, hamming, hanning: Add periodic window option (bug #43305)
Mike Miller <mtmiller@ieee.org>
parents:
19697
diff
changeset
|
25 ## If the optional argument @qcode{"periodic"} is given, the periodic form |
26fb4bfa4193
blackman, hamming, hanning: Add periodic window option (bug #43305)
Mike Miller <mtmiller@ieee.org>
parents:
19697
diff
changeset
|
26 ## of the window is returned. This is equivalent to the window of length |
26fb4bfa4193
blackman, hamming, hanning: Add periodic window option (bug #43305)
Mike Miller <mtmiller@ieee.org>
parents:
19697
diff
changeset
|
27 ## @var{m}+1 with the last coefficient removed. The optional argument |
26fb4bfa4193
blackman, hamming, hanning: Add periodic window option (bug #43305)
Mike Miller <mtmiller@ieee.org>
parents:
19697
diff
changeset
|
28 ## @qcode{"symmetric"} is equivalent to not specifying a second argument. |
26fb4bfa4193
blackman, hamming, hanning: Add periodic window option (bug #43305)
Mike Miller <mtmiller@ieee.org>
parents:
19697
diff
changeset
|
29 ## |
20165
f1d0f506ee78
doc: Update more docstrings to have one sentence summary as first line.
Rik <rik@octave.org>
parents:
20026
diff
changeset
|
30 ## For a definition of the Blackman window, see, e.g., |
19040
0850b5212619
doc: Add @nospell macro around proper names in documentation.
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
31 ## @nospell{A.V. Oppenheim & R. W. Schafer}, |
0850b5212619
doc: Add @nospell macro around proper names in documentation.
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
32 ## @cite{Discrete-Time Signal Processing}. |
3449 | 33 ## @end deftypefn |
3191 | 34 |
3457 | 35 ## Author: AW <Andreas.Weingessel@ci.tuwien.ac.at> |
36 ## Description: Coefficients of the Blackman window | |
3426 | 37 |
19826
26fb4bfa4193
blackman, hamming, hanning: Add periodic window option (bug #43305)
Mike Miller <mtmiller@ieee.org>
parents:
19697
diff
changeset
|
38 function c = blackman (m, opt) |
3426 | 39 |
19826
26fb4bfa4193
blackman, hamming, hanning: Add periodic window option (bug #43305)
Mike Miller <mtmiller@ieee.org>
parents:
19697
diff
changeset
|
40 if (nargin < 1 || nargin > 2) |
6046 | 41 print_usage (); |
3191 | 42 endif |
3426 | 43 |
13279
984359717d71
Use common code idiom for checking whether a double value is an integer.
Rik <octave@nomad.inbox5.com>
parents:
13068
diff
changeset
|
44 if (! (isscalar (m) && (m == fix (m)) && (m > 0))) |
19826
26fb4bfa4193
blackman, hamming, hanning: Add periodic window option (bug #43305)
Mike Miller <mtmiller@ieee.org>
parents:
19697
diff
changeset
|
45 error ("blackman: M must be a positive integer"); |
26fb4bfa4193
blackman, hamming, hanning: Add periodic window option (bug #43305)
Mike Miller <mtmiller@ieee.org>
parents:
19697
diff
changeset
|
46 endif |
26fb4bfa4193
blackman, hamming, hanning: Add periodic window option (bug #43305)
Mike Miller <mtmiller@ieee.org>
parents:
19697
diff
changeset
|
47 |
20026
230c1e2a678d
blackman, hamming, hanning: Simplify handling of periodic window option
Mike Miller <mtmiller@octave.org>
parents:
20025
diff
changeset
|
48 N = m - 1; |
19826
26fb4bfa4193
blackman, hamming, hanning: Add periodic window option (bug #43305)
Mike Miller <mtmiller@ieee.org>
parents:
19697
diff
changeset
|
49 if (nargin == 2) |
26fb4bfa4193
blackman, hamming, hanning: Add periodic window option (bug #43305)
Mike Miller <mtmiller@ieee.org>
parents:
19697
diff
changeset
|
50 switch (opt) |
26fb4bfa4193
blackman, hamming, hanning: Add periodic window option (bug #43305)
Mike Miller <mtmiller@ieee.org>
parents:
19697
diff
changeset
|
51 case "periodic" |
20026
230c1e2a678d
blackman, hamming, hanning: Simplify handling of periodic window option
Mike Miller <mtmiller@octave.org>
parents:
20025
diff
changeset
|
52 N = m; |
19826
26fb4bfa4193
blackman, hamming, hanning: Add periodic window option (bug #43305)
Mike Miller <mtmiller@ieee.org>
parents:
19697
diff
changeset
|
53 case "symmetric" |
26fb4bfa4193
blackman, hamming, hanning: Add periodic window option (bug #43305)
Mike Miller <mtmiller@ieee.org>
parents:
19697
diff
changeset
|
54 ## Default option, same as no option specified. |
26fb4bfa4193
blackman, hamming, hanning: Add periodic window option (bug #43305)
Mike Miller <mtmiller@ieee.org>
parents:
19697
diff
changeset
|
55 otherwise |
26fb4bfa4193
blackman, hamming, hanning: Add periodic window option (bug #43305)
Mike Miller <mtmiller@ieee.org>
parents:
19697
diff
changeset
|
56 error ('blackman: window type must be either "periodic" or "symmetric"'); |
26fb4bfa4193
blackman, hamming, hanning: Add periodic window option (bug #43305)
Mike Miller <mtmiller@ieee.org>
parents:
19697
diff
changeset
|
57 endswitch |
3191 | 58 endif |
3426 | 59 |
3191 | 60 if (m == 1) |
61 c = 1; | |
62 else | |
20231
83792dd9bcc1
Use in-place operators in m-files where possible.
Rik <rik@octave.org>
parents:
20165
diff
changeset
|
63 m -= 1; |
20026
230c1e2a678d
blackman, hamming, hanning: Simplify handling of periodic window option
Mike Miller <mtmiller@octave.org>
parents:
20025
diff
changeset
|
64 k = (0 : m)' / N; |
3191 | 65 c = 0.42 - 0.5 * cos (2 * pi * k) + 0.08 * cos (4 * pi * k); |
66 endif | |
3426 | 67 |
3191 | 68 endfunction |
13068
542891ebfcdb
codesprint: 9 tests for bartlett.m
Andriy Shinkarchuck <adriano32.gnu@gmail.com>
parents:
11523
diff
changeset
|
69 |
14363
f3d52523cde1
Use Octave coding conventions in all m-file %!test blocks
Rik <octave@nomad.inbox5.com>
parents:
14138
diff
changeset
|
70 |
f3d52523cde1
Use Octave coding conventions in all m-file %!test blocks
Rik <octave@nomad.inbox5.com>
parents:
14138
diff
changeset
|
71 %!assert (blackman (1), 1) |
f3d52523cde1
Use Octave coding conventions in all m-file %!test blocks
Rik <octave@nomad.inbox5.com>
parents:
14138
diff
changeset
|
72 %!assert (blackman (2), zeros (2,1), 1e-6) |
20025
d209fbae38ae
Fix tests for bartlett, blackman, hamming, and hanning
Mike Miller <mtmiller@octave.org>
parents:
19826
diff
changeset
|
73 %!assert (blackman (15), flip (blackman (15)), 5*eps) |
d209fbae38ae
Fix tests for bartlett, blackman, hamming, and hanning
Mike Miller <mtmiller@octave.org>
parents:
19826
diff
changeset
|
74 %!assert (blackman (16), flip (blackman (16)), 5*eps) |
13068
542891ebfcdb
codesprint: 9 tests for bartlett.m
Andriy Shinkarchuck <adriano32.gnu@gmail.com>
parents:
11523
diff
changeset
|
75 %!test |
542891ebfcdb
codesprint: 9 tests for bartlett.m
Andriy Shinkarchuck <adriano32.gnu@gmail.com>
parents:
11523
diff
changeset
|
76 %! N = 9; |
542891ebfcdb
codesprint: 9 tests for bartlett.m
Andriy Shinkarchuck <adriano32.gnu@gmail.com>
parents:
11523
diff
changeset
|
77 %! A = blackman (N); |
20025
d209fbae38ae
Fix tests for bartlett, blackman, hamming, and hanning
Mike Miller <mtmiller@octave.org>
parents:
19826
diff
changeset
|
78 %! assert (A(ceil (N/2)), 1, 1e-6); |
d209fbae38ae
Fix tests for bartlett, blackman, hamming, and hanning
Mike Miller <mtmiller@octave.org>
parents:
19826
diff
changeset
|
79 %! assert ([A(1), A(length (A))], zeros (1,2), 1e-6); |
13068
542891ebfcdb
codesprint: 9 tests for bartlett.m
Andriy Shinkarchuck <adriano32.gnu@gmail.com>
parents:
11523
diff
changeset
|
80 |
21317
a4faec57f4c8
maint: remove semicolon after %!assert tests to follow Octave conventions.
Rik <rik@octave.org>
parents:
20852
diff
changeset
|
81 %!assert (blackman (15), blackman (15, "symmetric")) |
a4faec57f4c8
maint: remove semicolon after %!assert tests to follow Octave conventions.
Rik <rik@octave.org>
parents:
20852
diff
changeset
|
82 %!assert (blackman (16)(1:15), blackman (15, "periodic")) |
19826
26fb4bfa4193
blackman, hamming, hanning: Add periodic window option (bug #43305)
Mike Miller <mtmiller@ieee.org>
parents:
19697
diff
changeset
|
83 %!test |
26fb4bfa4193
blackman, hamming, hanning: Add periodic window option (bug #43305)
Mike Miller <mtmiller@ieee.org>
parents:
19697
diff
changeset
|
84 %! N = 16; |
26fb4bfa4193
blackman, hamming, hanning: Add periodic window option (bug #43305)
Mike Miller <mtmiller@ieee.org>
parents:
19697
diff
changeset
|
85 %! A = blackman (N, "periodic"); |
20025
d209fbae38ae
Fix tests for bartlett, blackman, hamming, and hanning
Mike Miller <mtmiller@octave.org>
parents:
19826
diff
changeset
|
86 %! assert (A(N/2 + 1), 1, 1e-6); |
19826
26fb4bfa4193
blackman, hamming, hanning: Add periodic window option (bug #43305)
Mike Miller <mtmiller@ieee.org>
parents:
19697
diff
changeset
|
87 |
14363
f3d52523cde1
Use Octave coding conventions in all m-file %!test blocks
Rik <octave@nomad.inbox5.com>
parents:
14138
diff
changeset
|
88 %!error blackman () |
f3d52523cde1
Use Octave coding conventions in all m-file %!test blocks
Rik <octave@nomad.inbox5.com>
parents:
14138
diff
changeset
|
89 %!error blackman (0.5) |
f3d52523cde1
Use Octave coding conventions in all m-file %!test blocks
Rik <octave@nomad.inbox5.com>
parents:
14138
diff
changeset
|
90 %!error blackman (-1) |
f3d52523cde1
Use Octave coding conventions in all m-file %!test blocks
Rik <octave@nomad.inbox5.com>
parents:
14138
diff
changeset
|
91 %!error blackman (ones (1,4)) |
21581
6fab85c1538f
maint: Follow Octave conventions for use of semicolon in BIST tests.
Rik <rik@octave.org>
parents:
21317
diff
changeset
|
92 %!error blackman (1, "invalid") |