home *** CD-ROM | disk | FTP | other *** search
/ OS/2 Shareware BBS: 10 Tools / 10-Tools.zip / octa21eb.zip / octave / SCRIPTS.ZIP / scripts.fat / signal / rect_sw.m < prev    next >
Text File  |  1999-04-29  |  1KB  |  66 lines
  1. ## Copyright (C) 1995, 1996, 1997  Friedrich Leisch
  2. ## 
  3. ## This program is free software; you can redistribute it and/or modify
  4. ## it under the terms of the GNU General Public License as published by
  5. ## the Free Software Foundation; either version 2, or (at your option)
  6. ## any later version.
  7. ## 
  8. ## This program is distributed in the hope that it will be useful, but
  9. ## WITHOUT ANY WARRANTY; without even the implied warranty of
  10. ## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  11. ## General Public License for more details. 
  12. ## 
  13. ## You should have received a copy of the GNU General Public License
  14. ## along with this file.  If not, write to the Free Software Foundation,
  15. ## 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA.
  16.  
  17. ## usage:  retval = rect_sw (n, b)
  18. ##
  19. ## Rectangular spectral window. Subfunction used for spectral density
  20. ## estimation.
  21.   
  22. ## Author:  FL <Friedrich.Leisch@ci.tuwien.ac.at>
  23. ## Description:  Rectangular spectral window
  24.  
  25. function retval = rect_sw (n, b)
  26.   
  27.   retval = zeros (n, 1);
  28.   retval(1) = 2 / b + 1;
  29.  
  30.   l = (2:n)' - 1;
  31.   l = 2 * pi * l / n;
  32.  
  33.   retval(2:n) = sin( (2/b + 1) * l / 2 ) ./ sin (l / 2);
  34.     
  35. endfunction
  36.  
  37.  
  38.  
  39.  
  40.  
  41.  
  42.  
  43.  
  44.  
  45.  
  46.  
  47.  
  48.  
  49.  
  50.  
  51.  
  52.  
  53.  
  54.  
  55.  
  56.  
  57.  
  58.  
  59.  
  60.  
  61.  
  62.  
  63.  
  64.  
  65.  
  66.