octave_packages/m/signal/private/rectangle_sw.m

   1 ## Copyright (C) 1995-2012 Friedrich Leisch
   2 ##
   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
   7 ## the Free Software Foundation; either version 3 of the License, or (at
   8 ## your option) any later version.
   9 ##
  10 ## Octave is distributed in the hope that it will be useful, but
  11 ## WITHOUT ANY WARRANTY; without even the implied warranty of
  12 ## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  13 ## General Public License for more details.
  14 ##
  15 ## You should have received a copy of the GNU General Public License
  16 ## along with Octave; see the file COPYING.  If not, see
  17 ## <http://www.gnu.org/licenses/>.
  18
  19 ## -*- texinfo -*-
  20 ## @deftypefn {Function File} {} rectangle_sw (@var{n}, @var{b})
  21 ## Rectangular spectral window.  Subfunction used for spectral density
  22 ## estimation.
  23 ## @end deftypefn
  24
  25 ## Author: FL <Friedrich.Leisch@ci.tuwien.ac.at>
  26 ## Description: Rectangular spectral window
  27
  28 function retval = rectangle_sw (n, b)
  29
  30   if (nargin != 2)
  31     print_usage ();
  32   endif
  33
  34   retval = zeros (n, 1);
  35   retval(1) = 2 / b + 1;
  36
  37   l = (2:n)' - 1;
  38   l = 2 * pi * l / n;
  39
  40   retval(2:n) = sin ((2/b + 1) * l / 2) ./ sin (l / 2);
  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