Breach doc: Simulations

Home       Download and install       Systems       Simulation       Temporal Logic       Parameter sets       Graphical User Interfaces      

The first thing one wants to do with a dynamical system is computing its evolution over time given some values for its initial conditions and various parameters. In Breach, this is done using the ComputeTraj command. It takes three arguments:

  • A system Sys, obviously corresponding to the dynamical system under study
  • A Parameter set P, which is a structure assigning values to parameters of the system
  • A time specification time, which indicates the time instants for which we want to compute the state of system.

Computing the nominal trajectory

The system structure Sys contains nominal (or default) values for the system parameters. To compute a trajectory corresponding to these nominal parameters, we can obtain the corresponding nominal parameter set by using the CreateParamSet routine. The example below shows how to do it for the Lorentz example (in Examples/Lorentz folder):

%% Init system and parameter set (create the Sys structure)
CreateSystem; 

%% Create nominal parameter set 
P0 = CreateParamSet(Sys); 

%% Compute nominal trajectories and output result for t in [0 10] with time step .01 
time = 0:.01:10;  
Pf = ComputeTraj(Sys, P0, [0 10);

To plot the result, one can use two different routines. SplotTraj plots the result in the phase space (in this case, in the space file:.). The code

SplotTraj(Ptraj)

produces:

images/SplotTraj_Ex0.png

SplotVar plots each variable of the simulation against time. Below is the output of

SplotVar(Ptraj)

images/SplotVar_Ex0.png

Home       Download and install       Systems       Simulation       Temporal Logic       Parameter sets       Graphical User Interfaces