Commit e4d67bc7 authored by FAGES Francois's avatar FAGES Francois
Browse files

RobustMonitoringAsInMaster

parent 2bc70f4f
%% Cell type:code id: tags:
```
load(ball.bc).
```
%% Cell type:code id: tags:
```
list_model.
```
%%%% Output: execute_result
MA(c)for _=[v]=>x.
g/c for v=>_.
MA(D*c)for v=[v]=>_.
present(x,x0).
present(y,3.5).
parameter(
x0 = 8.725,
D = -0.05,
K = 0.75,
g = 9.8,
c = 1.0
).
add_event(x<=0, c = -K*c).
%% Cell type:code id: tags:
```
list_ode.
```
%%%% Output: execute_result
[0] d(v)/dt= -1*g/c-D*c*v^2
[1] d(x)/dt=c*v
[2] d(y)/dt=0
%% Cell type:code id: tags:
```
numerical_simulation.
plot.
```
%%%% Output: display_data
%% Cell type:code id: tags:
```
robustness(F(x < h /\ F(x > h)), [x0, K, D], [h -> 3.5]).
```
%%%% Output: execute_result
Robustness degree: 0.884135
%% Cell type:code id: tags:
```
validity_domain(F(x < h /\ F(x > h))).
```
%%%% Output: execute_result
h<3.96644/\h> -1.1e-5
%% Cell type:code id: tags:
```
satisfaction_degree(F(x < h /\ F(x > h)), [h -> 3.5]).
```
%%%% Output: execute_result
1.466400
%% Cell type:code id: tags:
```
search_parameters(F(x < h /\ F(x > h)), [5 <= x0 <= 10, 0.5 <= K <= 1.0 , -0.1 <= D <= 0.0], [h -> 3.5], cmaes_stop_fitness: -0.5).
```
%%%% Output: execute_result
Time: 2.449 s
Stopping reason: Fitness: function value -5.72e-01 <= stopFitness (-5.00e-01)
Best satisfaction degree: 2.339104
[0] parameter(x0=9.555198250463423)
[1] parameter(K=0.8014006370116868)
[2] parameter(D= -0.058468364144163905)
%% Cell type:code id: tags:
```
numerical_simulation.
plot.
```
%%%% Output: display_data
%% Cell type:code id: tags:
```
satisfaction_degree(F(x < h /\ F(x > h)), [h -> 3.5]).
```
%%%% Output: execute_result
2.339120
%% Cell type:code id: tags:
```
robustness(F(x < h /\ F(x > h)), [x0, K, D], [h -> 3.5]).
```
%%%% Output: execute_result
Robustness degree: 0.930547
%% Cell type:code id: tags:
```
```
This source diff could not be displayed because it is too large. You can view the blob instead.
Examples of robustness computation and optimization in BIOCHAM-4
http://lifeware.inria.fr/biocham4
Inria Saclay Ile de France
The .ipynb files are BIOCHAM-4 notebooks directly executable on our server
(use shift-return to execute a cell)
The .bc files are BIOCHAM-4 files loaded by the notebooks
The reference manual and tutorials of BIOCHAM-4 are available at http://lifeware.inria.fr/biocham4
%% Cell type:code id: tags:
```
load(library:examples/cell_cycle/Tyson_1991.bc).
```
%% Cell type:code id: tags:
```
option(time:150).
```
%% Cell type:code id: tags:
```
parameter(k8=100).
```
%% Cell type:code id: tags:
```
numerical_simulation.
plot.
```
%%%% Output: display_data
%% Cell type:code id: tags:
```
satisfaction_degree(exists(max, exists(min, F(Cdc2-Cyclin~{p1} > max /\ F(Cdc2-Cyclin~{p1} < min /\ F(Cdc2-Cyclin~{p1} > max /\ F(Cdc2-Cyclin~{p1} < min)))) /\ max - min > amp)), [amp -> 0.19]).
```
%%%% Output: execute_result
0.974449
%% Cell type:code id: tags:
```
robustness(exists(max, exists(min, F(Cdc2-Cyclin~{p1} > max /\ F(Cdc2-Cyclin~{p1} < min /\ F(Cdc2-Cyclin~{p1} > max /\ F(Cdc2-Cyclin~{p1} < min)))) /\ max - min > amp)), [k4, k6], [amp -> 0.19]).
```
%%%% Output: execute_result
Robustness degree: 0.958877
%% Cell type:code id: tags:
```
parameter(k4=550, k6=2).
```
%% Cell type:code id: tags:
```
numerical_simulation. plot.
```
%%%% Output: display_data
%% Cell type:code id: tags:
```
satisfaction_degree(exists(max, exists(min, F(Cdc2-Cyclin~{p1} > max /\ F(Cdc2-Cyclin~{p1} < min /\ F(Cdc2-Cyclin~{p1} > max /\ F(Cdc2-Cyclin~{p1} < min)))) /\ max - min > amp)), [amp -> 0.19]).
```
%%%% Output: execute_result
1.009980
%% Cell type:code id: tags:
```
robustness(exists(max, exists(min, F(Cdc2-Cyclin~{p1} > max /\ F(Cdc2-Cyclin~{p1} < min /\ F(Cdc2-Cyclin~{p1} > max /\ F(Cdc2-Cyclin~{p1} < min)))) /\ max - min > amp)), [k4, k6], [amp -> 0.19], robustness_coeff_var: 0.1).
```
%%%% Output: execute_result
Robustness degree: 0.99358
%% Cell type:code id: tags:
```
```
%% Cell type:code id: tags:
```
load(library:examples/cell_cycle/Tyson_1991.bc).
```
%% Cell type:code id: tags:
```
option(time:150).
```
%% Cell type:code id: tags:
```
parameter(k8=100).
```
%% Cell type:code id: tags:
```
numerical_simulation.
plot.
```
%%%% Output: display_data
%% Cell type:code id: tags:
```
satisfaction_degree(exists(max, exists(min, F(Cdc2-Cyclin~{p1} > max /\ F(Cdc2-Cyclin~{p1} < min /\ F(Cdc2-Cyclin~{p1} > max /\ F(Cdc2-Cyclin~{p1} < min)))) /\ max - min > amp)), [amp -> 0.19]).
```
%%%% Output: execute_result
0.974450
%% Cell type:code id: tags:
```
seed(0).
robustness(exists(max, exists(min, F(Cdc2-Cyclin~{p1} > max /\ F(Cdc2-Cyclin~{p1} < min /\ F(Cdc2-Cyclin~{p1} > max /\ F(Cdc2-Cyclin~{p1} < min)))) /\ max - min > amp)), [k4, k6], [amp -> 0.19]).
```
%%%% Output: execute_result
Time: 6.363 s
Time: 7.494 s
Robustness degree: 0.978333
%% Cell type:code id: tags:
```
parameter(k4=550, k6=2).
```
%% Cell type:code id: tags:
```
numerical_simulation. plot.
```
%%%% Output: display_data
%% Cell type:code id: tags:
```
satisfaction_degree(exists(max, exists(min, F(Cdc2-Cyclin~{p1} > max /\ F(Cdc2-Cyclin~{p1} < min /\ F(Cdc2-Cyclin~{p1} > max /\ F(Cdc2-Cyclin~{p1} < min)))) /\ max - min > amp)), [amp -> 0.19]).
```
%%%% Output: execute_result
1.009960
%% Cell type:code id: tags:
```
seed(0).
robustness(exists(max, exists(min, F(Cdc2-Cyclin~{p1} > max /\ F(Cdc2-Cyclin~{p1} < min /\ F(Cdc2-Cyclin~{p1} > max /\ F(Cdc2-Cyclin~{p1} < min)))) /\ max - min > amp)), [k4, k6], [amp -> 0.19], robustness_coeff_var: 0.1).
```
%%%% Output: execute_result
Time: 8.757 s
Time: 10.956 s
Robustness degree: 0.996442
%% Cell type:code id: tags:
```
```
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment