root/applications/mpdm/TR2244/merger_iter_debug.m @ 1305

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

eol-native

  • Property svn:eol-style set to native
Line 
1itload('iter_cond_debug.it');
2
3ndat = 1000; % check if true!
4
5XL = [-3 3];
6YL= XL;
7
8iters =[];% [0 2 39]
9
10figure(1);
11niters = length(iters);
12noi = 1;
13for it=iters
14%    figure(it+1)
15    subplot(5,niters,noi);
16    si = num2str(it);
17    eval(['contour_2(Smp' si '(1,:),Smp' si '(2,:),exp(Mpdf' si '''))']);
18    if it==iters(1),
19    ylabel('Proposal density ');
20    end
21    set(gca,'XLim',XL);
22    set(gca,'YLim',YL);
23       
24    subplot(5,niters,noi+niters);
25    eval(['contour_2(Smp' si '(1,:),Smp' si '(2,:),exp(lW' si '(1,:)))']);
26    if it==iters(1),
27    ylabel('First source');
28    end
29    set(gca,'XLim',XL);
30    set(gca,'YLim',YL);
31   
32    subplot(5,niters,noi+2*niters);
33    eval(['contour_2(Smp' si '(1,:),Smp' si '(2,:),exp(lW' si '(2,:)))']);
34    if it==iters(1),
35    ylabel('Second source');
36    end
37    set(gca,'XLim',XL);
38    set(gca,'YLim',YL);
39   
40   
41    subplot(5,niters,noi+3*niters);
42    hold off
43    eval(['contour_2(Smp' si '(1,:),Smp' si '(2,:),w' si ''')']);
44    if it==iters(1),
45    ylabel('Merged density');
46    end
47    set(gca,'XLim',XL);
48    set(gca,'YLim',YL);
49   
50    subplot(5,niters,noi+4*niters);
51    hold off
52    eval(['contour_2(Smp' si '(1,:),Smp' si '(2,:),exp(lw_m' si '''))']);
53    if it==iters(1),
54    ylabel('Importance density');
55    end
56    set(gca,'XLim',XL);
57    set(gca,'YLim',YL);
58    noi = noi+1;
59end
60
61prop_mean=zeros(2,40);
62prop_var=zeros(2,40);
63for it=0:39
64    si = num2str(it);
65    eval(['eff_w(it+1) = 1/(w' si '''*w' si ');']);
66    eval(['prop_mean(1:2,it+1) = Mpred_mean' si ';']);
67    eval(['prop_var(1:2,it+1) = Mpred_var' si '(1:2);']);
68end
69
70figure(2)
71subplot(1,2,1)
72hold off
73plot(eff_w)
74xlabel('iteration');
75ylabel('Effective sample size');
76%title('Non-normalized importance weights')
77
78subplot(1,2,2);
79hold off
80plot(prop_mean');
81hold on
82plot(prop_mean'+2*sqrt(prop_var'),'.');
83plot(prop_mean'-2*sqrt(prop_var'),'.');
84
85%plot(ones(size(prop_mean,2),1)*[3 2],'--')
86xlabel('iteration');
87title('Mean value of the joint merger $\tilde{f}(\bm{x})$');
88
89% itload('../merger_iter_test.it');
90% XG = reshape(Grid(1,:),Npoints,Npoints);
91% YG = reshape(Grid(2,:),Npoints,Npoints);
92%
93% figure(2);
94% M1 = reshape(exp(Res1),Npoints,Npoints);
95% contour(XG,YG,M1,7);
96%
97% figure(3);
98% hold off
99% mm = max(max(Res2));
100% for i=1:size(Res2,1)
101%     M2 = reshape(Res2(i,:),Npoints,Npoints);
102%     contour(XG,YG,M2,[0:mm/7:mm]);
103%     hold on
104% end
105
Note: See TracBrowser for help on using the browser.