Mercurial > hg > ltpda

comparison m-toolbox/classes/@specwin/win_hft196d.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_HFT196D returns HFT196D 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_hft196d.m,v 1.5 2011/05/25 07:30:09 mauro Exp $
8 %
9
10 function varargout = win_hft196d(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.979280420 * cos (z) + ...
22 1.710288951 * cos (2 * z) - 1.081629853 * cos (3 * z) + ...
23 0.448734314 * cos (4 * z) - 0.112376628 * cos (5 * z) + ...
24 0.015122992 * cos (6 * z) - 0.000871252 * cos (7 * z) + ...
25 0.000011896 * cos (8 * z));
26
27 case 'define'
28 % Make window struct
29 w.type = 'HFT196D';
30 w.alpha = 0;
31 w.psll = 196.2;
32 w.rov = 82.3;
33 w.nenbw = 5.1134;
34 w.w3db = 5.0308;
35 w.flatness = 0.0013;
36 w.len = N;
37 w.skip = 0;
38
39 varargout{1} = w;
40 end
41
42 % END