Mercurial > hg > ltpda

comparison m-toolbox/classes/@specwin/win_hft223d.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_HFT223D returns HFT223D 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_hft223d.m,v 1.5 2011/05/25 07:30:08 mauro Exp $
8 %
9
10 function varargout = win_hft223d(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.98298997309 * cos(z) + ...
22 1.75556083063 * cos (2 * z) - 1.19037717712 * cos (3 * z) + ...
23 0.56155440797 * cos (4 * z) - 0.17296769663 * cos (5 * z) + ...
24 0.03233247087 * cos (6 * z) - 0.00324954578 * cos (7 * z) + ...
25 0.00013801040 * cos (8 * z) - 0.00000132725 * cos (9 * z));
26 case 'define'
27 % Make window struct
28 w.type = 'HFT223D';
29 w.alpha = 0;
30 w.psll = 223;
31 w.rov = 83.3;
32 w.nenbw = 5.3888;
33 w.w3db = 5.3;
34 w.flatness = -0.0011;
35 w.len = N;
36 w.skip = 0;
37
38 varargout{1} = w;
39 end
40
41 % END