DTI / ISEA / LSM / Titolo principale della presentazione CACSD with Python Prof. Roberto Bucher SUPSI - ISEA – SMT Rapperswil 28.01.2014 SGA Meeting und Preisübergabe Förderpreis 2013 DTI / ISEA / LSM / Titolo principale della presentazione Outline • • • • • • • Introduction CACSD environment Python Examples Comparison Matlab – Scicoslab -Python Code generation Conclusions DTI / ISEA / LSM / Titolo principale della presentazione Introduction • Matlab / Simulink are “de facto” THE TOOLS in education • But... – What happens wenn our (ex-) students have to develop advanced control systems in little companies? – Mathworks licenses are too expensive for most little companies! DTI / ISEA / LSM / Titolo principale della presentazione CACSD - Matlab alternatives CACSD Tool Matrix/Algebra Hybrid Simulator Discrete Events Code generation for emb. Sys. Matlab ✔ Matlab ✔ Simulink ✔ Stateflow ✔ RTW Scicoslab/NSP ✔ Scilab ✔ Scicos ✔ SMCube1 Scilab ✔ Scilab ✔ xCos - (-)3 Octave ✔ Matlab like - - - 1 Commercial Free + Commercial (RTAICodeGen, E4Coder) 3 Modified RTAICodeGen, Scilab 6.x 2 ✔ 2 DTI / ISEA / LSM / Titolo principale della presentazione CACSD - Matlab alternatives CACSD Tool Matrix/Algebra Hybrid Simulator Discrete Events Code generation for emb. Sys. Matlab ✔ Matlab ✔ Simulink ✔ Stateflow ✔ RTW Scicoslab/NSP ✔ Scilab ✔ Scicos ✔ SMCube1 Scilab ✔ Scilab ✔ xCos - (-)3 Octave ✔ Matlab like - - - Python ✔ (✔) Simpy4 (✔) 1 Commercial Free + Commercial (RTAICodeGen, E4Coder) 3 Modified RTAICodeGen, Scilab 6.x 4 No graphical editor 2 ✔ 2 DTI / ISEA / LSM / Titolo principale della presentazione Python DTI / ISEA / LSM / Titolo principale della presentazione Python • Python is a general-purpose, high-level, object oriented programming language similar to Matlab, Scicoslab, Scilab and Octave. • Its implementation was started in December 1989 (Guido van Rossum at Centrum Wiskunde & Informatica in the Netherlands) • Python is often used as a scripting language, but is also used in a wide range of non-scripting contexts. • Distributed under the Python Software Foundation License (PSFL) similar to BSD license. • Very large community of developers! DTI / ISEA / LSM / Titolo principale della presentazione Python • Numpy numeric array • Scipy scientific functions – FFT, Integrate, Interpolate, Optimize, Signal, Sparse, Statistic • Matplotlib Plot functions • PyQt Integration with Qt Libraries • Python-Qwt5 Integration with Qt Widgets • ... • Python Control Systems Library (python-control) DTI / ISEA / LSM / Titolo principale della presentazione Python for Matlab users • The python scripting language is similar to Matlab – http://wiki.scipy.org/NumPy_for_Matlab_Users • Main differencies: – Indexing (Python is 0 based!) – Array/Matrix distinction (linear algebra operations!) – Importing of functions • from np.linalg import * • from control.matlab import ss, tf DTI / ISEA / LSM / Titolo principale della presentazione Python Control Systems Library (python-control) • • 2009-2011 – Richard Murray1 – Caltech • Control-0.3a – Continous Time systems – SISO Systems 2011 – Roberto Bucher (SUPSI) • yottalab.py – MIMO systems – Discrete Time systems (Field Tsamp!) – Set of functions to handle discrete time systems. – Design of complex discrete controllers was possible! – yottalab.py distributed with control toolbox as separate toolbox 1 Aström and Murray – Feedback Systems DTI / ISEA / LSM / Titolo principale della presentazione Python Control Systems Library (python-control) • Situation 2014 – Richard Murray (Caltech) • Control-0.6c – Some ideas and functions migrated from yottalab to python-control – Discrete Time systems (dt: sampling time) – Functions related to discrete time systems – Roberto Bucher (SUPSI) • Supsictrl – yottalab.py » Few functions not cleanly implemented by Richard Murray » Some additional functions to simplify the controller design (observer, reduced order observer, controllers in compact form, anti windup design, etc.) – RCPblk.py » Toolbox for Rapid Controller Prototyping, hybrid simulation and Code generation for RT systems DTI / ISEA / LSM / Titolo principale della presentazione Python for control – Other related projects • Ryan Krauss (Southern Illinois University Edwardsville) – controls.py • Some routines related to transfer functions (continous and discrete) only for education purposes. • Python control systems, with only one developper, Rafael G. Martins (only for education purposes) Stephen Boyd and Lieven Vandenberghe are now mostly developing optimization libraries, but a few applications specific to control systems can be found on their page • DTI / ISEA / LSM / Titolo principale della presentazione Python for Control • Files – Ipython – Numpy – Scipy – Matplotlib – Avventi Slycot – Control-0.6c – Yottalab.py – RCPblk.py – pyBlockEdit DTI / ISEA / LSM / Titolo principale della presentazione Python for Control • Files – Ipython – Numpy – Scipy – Matplotlib – Avventi Slycot – Control-0.6c – Yottalab.py – RCPblk.py – pyBlockEdit DTI / ISEA / LSM / Titolo principale della presentazione Python for Control • Files – Ipython – Numpy – Scipy – Matplotlib – Avventi Slycot – Control-0.6c – Yottalab.py – RCPblk.py – pyBlockEdit DTI / ISEA / LSM / Titolo principale della presentazione Python for Control • Files – Ipython – Numpy – Scipy – Matplotlib – Avventi Slycot – Control-0.6c – Yottalab.py – RCPblk.py – pyBlockEdit DTI / ISEA / LSM / Titolo principale della presentazione Comparison Matlab / Scicoslab / Python DTI / ISEA / LSM / Titolo principale della presentazione Matlab + Simulink DTI / ISEA / LSM / Titolo principale della presentazione Matlab + Simulink DTI / ISEA / LSM / Titolo principale della presentazione Scicoslab + Modelica DTI / ISEA / LSM / Titolo principale della presentazione Scicoslab + Modelica DTI / ISEA / LSM / Titolo principale della presentazione Scicoslab + Modelica DTI / ISEA / LSM / Titolo principale della presentazione Python Demo 1 DTI / ISEA / LSM / Titolo principale della presentazione Python + CodeGen + Simulation fname = 'disks' genCode(fname,4,0.01,[sq1, Plant, contr, prnt]) genMake(fname) DTI / ISEA / LSM / Titolo principale della presentazione Python + CodeGen + Simulation fname = 'disks' genCode(fname,4,0.01,[sq1, Plant, contr, prnt]) genMake(fname) DTI / ISEA / LSM / Titolo principale della presentazione Python + CodeGen + Simulation fname = 'disks' genCode(fname,4,0.01,[sq1, Plant, contr, prnt]) genMake(fname) DTI / ISEA / LSM / Titolo principale della presentazione Code generation • • Find the right block sequence Generate the C code with 3 functions: – Init • Initialization – ISR • Outputs update • States update – Termination DTI / ISEA / LSM / Titolo principale della presentazione Code generation Standard C-Code DsPIC board Linux ARM PC with Linux RTAI or PREEMPT RT DTI / ISEA / LSM / Titolo principale della presentazione Python Demo 2 DTI / ISEA / LSM / Titolo principale della presentazione Real-time execution • • • State feedback controller with integral part Code generated from Python as PREEMPT_RT executable 10 ms Sampling time DTI / ISEA / LSM / Titolo principale della presentazione Real-time execution • • • State feedback controller with integral part Code generated from Python as PREEMPT_RT executable 10 ms Sampling time DTI / ISEA / LSM / Titolo principale della presentazione Real-time execution • • • State feedback controller with integral part Code generated from Python as PREEMPT_RT executable 10 ms Sampling time DTI / ISEA / LSM / Titolo principale della presentazione Conclusions • Which alternatives to Matlab / Simulink / RTW? • Python for control is promising, but not enough mature yet • Python control is an education toolbox, but not ready for industrial development • But Freedom to develop what we want and need • Urgent – Rewrite some functions (ex. Pole placement for MIMO system gives bad results) – Graphical editor for block diagram! DTI / ISEA / LSM / Titolo principale della presentazione Questions?
© Copyright 2024 Paperzz