Mercurial > hg > ltpda

comparison m-toolbox/test/test_noisegen.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 function test_noisegen()
2
3 % TEST_NOISEGEN tests the Franklin noise-generator algorithm.
4 %
5 %
6 % $Id: test_noisegen.m,v 1.11 2008/07/18 15:53:58 anneke Exp $
7 %
8
9
10
11 %% poles and zeros
12 pole = [pz(5,1) pz(1)];
13 zero = pz(3);
14 gain = 1;
15 pzm = pzmodel(gain, pole, zero);%% convert to get a and b coefficients
16
17 % parameter list for noise generation
18 pl = plist('nsecs', 100000,'fs', 10,'pzmodel', pzm, 'ndigits',50);
19
20 %% calling the noisegenerator
21 b = ao(pl);
22
23 %% Make LPSD
24 % Window function
25 w = specwin('Kaiser', 10, 150);
26 % parameter list for lpsd
27 lpl = plist('Kdes', 1000,'Kmin', 2,'Jdes', 1000,'Win', w,'order',0,'scale','ASD');
28
29 % LPSD
30 noise = lpsd(b, lpl); % matlab
31
32 rpl = plist('f1',1e-5,'f2',10,'nf',5e3);
33 %% Plot
34 % iplot(b)
35 % iplot(analyt)
36 plpl = plist('lengend','off');
37 iplot(noise,abs(resp(pzm,rpl)),plpl)
38
39 % plot(noise.hist)