Cyber-Physical Systems Modeling and Simulation of Continuous Systems
|
|
- Erick Boyd
- 5 years ago
- Views:
Transcription
1 Cyber-Physical Systems Modeling and Simulation of Continuous Systems Matthias Althoff TU München 29. May 2015 Matthias Althoff Modeling and Simulation of Cont. Systems 29. May / 38
2 Ordinary Differential Equations Continuous Systems Continuous systems describe the dynamics of systems that can be described by a set of continuous states (e.g. position, temperature, etc.). Starting from an initial state vector x(t 0 ) and for a given input signal u(t), a continuous system creates a continuous output signal y(t). We always assume that the state, the input, and the output are vectors of proper dimension. Static systems are special cases, where the output directly follows from a function of the input: y(t) = g(u(t)). u(t) continuous system y(t) Train example from the first lecture: desired velocity v d velocity v t + controller t Matthias Althoff Modeling and Simulation of Cont. Systems 29. May / 38
3 Ordinary Differential Equations Continuous System in State Space Form Definition A continuous system in state space form C is a tuple (ordered set): where x(0) R n is the initial state, and C = (X,U c,y c,f,g,x(0)), X R n set of states U c R m set of continuous inputs Y c R o set of continuous outputs f : X U c X flow function ẋ(t) = f(x(t),u(t)) g : X U c Y c output function y(t) = g(x(t),u(t)) Matthias Althoff Modeling and Simulation of Cont. Systems 29. May / 38
4 Ordinary Differential Equations Example I: Mass Spring Damper (MSD) System mass m position s stiffness c damping coefficient d Matthias Althoff Modeling and Simulation of Cont. Systems 29. May / 38
5 Ordinary Differential Equations Equations of Each Component (MSD System) Variables: s (position), v = ṡ (velocity), a = v = s (acceleration), F e (external force); Parameters: m (mass), d (damping coefficient), c (stiffness), s 0 (position of relaxed spring) Mass F a s mass m F a = ma = m s. Damper F d s F d = d v = d ṡ. Spring F s s F s = c (s s 0 ). Sum of forces equals 0: m s +d ṡ +c(s s 0 )+F e = 0 Matthias Althoff Modeling and Simulation of Cont. Systems 29. May / 38
6 Ordinary Differential Equations Transformation to State Space Form (MSD System) The differential equation m s +d ṡ +c(s s 0 )+F e = 0 is of second order. To transform it in state space form we need a set of first order differential equations. This is achieved by introducing a new state variable for each higher derivative that has to be lowered: x 1 := s x 2 := v We also introduce the input u(t) = F e (t). The final set of differential equations is ẋ 1 = x 2, ẋ 2 = 1 m ( c(x 1 s 0 ) d x 2 u), which has the form ẋ(t) = f(x(t),u(t)). The output is the position of the mass: y = x 1. Matthias Althoff Modeling and Simulation of Cont. Systems 29. May / 38
7 Ordinary Differential Equations Example II: RLC Circuit inductance L current i e resistance R capacitance C voltage v Matthias Althoff Modeling and Simulation of Cont. Systems 29. May / 38
8 Ordinary Differential Equations Equations of each Component (RLC Circuit) Variables: i (current), v (voltage), i e (external current); Parameters: R (resistance), C (capacitance), L (inductance). Resistance resistance R i R v i R = 1 R v. Conductance i C v i C = C v. Inductance i L v i L = 1 L t t 0 v(τ)dτ +i L (t 0 ). Sum of currents equals 0: C v + 1 R v + 1 t L t 0 v(τ)dτ +i L (t 0 )+i e = 0. Differentiation with respect to time yields: C v + 1 R v + 1 L v + d dt i e = 0 Matthias Althoff Modeling and Simulation of Cont. Systems 29. May / 38
9 Ordinary Differential Equations Transformation to State Space Form (RLC Circuit) The differential equation C v + 1 R v + 1 L v + d dt i e = 0 is of second order. Analogously to the MSD system, we introduce state variables: x 1 := v x 2 := v We also introduce the input u(t) = d dt i e(t). The final set of differential equations is ẋ 1 = x 2, ẋ 2 = 1 C ( 1 L x 1 1 R x 2 u), which has the form ẋ(t) = f(x(t),u(t)). The output is the voltage: y = x 1. Matthias Althoff Modeling and Simulation of Cont. Systems 29. May / 38
10 Ordinary Differential Equations Comparison of Example I and II The mechanical and the electrical system have the same structure. Mechanical: Electrical: ẋ 1 = x 2, ẋ 2 = 1 m ( c(x 1 s 0 ) d x 2 u), y = x 1, ẋ 1 = x 2, ẋ 2 = 1 C ( 1 L x 1 1 R x 2 u), y = x 1, Their behavior is identical for m = C, c = 1 L, d = 1 R, and s 0 = 0. From now on, we will only use the abstract state space form, which represents all kinds of systems: mechanical, electrical, thermal, chemical, etc. Matthias Althoff Modeling and Simulation of Cont. Systems 29. May / 38
11 Analytical Solutions of Ordinary Differential Equations Analytical Solutions For a very limited class of ODEs, analytical solutions exist (i.e. solutions that can be expressed analytically in terms of a finite number of certain well-known functions). There exist formularies for known solutions. We will only look at the technique separation of variables. Given is the ODE Rearranging yields dυ dξ = g(ξ) f(υ). 1 dυ f(υ) dξ = g(ξ), and integrating both sides with respect to ξ results in 1 dυ f(υ) dξ dξ = g(ξ)dξ, or equivalently, 1 f(υ) dυ = g(ξ)dξ. Matthias Althoff Modeling and Simulation of Cont. Systems 29. May / 38
12 Analytical Solutions of Ordinary Differential Equations Example: Solution for Torricelli s Law (I) 2R h(t) The separation of variables makes it possible to solve the differential equation for a liquid that flows out of a tank according to Toricelli s law: 2r ḣ = k h, k = r R 2g, where h is the height of the fluid in the tank, r and R are the radii of the hole and the tank, and g is the gravity constant, see figure above. With υ = h, ξ = t, we obtain dυ dξ = g(ξ) f(υ), g(ξ) = 1, f(υ) = k υ = kυ 0.5. Matthias Althoff Modeling and Simulation of Cont. Systems 29. May / 38
13 Analytical Solutions of Ordinary Differential Equations Example: Solution for Torricelli s Law (II) Integration yields 1 f(υ) dυ = g(ξ)dξ = 1 k υ 0.5 dυ = 2 k υ0.5 +C 1 1dξ = ξ +C 2 Since 1 f(υ) dυ = g(ξ)dξ, one obtains 2 k υ0.5 = ξ +C υ 0.5 = k (ξ +C) 2 ( υ = k ) 2. 2 (ξ +C) The constant C is obtained from the initial state υ(0) = h 0 : υ(0) = ( k 2 C)2 = h 0 C = ± 2 k h0. Matthias Althoff Modeling and Simulation of Cont. Systems 29. May / 38
14 Analytical Solutions of Ordinary Differential Equations Example: Solution for Torricelli s Law (III) Inserting υ = h and ξ = t yields the solution ( h(t) = h0 k ) 2. 2 t Since the differential equation is no longer defined for h(t) = 0 (tank is empty), we obtain the final result ( h0 ) 2 k h(t) = 2 t for 0 t 2 h 0 k, 0 for 2 h 0 k < t. h(t) h 0 2 h 0 k t Matthias Althoff Modeling and Simulation of Cont. Systems 29. May / 38
15 Analytical Solutions of Ordinary Differential Equations Linear Time-Invariant Ordinary Differential Equations An important class are linear time invariant (LTI) ordinary differential equations since many technical systems can be described by LTI systems, LTI systems have an analytical solution, the solution for the initial state and the input can be obtained separately (superposition principle, see later). A one-dimensional LTI system is ẋ(t) = ax(t)+u(t), where x is the state and u is the input. We first consider the homogeneous solution (u(t) = 0): ẋ(t) = ax(t), for which separation of variables can be applied using υ = x, ξ = t. Matthias Althoff Modeling and Simulation of Cont. Systems 29. May / 38
16 Analytical Solutions of Ordinary Differential Equations Homogeneous Solution of Scalar LTI Systems Integration yields 1 f(υ) dυ = 1 aυ dυ = 1 a ln υ +C 1; Since 1 f(υ) dυ = g(ξ)dξ, one obtains 1 ln υ = ξ +C a υ = e a(ξ+c) υ = C e (aξ) g(ξ)dξ = The constant C is obtained from the initial state υ(t 0 ) = x 0 : υ(t 0 ) = C e (at0) = x 0 C = e ( at0) x 0. 1dξ = ξ +C 2 Homogeneous solution x (h) (t) = e a(t t0) x 0 Matthias Althoff Modeling and Simulation of Cont. Systems 29. May / 38
17 Analytical Solutions of Ordinary Differential Equations Particular Solution of Scalar LTI Systems We assume that variation of parameters (German: Variation der Konstanten ) yields the particular solution x (p) (t) for x 0 = 0, u(t) 0: x (p) (t) = c(t)x (h) (t). Re-arranging the LTI system and inserting x (p) (t) results in from which follows that c(t) = u(t) = ẋ (p) (t) ax (p) (t) t = ċ(t)x (h) (t)+c(t)ẋ (h) (t) c(t)ax (h) (t) }{{} =0 since ẋ (h) (t)=ax (h) (t) 1 u(τ)dτ +C = t 0 x (h) (τ) Particular solution (x 0 = 0 C = 0) x (p) (t) = c(t)x (h) (t) = t t t 0 e a(τ t0) u(τ)dτ +C t 0 e a(t τ) u(τ)dτ Matthias Althoff Modeling and Simulation of Cont. Systems 29. May / 38
18 Analytical Solutions of Ordinary Differential Equations Superposition Principle of LTI Systems Superposition principle If x (1) (t) and x (2) (t) are solutions of ẋ(t) = ax(t)+u(t) for x (1) 0, x(2) 0 and u (1) (t), u (2) (t), respectively, then x(t) = αx (1) (t)+βx (2) (t) is a solution for x 0 = αx (1) 0 +βx (2) 0 and u(t) = αu (1) (t)+βu (2) (t). Proof: x(t) = αx (1) (t)+βx (2) (t) ẋ(t) = αẋ (1) (t)+βẋ (2) (t) = α ( ax (1) (t)+u (1) (t) ) +β ( ax (2) (t)+u (2) (t) ) = a ( αx (1) (t)+βx (2) (t) ) +αu (1) (t)+βu (2) (t) = ax(t)+u(t). Matthias Althoff Modeling and Simulation of Cont. Systems 29. May / 38
19 Analytical Solutions of Ordinary Differential Equations Solution of Scalar LTI Systems Due to the superposition principle, the final solution is a combination of the homogeneous solution and the particulate solution: Solution for one-dimensional LTI systems x(t) = x (h) (t)+x (p) (t) = e a(t t 0) x 0 + t t 0 e a(t τ) u(τ)dτ Example: ẋ(t) = x(t)+u(t), x 0 = 1, u(t) = 2 (const) The solution is x (h) (t) = e t x (p) (t) = 2 t t 0 e (t τ) dτ = 2(1 e (t t 0) ) The particular solution always becomes dominant for stable systems. Why? x (h) (t) x (p) (t) x(t) t 5 10 Matthias Althoff Modeling and Simulation of Cont. Systems 29. May / 38
20 Analytical Solutions of Ordinary Differential Equations Linear Time-Invariant Systems (n-dimensional) We are now considering n-dimensional instead of one-dimensional LTI systems: ẋ(t) = Ax(t)+u(t), where x R n is the state, u(t) R n is the input, and A R n n is the system matrix. The solution is analogous to one-dimensional LTI systems: Solution for n-dimensional LTI systems x(t) = x (h) (t)+x (p) (t) = e A(t t0) x 0 + t t 0 e A(t τ) u(τ)dτ The exponential matrix is defined as the power series e X 1 = k! Xk = I +X + 1 2! X ! X k=0 where I is the identity matrix. From this follows that d dt eat = d dt I + d dt At + d 1 dt 2! (At) = A+A 1 1! (At)+A 1 2! (At) = Ae At Matthias Althoff Modeling and Simulation of Cont. Systems 29. May / 38
21 Analytical Solutions of Ordinary Differential Equations Proof of Solution for Linear Time-Invariant Systems To prove correctness, we require the Leibniz integral rule: d dx v(x) u(x) Differentiation of x(t) yields F(x,φ)dφ = F(x,v) dv v(x) dx F(x,u)du dx + u(x) x F(x,φ)dφ t ẋ(t) = d dt ea(t t0) x 0 + d e A(t τ) u(τ)dτ dt t } 0 {{} =Ae A(t t0) x 0 +u(t)+ ( t =A e A(t t0) x 0 + Leibniz t t 0 Ae A(t τ) u(τ)dτ e A(t τ) u(τ)dτ t 0 =Ax(t) + u(t) ) +u(t) Matthias Althoff Modeling and Simulation of Cont. Systems 29. May / 38
22 Analytical Solutions of Ordinary Differential Equations Example: Mass Spring Damper (MSD) System Reminder: ẋ 1 = x 2, ẋ 2 = 1 m ( c(x 1 s 0 ) d x 2 û). After introducing the state matrix A and the input vector u as [ ] [ ] A = c m d, u = c m m s 0 mû 1 The system dynamics can be written as ẋ(t) = Ax(t)+u(t). We use parameters from a quarter car model: m [kg] c [N/m] d [Ns/m] s 0 [m] Matthias Althoff Modeling and Simulation of Cont. Systems 29. May / 38
23 Analytical Solutions of Ordinary Differential Equations Exact Solution of MSD System x1(t) t 1 x2(t) 0 1 x2(t) t x 1(t) Matthias Althoff Modeling and Simulation of Cont. Systems 29. May / 38
24 Numerical Solutions of Ordinary Differential Equations Numerical Solutions Analytical solutions only exist for specific nonlinear ODEs and linear ODEs. All other systems have to be solved using numerical techniques. We will start with simple techniques and move towards more sophisticated techniques used in commercial tools (such as MATLAB/Simulink): Euler method Heun method Runge-Kutta method Runge-Kutta-Fehlberg method For simplicity, we will not consider inputs from now on. Inputs, however, can be easily integrated in each method. Matthias Althoff Modeling and Simulation of Cont. Systems 29. May / 38
25 Numerical Solutions of Ordinary Differential Equations Euler Method The Euler method is the simplest algorithm for approximating solutions of ODEs. The solution of the ODE ẋ(t) = f(x(t)) for a future point in time is x(t k+1 ) = x(t k )+ tk+1 t k f(x(τ))dτ. Since the integral cannot be exactly obtained in general, one can use: Euler method x(t k+1 ) x(t k )+(t k+1 t k )f(x(t k )). }{{} =:h The Euler method is an explicit method since it only uses values from the past. x(t) exact solution x(t k ) x(t k+1 ) h t k t k+1 Matthias Althoff Modeling and Simulation of Cont. Systems 29. May / 38 t
26 Numerical Solutions of Ordinary Differential Equations Heun Method Heun s method does not only use the slope at the previous time step, but the average slope from the previous and the next time step: x(t k+1 ) x(t k )+(t k+1 t }{{ k ) 1 } 2 =:h (f ( x(t k ) ) +f ( x(t k+1 ) )) The above scheme is also referred to as the trapezoidal rule, which is an implicit method since it requires future values. Solving for x(t k+1 ) requires solving a set of nonlinear equations, which is typically done using Newton s method. Heun s method approximates x(t k+1 ) by Euler s method: Heun s method x(t k+1 ) x(t k )+h 1 2 (K 1 +K 2 ) K 1 = f ( x(t k ) ) K 2 = f ( x(t k )+hk 1 ) Matthias Althoff Modeling and Simulation of Cont. Systems 29. May / 38
27 Numerical Solutions of Ordinary Differential Equations Illustration of Heun Method x(t) K 2 x(t k ) K 1 K 2 1 (K1 +K2) 2 x(t k+1 ) t k t k+1 t Matthias Althoff Modeling and Simulation of Cont. Systems 29. May / 38
28 Numerical Solutions of Ordinary Differential Equations Runge-Kutta Method The Runge-Kutta method is a very popular method to numerically solve ODEs. To further improve the accuracy, a new supporting point and four function evaluations are used: Runge-Kutta method x(t k+1 ) x(t k )+h 1 6 (K 1 +2K 2 +2K 3 +K 4 ) K 1 = f ( x(t k ) ) K 2 = f ( x(t k )+ 1 2 hk ) 1 K 3 = f ( x(t k )+ 1 2 hk ) 2 K 4 = f ( ) x(t k )+hk 3 Matthias Althoff Modeling and Simulation of Cont. Systems 29. May / 38
29 Numerical Solutions of Ordinary Differential Equations Illustration of Runge-Kutta Method K 4 x(t) K 3 x(t k+1 ) x(t k ) K 1 K 2 t k t k + h t k+1 2 t Matthias Althoff Modeling and Simulation of Cont. Systems 29. May / 38
30 Numerical Solutions of Ordinary Differential Equations Local Truncation Error The Heun method requires 2 function evaluations and the Runge-Kutta method even 4 function evaluations. Is the additionally effort justified? Yes! The accuracy is orders of magnitude better, making it possible to extend the step size (see exercise). Order of numerical methods The method has order p if there exist constants h and C such that x e (t k +h) x(t k +h) h C hp One can show that method: Euler Heun Runge-Kutta order: Matthias Althoff Modeling and Simulation of Cont. Systems 29. May / 38
31 Numerical Solutions of Ordinary Differential Equations Comparison of Methods for h=0.01 We use the MSD system from previous slides to compute the error for h = 0.01 to the exact solution (exact solution exists for linear systems) Euler Heun Runge Kutta 0.01 error of x t Matthias Althoff Modeling and Simulation of Cont. Systems 29. May / 38
32 Numerical Solutions of Ordinary Differential Equations Comparison of Methods for h=0.01 We use the MSD system from previous slides to compute the error for h = 0.01 to the exact solution (exact solution exists for linear systems) Euler Heun Runge Kutta error of x t Matthias Althoff Modeling and Simulation of Cont. Systems 29. May / 38
33 Numerical Solutions of Ordinary Differential Equations Comparison of Methods for h=0.1 We use the MSD system from previous slides to compute the error for h = 0.1 to the exact solution (exact solution exists for linear systems). 0.5 Euler Heun Runge Kutta error of x t Matthias Althoff Modeling and Simulation of Cont. Systems 29. May / 38
34 Numerical Solutions of Ordinary Differential Equations Comparison of Methods for h=0.1 We use the MSD system from previous slides to compute the error for h = 0.1 to the exact solution (exact solution exists for linear systems). 0.5 Euler Heun Runge Kutta error of x t Matthias Althoff Modeling and Simulation of Cont. Systems 29. May / 38
35 Numerical Solutions of Ordinary Differential Equations Discretization and Rounding Error Numerical techniques suffer from two kinds of errors: Discretization error Refers to the error due to time discretization. Smaller step sizes h decrease the discretization error. Rounding error Refers to the error due to the finite number of digits to represent numbers in computers. Larger step sizes h increase the number of time steps and thus the number of computations with rounding errors. error There exists an optimal step size for the total error: discretization error total error rounding error h Matthias Althoff Modeling and Simulation of Cont. Systems 29. May / 38
36 Numerical Solutions of Ordinary Differential Equations Comparison of Step Sizes Error using Runge-Kutta for different step sizes: x h=1e 3 6 x h=1e 4 error of x error of x t t 5 x h=1e 5 error of x The best step size is h = 1e t Matthias Althoff Modeling and Simulation of Cont. Systems 29. May / 38
37 Numerical Solutions of Ordinary Differential Equations Adaptive Step Size The discretization error varies depending on the characteristics of the ODE. In some areas, f(x) is sensitive with respect to changes of x, i.e. f(x) is large, while it is small elsewhere. x A popoular method that is based on the Runge-Kutta family of methods is the Runge-Kutta-Fehlberg method. It is implemented in MATLAB/Simulink as ode45. Basic idea of the Runge-Kutta-Fehlberg method 1 Compute one-step solution with lower order Runge-Kutta method (order 4). 2 Compute one-step solution with higher order Runge-Kutta method (order 5). 3 If the error is large, the step size is shortened, otherwise the step size is enlarged. Matthias Althoff Modeling and Simulation of Cont. Systems 29. May / 38
38 Numerical Solutions of Ordinary Differential Equations Runge-Kutta-Fehlberg Method K 1 = f ( x(t k ) ), K 2 = f ( x(t k )+ 1 4 hk ) 1 K 3 = f ( x(t k ) hk hk ) 2 K 4 = f ( x(t k ) hk hk hk 3 K 5 = f ( x(t k ) hk 1 8hK hk hk 4 K 6 = f ( x(t k ) 8 27 hk 1 +2hK hk hk hk 5 The solution of the 4 th order method is ( 25 x(t k+1 ) = x(t k )+h 216 K K K 4 1 ) 5 K 5 The solution of the 5 th order method is ( 16 x(t k+1 ) = x(t k )+h 135 K K K K ) 55 K 6 The new step size for a specified error tolerance tol is ( tolh h := 2 x(t k+1 ) x(t k+1 ) )1 4 h. ) ) ) (K 2 not used). (K 2 not used). Matthias Althoff Modeling and Simulation of Cont. Systems 29. May / 38
39 Numerical Solutions of Ordinary Differential Equations Adaptive Step size for MSD System We use the MSD system from previous slides to compute the varying step size h for tol = 1e 4 and the initial step size h = h step nr Matthias Althoff Modeling and Simulation of Cont. Systems 29. May / 38
40 Numerical Solutions of Ordinary Differential Equations Adaptive Step size for MSD System We use the MSD system from previous slides to compute the varying step size h for tol = 1e 4 and the initial step size h = Runge Kutta Fehlberg exact 0.1 x t The solution is very accurate for this system. However: One should always be skeptical with numerically obtained results. Matthias Althoff Modeling and Simulation of Cont. Systems 29. May / 38
41 Numerical Solutions of Ordinary Differential Equations Conclusions The main modeling formalism for dynamical continuous systems are ordinary differential equations (ODEs). There exist only a very limited amount of analytical solutions for OSEs. If they exist, one should prefer them over numerical solutions, of course. There exist analytical solutions for linear time-invariant systems. Numerical solutions of ODEs are best obtained with higher-order methods due to their preferred ratio of accuracy and computation time. Adaptive step size typically provides better results than fixed step size. One should always be skeptical with numerically obtained results. Matthias Althoff Modeling and Simulation of Cont. Systems 29. May / 38
Lecture 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 informationOrdinary Differential Equations II
Ordinary Differential Equations II CS 205A: Mathematical Methods for Robotics, Vision, and Graphics Justin Solomon CS 205A: Mathematical Methods Ordinary Differential Equations II 1 / 29 Almost Done! No
More informationECEN 420 LINEAR CONTROL SYSTEMS. Lecture 6 Mathematical Representation of Physical Systems II 1/67
1/67 ECEN 420 LINEAR CONTROL SYSTEMS Lecture 6 Mathematical Representation of Physical Systems II State Variable Models for Dynamic Systems u 1 u 2 u ṙ. Internal Variables x 1, x 2 x n y 1 y 2. y m Figure
More informationOrdinary Differential Equations II
Ordinary Differential Equations II CS 205A: Mathematical Methods for Robotics, Vision, and Graphics Justin Solomon CS 205A: Mathematical Methods Ordinary Differential Equations II 1 / 33 Almost Done! Last
More informationDiscrete and continuous dynamic systems
Discrete and continuous dynamic systems Bounded input bounded output (BIBO) and asymptotic stability Continuous and discrete time linear time-invariant systems Katalin Hangos University of Pannonia Faculty
More informationControl Systems. Time response
Control Systems Time response L. Lanari outline zero-state solution matrix exponential total response (sum of zero-state and zero-input responses) Dirac impulse impulse response change of coordinates (state)
More informationDynamical systems: basic concepts
Dynamical systems: basic concepts Daniele Carnevale Dipartimento di Ing. Civile ed Ing. Informatica (DICII), University of Rome Tor Vergata Fondamenti di Automatica e Controlli Automatici A.A. 2014-2015
More informationNumerical Methods for the Solution of Differential Equations
Numerical Methods for the Solution of Differential Equations Markus Grasmair Vienna, winter term 2011 2012 Analytical Solutions of Ordinary Differential Equations 1. Find the general solution of the differential
More informationDifferential and Difference LTI systems
Signals and Systems Lecture: 6 Differential and Difference LTI systems Differential and difference linear time-invariant (LTI) systems constitute an extremely important class of systems in engineering.
More informationIntroduction to Modern Control MT 2016
CDT Autonomous and Intelligent Machines & Systems Introduction to Modern Control MT 2016 Alessandro Abate Lecture 2 First-order ordinary differential equations (ODE) Solution of a linear ODE Hints to nonlinear
More informationModeling and Experimentation: Compound Pendulum
Modeling and Experimentation: Compound Pendulum Prof. R.G. Longoria Department of Mechanical Engineering The University of Texas at Austin Fall 2014 Overview This lab focuses on developing a mathematical
More informationSolving Ordinary Differential equations
Solving Ordinary Differential equations Taylor methods can be used to build explicit methods with higher order of convergence than Euler s method. The main difficult of these methods is the computation
More informationLinear Systems Theory
ME 3253 Linear Systems Theory Review Class Overview and Introduction 1. How to build dynamic system model for physical system? 2. How to analyze the dynamic system? -- Time domain -- Frequency domain (Laplace
More informationOrdinary differential equations. Phys 420/580 Lecture 8
Ordinary differential equations Phys 420/580 Lecture 8 Most physical laws are expressed as differential equations These come in three flavours: initial-value problems boundary-value problems eigenvalue
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 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 informationLecture 17: Ordinary Differential Equation II. First Order (continued)
Lecture 17: Ordinary Differential Equation II. First Order (continued) 1. Key points Maple commands dsolve dsolve[interactive] dsolve(numeric) 2. Linear first order ODE: y' = q(x) - p(x) y In general,
More informationModeling and Simulation Revision IV D R. T A R E K A. T U T U N J I P H I L A D E L P H I A U N I V E R S I T Y, J O R D A N
Modeling and Simulation Revision IV D R. T A R E K A. T U T U N J I P H I L A D E L P H I A U N I V E R S I T Y, J O R D A N 2 0 1 7 Modeling Modeling is the process of representing the behavior of a real
More informationSchool of Engineering Faculty of Built Environment, Engineering, Technology & Design
Module Name and Code : ENG60803 Real Time Instrumentation Semester and Year : Semester 5/6, Year 3 Lecture Number/ Week : Lecture 3, Week 3 Learning Outcome (s) : LO5 Module Co-ordinator/Tutor : Dr. Phang
More informationLecture Notes of EE 714
Lecture Notes of EE 714 Lecture 1 Motivation Systems theory that we have studied so far deals with the notion of specified input and output spaces. But there are systems which do not have a clear demarcation
More informationChapter 1 Fundamental Concepts
Chapter 1 Fundamental Concepts 1 Signals A signal is a pattern of variation of a physical quantity, often as a function of time (but also space, distance, position, etc). These quantities are usually the
More informationControl Systems. Time response. L. Lanari
Control Systems Time response L. Lanari outline zero-state solution matrix exponential total response (sum of zero-state and zero-input responses) Dirac impulse impulse response change of coordinates (state)
More informationLogistic Map, Euler & Runge-Kutta Method and Lotka-Volterra Equations
Logistic Map, Euler & Runge-Kutta Method and Lotka-Volterra Equations S. Y. Ha and J. Park Department of Mathematical Sciences Seoul National University Sep 23, 2013 Contents 1 Logistic Map 2 Euler and
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 informationNumerical Methods for ODEs. Lectures for PSU Summer Programs Xiantao Li
Numerical Methods for ODEs Lectures for PSU Summer Programs Xiantao Li Outline Introduction Some Challenges Numerical methods for ODEs Stiff ODEs Accuracy Constrained dynamics Stability Coarse-graining
More informationChapter 1 Fundamental Concepts
Chapter 1 Fundamental Concepts Signals A signal is a pattern of variation of a physical quantity as a function of time, space, distance, position, temperature, pressure, etc. These quantities are usually
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 informationManifesto on Numerical Integration of Equations of Motion Using Matlab
Manifesto on Numerical Integration of Equations of Motion Using Matlab C. Hall April 11, 2002 This handout is intended to help you understand numerical integration and to put it into practice using Matlab
More informationDifferential Equation (DE): An equation relating an unknown function and one or more of its derivatives.
Lexicon Differential Equation (DE): An equation relating an unknown function and one or more of its derivatives. Ordinary Differential Equation (ODE): A differential equation that contains only ordinary
More informationLecture 1 From Continuous-Time to Discrete-Time
Lecture From Continuous-Time to Discrete-Time Outline. Continuous and Discrete-Time Signals and Systems................. What is a signal?................................2 What is a system?.............................
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 informationContinuous Dynamics Solving LTI state-space equations גרא וייס המחלקה למדעי המחשב אוניברסיטת בן-גוריון
Continuous Dynamics Solving LTI state-space equations גרא וייס המחלקה למדעי המחשב אוניברסיטת בן-גוריון 2 State Space Models For a causal system with m inputs u t R m and p outputs y t R p, an nth-order
More informationProblem 1 Cost of an Infinite Horizon LQR
THE UNIVERSITY OF TEXAS AT SAN ANTONIO EE 5243 INTRODUCTION TO CYBER-PHYSICAL SYSTEMS H O M E W O R K # 5 Ahmad F. Taha October 12, 215 Homework Instructions: 1. Type your solutions in the LATEX homework
More informationNumerical Analysis of Differential Equations Numerical Solution of Parabolic Equations
Numerical Analysis of Differential Equations 215 6 Numerical Solution of Parabolic Equations 6 Numerical Solution of Parabolic Equations TU Bergakademie Freiberg, SS 2012 Numerical Analysis of Differential
More informationChapter 6 - Ordinary Differential Equations
Chapter 6 - Ordinary Differential Equations 7.1 Solving Initial-Value Problems In this chapter, we will be interested in the solution of ordinary differential equations. Ordinary differential equations
More informationACM/CMS 107 Linear Analysis & Applications Fall 2016 Assignment 4: Linear ODEs and Control Theory Due: 5th December 2016
ACM/CMS 17 Linear Analysis & Applications Fall 216 Assignment 4: Linear ODEs and Control Theory Due: 5th December 216 Introduction Systems of ordinary differential equations (ODEs) can be used to describe
More informationAC&ST AUTOMATIC CONTROL AND SYSTEM THEORY SYSTEMS AND MODELS. Claudio Melchiorri
C. Melchiorri (DEI) Automatic Control & System Theory 1 AUTOMATIC CONTROL AND SYSTEM THEORY SYSTEMS AND MODELS Claudio Melchiorri Dipartimento di Ingegneria dell Energia Elettrica e dell Informazione (DEI)
More informationIntroduction to First Order Equations Sections
A B I L E N E C H R I S T I A N U N I V E R S I T Y Department of Mathematics Introduction to First Order Equations Sections 2.1-2.3 Dr. John Ehrke Department of Mathematics Fall 2012 Course Goals The
More information1 Controllability and Observability
1 Controllability and Observability 1.1 Linear Time-Invariant (LTI) Systems State-space: Dimensions: Notation Transfer function: ẋ = Ax+Bu, x() = x, y = Cx+Du. x R n, u R m, y R p. Note that H(s) is always
More informationNonlinear Control Lecture # 14 Input-Output Stability. Nonlinear Control
Nonlinear Control Lecture # 14 Input-Output Stability L Stability Input-Output Models: y = Hu u(t) is a piecewise continuous function of t and belongs to a linear space of signals The space of bounded
More informationChap. 20: Initial-Value Problems
Chap. 20: Initial-Value Problems Ordinary Differential Equations Goal: to solve differential equations of the form: dy dt f t, y The methods in this chapter are all one-step methods and have the general
More informationStreamline calculations. Lecture note 2
Streamline calculations. Lecture note 2 February 26, 2007 1 Recapitulation from previous lecture Definition of a streamline x(τ) = s(τ), dx(τ) dτ = v(x,t), x(0) = x 0 (1) Divergence free, irrotational
More informationNUMERICAL SOLUTION OF ODE IVPs. Overview
NUMERICAL SOLUTION OF ODE IVPs 1 Quick review of direction fields Overview 2 A reminder about and 3 Important test: Is the ODE initial value problem? 4 Fundamental concepts: Euler s Method 5 Fundamental
More informationENGI 3424 First Order ODEs Page 1-01
ENGI 344 First Order ODEs Page 1-01 1. Ordinary Differential Equations Equations involving only one independent variable and one or more dependent variables, together with their derivatives with respect
More informationSystems of Ordinary Differential Equations
Systems of Ordinary Differential Equations MATH 365 Ordinary Differential Equations J Robert Buchanan Department of Mathematics Fall 2018 Objectives Many physical problems involve a number of separate
More informationCIS 4930/6930: Principles of Cyber-Physical Systems
CIS 4930/6930: Principles of Cyber-Physical Systems Chapter 2: Continuous Dynamics Hao Zheng Department of Computer Science and Engineering University of South Florida H. Zheng (CSE USF) CIS 4930/6930:
More informationMath 266: Ordinary Differential Equations
Math 266: Ordinary Differential Equations Long Jin Purdue University, Spring 2018 Basic information Lectures: MWF 8:30-9:20(111)/9:30-10:20(121), UNIV 103 Instructor: Long Jin (long249@purdue.edu) Office
More informationPowerPoints organized by Dr. Michael R. Gustafson II, Duke University
Part 6 Chapter 20 Initial-Value Problems PowerPoints organized by Dr. Michael R. Gustafson II, Duke University All images copyright The McGraw-Hill Companies, Inc. Permission required for reproduction
More informationChapter 3 Second Order Linear Equations
Partial Differential Equations (Math 3303) A Ë@ Õæ Aë áöß @. X. @ 2015-2014 ú GA JË@ É Ë@ Chapter 3 Second Order Linear Equations Second-order partial differential equations for an known function u(x,
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 informationNumerical Integration of Ordinary Differential Equations for Initial Value Problems
Numerical Integration of Ordinary Differential Equations for Initial Value Problems Gerald Recktenwald Portland State University Department of Mechanical Engineering gerry@me.pdx.edu These slides are a
More informationSolution: (a) Before opening the parachute, the differential equation is given by: dv dt. = v. v(0) = 0
Math 2250 Lab 4 Name/Unid: 1. (25 points) A man bails out of an airplane at the altitute of 12,000 ft, falls freely for 20 s, then opens his parachute. Assuming linear air resistance ρv ft/s 2, taking
More informationSection 3.4. Second Order Nonhomogeneous. The corresponding homogeneous equation. is called the reduced equation of (N).
Section 3.4. Second Order Nonhomogeneous Equations y + p(x)y + q(x)y = f(x) (N) The corresponding homogeneous equation y + p(x)y + q(x)y = 0 (H) is called the reduced equation of (N). 1 General Results
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 informationDefinition of differential equations and their classification. Methods of solution of first-order differential equations
Introduction to differential equations: overview Definition of differential equations and their classification Solutions of differential equations Initial value problems Existence and uniqueness Mathematical
More informationAnalog Signals and Systems and their properties
Analog Signals and Systems and their properties Main Course Objective: Recall course objectives Understand the fundamentals of systems/signals interaction (know how systems can transform or filter signals)
More informationsc Control Systems Design Q.1, Sem.1, Ac. Yr. 2010/11
sc46 - Control Systems Design Q Sem Ac Yr / Mock Exam originally given November 5 9 Notes: Please be reminded that only an A4 paper with formulas may be used during the exam no other material is to be
More informationLecture 2 and 3: Controllability of DT-LTI systems
1 Lecture 2 and 3: Controllability of DT-LTI systems Spring 2013 - EE 194, Advanced Control (Prof Khan) January 23 (Wed) and 28 (Mon), 2013 I LTI SYSTEMS Recall that continuous-time LTI systems can be
More informationSection 3.4. Second Order Nonhomogeneous. The corresponding homogeneous equation
Section 3.4. Second Order Nonhomogeneous Equations y + p(x)y + q(x)y = f(x) (N) The corresponding homogeneous equation y + p(x)y + q(x)y = 0 (H) is called the reduced equation of (N). 1 General Results
More informationBasic Theory of Differential Equations
page 104 104 CHAPTER 1 First-Order Differential Equations 16. The following initial-value problem arises in the analysis of a cable suspended between two fixed points y = 1 a 1 + (y ) 2, y(0) = a, y (0)
More informationObjective. 1 Specification/modeling of the controlled system. 2 Specification of a performance criterion
Optimal Control Problem Formulation Optimal Control Lectures 17-18: Problem Formulation Benoît Chachuat Department of Chemical Engineering Spring 2009 Objective Determine the control
More informationMathematics for Engineers II. lectures. Differential Equations
Differential Equations Examples for differential equations Newton s second law for a point mass Consider a particle of mass m subject to net force a F. Newton s second law states that the vector acceleration
More informationECE 602 Solution to Homework Assignment 1
ECE 6 Solution to Assignment 1 December 1, 7 1 ECE 6 Solution to Homework Assignment 1 1. For a function to be linear, it must satisfy y = f() = in addition to y = ax, so only the first graph represents
More informationOrdinary Differential Equations
Chapter 13 Ordinary Differential Equations We motivated the problem of interpolation in Chapter 11 by transitioning from analzying to finding functions. That is, in problems like interpolation and regression,
More informationDefinition of Dynamic System
Capitolo 0. INTRODUCTION 1.1 Definition of Dynamic System There are various types of dynamic systems: continuous-time, discrete-time, linear, non-linear, lumped systems, distributed systems, finite states,
More informationElementary Differential Equations
Elementary Differential Equations George Voutsadakis 1 1 Mathematics and Computer Science Lake Superior State University LSSU Math 310 George Voutsadakis (LSSU) Differential Equations January 2014 1 /
More informationLinearization problem. The simplest example
Linear Systems Lecture 3 1 problem Consider a non-linear time-invariant system of the form ( ẋ(t f x(t u(t y(t g ( x(t u(t (1 such that x R n u R m y R p and Slide 1 A: f(xu f(xu g(xu and g(xu exist and
More informationMATH 1231 MATHEMATICS 1B Calculus Section 3A: - First order ODEs.
MATH 1231 MATHEMATICS 1B 2010. For use in Dr Chris Tisdell s lectures. Calculus Section 3A: - First order ODEs. Created and compiled by Chris Tisdell S1: What is an ODE? S2: Motivation S3: Types and orders
More informationLecture 2 ELE 301: Signals and Systems
Lecture 2 ELE 301: Signals and Systems Prof. Paul Cuff Princeton University Fall 2011-12 Cuff (Lecture 2) ELE 301: Signals and Systems Fall 2011-12 1 / 70 Models of Continuous Time Signals Today s topics:
More informationConsistency and Convergence
Jim Lambers MAT 77 Fall Semester 010-11 Lecture 0 Notes These notes correspond to Sections 1.3, 1.4 and 1.5 in the text. Consistency and Convergence We have learned that the numerical solution obtained
More informationModeling and Simulation Revision III D R. T A R E K A. T U T U N J I P H I L A D E L P H I A U N I V E R S I T Y, J O R D A N
Modeling and Simulation Revision III D R. T A R E K A. T U T U N J I P H I L A D E L P H I A U N I V E R S I T Y, J O R D A N 0 1 4 Block Diagrams Block diagram models consist of two fundamental objects:
More informationAnalysis and Control of Multi-Robot Systems. Elements of Port-Hamiltonian Modeling
Elective in Robotics 2014/2015 Analysis and Control of Multi-Robot Systems Elements of Port-Hamiltonian Modeling Dr. Paolo Robuffo Giordano CNRS, Irisa/Inria! Rennes, France Introduction to Port-Hamiltonian
More informationInitial value problems for ordinary differential equations
Initial value problems for ordinary differential equations Xiaojing Ye, Math & Stat, Georgia State University Spring 2019 Numerical Analysis II Xiaojing Ye, Math & Stat, Georgia State University 1 IVP
More informationRobotics. Dynamics. Marc Toussaint U Stuttgart
Robotics Dynamics 1D point mass, damping & oscillation, PID, dynamics of mechanical systems, Euler-Lagrange equation, Newton-Euler recursion, general robot dynamics, joint space control, reference trajectory
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 information[ ] is a vector of size p.
Lecture 11 Copyright by Hongyun Wang, UCSC Recap: General form of explicit Runger-Kutta methods for solving y = F( y, t) k i = hfy n + i1 j =1 c ij k j, t n + d i h, i = 1,, p + p j =1 b j k j A Runge-Kutta
More informationIndex. Index. More information. in this web service Cambridge University Press
A-type elements, 4 7, 18, 31, 168, 198, 202, 219, 220, 222, 225 A-type variables. See Across variable ac current, 172, 251 ac induction motor, 251 Acceleration rotational, 30 translational, 16 Accumulator,
More informationMusic 420 Lecture Elementary Finite Different Schemes
Music 420 Lecture Elementary Finite Different Schemes Julius O. Smith III (jos@ccrma.stanford.edu) Center for Computer Research in Music and Acoustics (CCRMA) Department of Music, Stanford University Stanford,
More informationEE292: Fundamentals of ECE
EE292: Fundamentals of ECE Fall 2012 TTh 10:00-11:15 SEB 1242 Lecture 14 121011 http://www.ee.unlv.edu/~b1morris/ee292/ 2 Outline Review Steady-State Analysis RC Circuits RL Circuits 3 DC Steady-State
More informationNumerical Differential Equations: IVP
Chapter 11 Numerical Differential Equations: IVP **** 4/16/13 EC (Incomplete) 11.1 Initial Value Problem for Ordinary Differential Equations We consider the problem of numerically solving a differential
More informationBoyce/DiPrima/Meade 11 th ed, Ch 1.1: Basic Mathematical Models; Direction Fields
Boyce/DiPrima/Meade 11 th ed, Ch 1.1: Basic Mathematical Models; Direction Fields Elementary Differential Equations and Boundary Value Problems, 11 th edition, by William E. Boyce, Richard C. DiPrima,
More informationLecture A1 : Systems and system models
Lecture A1 : Systems and system models Jan Swevers July 2006 Aim of this lecture : Understand the process of system modelling (different steps). Define the class of systems that will be considered in this
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 informationComputational Physics (6810): Session 8
Computational Physics (6810): Session 8 Dick Furnstahl Nuclear Theory Group OSU Physics Department February 24, 2014 Differential equation solving Session 7 Preview Session 8 Stuff Solving differential
More informationFDM for parabolic equations
FDM for parabolic equations Consider the heat equation where Well-posed problem Existence & Uniqueness Mass & Energy decreasing FDM for parabolic equations CNFD Crank-Nicolson + 2 nd order finite difference
More informationENGI9496 Lecture Notes State-Space Equation Generation
ENGI9496 Lecture Notes State-Space Equation Generation. State Equations and Variables - Definitions The end goal of model formulation is to simulate a system s behaviour on a computer. A set of coherent
More informationOrdinary Differential Equations
Chapter 7 Ordinary Differential Equations Differential equations are an extremely important tool in various science and engineering disciplines. Laws of nature are most often expressed as different equations.
More informationAMJAD HASOON Process Control Lec4.
Multiple Inputs Control systems often have more than one input. For example, there can be the input signal indicating the required value of the controlled variable and also an input or inputs due to disturbances
More informationLesson 14: Van der Pol Circuit and ode23s
Lesson 4: Van der Pol Circuit and ode3s 4. Applied Problem. A series LRC circuit when coupled via mutual inductance with a triode circuit can generate a sequence of pulsing currents that have very rapid
More informationFourth Order RK-Method
Fourth Order RK-Method The most commonly used method is Runge-Kutta fourth order method. The fourth order RK-method is y i+1 = y i + 1 6 (k 1 + 2k 2 + 2k 3 + k 4 ), Ordinary Differential Equations (ODE)
More informationIntroduction to multiscale modeling and simulation. Explicit methods for ODEs : forward Euler. y n+1 = y n + tf(y n ) dy dt = f(y), y(0) = y 0
Introduction to multiscale modeling and simulation Lecture 5 Numerical methods for ODEs, SDEs and PDEs The need for multiscale methods Two generic frameworks for multiscale computation Explicit methods
More informationCompare the measured behavior with desired behavior, and
Chapter Introduction The objective of this text is to introduce methods for design and analysis of classical feedback control systems. The basic principle of feedback is to: Use a sensor to measure the
More informationMATH 23 Exam 2 Review Solutions
MATH 23 Exam 2 Review Solutions Problem 1. Use the method of reduction of order to find a second solution of the given differential equation x 2 y (x 0.1875)y = 0, x > 0, y 1 (x) = x 1/4 e 2 x Solution
More informationIntroduction to the Numerical Solution of IVP for ODE
Introduction to the Numerical Solution of IVP for ODE 45 Introduction to the Numerical Solution of IVP for ODE Consider the IVP: DE x = f(t, x), IC x(a) = x a. For simplicity, we will assume here that
More informationFinal Exam Sample Problems, Math 246, Spring 2018
Final Exam Sample Problems, Math 246, Spring 2018 1) Consider the differential equation dy dt = 9 y2 )y 2. a) Find all of its stationary points and classify their stability. b) Sketch its phase-line portrait
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 informationLecture 1: Pragmatic Introduction to Stochastic Differential Equations
Lecture 1: Pragmatic Introduction to Stochastic Differential Equations Simo Särkkä Aalto University, Finland (visiting at Oxford University, UK) November 13, 2013 Simo Särkkä (Aalto) Lecture 1: Pragmatic
More informationLinear Systems. Chapter Basic Definitions
Chapter 5 Linear Systems Few physical elements display truly linear characteristics. For example the relation between force on a spring and displacement of the spring is always nonlinear to some degree.
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 informationProblem Set 3: Solution Due on Mon. 7 th Oct. in class. Fall 2013
EE 56: Digital Control Systems Problem Set 3: Solution Due on Mon 7 th Oct in class Fall 23 Problem For the causal LTI system described by the difference equation y k + 2 y k = x k, () (a) By first finding
More information