## [d4a496]: inst / nuttallwin.m  Maximize  Restore  History

### 49 lines (41 with data), 1.4 kB

 ``` 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48``` ```## Copyright (C) 2007 Sylvain Pelissier ## ## This program is free software; you can redistribute it and/or modify it under ## the terms of the GNU General Public License as published by the Free Software ## Foundation; either version 3 of the License, or (at your option) any later ## version. ## ## This program is distributed in the hope that it will be useful, but WITHOUT ## ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or ## FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more ## details. ## ## You should have received a copy of the GNU General Public License along with ## this program; if not, see . ## -*- texinfo -*- ## @deftypefn {Function File} {[@var{w}] =} nuttallwin (@var{m}) ## Compute the Blackman-Harris window defined by Nuttall of length @var{m}. ## @seealso{blackman, blackmanharris} ## @end deftypefn function [w] = nuttallwin(m) if (nargin != 1); print_usage; endif if(m < 0) error('M must be positive'); endif if(m ~= floor(m)) m = round(m); warning('M rounded to the nearest integer.'); endif if (m == 1) w = 1; else N = m - 1; a0 = 0.355768; a1 = 0.487396; a2 = 0.144232; a3 = 0.012604; n = -N/2:N/2; w = a0 + a1.*cos(2.*pi.*n./N) + a2.*cos(4.*pi.*n./N) + a3.*cos(6.*pi.*n./N); w = w'; endif endfunction ```