Mercurial > hg > octave-nkf
annotate scripts/signal/private/triangle_sw.m @ 14138:72c96de7a403 stable
maint: update copyright notices for 2012
| author | John W. Eaton <jwe@octave.org> |
|---|---|
| date | Mon, 02 Jan 2012 14:25:41 -0500 |
| parents | 0fed4935de94 |
| children | 5d3a684236b0 |
| rev | line source |
|---|---|
|
14138
72c96de7a403
maint: update copyright notices for 2012
John W. Eaton <jwe@octave.org>
parents:
12579
diff
changeset
|
1 ## Copyright (C) 1995-2012 Friedrich Leisch |
| 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 -*- |
| 20 ## @deftypefn {Function File} {} triangle_sw (@var{n}, @var{b}) | |
| 3191 | 21 ## Triangular spectral window. Subfunction used for spectral density |
| 22 ## estimation. | |
| 3449 | 23 ## @end deftypefn |
| 3426 | 24 |
| 3457 | 25 ## Author: FL <Friedrich.Leisch@ci.tuwien.ac.at> |
| 26 ## Description: Triangular spectral window | |
| 3191 | 27 |
| 28 function retval = triangle_sw (n, b) | |
| 3426 | 29 |
| 7125 | 30 if (nargin != 2) |
| 31 print_usage (); | |
| 32 endif | |
| 33 | |
| 3191 | 34 retval = zeros(n,1); |
| 35 retval(1) = 1 / b; | |
| 36 | |
| 37 l = (2:n)' - 1; | |
| 38 l = 2 * pi * l / n; | |
| 3426 | 39 |
| 3457 | 40 retval(2:n) = b * (sin (l / (2*b)) ./ sin (l / 2)).^2; |
| 3191 | 41 |
| 42 endfunction | |
| 43 | |
| 44 | |
| 45 | |
| 46 | |
| 47 | |
| 48 | |
| 49 | |
| 50 | |
| 51 | |
| 52 | |
| 53 | |
| 54 | |
| 55 | |
| 56 | |
| 57 | |
| 58 | |
| 59 | |
| 60 | |
| 61 | |
| 62 | |
| 63 | |
| 64 | |
| 65 | |
| 66 | |
| 67 | |
| 68 | |
| 69 | |
| 70 | |
| 71 | |
| 72 |