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
|
s_system_name = sprintf("s%s", system_name); # Name of sensitivity system
## Set up parameters
sim = eval([s_system_name, "_simpar;"]); # Simulation parameter
sym = eval([s_system_name, "_sympar;"]); # Parameter names
x0 = eval([s_system_name, "_state(par_0);"]); # Initial state
## Set up the free parameter list
free = [];
[n,m] = size(par_names);
for i = 1:n
p_name = deblank(par_names(i,:));
s_name = sprintf("%ss", p_name);
if struct_contains(sym, p_name)
i_p = eval(sprintf("sym.%s;", p_name));
if struct_contains(sym, s_name)
i_s = eval(sprintf("sym.%s;", s_name));
free_i = eval(sprintf("[%i,%i];", i_p, i_s));
free = [free; free_i];
else
printf("Sensitivity parameter %s does not exist: ignoring \
parameter %s\n", s_name, p_name);
endif
else
printf("Parameter %s does not exist: ignoring\n", p_name)
endif
endfor
[par,Par,Error,Y] = ppp_optimise(s_system_name,x0,par_0,sim,u,y_0,free,extras);
endfunction
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
>
|
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
|
s_system_name = sprintf("s%s", system_name); # Name of sensitivity system
## Set up parameters
sim = eval([s_system_name, "_simpar;"]); # Simulation parameter
sym = eval([s_system_name, "_sympar;"]); # Parameter names
x0 = eval([s_system_name, "_state(par_0);"]); # Initial state
# ## Set up the free parameter list
# free = [];
# [n,m] = size(par_names);
# for i = 1:n
# p_name = deblank(par_names(i,:));
# s_name = sprintf("%ss", p_name);
# if struct_contains(sym, p_name)
# i_p = eval(sprintf("sym.%s;", p_name));
# if struct_contains(sym, s_name)
# i_s = eval(sprintf("sym.%s;", s_name));
# free_i = eval(sprintf("[%i,%i];", i_p, i_s));
# free = [free; free_i];
# else
# printf("Sensitivity parameter %s does not exist: ignoring \
# parameter %s\n", s_name, p_name);
# endif
# else
# printf("Parameter %s does not exist: ignoring\n", p_name)
# endif
# endfor
# free
free = ppp_indices(par_names,sym);
[par,Par,Error,Y] = ppp_optimise(s_system_name,x0,par_0,sim,u,y_0,free,extras);
endfunction
|