octave_packages/signal-1.1.3/bohmanwin.m

   1 ## Copyright (C) 2007 Sylvain Pelissier <sylvain.pelissier@gmail.com>
   2 ##
   3 ## This program is free software; you can redistribute it and/or modify it under
   4 ## the terms of the GNU General Public License as published by the Free Software
   5 ## Foundation; either version 3 of the License, or (at your option) any later
   6 ## version.
   7 ##
   8 ## This program is distributed in the hope that it will be useful, but WITHOUT
   9 ## ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
  10 ## FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more
  11 ## details.
  12 ##
  13 ## You should have received a copy of the GNU General Public License along with
  14 ## this program; if not, see <http://www.gnu.org/licenses/>.
  15
  16 ## -*- texinfo -*-
  17 ## @deftypefn {Function File} {[@var{w}] =} bohmanwin(@var{L})
  18 ## Compute the Bohman window of lenght L.
  19 ## @seealso{rectwin,  bartlett}
  20 ## @end deftypefn
  21
  22 function [w] = bohmanwin(L)
  23   if (nargin < 1)
  24     print_usage
  25   elseif(! isscalar(L))
  26     error("L must be a number");
  27   elseif(L < 0)
  28     error('L must be positive');
  29   end
  30
  31   if(L ~= floor(L))
  32     L = round(L);
  33     warning('L rounded to the nearest integer.');
  34   end
  35
  36   if(L == 0)
  37     w = [];
  38
  39   elseif(L == 1)
  40     w = 1;
  41
  42   else
  43     N = L-1;
  44     n = -N/2:N/2;
  45
  46     w = (1-2.*abs(n)./N).*cos(2.*pi.*abs(n)./N) + (1./pi).*sin(2.*pi.*abs(n)./N);
  47     w(1) = 0;
  48     w(length(w))=0;
  49     w = w';
  50   end
  51 endfunction