forked from LionelBenoit/Local-rainfall-model
-
Notifications
You must be signed in to change notification settings - Fork 0
/
plot_result_metropolis.m
76 lines (50 loc) · 1.1 KB
/
plot_result_metropolis.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
function [] = plot_result_metropolis( V_M, num_fig )
[sx,sy]=size(V_M);
figure(num_fig)
clf;
subplot(4,6,1)
hist(V_M(:,1),50)
subplot(4,6,7)
plot(V_M(:,1),1:length(V_M(:,1)))
subplot(4,6,2)
hist(V_M(:,2),50)
subplot(4,6,8)
plot(V_M(:,2),1:length(V_M(:,2)))
subplot(4,6,3)
hist(V_M(:,3),50)
subplot(4,6,9)
plot(V_M(:,3),1:length(V_M(:,3)))
subplot(4,6,4)
hist(V_M(:,4),50)
subplot(4,6,10)
plot(V_M(:,4),1:length(V_M(:,4)))
subplot(4,6,5)
hist(V_M(:,5),50)
subplot(4,6,11)
plot(V_M(:,5),1:length(V_M(:,5)))
subplot(4,6,6)
hist(V_M(:,6),50)
subplot(4,6,12)
plot(V_M(:,6),1:length(V_M(:,6)))
%---
subplot(4,6,13)
hist(V_M(:,7),50)
subplot(4,6,19)
plot(V_M(:,7),1:length(V_M(:,7)))
subplot(4,6,14)
hist(V_M(:,8),50)
subplot(4,6,20)
plot(V_M(:,8),1:length(V_M(:,8)))
subplot(4,6,15)
hist(V_M(:,9),50)
subplot(4,6,21)
plot(V_M(:,9),1:length(V_M(:,9)))
subplot(4,6,16)
hist(V_M(:,10),50)
subplot(4,6,22)
plot(V_M(:,10),1:length(V_M(:,10)))
subplot(4,6,17)
hist(V_M(:,11),50)
subplot(4,6,23)
plot(V_M(:,11),1:length(V_M(:,11)))
end