Mercurial > hg > ltpda

comparison m-toolbox/classes/@specwin/win_hft248d.m @ 0:f0afece42f48

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Import.
author Daniele Nicolodi <nicolodi@science.unitn.it>
date Wed, 23 Nov 2011 19:22:13 +0100
parents
children
comparison
equal deleted inserted replaced
-1:000000000000 0:f0afece42f48
1 % WIN_HFT248D returns HFT248D window, with N points.
2 % If mode == 'define', the window values will be empty and all the features
3 % If mode == 'build', the window values will be calculated
4 %
5 % M Hewitson 18-05-07
6 %
7 % $Id: win_hft248d.m,v 1.5 2011/05/25 07:30:08 mauro Exp $
8 %
9
10 function varargout = win_hft248d(w, mode, N)
11
12 switch lower(mode)
13
14 case 'build'
15 n_args = nargin;
16 if n_args < 3
17 N = w.len;
18 end
19 % Calculate the values of the window
20 z = [0:N-1]./N * 2 * pi;
21 varargout{1} = (1 - 1.985844164102 * cos(z) + ...
22 1.791176438506 * cos (2 * z) - 1.282075284005 * cos (3 * z) + ...
23 0.667777530266 * cos (4 * z) - 0.240160796576 * cos (5 * z) + ...
24 0.056656381764 * cos (6 * z) - 0.008134974479 * cos (7 * z) + ...
25 0.000624544650 * cos (8 * z) - 0.000019808998 * cos (9 * z) + ...
26 0.000000132974 * cos (10 * z));
27
28 case 'define'
29 % Make window struct
30 w.type = 'HFT248D';
31 w.alpha = 0;
32 w.psll = 248.4;
33 w.rov = 84.1;
34 w.nenbw = 5.6512;
35 w.w3db = 5.5567;
36 w.flatness = 0.0009;
37 w.len = N;
38 w.skip = 0;
39
40 varargout{1} = w;
41 end
42
43 % END