Mercurial > hg > ltpda

comparison m-toolbox/test/lisovfit/test_lisovfit.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 mc
2
3 %%
4 % make ao of desired response
5 apl = plist('filename', 'mdc2_1_tf11_np.dat', 'type', 'fsdata')
6 a = ao(apl);
7
8 % now make LISO fit file
9
10 % start with a guess pzmodel
11
12 z0(1) = pz(1e-6);
13 zl(1) = pz(1e-07);
14 zu(1) = pz(1e-05);
15
16 p0(1) = pz(400e-9);
17 pl(1) = pz(1e-07);
18 pu(1) = pz(1e-05);
19
20 p0(2) = pz(1e-6);
21 pl(2) = pz(1e-07);
22 pu(2) = pz(1e-05);
23
24 p0(3) = pz(30e-6);
25 pl(3) = pz(1e-06);
26 pu(3) = pz(1e-4);
27
28 z0(2) = pz(50e-6);
29 zl(2) = pz(1e-06);
30 zu(2) = pz(1e-4);
31
32 z0(3) = pz(50e-6);
33 zl(3) = pz(1e-06);
34 zu(3) = pz(1e-4);
35
36 z0(4) = pz(50e-6);
37 zl(4) = pz(1e-06);
38 zu(4) = pz(1e-4);
39
40 z0(5) = pz(50e-6);
41 zl(5) = pz(1e-06);
42 zu(5) = pz(1e-4);
43
44 p0(4) = pz(7e-3);
45 pl(4) = pz(1e-4);
46 pu(4) = pz(1e-1);
47
48 p0(5) = pz(7e-3);
49 pl(5) = pz(1e-4);
50 pu(5) = pz(1e-1);
51
52 p0(6) = pz(40e-3);
53 pl(6) = pz(1e-4);
54 pu(6) = pz(1);
55
56 p0(7) = pz(40e-3);
57 pl(7) = pz(1e-4);
58 pu(7) = pz(1);
59
60 z0(6) = pz(0.5);
61 zl(6) = pz(0.01);
62 zu(6) = pz(1);
63
64 p0(8) = pz(0.1);
65 pl(8) = pz(0.01);
66 pu(8) = pz(10);
67
68 z0(6) = pz(0.5);
69 zl(6) = pz(0.01);
70 zu(6) = pz(1);
71
72
73 pzm0 = pzmodel(1e-10, p0, z0);
74 pzml = pzmodel(1e-12, pl, zl);
75 pzmu = pzmodel(1e-8, pu, zu);
76
77 resp(pzm0,pzml,pzmu)
78
79
80 %%
81 fpl = plist('PZM0', pzm0, 'PZML', pzml, 'PZMU', pzmu, 'delay', [0 1 10], 'f1', 1e-4, 'f2', 1, 'nf', 100);
82
83 pzmfit = lisovfit(a, fpl);
84
85 %%
86
87 f0 = resp(pzm0);
88 rfit = resp(pzmfit);
89
90 iplot(f0, rfit)
91