root/applications/mpdm/TR2244/merger_iter_debug3.m @ 833

Revision 706, 1.1 kB (checked in by smidl, 15 years ago)

eol-native

  • Property svn:eol-style set to native
Line 
1itload('iter_cond_debug3.it');
2
3ndat = 1000; % check if true!
4
5prop_mean=zeros(3,40);
6prop_var=zeros(3,40);
7for it=0:39
8    si = num2str(it);
9    eval(['eff_w(it+1) = 1/(w' si '''*w' si ');']);
10    eval(['prop_mean(1:3,it+1) = Mpred_mean' si ';']);
11    eval(['prop_var(1:3,it+1) = Mpred_var' si '(1:3);']);
12end
13
14figure(2)
15subplot(1,2,1)
16hold off
17plot(eff_w)
18xlabel('iteration');
19title('Effective sample size');
20%title('Non-normalized importance weights')
21
22subplot(1,2,2);
23hold off
24plot(prop_mean');
25hold on
26plot(prop_mean'+2*sqrt(prop_var'),'.');
27plot(prop_mean'-2*sqrt(prop_var'),'.');
28
29%plot(ones(size(prop_mean,2),1)*[3 2],'--')
30xlabel('iteration');
31title('Mean value of the joint merger $\tilde{f}(\bm{x})$');
32
33% itload('../merger_iter_test.it');
34% XG = reshape(Grid(1,:),Npoints,Npoints);
35% YG = reshape(Grid(2,:),Npoints,Npoints);
36%
37% figure(2);
38% M1 = reshape(exp(Res1),Npoints,Npoints);
39% contour(XG,YG,M1,7);
40%
41% figure(3);
42% hold off
43% mm = max(max(Res2));
44% for i=1:size(Res2,1)
45%     M2 = reshape(Res2(i,:),Npoints,Npoints);
46%     contour(XG,YG,M2,[0:mm/7:mm]);
47%     hold on
48% end
49
Note: See TracBrowser for help on using the browser.