Numerical Analysis and Reachability
|
|
- Derek Stokes
- 5 years ago
- Views:
Transcription
1 Numerical Analysis and Reachability Hybrid Systems: Reachability Goran Frehse February 22, 2017 Ph.D. Course: Hybrid Systems Ph.D. in Information Technology Dipartimento di Elettronica, Informazione e Bioingegneria (DEIB) Politecnico di Milano
2 Overview Numerical Analysis Running Example Solving ODEs Computing Trajectories and Jumps Reachability with Piecewise Affine Dynamics Reachability with Nonlinear Dynamics 2
3 From Numbers to Sets Symbolic state-space exploration can be seen as a generalization of numerical analysis from numbers to sets. Look at numerical analysis to understand the principles see the difficulties 3
4 Running Example: Ball on String L L x 0 m F s x 0 m F g F g (a) extension (b) freefall 4
5 Equations of Motion dynamics in freefall when x 0, with mass m, mẍ = F g = mg. dynamics in extension when x 0, with spring constant k, damping factor d, mẍ = F g + F s = mg kx dẋ. transition when x = L, collision factor c [0, 1], ẋ = cẋ. 5
6 Equations of Motion as 1st order ODEs auxiliary variable v = ẋ, so v = ẍ. dynamics in freefall when x 0, ) ( ) ( ) ( ) (ẋ 0 1 x 0 = + v 0 0 v g dynamics in extension when x 0, ) ( ) ( ) ( ) (ẋ 0 1 x 0 = + v v g k m d m 6
7 Jumps as 1st order Difference Equations transition when x = L, collision factor c [0, 1], ẋ = cẋ. difference equation using auxiliary variable v = ẋ, denoting by x the target value of x: ( ) ( ) ( ) x 1 0 x = v 0 c v 7
8 Hybrid automaton model clip from SpaceEx Model Editor 1 1 G. Frehse, C. L. Guernic, A. Donzé, R. Ray, O. Lebeltel, R. Ripado, A. Girard, T. Dang, and O. Maler, SpaceEx: Scalable verification of hybrid systems, in CAV 11, ser. LNCS, Springer,
9 Behavior over Time 1 x 2 position x 0 x 1 x 3 x 4 x 5 1 x t velocity v 5 0 v 0 v 1 v 2 v 2 v 4 v 5 5 v t 9
10 Behavior in the State Space ξ 0 (δ 0 ) = x 1 5 ξ 3 (δ 3 ) = x 4 ξ 1 (δ 1 ) velocity v 0 x 0 ξ 4 (δ 4 ) = x 5 x 2 5 ξ 2 (δ 2 ) = x position x 10
11 Overview Numerical Analysis Running Example Solving ODEs Computing Trajectories and Jumps Reachability with Piecewise Affine Dynamics Reachability with Nonlinear Dynamics 11
12 Solving ODEs Given an ordinary differential equation (ODE) ẋ = f(x), with initial value x 0, find ξ(t) with ξ(0) = x 0 and ξ(t) = f(ξ(t)) for all t 0. Numerical solution by computing x 0,..., x N such that x i ξ(t i ) at time points t 0,..., t N. 2 Using fixed time step h: t i = ih. 2 R. P. Canale and S. C. Chapra, Numerical methods for engineers, Mc Graw Hill, New York,
13 Euler s Method Compute x 0,..., x N with the sequence x i+1 = x i + f(x i )h. Comparing to Taylor series around x i, (n 1) x i+1 = x i + ẋ i h + ẍi x i 2! h n! hn +, obtain estimate of local error ε a = O(h 2 ). global error ε g = O(h) first-order method accuracy limited by numerical roundoff error O(1/h) 13
14 Ball on String in Extension: Euler s Method 1 h = 0.05 position x 0 h = h = t 14
15 Ball on String in Extension: Euler s Method x 2 x 3 velocity v 5 0 x 1 x h = h = x 15
16 Stability The linear ODE ẋ = ax, converges to zero iff a < 0. Euler s method x i+1 = x i + f(x i )h = x i + ax i h = (1 + ah)x i converges to zero iff 1 + ah < 1 conditionally stable. 16
17 Backwards Euler Method Compute x 0,..., x N with the sequence x i+1 = x i + f(x i+1 )h, solved for x i+1 at each i using root-finding (Newton s method). implicit method Backwards Euler for ẋ = ax, x i+1 = x i + ax i+1 h = 1 1 ah x i converges for all a < 0, h > 0 unconditionally stable. 17
18 Runge-Kutta Methods Explicit Runge-Kutta methods compute the sequence x i+1 = x i + ϕ(x i, h)h, ϕ(x i, h) = a 1 k 1 + a 2 k a n k n, weights a i,q ij and derivative k j = f(ˆx i j ) at intermediate states ˆx 1 i = x i, ˆx 2 i = x i + q 11 k 1 h, ˆx 3 i = x i + q 21 k 1 h + q 21 k 2 h,. ˆx n i = x i + q (n 1)1 k 1 h + q (n 1)2 k 2 h + + q (n 1)(n 1) k n 1 h 18
19 (Kutta, 1901) Runge-Kutta Methods Runge-Kutta method defined by n and parameters a i,q ij chosen to match first n terms of Taylor series. Remaining degrees of freedom used to optimize, e.g., truncation error O(h n+1 ) and global error O(h n ) for n = 2,..., 5. Ralston s method has the smallest truncation error for n = 2: k 1 = f(ˆx 1 i ), ˆx 1 i = x i, k 2 = f(ˆx 2 i ), ˆx 2 i = x i + 3/4k 1 h, ϕ(x i, h) = 1k k
20 Ball on String in Extension: Ralston s Method position x x 2 ˆx 2 2 x 3 ˆx 2 3 h = 0.05 h = ˆx 2 x t 20
21 Ball on String in Extension: Ralston s Method ˆx 2 1 x 1 ˆx h = 0.05 x 2 h = 0.1 velocity v 0 x 0 x 3 ˆx ˆx 2 4 x x 21
22 Error Estimation Error estimation using difference between results for time steps h/2 and h, results for (n 1)th and nth order solver. Runge-Kutta-Fehlberg (RKF) methods (Fehlberg, 1970) compare (n 1)th and nth order RK, reuse intermediate results, no more evaluations than nth order RK. Popular: RKF 2(3) and RKF 4(5), aka and. 22
23 Adaptive Time Steps Adapt time step using error estimation ε a and desired error ϵ d. Heuristics 3 : h h ϵd /ε a 0.25 if εa < ϵ d h h ϵ d /ε a 0.2 if εa ϵ d 3 W. H. Press, Numerical recipes 3rd edition: The art of scientific computing. Cambridge University Press,
24 Stiff Systems Stiff ODEs have time constants (Eigenvalues) differing by a factor of 1000 or more. Solvers take tiny time steps throughout the entire time horizon Special solvers are available for stiff ODEs, using implicit methods to achieve stability at larger time steps. Ball on String example is stiff for small mass. 24
25 Ball on String: Stiff for Small Mass 0 h = position x velocity v x 0 h = { , t , t > 0.4 h = h = { , t , t > 0.4 h = h = 5 10 v t t 25
26 Overview Numerical Analysis Running Example Solving ODEs Computing Trajectories and Jumps Reachability with Piecewise Affine Dynamics Reachability with Nonlinear Dynamics 26
27 Computing Trajectories and Jumps Numerical simulation of hybrid automata: use ODE solver to approximate trajectories, detect when trajectory enters guard, detect when trajectory leaves invariant. ODE solver offer zero crossing detection using root-finding algorithms. Detecting guards/invariants using root functions is computationally expensive and potentially inaccurate. 27
28 Shortcomings 4 Missed roots violations of invariant or entering guard go undetected. Increased cost ODE solvers reuse intermediate states for increasing time sequence. Lost through back-and-forth of root-finding. Spurious behavior Numerically approximated state may lie slightly outside the guard or invariant, so constraints are relaxed 4 F. Zhang, M. Yeddanapudi, and P. Mosterman, Zero-crossing location and detection algorithms for hybrid system simulation, in IFAC World Congress, 2008, pp
29 Zeno Behavior Zeno behavior occurs if infinitely many events occur in a bounded time interval. 5 Chattering Zeno: zero-time events Genuine Zeno: event times converge towards a fixed point Simulator seems to get stuck as switching times converge. Ball on String example zeno if upside-down (negative gravity) bouncing ball. 5 A. D. Ames and S. Sastry, Characterization of zeno behavior in hybrid systems using homological methods, in ACC 05,
30 Ball on String: Zeno Behavior 1 x 0 position x t velocity v v t 30
31 Accounting for Nondeterminism The biggest challenge is nondeterminism: select initial state and successor states in jump relation; choose between transitions if guards overlap; choose jump time from interval of time; differential inclusions, such as ẋ [ 1, 1], require to pick derivative for each time step; Number of runs increases exponentially with each choice. Simulators like Simulink, Modelica, or Ptolemy use purely deterministic models that jump as soon as possible. 31
32 Overview Numerical Analysis Reachability with Piecewise Affine Dynamics Set Representations SpaceEx (advertisement) EMB Case Study Reachability with Nonlinear Dynamics 32
33 Piecewise Affine Dynamics Hybrid automata with piecewise affine dynamics (PWA) initial states and invariants are polyhedra, flows are affine ODEs ẋ = Ax + Bu, u U, jumps have a guard set and assignments x = Cx + w, w W. 33
34 Exercise: When will Lake Mead go dry? 6 by Serge Melki from Indianapolis, USA [CC BY 2.0], via Wikimedia Commons Lake Mead is fed by the Colorado River r in 15 million acre feet per year (maf/yr). The feed rate is expected to decrease by up to 30% over t dec = 50 years. Water is taken out of the lake at a rate of 14 maf/yr, plus 1.7 maf/yr due to evaporation and infiltration, resulting in a constant loss of r out 15.7 maf/yr. 6 Barnett and Pierce. When will Lake Mead run dry? J. Water Resources Research,
35 Exercise: When will Lake Mead go dry? Switched Affine Model: In the rainy season, April to June, r in,rainy = 30 maf/yr, while in the rest of the year, dry season, r in,dry = 10 maf/yr. a) Model that r in decreases linearly by 30% in t dec years. b) In what time frame could Lake Mead be dry? 35
36 Reachability Discrete Successors Edge e = (l, α, l ) with guard x G and nondeterministic assignment x = Cx + w, w W, post e (l P) = l ( C(P G) W ) Inv(l ). 36
37 Reachability Continuous successors ẋ = Ax + Bu, u U, trajectory ξ(t) from ξ(0) = x 0 for given input signal ζ(t) U: ξ x0,ζ(t) = e At x 0 + t 0 e A(t s) Bζ(s)ds. reachable states from set X 0 for any input signal: Y t = t 0 X t = e At X 0 Y t, t/δ e As BUds = lim e Aδk δbu. δ 0 k=0 37
38 Computing a Convex Cover Ω 2 X 2δ X δ Ω 1 X 0 Ω 0 Compute Ω 0, Ω 1,... such that X t Ω 0 Ω t T 38
39 Time Discretization Ω 2 X 2δ X δ Ω 1 X 0 Semi-group property: (X kδ ) δ = X (k+1)δ Time discretization: X (k+1)δ = e Aδ X kδ Y δ. Ω 0 Given initial approximations Ω 0 and Ψ δ such that X t Ω 0, Y δ Ψ δ, 0 t δ X t is covered by the sequence Ω k+1 = e Aδ Ω k Ψ δ. 39
40 Ball on String: Reachable States (clip from SpaceEx output) 40
41 Initial Approximations X δ Ω 0 X 0 (a) convex hull and pushing facets (b) convex hull and bloating 41
42 Initial Approximations Forward Bloating Bloating based on norms: 7 Ω 0 = chull(x 0 e Aδ X 0 ) (α δ + β δ )B, Ψ δ = β δ B, α δ = µ(x 0 ) (e A δ 1 A δ), β δ = 1 µ(bu) A (e A δ 1), with radius µ(x ) = max x X x and unit ball B. 7 A. Girard, Reachability of uncertain linear systems using zonotopes, in HSCC, 2005, pp
43 Initial Approximations Forward Bloating X 2δ X δ ˆΩ 2 ˆΩ 1 X 0 Forward bloating is tight on X 0 and bloated on X δ. Improvements: intersect forward bloating with backward bloating bloat based on interpolation error (shown before) ˆΩ 0 43
44 Wrapping Effect Appr(e Aδ Appr(e Aδ X 0 )) Appr(e Aδ X 0 ) Appr(e A2δ X 0 ) Appr(e Aδ X 0 ) e Aδ X 0 e Aδ X 0 X 0 X 0 (a) with wrapping effect (b) using a wrapping-free algorithm avoid increasing complexity through approximation ˆΩ k+1 = Appr(e Aδ ˆΩk Ψ δ ). wrapping effect: error accumulation 44
45 Wrapping Effect Solution: Split sequence 8 ˆΨ k+1 = Appr(e Akδ Ψ δ ) ˆΨ k, with ˆΨ 0 = {0}, ˆΩ k = Appr(e Akδ Ω 0 ) ˆΨ k. satisfies ˆΩ k = Appr(Ω k ) (wrapping-free) if e.g., bounding box. Appr(P Q) = Appr(P) Appr(Q), 8 A. Girard, C. L. Guernic, and O. Maler, Efficient computation of reachable sets of linear time-invariant systems with inputs, in HSCC, 2006, pp
46 Polyhedra Ω 2 Ω 1 X 2δ Ω 0 X δ X 0 #ops polyhedra m constr. k gen. convex hull exp 1 Minkowski sum exp nk 2 linear map n 2 m / exp n 2 k intersection 1 exp 46
47 Ellipsoids 9 X 2δ X δ ˆΩ 2 ˆΩ 1 X 0 ˆΩ 0 #ops polyhedra ellipsoids m constr. k gen. n n matrix convex hull exp 1 approx Minkowski sum exp nk 2 approx linear map n 2 m / exp n 2 k n 3 intersection 1 exp approx 9 A. B. Kurzhanski and P. Varaiya, Dynamics and control of trajectory tubes. Springer,
48 Zonotopes v 2 v 1 v 4 c v 3 Zonotope with center c R n and generators v 1,..., v k R n P = { c + k i=1 } α i v i α i [ 1, 1]. subclass of centrally symmetric polytopes 48
49 Zonotopes alternative representation: affine image of hypercube P = c MB k. generator matrix: M = [v 1... v k ] unit hypercube in k dimensions: { } B k = x R k x 1. 49
50 Zonotopes convex hull: approximation 10 CH(P e Aδ P) 1 2 { c + e Aδ c + k P i=1 α i (v i + e Aδ v i ) + 2k P i=k P +1 α i (v i e Aδ v i ) } + α 2kP +1(c e Aδ c) α i [ 1, 1]. binary operations: 2n 2 (k + 1) + 2n(k + 2) 10 A. Girard, Reachability of uncertain linear systems using zonotopes, in HSCC, 2005, pp
51 Zonotopes Minkowski sum: add centers and join generators P Q = { c P + c Q + k P i=1 α i v P i + k P +k Q i=k P +1 α i v Q i } α i [ 1, 1]. binary operations: n 51
52 Zonotopes linear map: map center and generators MP = { Mc + k i=1 binary operations: n 2 (k + 1) } α i Mv i α i [ 1, 1]. 52
53 Zonotopes in Reachability reachable states from set X 0 for any input signal: approximation: t/δ X t = e At X 0 lim e Aδk δu. δ 0 k=0 t/δ ˆX t = e At X 0 e Aδk (δu E U,δ ). k=0 with X 0, U zonotopes, ˆX t is a zonotope for any δ > 0 53
54 Zonotopes 11 X 2δ X δ ˆΩ 2 ˆΩ 1 X 0 ˆΩ 0 #ops polyhedra ellipsoids zonotopes m constr. k gen. n n matrix k generators convex hull exp 1 approx n 2 k approx Minkowski sum exp nk 2 approx n linear map n 2 m / exp n 2 k n 3 n 2 k intersection 1 exp approx approx 11 A. Girard, Reachability of uncertain linear systems using zonotopes, in HSCC, 2005, pp
55 Support Functions d 2 ρ P (d) d d 3 d 1 0 d 4 (a) support function in direction d (b) outer approximation support function = linear optimization (efficient!) ρ P (d) = max{d T x x P}. computed values define polyhedral outer approximation P D = d D{ d T x ρ P (d) }. 55
56 Support Functions d 2 ρ P (d) d d 3 d 1 0 d 4 (a) support function in direction d (b) outer approximation linear map: ρ MX (l) = ρ X (M T l), O(mn), convex hull: ρ chull(p Q) (l) = max{ρ P (l), ρ Q (l)}, O(1), Minkowski sum: ρ X Y (l) = ρ X (l) + ρ Y (l), O(1). 56
57 (Le Guernic, Girard, 09)[12] Support Functions X 2δ X 2δ X δ X δ ˆΩ 2 ˆΩ 2 ˆΩ 1 X 0 ˆΩ 1 X 0 ˆΩ 0 ˆΩ 0 support functions: lazy approximation on demand #ops polyhedra ellipsoids zonotopes support f. m constr. k gen. n n matrix k generators convex hull exp 1 approx n 2 k approx 2#ops Minkowski sum exp nk 2 approx n 2#ops linear map n 2 m / exp n 2 k n 3 n 2 k n 2 + #ops intersection 1 exp approx approx opt. / approx 57
58 Example: Switched Oscillator Switched oscillator 2 continuous variables 4 discrete states similar to many circuits (Buck converters, ) plus linear filter m continuous variables dampens output signal affine dynamics total 2 + m continuous variables 58
59 Example: Switched Oscillator Low number of directions sufficient? here: 6 state variables 12 box constraints (axis directions) 72 octagonal constraints (± x i ± x j ) 59
60 Example: Switched Oscillator Scalability Measurements: fixpoint reached in O(nm 2 ) time box constraints: O(n 3 ) octagonal constraints: O(n 5 ) runtime [s] number of variables n 60
61 Example: Controlled Helicopter Photo by Andrew P Clarke 28-dim model of a Westland Lynx helicopter 8-dim model of flight dynamics 20-dim continuous H controller for disturbance rejection stiff, highly coupled dynamics 61
62 28 state variables + clock Example: Helicopter CAV 11: 1440 sets in 5.9s 1440 time steps 62
63 Example: Helicopter 28 state variables + clock convex in 29 dimensions! HSCC 13: 32 sets in 15.2s (4.8s clustering) time steps, median
64 Example: Chaotic Circuit piecewise linear Rössler-like circuit Pisarchik, Jaimes-Reátegui. ICCSDS 05 added nondet. disturbances 3 variables, hard! 64
65 Overview Numerical Analysis Reachability with Piecewise Affine Dynamics Set Representations SpaceEx (advertisement) EMB Case Study Reachability with Nonlinear Dynamics 65
66 SpaceEx Verification Platform Browser-based GUI 2D/3D output runs remotely 66
67 SpaceEx Model Editor Components = Hybrid Automata real-values variables ODE, linear DAE 67
68 SpaceEx Model Editor Block diagrams connect components templates, nesting 68
69 SpaceEx Reachability Algorithms PHAVer constant dynamics (LHA) formally sound and exact Support Function Algo many continuous variables low discrete complexity Simulation nonlinear dynamics based on CVODE spaceex.imag.fr 69
70 Overview Numerical Analysis Reachability with Piecewise Affine Dynamics Set Representations SpaceEx (advertisement) EMB Case Study Reachability with Nonlinear Dynamics 70
71 Case Study: Electro-Mechanical Brake 71
72 Case Study: Electro-Mechanical Brake Controller Implementation discrete time fixed-point arithmetic multi-tasking processor: scheduling with uncertain frequency worst-case analysis too conservative 72
73 Case Study: Electro-Mechanical Brake 73
74 Case Study: Electro-Mechanical Brake Typical Worst-Case Execution Time limit missed schedules per time interval 74
75 Case Study: Electro-Mechanical Brake caliper position time 75
76 Case Study: Electro-Mechanical Brake caliper position only failure hard to detect artificial failure case (inconsistent with classical theory) time 76
77 Case Study: Electro-Mechanical Brake current physical properties: maximum impulse on contact (measured via current) time 77
78 Overview Numerical Analysis Reachability with Piecewise Affine Dynamics Reachability with Nonlinear Dynamics 78
79 Nonlinear Dynamics Linearization ẋ = f(x), with f globally Lipschitz continuous. Linearization: choose domain S (partition, sliding window) overapproximate in S with ẋ = Ax + u, u U linearizing f(x) around x 0 S gives a ij = f i x j and b = f(x 0 ) Ax 0. x=x0 U = Appr {f(x) (Ax + b), x S} b. 79
80 Example: Van der Pol Oscillator 12 ẋ = y ẏ = y(1 x 2 ) x hybridization: here triangular partition of size 0.05 partitioning generally doesn t scale well 12 E. Asarin, T. Dang, and A. Girard, Hybridization methods for the analysis of nonlinear systems, Acta Inf., vol. 43, no. 7, pp ,
81 Nonlinear Dynamics Polynomial Approximations Bernstein polynomials for polynomial f(x) polyhedral approximation of successors 13 Taylor models polynomial approximations of Taylor expansion represent sets with polynomials Flow* verification tool [15] 13 T. Dang and R. Testylier, Reachability analysis for polynomial dynamical systems using the bernstein expansion, Reliable Computing, vol. 17, no. 2, pp ,
82 References [2] R. Alur, C. Courcoubetis, N. Halbwachs, T. Henzinger, P.-H. Ho, X. Nicollin, A. Olivero, J. Sifakis, and S. Yovine, The algorithmic analysis of hybrid systems, Theoretical Computer Science, vol. 138, pp. 3 34, [3] T. A. Henzinger, The theory of hybrid automata., in LICS, Los Alamitos: IEEE Computer Society, 1996, pp [8] T. A. Henzinger, P. W. Kopke, A. Puri, and P. Varaiya, What s decidable about hybrid automata? Journal of Computer and System Sciences, vol. 57, pp , [12] C. Le Guernic and A. Girard, Reachability analysis of linear systems using support functions, Nonlinear Analysis: Hybrid Systems, vol. 4, no. 2, pp , [15] X. Chen, E. Ábrahám, and S. Sankaranarayanan, Taylor model flowpipe construction for non-linear hybrid systems, in RTSS, IEEE Computer Society, 2012, pp , ISBN:
An Introduction to Hybrid Automata, Numerical Simulation and Reachability Analysis
An Introduction to Hybrid Automata, Numerical Simulation and Reachability Analysis Goran Frehse SyDe Summer School, September 10, 2015 Univ. Grenoble Alpes Verimag, 2 avenue de Vignate, Centre Equation,
More informationModel Checking of Hybrid Systems
Model Checking of Hybrid Systems Goran Frehse AVACS Autumn School, October 1, 2015 Univ. Grenoble Alpes Verimag, 2 avenue de Vignate, Centre Equation, 38610 Gières, France, Overview Hybrid Automata Set-Based
More informationModeling and Analysis of Hybrid Systems
Modeling and Analysis of Hybrid Systems 7. Linear hybrid automata II Prof. Dr. Erika Ábrahám Informatik 2 - LuFG Theory of Hybrid Systems RWTH Aachen University Szeged, Hungary, 27 September - 6 October
More informationEECS 144/244: System Modeling, Analysis, and Optimization
EECS 144/244: System Modeling, Analysis, and Optimization Continuous Systems Lecture: Hybrid Systems Alexandre Donzé University of California, Berkeley April 5, 2013 Alexandre Donzé: EECS 144/244 Hybrid
More informationReachability Analysis of Hybrid Systems using Support Functions
Reachability Analysis of Hybrid Systems using Support Functions Colas Le Guernic 1 and Antoine Girard 2 1 Verimag, Université de Grenoble 2 Laboratoire Jean Kuntzmann, Université de Grenoble {Colas.Le-Guernic,Antoine.Girard}@imag.fr
More informationModeling and Analysis of Hybrid Systems
Modeling and Analysis of Hybrid Systems Linear hybrid automata II: Approximation of reachable state sets Prof. Dr. Erika Ábrahám Informatik 2 - Theory of Hybrid Systems RWTH Aachen University SS 2015 Ábrahám
More informationWork in Progress: Reachability Analysis for Time-triggered Hybrid Systems, The Platoon Benchmark
Work in Progress: Reachability Analysis for Time-triggered Hybrid Systems, The Platoon Benchmark François Bidet LIX, École polytechnique, CNRS Université Paris-Saclay 91128 Palaiseau, France francois.bidet@polytechnique.edu
More informationVerification of Hybrid Systems
Verification of Hybrid Systems Goran Frehse Universite Grenoble 1, Verimag - with work from Thao Dang, Antoine Girard and Colas Le Guernic - MOVEP 08, June 25, 2008 1 Boeing & Tupolew Collision B757-200
More informationTime Domain Verification of Oscillator Circuit Properties
Electronic Notes in Theoretical Computer Science 153 (2006) 9 22 www.elsevier.com/locate/entcs Time Domain Verification of Oscillator Circuit Properties Goran Frehse, Bruce H. Krogh, Rob A. Rutenbar 1
More informationTime-Aware Abstractions in HybridSal
Time-Aware Abstractions in HybridSal Ashish Tiwari SRI International, Menlo Park, CA Abstract. HybridSal is a tool for enabling verification of hybrid systems using infinite bounded model checking and
More informationVerifying Safety Properties of Hybrid Systems.
Verifying Safety Properties of Hybrid Systems. Sriram Sankaranarayanan University of Colorado, Boulder, CO. October 22, 2010. Talk Outline 1. Formal Verification 2. Hybrid Systems 3. Invariant Synthesis
More informationCyber-Physical Systems Modeling and Simulation of Hybrid Systems
Cyber-Physical Systems Modeling and Simulation of Hybrid Systems Matthias Althoff TU München 05. June 2015 Matthias Althoff Modeling and Simulation of Hybrid Systems 05. June 2015 1 / 28 Overview Overview
More informationApproximately Bisimilar Finite Abstractions of Stable Linear Systems
Approximately Bisimilar Finite Abstractions of Stable Linear Systems Antoine Girard Université Joseph Fourier Laboratoire de Modélisation et Calcul B.P. 53, 38041 Grenoble, France Antoine.Girard@imag.fr
More informationHybrid systems and computer science a short tutorial
Hybrid systems and computer science a short tutorial Eugene Asarin Université Paris 7 - LIAFA SFM 04 - RT, Bertinoro p. 1/4 Introductory equations Hybrid Systems = Discrete+Continuous SFM 04 - RT, Bertinoro
More informationVerifying Global Convergence for a Digital Phase-Locked Loop
Verifying Global Convergence for a Digital Phase-Locked Loop Jijie Wei & Yan Peng & Mark Greenstreet & Grace Yu University of British Columbia Vancouver, Canada October 22, 2013 Wei & Peng & Greenstreet
More informationRigorous Simulation-Based Analysis of Linear Hybrid Systems
Rigorous Simulation-Based Analysis of Linear Hybrid Systems Stanley Bak 1 and Parasara Sridhar Duggirala 2 1 Air Force Research Laboratory 2 University of Connecticut Abstract. Design analysis of Cyber-Physical
More informationUser s Manual of Flow* Version 2.0.0
User s Manual of Flow* Version 2.0.0 Xin Chen University of Colorado, Boulder 1 Introduction Flow* is a tool for safety verification of hybrid systems. Given a real-valued interval, a natural number m,
More informationAlgorithmic Verification of Stability of Hybrid Systems
Algorithmic Verification of Stability of Hybrid Systems Pavithra Prabhakar Kansas State University University of Kansas February 24, 2017 1 Cyber-Physical Systems (CPS) Systems in which software "cyber"
More informationAnalysis of a Boost Converter Circuit Using Linear Hybrid Automata
Analysis of a Boost Converter Circuit Using Linear Hybrid Automata Ulrich Kühne LSV ENS de Cachan, 94235 Cachan Cedex, France, kuehne@lsv.ens-cachan.fr 1 Introduction Boost converter circuits are an important
More informationFormally Correct Monitors for Hybrid Automata. Verimag Research Report n o TR
Formally Correct Monitors for Hybrid Automata Goran Frehse, Nikolaos Kekatos, Dejan Nickovic Verimag Research Report n o TR-2017-5 September 20, 2017 Verimag, University of Grenoble Alpes, Grenoble, France.
More informationReachability Analysis: State of the Art for Various System Classes
Reachability Analysis: State of the Art for Various System Classes Matthias Althoff Carnegie Mellon University October 19, 2011 Matthias Althoff (CMU) Reachability Analysis October 19, 2011 1 / 16 Introduction
More informationControl of Industrial Systems to Avoid Failures: Application to Electrical System
Control of Industrial Systems to Avoid Failures: Application to Electrical System Yamen EL TOUATI Department of Computer Science Faculty of Computing and Information Technology Northern Border University
More informationECE257 Numerical Methods and Scientific Computing. Ordinary Differential Equations
ECE257 Numerical Methods and Scientific Computing Ordinary Differential Equations Today s s class: Stiffness Multistep Methods Stiff Equations Stiffness occurs in a problem where two or more independent
More informationReachability Analysis of Hybrid Systems Using Support Functions
Reachability Analysis of Hybrid Systems Using Support Functions Colas Le Guernic 1 and Antoine Girard 2 1 Verimag, Université de Grenoble 2 Laboratoire Jean Kuntzmann, Université de Grenoble {Colas.Le-Guernic,Antoine.Girard}@imag.fr
More informationLarge-Scale Linear Systems from Order-Reduction (Benchmark Proposal)
EPiC Series in Computing Volume 43, 2017, Pages 60 67 ARCH16. 3rd International Workshop on Applied Verification for Continuous and Hybrid Systems Large-Scale Linear Systems from Order-Reduction (Benchmark
More informationReachability Analysis of Nonlinear and Hybrid Systems using Zonotopes May 7, / 56
Reachability Analysis of Nonlinear and Hybrid Systems using Zonotopes Matthias Althoff Carnegie Mellon Univ. May 7, 2010 Reachability Analysis of Nonlinear and Hybrid Systems using Zonotopes May 7, 2010
More informationVerifying Analog Oscillator Circuits Using Forward/Backward Abstraction Refinement
Verifying Analog Oscillator Circuits Using Forward/Backward Abstraction Refinement Goran Frehse VERIMAG 2, av. de Vignate 38610 Gières, France goran.frehse@imag.fr Bruce H. Krogh, Rob A. Rutenbar Carnegie
More informationHyLAA: A Tool for Computing Simulation-Equivalent Reachability for Linear Systems
HyLAA: A Tool for Computing Simulation-Equivalent Reachability for Linear Systems ABSTRACT Stanley Bak Air Force Research Laboratory Aerospace Systems Directorate stanleybak@gmail.com Simulations are a
More informationVerification of analog and mixed-signal circuits using hybrid systems techniques
FMCAD, November 2004, Austin Verification of analog and mixed-signal circuits using hybrid systems techniques Thao Dang, Alexandre Donze, Oded Maler VERIMAG Grenoble, France Plan 1. Introduction 2. Verification
More informationLarge-Scale Linear Systems from Order-Reduction (Benchmark Proposal)
Large-Scale Linear Systems from Order-Reduction (Benchmark Proposal) Hoang-Dung Tran, Luan Viet Nguyen, and Taylor T. Johnson v0.1, 2016-02-18 Abstract This benchmark suite is composed of nine examples
More informationHybrid Automata and ɛ-analysis on a Neural Oscillator
Hybrid Automata and ɛ-analysis on a Neural Oscillator A. Casagrande 1 T. Dreossi 2 C. Piazza 2 1 DMG, University of Trieste, Italy 2 DIMI, University of Udine, Italy Intuitively... Motivations: Reachability
More informationApproximating Continuous Systems by Timed Automata
Approximating Continuous Systems by Timed Automata Oded Maler 1 and Grégory Batt 2 1 Verimag-UJF-CNRS, 2 Av. de Vignate, 38610 Gières, France Oded.Maler@imag.fr 2 INRIA Rocquencourt, 78153 Le Chesnay,
More informationParametric Verification and Test Coverage for Hybrid Automata Using the Inverse Method
Parametric Verification and Test Coverage for Hybrid Automata Using the Inverse Method Laurent Fribourg and Ulrich Kühne LSV ENS de Cachan, 94235 Cachan, France {kuehne,fribourg}@lsv.ens-cachan.fr Abstract.
More informationAPPROXIMATE SIMULATION RELATIONS FOR HYBRID SYSTEMS 1. Antoine Girard A. Agung Julius George J. Pappas
APPROXIMATE SIMULATION RELATIONS FOR HYBRID SYSTEMS 1 Antoine Girard A. Agung Julius George J. Pappas Department of Electrical and Systems Engineering University of Pennsylvania Philadelphia, PA 1914 {agirard,agung,pappasg}@seas.upenn.edu
More informationLecture 6: Reachability Analysis of Timed and Hybrid Automata
University of Illinois at Urbana-Champaign Lecture 6: Reachability Analysis of Timed and Hybrid Automata Sayan Mitra Special Classes of Hybrid Automata Timed Automata ß Rectangular Initialized HA Rectangular
More informationSynthesizing Switching Logic using Constraint Solving
Synthesizing Switching Logic using Constraint Solving Ankur Taly 1, Sumit Gulwani 2, and Ashish Tiwari 3 1 Computer Science Dept., Stanford University ataly@stanford.edu 2 Microsoft Research, Redmond,
More informationDifferential Equations
Differential Equations Overview of differential equation! Initial value problem! Explicit numeric methods! Implicit numeric methods! Modular implementation Physics-based simulation An algorithm that
More informationc 2011 Kyoung-Dae Kim
c 2011 Kyoung-Dae Kim MIDDLEWARE AND CONTROL OF CYBER-PHYSICAL SYSTEMS: TEMPORAL GUARANTEES AND HYBRID SYSTEM ANALYSIS BY KYOUNG-DAE KIM DISSERTATION Submitted in partial fulfillment of the requirements
More informationLecture 4: Numerical solution of ordinary differential equations
Lecture 4: Numerical solution of ordinary differential equations Department of Mathematics, ETH Zürich General explicit one-step method: Consistency; Stability; Convergence. High-order methods: Taylor
More informationAutomata-theoretic analysis of hybrid systems
Automata-theoretic analysis of hybrid systems Madhavan Mukund SPIC Mathematical Institute 92, G N Chetty Road Chennai 600 017, India Email: madhavan@smi.ernet.in URL: http://www.smi.ernet.in/~madhavan
More informationNonlinear Real Arithmetic and δ-satisfiability. Paolo Zuliani
Nonlinear Real Arithmetic and δ-satisfiability Paolo Zuliani School of Computing Science Newcastle University, UK (Slides courtesy of Sicun Gao, UCSD) 1 / 27 Introduction We use hybrid systems for modelling
More informationThe Controlled Composition Analysis of Hybrid Automata
The Controlled Composition Analysis of Hybrid Automata Ying Shang Michael D. Lemmon Department of Electrical Engineering University of Notre Dame Notre Dame IN 46556 USA Abstract A controlled hybrid automaton
More informationApplied Math for Engineers
Applied Math for Engineers Ming Zhong Lecture 15 March 28, 2018 Ming Zhong (JHU) AMS Spring 2018 1 / 28 Recap Table of Contents 1 Recap 2 Numerical ODEs: Single Step Methods 3 Multistep Methods 4 Method
More informationFrom Non-Zenoness Verification to Termination
From Non-Zenoness Verification to Termination Pierre Ganty 1, Samir Genaim 2, Ratan Lal 1, and Pavithra Prabhakar 1 1 IMDEA Software Institute, Madrid, Spain 2 Universidad Complutense de Madrid, Spain
More informationHybrid Systems - Lecture n. 3 Lyapunov stability
OUTLINE Focus: stability of equilibrium point Hybrid Systems - Lecture n. 3 Lyapunov stability Maria Prandini DEI - Politecnico di Milano E-mail: prandini@elet.polimi.it continuous systems decribed by
More informationMCE693/793: Analysis and Control of Nonlinear Systems
MCE693/793: Analysis and Control of Nonlinear Systems Systems of Differential Equations Phase Plane Analysis Hanz Richter Mechanical Engineering Department Cleveland State University Systems of Nonlinear
More informationReachability Analysis for One Dimensional Linear Parabolic Equations
Reachability Analysis for One Dimensional Linear Parabolic Equations Hoang-Dung Tran Weiming Xiang Stanley Bak Taylor T. Johnson Vanderbilt University, TN 37023, USA. Air Force Research Laboratory, USA.
More informationThe Finite Element Method for the Analysis of Non-Linear and Dynamic Systems: Non-Linear Dynamics Part I
The Finite Element Method for the Analysis of Non-Linear and Dynamic Systems: Non-Linear Dynamics Part I Prof. Dr. Eleni Chatzi Dr. Giuseppe Abbiati, Dr. Konstantinos Agathos Lecture 5/Part A - 23 November,
More informationHybridSAL Relational Abstracter
HybridSAL Relational Abstracter Ashish Tiwari SRI International, Menlo Park, CA. ashish.tiwari@sri.com Abstract. In this paper, we present the HybridSAL relational abstracter a tool for verifying continuous
More informationPiecewise-Affine Approximations for a Powertrain Control Verification Benchmark
Piecewise-Affine Approximations for a Powertrain Control Verification Benchmark Jyotirmoy V. Deshmukh, Hisahiro Ito, Xiaoqing Jin, James Kapinski, Ken Butts, Jürgen Gerhard, Behzad Samadi, Kevin Walker,
More informationModel Predictive Real-Time Monitoring of Linear Systems
Model Predictive Real-Time Monitoring of Linear Systems Xin Chen and Sriram Sankaranarayanan University of Colorado, Boulder, CO Email: {xinchen,srirams}@colorado.edu Easy to Reuse * Consistent * Well
More informationStatic-Dynamic Analysis of Security Metrics
Static-Dynamic Analysis of Security Metrics for Cyber-Physical Systems Sayan Mitra (PI), Geir Dullerud (co-pi), Swarat Chaudhuri (co-pi) University of Illinois at Urbana Champaign NSA SoS Quarterly meeting,
More informationVerification of temporal properties on hybrid automata by simulation relations
Proceedings of the 46th IEEE Conference on Decision and Control New Orleans, LA, USA, Dec. 2-4, 2007 Verification of temporal properties on hybrid automata by simulation relations A. D Innocenzo, A.A.
More informationComputing Reachable States for Nonlinear Biological Models
Computing Reachable States for Nonlinear Biological Models Thao Dang 1, Colas Le Guernic 1, Oded Maler 1 cnrs-verimag, 2, av. de Vignate, 38610 Gieres, France Email addresses: Thao.Dang@imag.fr (Thao Dang),
More informationSynthesizing Switching Logic using Constraint Solving
Synthesizing Switching Logic using Constraint Solving Ankur Taly 1, Sumit Gulwani 2, and Ashish Tiwari 3 1 Computer Science Dept., Stanford University ataly@stanford.edu 2 Microsoft Research, Redmond,
More informationTutorial: Software Tools for Hybrid Systems Verification, Transformation, and Synthesis: C2E2, HyST, and TuLiP
2016 IEEE Conference on Control Applications (CCA) Part of 2016 IEEE Multi-Conference on Systems and Control September 19-22, 2016. Buenos Aires, Argentina Tutorial: Software Tools for Hybrid Systems Verification,
More informationTesting System Conformance for Cyber-Physical Systems
Testing System Conformance for Cyber-Physical Systems Testing systems by walking the dog Rupak Majumdar Max Planck Institute for Software Systems Joint work with Vinayak Prabhu (MPI-SWS) and Jyo Deshmukh
More informationONR MURI AIRFOILS: Animal Inspired Robust Flight with Outer and Inner Loop Strategies. Calin Belta
ONR MURI AIRFOILS: Animal Inspired Robust Flight with Outer and Inner Loop Strategies Provable safety for animal inspired agile flight Calin Belta Hybrid and Networked Systems (HyNeSs) Lab Department of
More informationIdentification Methods for Structural Systems. Prof. Dr. Eleni Chatzi Lecture March, 2016
Prof. Dr. Eleni Chatzi Lecture 4-09. March, 2016 Fundamentals Overview Multiple DOF Systems State-space Formulation Eigenvalue Analysis The Mode Superposition Method The effect of Damping on Structural
More informationCompositionally Analyzing a Proportional-Integral Controller Family
2 5th IEEE Conference on Decision and Control and European Control Conference (CDC-ECC) Orlando, FL, USA, December 2-5, 2 Compositionally Analyzing a Proportional-Integral Controller Family Ashish Tiwari
More informationGeometric Programming Relaxations for Linear System Reachability
Geometric Programg Relaxations for Linear System Reachability Hakan Yazarel and George J. Pappas Abstract One of the main obstacles in the safety analysis of continuous and hybrid systems has been the
More informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 9 Initial Value Problems for Ordinary Differential Equations Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign
More informationEmbedded Systems 5. Synchronous Composition. Lee/Seshia Section 6.2
Embedded Systems 5-1 - Synchronous Composition Lee/Seshia Section 6.2 Important semantic model for concurrent composition Here: composition of actors Foundation of Statecharts, Simulink, synchronous programming
More informationas support functions [18] and polynomials [34].
Decomposed Reachability Analysis for Nonlinear Systems Xin Chen University of Colorado, Boulder, CO xinchen@colorado.edu Sriram Sankaranarayanan University of Colorado, Boulder, CO srirams@colorado.edu
More informationCHAPTER 10: Numerical Methods for DAEs
CHAPTER 10: Numerical Methods for DAEs Numerical approaches for the solution of DAEs divide roughly into two classes: 1. direct discretization 2. reformulation (index reduction) plus discretization Direct
More informationOperational models of piecewise-smooth systems
Operational models of piecewise-smooth systems Andrew Sogokon Carnegie Mellon University Pittsburgh, PA, USA asogokon@cs.cmu.edu Khalil Ghorbal INRIA Rennes, France khalil.ghorbal@inria.fr Taylor T. Johnson
More informationAbstraction-Based Guided Search for Hybrid Systems
Abstraction-Based Guided Search for Hybrid Systems Sergiy Bogomolov 1, Alexandre Donzé 2, Goran Frehse 3, Radu Grosu 4, Taylor T. Johnson 5, Hamed Ladan 1, Andreas Podelski 1, and Martin Wehrle 6 1 University
More informationScalable Static Hybridization Methods for Analysis of Nonlinear Systems
Scalable Static Hybridization Methods for Analysis of Nonlinear Systems Stanley Bak Air Force Research Laboratory Information Directorate, USA Taylor T. Johnson University of Texas at Arlington, USA Sergiy
More informationAutomatic reachability analysis for nonlinear hybrid models with C2E2
Automatic reachability analysis for nonlinear hybrid models with C2E2 Chuchu Fan 1, Bolun Qi 1, Sayan Mitra 1, Mahesh Viswanathan 1, and Parasara Sridhar Duggirala 2 1 University of Illinois, Urbana-Champaign
More informationNumerical Methods - Initial Value Problems for ODEs
Numerical Methods - Initial Value Problems for ODEs Y. K. Goh Universiti Tunku Abdul Rahman 2013 Y. K. Goh (UTAR) Numerical Methods - Initial Value Problems for ODEs 2013 1 / 43 Outline 1 Initial Value
More informationAbstractions of hybrid systems: formal languages to describe dynamical behaviour
Abstractions of hybrid systems: formal languages to describe dynamical behaviour Rebekah Carter, Eva M. Navarro-López School of Computer Science, The University of Manchester Oxford Road, Manchester, M13
More informationOrdinary differential equations - Initial value problems
Education has produced a vast population able to read but unable to distinguish what is worth reading. G.M. TREVELYAN Chapter 6 Ordinary differential equations - Initial value problems In this chapter
More informationSafety Verification of Hybrid Systems Using Barrier Certificates
Safety Verification of Hybrid Systems Using Barrier Certificates Stephen Prajna 1 and Ali Jadbabaie 2 1 Control and Dynamical Systems, California Institute of Technology, Pasadena, CA 91125 - USA, prajna@cds.caltech.edu
More informationAn Introduction to Hybrid Systems Modeling
CS620, IIT BOMBAY An Introduction to Hybrid Systems Modeling Ashutosh Trivedi Department of Computer Science and Engineering, IIT Bombay CS620: New Trends in IT: Modeling and Verification of Cyber-Physical
More informationLecture IV: Time Discretization
Lecture IV: Time Discretization Motivation Kinematics: continuous motion in continuous time Computer simulation: Discrete time steps t Discrete Space (mesh particles) Updating Position Force induces acceleration.
More informationDiscrete-state Abstractions of Nonlinear Systems Using Multi-resolution Quantizer
Discrete-state Abstractions of Nonlinear Systems Using Multi-resolution Quantizer Yuichi Tazaki and Jun-ichi Imura Tokyo Institute of Technology, Ōokayama 2-12-1, Meguro, Tokyo, Japan {tazaki,imura}@cyb.mei.titech.ac.jp
More informationCS 450 Numerical Analysis. Chapter 9: Initial Value Problems for Ordinary Differential Equations
Lecture slides based on the textbook Scientific Computing: An Introductory Survey by Michael T. Heath, copyright c 2018 by the Society for Industrial and Applied Mathematics. http://www.siam.org/books/cl80
More informationLimited-Information Control of Hybrid Systems via Reachable Set Propagation
Limited-Information Control of Hybrid Systems via Reachable Set Propagation Daniel Liberzon Coordinated Science Laboratory University of Illinois at Urbana-Champaign Urbana, IL 61821, U.S.A. liberzon@uiuc.edu
More informationLecture V: The game-engine loop & Time Integration
Lecture V: The game-engine loop & Time Integration The Basic Game-Engine Loop Previous state: " #, %(#) ( #, )(#) Forces -(#) Integrate velocities and positions Resolve Interpenetrations Per-body change
More informationEuler s Method applied to the control of switched systems
Euler s Method applied to the control of switched systems FORMATS 2017 - Berlin Laurent Fribourg 1 September 6, 2017 1 LSV - CNRS & ENS Cachan L. Fribourg Euler s method and switched systems September
More informationStep Simulation Based Verification of Nonlinear Deterministic Hybrid System
Step Simulation Based Verification of Nonlinear Deterministic Hybrid System Ratnesh Kumar, Professor, IEEE Fellow PhD Student: Hao Ren Electrical and Computer Engineering Iowa State University Verification
More informationScalable Underapproximative Verification of Stochastic LTI Systems using Convexity and Compactness. HSCC 2018 Porto, Portugal
Scalable Underapproximative Verification of Stochastic LTI Systems using ity and Compactness Abraham Vinod and Meeko Oishi Electrical and Computer Engineering, University of New Mexico April 11, 2018 HSCC
More informationCourse on Hybrid Systems
Course on Hybrid Systems Maria Prandini Politecnico di Milano, Italy Organizer and lecturer: Maria Prandini Politecnico di Milano, Italy maria.prandini@polimi.it Additional lecturers: CONTACT INFO Goran
More informationSolving PDEs with PGI CUDA Fortran Part 4: Initial value problems for ordinary differential equations
Solving PDEs with PGI CUDA Fortran Part 4: Initial value problems for ordinary differential equations Outline ODEs and initial conditions. Explicit and implicit Euler methods. Runge-Kutta methods. Multistep
More informationFinite Difference and Finite Element Methods
Finite Difference and Finite Element Methods Georgy Gimel farb COMPSCI 369 Computational Science 1 / 39 1 Finite Differences Difference Equations 3 Finite Difference Methods: Euler FDMs 4 Finite Element
More informationPhysically Based Modeling Differential Equation Basics
Physically Based Modeling Differential Equation Basics Andrew Witkin and David Baraff Pixar Animation Studios Please note: This document is 2001 by Andrew Witkin and David Baraff. This chapter may be freely
More informationarxiv: v1 [cs.fl] 25 Nov 2018
Real-Time Systems Modeling and Analysis Lakhan Shiva Kamireddy [0000 0001 6007 5408] University of Colorado, Boulder CO 80302, USA {lakhan.kamireddy}@colorado.edu arxiv:1811.10083v1 [cs.fl] 25 Nov 2018
More informationLecture Notes to Accompany. Scientific Computing An Introductory Survey. by Michael T. Heath. Chapter 9
Lecture Notes to Accompany Scientific Computing An Introductory Survey Second Edition by Michael T. Heath Chapter 9 Initial Value Problems for Ordinary Differential Equations Copyright c 2001. Reproduction
More informationPiecewise-Affine Approximations for a Powertrain Control Verification Benchmark
EPiC Series in Computer Science Volume 34, 2015, Pages 98 112 ARCH14-15. 1st and 2nd International Workshop on Applied verification for Continuous and Hybrid Systems Piecewise-Affine Approximations for
More informationOn Reachability for Hybrid Automata over Bounded Time
On Reachability for Hybrid Automata over Bounded Time Thomas Brihaye, Laurent Doyen 2, Gilles Geeraerts 3, Joël Ouaknine 4, Jean-François Raskin 3, and James Worrell 4 Université de Mons, Belgium 2 LSV,
More informationIntegration of Ordinary Differential Equations
Integration of Ordinary Differential Equations Com S 477/577 Nov 7, 00 1 Introduction The solution of differential equations is an important problem that arises in a host of areas. Many differential equations
More informationECE 422/522 Power System Operations & Planning/Power Systems Analysis II : 7 - Transient Stability
ECE 4/5 Power System Operations & Planning/Power Systems Analysis II : 7 - Transient Stability Spring 014 Instructor: Kai Sun 1 Transient Stability The ability of the power system to maintain synchronism
More informationBounded Model Checking with SAT/SMT. Edmund M. Clarke School of Computer Science Carnegie Mellon University 1/39
Bounded Model Checking with SAT/SMT Edmund M. Clarke School of Computer Science Carnegie Mellon University 1/39 Recap: Symbolic Model Checking with BDDs Method used by most industrial strength model checkers:
More informationGraded Project #1. Part 1. Explicit Runge Kutta methods. Goals Differential Equations FMN130 Gustaf Söderlind and Carmen Arévalo
2008-11-07 Graded Project #1 Differential Equations FMN130 Gustaf Söderlind and Carmen Arévalo This homework is due to be handed in on Wednesday 12 November 2008 before 13:00 in the post box of the numerical
More information7. F.Balarin and A.Sangiovanni-Vincentelli, A Verication Strategy for Timing-
7. F.Balarin and A.Sangiovanni-Vincentelli, A Verication Strategy for Timing- Constrained Systems, Proc. 4th Workshop Computer-Aided Verication, Lecture Notes in Computer Science 663, Springer-Verlag,
More informationPart IB Numerical Analysis
Part IB Numerical Analysis Definitions Based on lectures by G. Moore Notes taken by Dexter Chua Lent 206 These notes are not endorsed by the lecturers, and I have modified them (often significantly) after
More informationScientific Computing II
Scientific Computing II Molecular Dynamics Numerics Michael Bader SCCS Technical University of Munich Summer 018 Recall: Molecular Dynamics System of ODEs resulting force acting on a molecule: F i = j
More informationAnnouncements. Review. Announcements. Piecewise Affine Quadratic Lyapunov Theory. EECE 571M/491M, Spring 2007 Lecture 9
EECE 571M/491M, Spring 2007 Lecture 9 Piecewise Affine Quadratic Lyapunov Theory Meeko Oishi, Ph.D. Electrical and Computer Engineering University of British Columbia, BC Announcements Lecture review examples
More informationPhysically Based Modeling: Principles and Practice Differential Equation Basics
Physically Based Modeling: Principles and Practice Differential Equation Basics Andrew Witkin and David Baraff Robotics Institute Carnegie Mellon University Please note: This document is 1997 by Andrew
More informationTimed Automata with Observers under Energy Constraints
Timed Automata with Observers under Energy Constraints Patricia Bouyer-Decitre Uli Fahrenberg Kim G. Larsen Nicolas Markey LSV, CNRS & ENS Cachan, France Aalborg Universitet, Danmark /9 Introduction The
More information