root/applications/mpdm/matlab/merger_iter_debug.m @ 1414

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

eol-native

  • Property svn:eol-style set to native
Line 
1itload('../merger_debug.it');
2
3ndat = 1000; % check if true!
4
5XL = [-1 4];
6YL= XL;
7
8iters = 0:9
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,:),w_is_' 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
61for it=0:25
62    si = num2str(it);
63    eval(['std_w(it+1) = std(w_is_' si '*200);']);
64    eval(['prop_mean(1:2,it+1) = Mpred_mean' si ';']);
65end
66
67figure(2)
68subplot(1,3,1)
69hold off
70plot(std_w)
71xlabel('iteration');
72ylabel('standard deviation');
73title('Non-normalized importance weights')
74
75subplot(1,3,2);
76hist (w_is_15*200);
77title('Histogram of importance weights');
78
79subplot(1,3,3);
80hold off
81plot(prop_mean');
82hold on
83plot(ones(size(prop_mean,2),1)*[3 2],'--')
84xlabel('iteration');
85title('Mean value of the merged density');
86
87% itload('../merger_iter_test.it');
88% XG = reshape(Grid(1,:),Npoints,Npoints);
89% YG = reshape(Grid(2,:),Npoints,Npoints);
90%
91% figure(2);
92% M1 = reshape(exp(Res1),Npoints,Npoints);
93% contour(XG,YG,M1,7);
94%
95% figure(3);
96% hold off
97% mm = max(max(Res2));
98% for i=1:size(Res2,1)
99%     M2 = reshape(Res2(i,:),Npoints,Npoints);
100%     contour(XG,YG,M2,[0:mm/7:mm]);
101%     hold on
102% end
103
Note: See TracBrowser for help on using the browser.