Cyber-Physical Systems Modeling and Simulation of Continuous Systems

Size: px
Start display at page:

Download "Cyber-Physical Systems Modeling and Simulation of Continuous Systems"

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 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 information

Ordinary Differential Equations II

Ordinary 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 information

ECEN 420 LINEAR CONTROL SYSTEMS. Lecture 6 Mathematical Representation of Physical Systems II 1/67

ECEN 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 information

Ordinary Differential Equations II

Ordinary 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 information

Discrete and continuous dynamic systems

Discrete 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 information

Control Systems. Time response

Control 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 information

Dynamical systems: basic concepts

Dynamical 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 information

Numerical Methods for the Solution of Differential Equations

Numerical 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 information

Differential and Difference LTI systems

Differential 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 information

Introduction to Modern Control MT 2016

Introduction 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 information

Modeling and Experimentation: Compound Pendulum

Modeling 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 information

Solving Ordinary Differential equations

Solving 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 information

Linear Systems Theory

Linear 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 information

Ordinary differential equations. Phys 420/580 Lecture 8

Ordinary 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 information

Lecture Notes to Accompany. Scientific Computing An Introductory Survey. by Michael T. Heath. Chapter 9

Lecture 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 information

Identification Methods for Structural Systems. Prof. Dr. Eleni Chatzi Lecture March, 2016

Identification 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 information

Lecture 17: Ordinary Differential Equation II. First Order (continued)

Lecture 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 information

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

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 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 information

School of Engineering Faculty of Built Environment, Engineering, Technology & Design

School 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 information

Lecture Notes of EE 714

Lecture 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 information

Chapter 1 Fundamental Concepts

Chapter 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 information

Control Systems. Time response. L. Lanari

Control 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 information

Logistic Map, Euler & Runge-Kutta Method and Lotka-Volterra Equations

Logistic 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 information

CS 450 Numerical Analysis. Chapter 9: Initial Value Problems for Ordinary Differential Equations

CS 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 information

Numerical Methods for ODEs. Lectures for PSU Summer Programs Xiantao Li

Numerical 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 information

Chapter 1 Fundamental Concepts

Chapter 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 information

Differential Equations

Differential 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 information

Manifesto on Numerical Integration of Equations of Motion Using Matlab

Manifesto 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 information

Differential Equation (DE): An equation relating an unknown function and one or more of its derivatives.

Differential 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 information

Lecture 1 From Continuous-Time to Discrete-Time

Lecture 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 information

Applied Math for Engineers

Applied 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 information

Continuous Dynamics Solving LTI state-space equations גרא וייס המחלקה למדעי המחשב אוניברסיטת בן-גוריון

Continuous 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 information

Problem 1 Cost of an Infinite Horizon LQR

Problem 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 information

Numerical Analysis of Differential Equations Numerical Solution of Parabolic Equations

Numerical 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 information

Chapter 6 - Ordinary Differential Equations

Chapter 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 information

ACM/CMS 107 Linear Analysis & Applications Fall 2016 Assignment 4: Linear ODEs and Control Theory Due: 5th December 2016

ACM/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 information

AC&ST AUTOMATIC CONTROL AND SYSTEM THEORY SYSTEMS AND MODELS. Claudio Melchiorri

AC&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 information

Introduction to First Order Equations Sections

Introduction 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 information

1 Controllability and Observability

1 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 information

Nonlinear Control Lecture # 14 Input-Output Stability. Nonlinear Control

Nonlinear 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 information

Chap. 20: Initial-Value Problems

Chap. 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 information

Streamline calculations. Lecture note 2

Streamline 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 information

NUMERICAL SOLUTION OF ODE IVPs. Overview

NUMERICAL 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 information

ENGI 3424 First Order ODEs Page 1-01

ENGI 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 information

Systems of Ordinary Differential Equations

Systems 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 information

CIS 4930/6930: Principles of Cyber-Physical Systems

CIS 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 information

Math 266: Ordinary Differential Equations

Math 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 information

PowerPoints organized by Dr. Michael R. Gustafson II, Duke University

PowerPoints 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 information

Chapter 3 Second Order Linear Equations

Chapter 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 information

Finite Difference and Finite Element Methods

Finite 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 information

Numerical Integration of Ordinary Differential Equations for Initial Value Problems

Numerical 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 information

Solution: (a) Before opening the parachute, the differential equation is given by: dv dt. = v. v(0) = 0

Solution: (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 information

Section 3.4. Second Order Nonhomogeneous. The corresponding homogeneous equation. is called the reduced equation of (N).

Section 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 information

Scientific Computing: An Introductory Survey

Scientific 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 information

Definition of differential equations and their classification. Methods of solution of first-order differential equations

Definition 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 information

Analog Signals and Systems and their properties

Analog 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 information

sc Control Systems Design Q.1, Sem.1, Ac. Yr. 2010/11

sc 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 information

Lecture 2 and 3: Controllability of DT-LTI systems

Lecture 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 information

Section 3.4. Second Order Nonhomogeneous. The corresponding homogeneous equation

Section 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 information

Basic Theory of Differential Equations

Basic 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 information

Objective. 1 Specification/modeling of the controlled system. 2 Specification of a performance criterion

Objective. 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 information

Mathematics for Engineers II. lectures. Differential Equations

Mathematics 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 information

ECE 602 Solution to Homework Assignment 1

ECE 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 information

Ordinary Differential Equations

Ordinary 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 information

Definition of Dynamic System

Definition 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 information

Elementary Differential Equations

Elementary 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 information

Linearization problem. The simplest example

Linearization 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 information

MATH 1231 MATHEMATICS 1B Calculus Section 3A: - First order ODEs.

MATH 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 information

Lecture 2 ELE 301: Signals and Systems

Lecture 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 information

Consistency and Convergence

Consistency 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 information

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

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 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 information

Analysis and Control of Multi-Robot Systems. Elements of Port-Hamiltonian Modeling

Analysis 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 information

Initial value problems for ordinary differential equations

Initial 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 information

Robotics. Dynamics. Marc Toussaint U Stuttgart

Robotics. 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 information

The 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 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.

[ ] 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 information

Index. Index. More information. in this web service Cambridge University Press

Index. 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 information

Music 420 Lecture Elementary Finite Different Schemes

Music 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 information

EE292: Fundamentals of ECE

EE292: 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 information

Numerical Differential Equations: IVP

Numerical 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 information

Boyce/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 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 information

Lecture A1 : Systems and system models

Lecture 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 information

ECE 422/522 Power System Operations & Planning/Power Systems Analysis II : 7 - Transient Stability

ECE 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 information

Computational Physics (6810): Session 8

Computational 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 information

FDM for parabolic equations

FDM 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 information

ENGI9496 Lecture Notes State-Space Equation Generation

ENGI9496 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 information

Ordinary Differential Equations

Ordinary 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 information

AMJAD HASOON Process Control Lec4.

AMJAD 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 information

Lesson 14: Van der Pol Circuit and ode23s

Lesson 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 information

Fourth Order RK-Method

Fourth 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 information

Introduction 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. 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 information

Compare the measured behavior with desired behavior, and

Compare 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 information

MATH 23 Exam 2 Review Solutions

MATH 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 information

Introduction to the Numerical Solution of IVP for ODE

Introduction 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 information

Final Exam Sample Problems, Math 246, Spring 2018

Final 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 information

Graded Project #1. Part 1. Explicit Runge Kutta methods. Goals Differential Equations FMN130 Gustaf Söderlind and Carmen Arévalo

Graded 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 information

Lecture 1: Pragmatic Introduction to Stochastic Differential Equations

Lecture 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 information

Linear Systems. Chapter Basic Definitions

Linear 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 information

Lecture V: The game-engine loop & Time Integration

Lecture 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 information

Problem Set 3: Solution Due on Mon. 7 th Oct. in class. Fall 2013

Problem 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