Commit 319abb5c authored by SOLIMAN Sylvain's avatar SOLIMAN Sylvain
Browse files

even more new versions

parent c5991337
%% Cell type:code id: tags: %% Cell type:code id: tags:
   
``` ```
load(ball.bc). load(ball.bc).
``` ```
   
%% Cell type:code id: tags: %% Cell type:code id: tags:
   
``` ```
list_model. list_model.
``` ```
   
%%%% Output: execute_result %%%% Output: execute_result
   
MA(c)for _=[v]=>x. MA(c)for _=[v]=>x.
g/c for v=>_. g/c for v=>_.
MA(D*c)for v=[v]=>_. MA(D*c)for v=[v]=>_.
present(x,x0). present(x,x0).
present(y,3.5). present(y,3.5).
parameter( parameter(
x0 = 8.725, x0 = 8.725,
D = -0.05, D = -0.05,
K = 0.75, K = 0.75,
g = 9.8, g = 9.8,
c = 1.0 c = 1.0
). ).
add_event(x<=0, c = -K*c). add_event(x<=0, c = -K*c).
   
%% Cell type:code id: tags: %% Cell type:code id: tags:
   
``` ```
list_ode. list_ode.
``` ```
   
%%%% Output: execute_result %%%% Output: execute_result
   
[0] d(v)/dt= -1*g/c-D*c*v^2 [0] d(v)/dt= -1*g/c-D*c*v^2
[1] d(x)/dt=c*v [1] d(x)/dt=c*v
[2] d(y)/dt=0 [2] d(y)/dt=0
   
%% Cell type:code id: tags: %% Cell type:code id: tags:
   
``` ```
numerical_simulation. numerical_simulation.
plot. plot.
``` ```
   
%%%% Output: display_data %%%% Output: display_data
   
   
%% Cell type:code id: tags: %% Cell type:code id: tags:
   
``` ```
%timeout 180
```
%% Cell type:code id: tags:
```
seed(0). seed(0).
robustness(F(x < h /\ F(x > h)), [x0, K, D], [h -> 3.5]). robustness(F(x < h /\ F(x > h)), [x0, K, D], [h -> 3.5]).
``` ```
   
%%%% Output: execute_result %%%% Output: execute_result
   
Time: 129.124 s Time: 129.124 s
Robustness degree: 0.88414 Robustness degree: 0.88414
   
%% Cell type:code id: tags: %% Cell type:code id: tags:
   
``` ```
validity_domain(F(x < h /\ F(x > h))). validity_domain(F(x < h /\ F(x > h))).
``` ```
   
%%%% Output: execute_result %%%% Output: execute_result
   
h<3.96646/\h> -1.0e-6 h<3.96646/\h> -1.0e-6
   
%% Cell type:code id: tags: %% Cell type:code id: tags:
   
``` ```
satisfaction_degree(F(x < h /\ F(x > h)), [h -> 3.5]). satisfaction_degree(F(x < h /\ F(x > h)), [h -> 3.5]).
``` ```
   
%%%% Output: execute_result %%%% Output: execute_result
   
1.466440 1.466440
   
%% Cell type:code id: tags: %% Cell type:code id: tags:
   
``` ```
seed(0). seed(0).
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). 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 %%%% Output: execute_result
   
Time: 22.517 s Time: 22.517 s
Stopping reason: Fitness: function value -5.72e-01 <= stopFitness (-5.00e-01) Stopping reason: Fitness: function value -5.72e-01 <= stopFitness (-5.00e-01)
Best satisfaction degree: 2.339136 Best satisfaction degree: 2.339136
[0] parameter(x0=9.555198250463423) [0] parameter(x0=9.555198250463423)
[1] parameter(K=0.8014006370116868) [1] parameter(K=0.8014006370116868)
[2] parameter(D= -0.058468364144163905) [2] parameter(D= -0.058468364144163905)
   
%% Cell type:code id: tags: %% Cell type:code id: tags:
   
``` ```
numerical_simulation. numerical_simulation.
plot. plot.
``` ```
   
%%%% Output: display_data %%%% Output: display_data
   
   
%% Cell type:code id: tags: %% Cell type:code id: tags:
   
``` ```
satisfaction_degree(F(x < h /\ F(x > h)), [h -> 3.5]). satisfaction_degree(F(x < h /\ F(x > h)), [h -> 3.5]).
``` ```
   
%%%% Output: execute_result %%%% Output: execute_result
   
2.339120 2.339120
   
%% Cell type:code id: tags: %% Cell type:code id: tags:
   
``` ```
seed(0). seed(0).
robustness(F(x < h /\ F(x > h)), [x0, K, D], [h -> 3.5]). robustness(F(x < h /\ F(x > h)), [x0, K, D], [h -> 3.5]).
``` ```
   
%%%% Output: execute_result %%%% Output: execute_result
   
Time: 82.057 s Time: 82.057 s
Robustness degree: 0.930547 Robustness degree: 0.930547
   
%% Cell type:code id: tags: %% Cell type:code id: tags:
   
``` ```
``` ```
This source diff could not be displayed because it is too large. You can view the blob instead.
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