-
Notifications
You must be signed in to change notification settings - Fork 1
/
Copy pathcreateParamSet_16Apr20.m
119 lines (115 loc) · 2.4 KB
/
createParamSet_16Apr20.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
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
function [] = createParamSet_16Apr20(tstep_abc , date_abc)
t_curr = tstep_abc;
date = date_abc;
%% Load parameters
paramDir = [pwd ,'/Params/'];
%paramSetMatrix = load([paramDir,'paramSets_calib_' , date , '_' , num2str(t_curr) , '.dat']);
pIdx = load([paramDir,'pIdx_calib_' , date , '_0.dat']);
% paramSetNew = [0.80; ... % t=0, test1
% 12.2; ...
% 25.0; ...
% 25.0; ...
% 12.0; ...
% 0.22; 0.50; ...
% 9.0; ...
% 26.0; ...
% 19.0; ...
% 14.0; ...
% 10.0; ...
% 0.58; 0.3; ...
% 0.25; 0.4921; ...
% 0.8; 1.1; ...
% 79.6; ...
% 78.6; ...
% 56.2; ...
% 54.6; ...
% 52.7; ...
% 0.01; 0.5173; ...
% 0.8; 2.0; ...
% 1.4; 2.0; ...
% 1.3; 2.0; ...
% 1.0; 5.0; ...
% 2.60; 1.41; ...
% 1.8; 0.97; ...
% 1.4; 0.2; ...
% 1.0; 0.2; ...
% 0.000575; ...
% 3.5];
% paramSetNew = [0.80; ... % t=0, test2
% 12.2; ...
% 25.0; ...
% 25.0; ...
% 12.0; ...
% 0.22; 0.50; ...
% 9.0; ...
% 26.0; ...
% 19.0; ...
% 14.0; ...
% 10.0; ...
% 0.58; 0.3; ...
% 0.25; 0.4921; ...
% 0.8; 1.1; ...
% 79.6*0.89; ...
% 78.6*0.89; ...
% 56.2; ...
% 54.6; ...
% 52.7; ...
% 0.01; 0.5173; ...
% 0.8; 2.0; ...
% 1.4; 2.0; ...
% 1.3; 2.0; ...
% 1.0; 5.0; ...
% 2.60; 1.41; ...
% 1.8; 0.97; ...
% 1.4; 0.2; ...
% 1.0; 0.2; ...
% 0.000575; ...
% 3.5];
% paramSetNew = [0.80; ... % t=0, test3 (round 1 calib)
% 12.2; ...
% 25.0; ...
% 25.0; ...
% 12.0; ...
% 0.22; 0.50; ...
% 9.0; ...
% 26.0; ...
% 19.0; ...
% 14.0; ...
% 10.0; ...
% 0.58; 0.3; ...
% 0.25; 0.4921; ...
% 70.844; ...
% 69.954; ...
% 56.2; ...
% 54.6; ...
% 52.7; ...
% 0.01; 0.5173; ...
% 2.60; 1.41; ...
% 0.000575; ...
% 3.5];
paramSetNew = [1.0; ... % 22Apr20, t=18 sample 1
12.0; ...
30.0; ...
9.0; ...
9.0; ...
0.1; 0.2; ...
11.0; ...
35.0; ...
12.0; ...
7.0; ...
8.0; ...
0.5; 0.3; ...
0.25; 0.5; ...
70.0; ...
80.0; ...
58.0; ...
45.0; ...
53; ...
0.08; 0.7; ...
1.7; 1.2; ...
0.0006; ...
2.5];
%% Save pIdx and parameter set
file = ['paramSets_calib_' , date , '_' , num2str(t_curr) , '_22Apr20_t18_s1' , '.dat'];
paramDir = [pwd , '/Params/'];
csvwrite([paramDir, file] , paramSetNew)