Mercurial > hg > ltpda

comparison m-toolbox/test/test_ao_mcmc.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 %% Simple test:
2
3 x = 1*ones(1000,1);
4 y = 3*x;
5
6 x = ao(x);
7 y = ao(y);
8
9 x0 = 3.3;
10 model = smodel(plist('expression','A*x','xvar','x','params','A','xvals',[x.y]));
11
12 noise = ao(plist('waveform', 'noise', 'fs', 10, 'nsecs', 100,'sigma',0.1));
13
14 pl = plist('N',1000,...
15 'J',1,...
16 'cov',0.1,... % [4e-5 1e-5 1e-6 1e-8 1e-9]
17 'range',{[0.5 5.5]},...
18 'param',{'A'},...
19 'model',model,...
20 'noise',noise,...
21 'search',false,...
22 'x0',[x0],...
23 'simplex',false,...
24 'plot',1);
25
26 b = mcmc(x,y,pl);
27
28 plot(b.chain)