REGLERTEKNIK AUTOMATIC CONTROL LINKÖPING
|
|
- Kevin York
- 6 years ago
- Views:
Transcription
1 Generating state space equations from a bond graph with dependent storage elements using singular perturbation theory. Krister Edstrom Department of Electrical Engineering Linkoping University, S Linkoping, Sweden WWW: edstrom@isy.liu.se August 25, 998 REGLERTEKNIK AUTOMATIC CONTROL LINKÖPING Report no.: LiTH-ISY-R-25 Submitted to ICBGM'99 Technical reports from the Automatic Control group in Linkoping are available by anonymous ftp at the address ftp.control.isy.liu.se. This report is contained in the compressed postscript le 25.ps.Z.
2 Generating state space equations from a bond graph with dependent storage elements using singular perturbation theory. Krister Edstrom (edstrom@isy.liu.se) Dept. of Electrical Engineering Linkopings Universitet SE Linkoping, Sweden Introduction It is well known that when allowing derivative causality in a bond graph, the generated mathematical description will be a DAE [5] (Dierential Algebraic Equation) as a mathematical description. For simulation purposes this is in most cases a suitable description. There is however a well developed theory for ODE:s (Ordinary Dierential Equations) in state space form. The ODE:s can be used to, e.g., analyze, or construct a controller for, the system. When generating a state space description from a bond graph without any derivative causality [5], algebraic loops or causal loops, the states have clear physical interpretations. For an I-element in the mechanical domain, the state is the momentum of the mass, and for a C-element in the electrical domain, the state is the charge of the capacitor. When generating a state space description from a bond graph with derivative causality, there should still be a physical interpretation. When the derivative causality is generated by two masses that are stiy connected, the state variable in a generated state space description should be the total momentum of the two masses. In this paper we will consider the derivative causality as a consequence of simplifying the model by neglecting energy dissipation. If the energy dissipation is small enough, it is a natural model simplication to remove it from the model. To allow generation of a state space description from a bond graph with connected storage elements, we will rst add the energy dissipation in the bond graph using a linear R-element. The dissipation will then be removed from the mathematical description by letting the value of the parameter in the constitutive relation of the R-element tend to zero. There exist methods for generating a state space description from a bond graph [4]. The method presented here have two properties that make it interesting: First there is a physical interpretation of the generated states, and second it relies on an interpretation of the dependency of the storage elements as the connection between the storage elements containing energy dissipation that tends to zero.
3 2 Singular perturbations In a singular perturbation framework [3] it is possible to study a system where the dynamics have two dierent time scales. A linear systems that is singularly perturbed can be expressed in the following form: _(t) _x(t) A () A 2 () A 2 () A 22 () (t) x(t) () When the value of the parameter becomes very small, the states in will have a signicantly smaller time constant than the states in x. By letting go to zero, the time constant of will go to zero. When looking at the step response of the system, will reach its steady state before x changes value, assuming that is very small. The singular perturbed model is in standard form if det A () 6. If the system is in standard form it can for be rewritten as: A () A 2 () _x(t) A 2 () A 22 () x(t) ) (2) (t) (t) _x(t) x(t)?a? ()A 2() A 22 ()? A 2 ()A? ()A 2() By assigning, the state space is reduced, and we get a DAE description of the system. As seen in Equation (2) it is possible to rewrite the equations in state space form with as an output variable if det A () 6. 3 Generating the state space description The generation of a state space description is done in ve steps as shown below. The described procedure handles one C-element with derivative causality, denoted by s, but the method can be extended to hold for I-elements, and for multiple elements with derivative causality. The bond graph is assumed to be nice, meaning that it should include no other storage elements with derivative causality, no algebraic loops and no causal loops.. Start with the bond graph with dependent storage elements. Add an R-element to relax the derivative causality. If s is connected to a -junction it is necessary to add a -junction as well. When the derivative causality is relaxed, the number of states increases with one. Let the value of the parameter in the linear constitutive relation of the R-element be. 2. Derive a state space description using SCAP [2]. _x Ax + Bu This is possible due to our assumption about a nice bond graph, since this assumptions implies that there are no zero-order causal paths (paths between storage elements or between R-elements or causal loops) or elements with derivative causality in the bond graph [5]. 3. Generate a transformation matrix T from the bond graph with derivative causality: T x T is found using a standard bond graph equation generation algorithm along the causal paths between the storage element with derivative causality and other storage elements. The generation of the transformation matrix is described more in detail in Section 4.. 2
4 4. Perform the change of basis of the state space description: where _ A + Bu A T AT? ; B T B Multiply the rst row by. The resulting state space description will have the following form where A; is invertible, see Section 4: A; + A ; A; A ; A; + B B u 5. Let!. _ A; A ; A; + B u Since A; is invertible we can according to Equation (2) replace with in the last n? equations: _ A; + B u (3) The resulting state space equations describe the dynamics of the bond graph with derivative causality. 4 Details and proofs In this section some details in the algorithm are explored. A theorem and a lemma stating some the results are formulated, but the proofs are left out and can be found in []. 4. Generation of a transformation matrix Here we will take a closer look at step 3 in the algorithm in Section 3. To generate the T matrix, we will use the bond graph in which the dependent storage elements are, i.e., the bond graph without the additional R-element. We will assume that the number of storage elements in the bond graph is n. An algorithm for generating T will look like this:. Propagate causality in the bond graph. 2. Find all storage elements that are connected to the storage element with derivative causality, s, by causal paths. Assume that there are l such storage elements, s included. Consider only the part of the bond graph connecting these elements. 3. Derive an equation x a 2 x 2 + : : : + a l x l, where x is the state variable in the constitutive relation for s and x 2 ; : : : ; x l are state variables belonging to constitutive relations for the other storage elements. Use a standard bond graph equation generation algorithm. Rewrite the equation as: x? a 2 x 2? : : :? a l x l (4) 3
5 4. Derive l? equations relating the rates of the state variables: _x ; : : : ; _x l : _x j j _x ; j 2; : : : ; l (5) This is also done with a standard bond graph equation generation algorithm. Rewrite the equations as _x j? j _x ; j 2; : : : ; l (6) and integrate: x j? j x ; j 2; : : : ; l (7) 5. Let the coecients in left hand sides of Equations (4) and (8) form a matrix T : T a 2 a 3 : : : a l? : : :? 2 : : : ? l : : : C A (8) 6. Extend T with an identity matrix to get a matrix T of dimension n, the number of storage elements in the bond graph: T T (9) I To show that T is of full rank we have the following lemma: Lemma If the bond graph has no algebraic loops, no causal loops and one storage element with derivative causality that is dependent of one or several other storage elements, then the following statement holds: T is invertible if we assume that at each C-element the energy half-arrow points towards the C-element. Proof: See []. 4.2 A state space form after the change of basis After the change of basis with the transformation matrix it will be clear that the system is a singular perturbed system in standard form. Theorem Let the bond graph have no algebraic loops, no causal loops and one storage element with derivative causality that is dependent of one or several other storage elements. Derive the state space description from the bond graph with the added R-element. _x Ax + Bu () Derive the matrix T from the bond graph without the R-element. Perform the change of basis to get a new state space description: _ A + Bu; where T x () 4
6 Hence A T AT? ; and B T B (2) If multiplying the rst row in Equation (), the structure of the matrices A; + A A ; A; A ; A; and B A and B will be the following: B B u (3) where A; 6. Proof: See []. Since A; 6 it is clear that the system is singular perturbed and in standard form. 5 Example We will derive a state space description from the bond graph in Figure. R C Se C 2 Figure : A bond graph example. The rst step is to introduce an R-element according to Figure 2. Its constitutive relation is assumed to be linear with parameter. R C Se C 2 R 2 Figure 2: The bond graph with causality relaxed by an R-element. Derive a state space description from the bond graph in Figure 2 using SCAP: _x _x 2? +r rc? c x x 2 + r x and x 2 are the states in the constitutive relations for the two C-elements and u is the eort given by the Se -element. u 5
7 Derive equations in both directions along the causal path between the two C-elements in Figure. x 2 c 2 c x _x? _x 2 Introduce two new state variables using the derived equations. c 2 c x? x 2 2 x + x 2 The transformation matrix is then: T c? Change the basis, and multiply the rst row with : A? c r+2cr+ r c(c+)r c(c+)r?? (c+)r (c+)r! ; B cr r Let! : _ 2? c+ c?? (c+)r (c+)r + u 2 r The state space description is thus _ 2? 2 + u (c + c 2 )r r There is a physical interpretation of the new state variable. If the bond graph models an electrical circuit, the causal path between the two C-elements shows that there are two capacitors connected in parallel. The two capacitors will behave as one with capacitance c + c 2 and with a state variable that describes the sum of the charges. This is exactly the structure of the state space equation above. References [] K. Edstrom. Singular perturbation analysis of a mode initialization algorithm for simulating mode switching system, long version. Technical report, Linkopings Universitet, 998. [2] D.C. Karnopp, D.L. Margolis, and R.C. Rosenberg. System dynamics, A unied approach. Wiley Interscience, 99. [3] P. Kokotovic, H.K. Khalil, and J. O'Reilly. Singular perturbation methods in control: Analysis and design. Academic Press, 986. [4] A. Rahmani and G. Dauphin-Tangui. Symbolic determination of state matrices from bond graph model with derivative causality. In Proc. of Computational Engineering in System Applications. IMACS, 998. [5] J. van Dijk. On the role of bond graph causality in modelling mechatronic systems. PhD thesis, University of Twente,
Modeling of Dynamic Systems: Notes on Bond Graphs Version 1.0 Copyright Diane L. Peters, Ph.D., P.E.
Modeling of Dynamic Systems: Notes on Bond Graphs Version 1.0 Copyright 2015 Diane L. Peters, Ph.D., P.E. Spring 2015 2 Contents 1 Overview of Dynamic Modeling 5 2 Bond Graph Basics 7 2.1 Causality.............................
More informationENGI9496 Modeling and Simulation of Dynamic Systems Bond Graphs
ENGI9496 Modeling and Simulation of Dynamic Systems Bond Graphs Topics covered so far: Analogies between mechanical (translation and rotation), fluid, and electrical systems o Review of domain-specific
More informationElectrical and Magnetic Modelling of a Power Transformer: A Bond Graph Approach
Vol:6, No:9, Electrical and Magnetic Modelling of a Power Transformer: A Bond Graph Approach Gilberto Gonzalez-A, Dunia Nuñez-P International Science Index, Electrical and Computer Engineering Vol:6, No:9,
More informationREGLERTEKNIK AUTOMATIC CONTROL LINKÖPING
Invariant Sets for a Class of Hybrid Systems Mats Jirstrand Department of Electrical Engineering Linkoping University, S-581 83 Linkoping, Sweden WWW: http://www.control.isy.liu.se Email: matsj@isy.liu.se
More informationResidual Bond Graph Sinks for Numerical Evaluation of Analytical Redundancy Relations in Model Based Single Fault Detection and Isolation
Residual Bond Graph Sinks for Numerical Evaluation of Analytical Redundancy Relations in Model Based Single Fault Detection and Isolation W. Borutzky Bonn-Rhein-Sieg University of Applied Sciences, D-53754
More informationAbstract. 1 Introduction. 2 Tearing TEARING IN BOND GRAPHS WITH DEPENDENT STORAGE ELEMENTS
CSA 96, IMACS Multiconference, Lille, France, July 9 12, 1996, pp. 1113 1119 1 TARING IN BOND GRAPHS WITH DPNDNT STORAG LMNTS Wolfgang Borutzky Dept. of Computer Science Cologne Polytechnic D-51643 Gummersbach
More informationAlternative Causality Assignment Procedures in Bond Graph for Mechanical Systems
Wilfrid Marquis-Favre e-mail: Wilfrid.Favre@lai.insa-lyon.fr Serge Scavarda Laboratoire d Automatique Industrielle, Institut National des Sciences Appliquées de Lyon, Bâtiment Antoine de Saint Exupéry,
More informationA relaxation of the strangeness index
echnical report from Automatic Control at Linköpings universitet A relaxation of the strangeness index Henrik idefelt, orkel Glad Division of Automatic Control E-mail: tidefelt@isy.liu.se, torkel@isy.liu.se
More informationEFFICIENT SIMULATION OF HYBRID SYSTEMS: AN APPLICATION TO ELECTRICAL POWER DISTRIBUTION SYSTEMS
EFFICIENT SIMULATION OF HYBRID SYSTEMS: AN APPLICATION TO ELECTRICAL POWER DISTRIBUTION SYSTEMS Indranil Roychoudhury, Matthew Daigle, Gautam Biswas, and Xenofon Koutsoukos Institute for Software Integrated
More informationMAS601 Design, Modeling & Simulation. Day 2
MAS601 Design, Modelling and Simulation of Mechatronic Systems, Semester 2, 2017. Page: 1 MAS601 Design, Modeling & Simulation Day 2 Analysis of causality and handling of algebraic loops to improve simulation
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 informationREGLERTEKNIK AUTOMATIC CONTROL LINKÖPING
Time-domain Identication of Dynamic Errors-in-variables Systems Using Periodic Excitation Signals Urban Forssell, Fredrik Gustafsson, Tomas McKelvey Department of Electrical Engineering Linkping University,
More informationCONSTITUTIVE HYBRID PROCESSES
CONSTITUTIVE HYBRID PROCESSES P.J.L. Cuijpers Technical University Eindhoven Eindhoven, Netherlands P.J.L.Cuijpers@tue.nl P.J. Mosterman The MathWorks, Inc. Natick, MA Pieter.Mosterman@mathworks.com J.F.
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 informationECE2262 Electric Circuit
ECE2262 Electric Circuit Chapter 7: FIRST AND SECOND-ORDER RL AND RC CIRCUITS Response to First-Order RL and RC Circuits Response to Second-Order RL and RC Circuits 1 2 7.1. Introduction 3 4 In dc steady
More information7.3 State Space Averaging!
7.3 State Space Averaging! A formal method for deriving the small-signal ac equations of a switching converter! Equivalent to the modeling method of the previous sections! Uses the state-space matrix description
More informationLS.1 Review of Linear Algebra
LS. LINEAR SYSTEMS LS.1 Review of Linear Algebra In these notes, we will investigate a way of handling a linear system of ODE s directly, instead of using elimination to reduce it to a single higher-order
More informationAircraft Pitch Attitude Control Using Backstepping
Aircraft Pitch Attitude Control Using Backstepping Ola Härkegård and Torkel Glad Department of Electrical Engineering Linköping University, SE-581 83 Linköping, Sweden WWW: http://www.control.isy.liu.se
More informationIn: Proc. BENELEARN-98, 8th Belgian-Dutch Conference on Machine Learning, pp 9-46, 998 Linear Quadratic Regulation using Reinforcement Learning Stephan ten Hagen? and Ben Krose Department of Mathematics,
More information3.4 Elementary Matrices and Matrix Inverse
Math 220: Summer 2015 3.4 Elementary Matrices and Matrix Inverse A n n elementary matrix is a matrix which is obtained from the n n identity matrix I n n by a single elementary row operation. Elementary
More informationAn efficient implementation of gradient and Hessian calculations of the coefficients of the characteristic polynomial of I XY
Technical report from Automatic Control at Linköpings universitet An efficient implementation of gradient and Hessian calculations of the coefficients of the characteristic polynomial of I XY Daniel Ankelhed
More informationBasic RL and RC Circuits R-L TRANSIENTS: STORAGE CYCLE. Engineering Collage Electrical Engineering Dep. Dr. Ibrahim Aljubouri
st Class Basic RL and RC Circuits The RL circuit with D.C (steady state) The inductor is short time at Calculate the inductor current for circuits shown below. I L E R A I L E R R 3 R R 3 I L I L R 3 R
More information7.6 The Inverse of a Square Matrix
7.6 The Inverse of a Square Matrix Copyright Cengage Learning. All rights reserved. What You Should Learn Verify that two matrices are inverses of each other. Use Gauss-Jordan elimination to find inverses
More informationChapter 2 Notes, Linear Algebra 5e Lay
Contents.1 Operations with Matrices..................................1.1 Addition and Subtraction.............................1. Multiplication by a scalar............................ 3.1.3 Multiplication
More information2005 AP PHYSICS C: ELECTRICITY AND MAGNETISM FREE-RESPONSE QUESTIONS
2005 AP PHYSICS C: ELECTRICITY AND MAGNETISM In the circuit shown above, resistors 1 and 2 of resistance R 1 and R 2, respectively, and an inductor of inductance L are connected to a battery of emf e and
More informationQualitative Dynamic Behavior of Physical System Models With Algebraic Loops
Presented at the Eleventh International Workshop on Principles of Diagnosis, Mexico. June 2000 Qualitative Dynamic Behavior of Physical System Models With Algebraic Loops Pieter J. Mosterman Institute
More informationModel reduction in the physical domain
48 Model reduction in the physical domain A Y Orbak, O S Turkay2*,**, E Eskinat and K Youcef-Toumi3 Department of Industrial Engineering, Uludag University, Bursa, Turkey 2Department of Mechanical Engineering,
More informationMAT 1332: CALCULUS FOR LIFE SCIENCES. Contents. 1. Review: Linear Algebra II Vectors and matrices Definition. 1.2.
MAT 1332: CALCULUS FOR LIFE SCIENCES JING LI Contents 1 Review: Linear Algebra II Vectors and matrices 1 11 Definition 1 12 Operations 1 2 Linear Algebra III Inverses and Determinants 1 21 Inverse Matrices
More informationREGLERTEKNIK AUTOMATIC CONTROL LINKÖPING
Expectation Maximization Segmentation Niclas Bergman Department of Electrical Engineering Linkoping University, S-581 83 Linkoping, Sweden WWW: http://www.control.isy.liu.se Email: niclas@isy.liu.se October
More informationLecture 9: Elementary Matrices
Lecture 9: Elementary Matrices Review of Row Reduced Echelon Form Consider the matrix A and the vector b defined as follows: 1 2 1 A b 3 8 5 A common technique to solve linear equations of the form Ax
More informationMatrix Algebra. Learning Objectives. Size of Matrix
Matrix Algebra 1 Learning Objectives 1. Find the sum and difference of two matrices 2. Find scalar multiples of a matrix 3. Find the product of two matrices 4. Find the inverse of a matrix 5. Solve a system
More informationPARAMETERIZATION OF STATE FEEDBACK GAINS FOR POLE PLACEMENT
PARAMETERIZATION OF STATE FEEDBACK GAINS FOR POLE PLACEMENT Hans Norlander Systems and Control, Department of Information Technology Uppsala University P O Box 337 SE 75105 UPPSALA, Sweden HansNorlander@ituuse
More informationMasteringPhysics: Assignment Print View. Problem 30.50
Page 1 of 15 Assignment Display Mode: View Printable Answers phy260s08 homework 13 Due at 11:00pm on Wednesday, May 14, 2008 View Grading Details Problem 3050 Description: A 15-cm-long nichrome wire is
More informationLinear Regression and Its Applications
Linear Regression and Its Applications Predrag Radivojac October 13, 2014 Given a data set D = {(x i, y i )} n the objective is to learn the relationship between features and the target. We usually start
More information(amperes) = (coulombs) (3.1) (seconds) Time varying current. (volts) =
3 Electrical Circuits 3. Basic Concepts Electric charge coulomb of negative change contains 624 0 8 electrons. Current ampere is a steady flow of coulomb of change pass a given point in a conductor in
More informationFrom Electrical Switched Networks to Hybrid Automata. Alessandro Cimatti 1, Sergio Mover 2, Mirko Sessa 1,3
1 2 3 From Electrical Switched Networks to Hybrid Automata Alessandro Cimatti 1, Sergio Mover 2, Mirko Sessa 1,3 Multidomain physical systems Multiple physical domains: electrical hydraulic mechanical
More informationMA 138 Calculus 2 with Life Science Applications Matrices (Section 9.2)
MA 38 Calculus 2 with Life Science Applications Matrices (Section 92) Alberto Corso albertocorso@ukyedu Department of Mathematics University of Kentucky Friday, March 3, 207 Identity Matrix and Inverse
More informationHybrid Simulation of Qualitative Bond Graph Model
Hybrid Simulation of Qualitative Bond Graph Model C.H. LO, Y.K. WONG and A.B. RAD Department of Electrical Engineering The Hong Kong Polytechnic University Hung Hom, Kowloon, Hong Kong Abstract: - Qualitative
More informationConceptual Questions for Review
Conceptual Questions for Review Chapter 1 1.1 Which vectors are linear combinations of v = (3, 1) and w = (4, 3)? 1.2 Compare the dot product of v = (3, 1) and w = (4, 3) to the product of their lengths.
More informationChapter 27. Circuits
Chapter 27 Circuits 1 1. Pumping Chagres We need to establish a potential difference between the ends of a device to make charge carriers follow through the device. To generate a steady flow of charges,
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 informationA Software of Generating a Symbolic Circuit Model with Computers for Wireless Power Transmission System
A Software of Generating a Symbolic Circuit Model with Computers for Wireless Power Transmission System Takuya Hirata Department of Materials and Informatics Interdisciplinary Graduate School of Agriculture
More informationBackstepping Designs for Aircraft Control What is there to gain?
Backstepping Designs for Aircraft Control What is there to gain? Ola Härkegård Division of Automatic Control Department of Electrical Engineering Linköpings universitet, SE-581 83 Linköping, Sweden WWW:
More informationAnalysis of the regularity, pointwise completeness and pointwise generacy of descriptor linear electrical circuits
Computer Applications in Electrical Engineering Vol. 4 Analysis o the regularity pointwise completeness pointwise generacy o descriptor linear electrical circuits Tadeusz Kaczorek Białystok University
More informationModule 6 : Solving Ordinary Differential Equations - Initial Value Problems (ODE-IVPs) Section 3 : Analytical Solutions of Linear ODE-IVPs
Module 6 : Solving Ordinary Differential Equations - Initial Value Problems (ODE-IVPs) Section 3 : Analytical Solutions of Linear ODE-IVPs 3 Analytical Solutions of Linear ODE-IVPs Before developing numerical
More informationAP Physics C. Electric Circuits III.C
AP Physics C Electric Circuits III.C III.C.1 Current, Resistance and Power The direction of conventional current Suppose the cross-sectional area of the conductor changes. If a conductor has no current,
More informationMatrix Arithmetic. (Multidimensional Math) Institute for Advanced Analytics. Shaina Race, PhD
Matrix Arithmetic Multidimensional Math) Shaina Race, PhD Institute for Advanced Analytics. TABLE OF CONTENTS Element-wise Operations Linear Combinations of Matrices and Vectors. Vector Multiplication
More informationoutput H = 2*H+P H=2*(H-P)
Ecient Algorithms for Multiplication on Elliptic Curves by Volker Muller TI-9/97 22. April 997 Institut fur theoretische Informatik Ecient Algorithms for Multiplication on Elliptic Curves Volker Muller
More informationCircuit theory for unusual inductor behaviour
Circuit theory for unusual inductor behaviour Horst Eckardt, Douglas W. Lindstrom Alpha Institute for Advanced Studies AIAS and Unied Physics Institute of Technology UPITEC July 9, 2015 Abstract The experiments
More informationKnowledge Discovery and Data Mining 1 (VO) ( )
Knowledge Discovery and Data Mining 1 (VO) (707.003) Review of Linear Algebra Denis Helic KTI, TU Graz Oct 9, 2014 Denis Helic (KTI, TU Graz) KDDM1 Oct 9, 2014 1 / 74 Big picture: KDDM Probability Theory
More informationREGLERTEKNIK AUTOMATIC CONTROL LINKÖPING
Recursive Least Squares and Accelerated Convergence in Stochastic Approximation Schemes Lennart Ljung Department of Electrical Engineering Linkping University, S-581 83 Linkping, Sweden WWW: http://www.control.isy.liu.se
More informationREGLERTEKNIK AUTOMATIC CONTROL LINKÖPING
An Alternative Motivation for the Indirect Approach to Closed-loop Identication Lennart Ljung and Urban Forssell Department of Electrical Engineering Linkping University, S-581 83 Linkping, Sweden WWW:
More informationBond Graphs. Peter Breedveld University of Twente, EWI/CE, P.O. Box 217, 7500 AE Enschede, Netherlands,
Bond Graphs Peter Breedveld University of Twente, EWI/CE, P.O. Box 217, 7500 AE Enschede, Netherlands, p.c.breedveld@utwente.nl 1 Introduction The topic area that has become commonly known as bond graph
More informationFactorization of singular integer matrices
Factorization of singular integer matrices Patrick Lenders School of Mathematics, Statistics and Computer Science, University of New England, Armidale, NSW 2351, Australia Jingling Xue School of Computer
More informationComputation of a canonical form for linear differential-algebraic equations
Computation of a canonical form for linear differential-algebraic equations Markus Gerdin Division of Automatic Control Department of Electrical Engineering Linköpings universitet, SE-581 83 Linköping,
More informationRepresentation of a general composition of Dirac structures
Representation of a general composition of Dirac structures Carles Batlle, Imma Massana and Ester Simó Abstract We provide explicit representations for the Dirac structure obtained from an arbitrary number
More informationTHERE are numerous number of products powered by
Symbolic Computer-Aided Design for Wireless Power Transmission Takuya Hirata, Univ. Miyazaki, Kazuya Yamaguchi, Univ. Miyazaki, Yuta Yamamoto, Univ. Miyazaki and Ichijo Hodaka, Univ. Miyazaki Abstract
More informationComputing Phase Noise Eigenfunctions Directly from Steady-State Jacobian Matrices
Computing Phase Noise Eigenfunctions Directly from Steady-State Jacobian Matrices Alper Demir David Long Jaijeet Roychowdhury Bell Laboratories Murray Hill New Jersey USA Abstract The main effort in oscillator
More informationPhysics 102 Spring 2006: Final Exam Multiple-Choice Questions
Last Name: First Name: Physics 102 Spring 2006: Final Exam Multiple-Choice Questions For questions 1 and 2, refer to the graph below, depicting the potential on the x-axis as a function of x V x 60 40
More informationEects of small delays on stability of singularly perturbed systems
Automatica 38 (2002) 897 902 www.elsevier.com/locate/automatica Technical Communique Eects of small delays on stability of singularly perturbed systems Emilia Fridman Department of Electrical Engineering
More informationLectures on Linear Algebra for IT
Lectures on Linear Algebra for IT by Mgr Tereza Kovářová, PhD following content of lectures by Ing Petr Beremlijski, PhD Department of Applied Mathematics, VSB - TU Ostrava Czech Republic 3 Inverse Matrix
More informationMath 240 Calculus III
The Calculus III Summer 2015, Session II Wednesday, July 8, 2015 Agenda 1. of the determinant 2. determinants 3. of determinants What is the determinant? Yesterday: Ax = b has a unique solution when A
More informationFirst-order transient
EIE209 Basic Electronics First-order transient Contents Inductor and capacitor Simple RC and RL circuits Transient solutions Constitutive relation An electrical element is defined by its relationship between
More informationMath Camp Notes: Linear Algebra I
Math Camp Notes: Linear Algebra I Basic Matrix Operations and Properties Consider two n m matrices: a a m A = a n a nm Then the basic matrix operations are as follows: a + b a m + b m A + B = a n + b n
More informationEC Control Systems- Question bank
MODULE I Topic Question mark Automatic control & modeling, Transfer function Write the merits and demerits of open loop and closed loop Month &Year May 12 Regula tion Compare open loop system with closed
More informationAC vs. DC Circuits. Constant voltage circuits. The voltage from an outlet is alternating voltage
Circuits AC vs. DC Circuits Constant voltage circuits Typically referred to as direct current or DC Computers, logic circuits, and battery operated devices are examples of DC circuits The voltage from
More information12. Perturbed Matrices
MAT334 : Applied Linear Algebra Mike Newman, winter 208 2. Perturbed Matrices motivation We want to solve a system Ax = b in a context where A and b are not known exactly. There might be experimental errors,
More informationKevin James. MTHSC 3110 Section 2.2 Inverses of Matrices
MTHSC 3110 Section 2.2 Inverses of Matrices Definition Suppose that T : R n R m is linear. We will say that T is invertible if for every b R m there is exactly one x R n so that T ( x) = b. Note If T is
More informationINRIA Rocquencourt, Le Chesnay Cedex (France) y Dept. of Mathematics, North Carolina State University, Raleigh NC USA
Nonlinear Observer Design using Implicit System Descriptions D. von Wissel, R. Nikoukhah, S. L. Campbell y and F. Delebecque INRIA Rocquencourt, 78 Le Chesnay Cedex (France) y Dept. of Mathematics, North
More informationNumerical Methods Lecture 2 Simultaneous Equations
CGN 42 - Computer Methods Numerical Methods Lecture 2 Simultaneous Equations Topics: matrix operations solving systems of equations Matrix operations: Adding / subtracting Transpose Multiplication Adding
More informationThe norms can also be characterized in terms of Riccati inequalities.
9 Analysis of stability and H norms Consider the causal, linear, time-invariant system ẋ(t = Ax(t + Bu(t y(t = Cx(t Denote the transfer function G(s := C (si A 1 B. Theorem 85 The following statements
More informationLesson 3. Inverse of Matrices by Determinants and Gauss-Jordan Method
Module 1: Matrices and Linear Algebra Lesson 3 Inverse of Matrices by Determinants and Gauss-Jordan Method 3.1 Introduction In lecture 1 we have seen addition and multiplication of matrices. Here we shall
More informationSwitched Mode Power Conversion Prof. L. Umanand Department of Electronics Systems Engineering Indian Institute of Science, Bangalore
Switched Mode Power Conversion Prof. L. Umanand Department of Electronics Systems Engineering Indian Institute of Science, Bangalore Lecture - 19 Modeling DC-DC convertors Good day to all of you. Today,
More informationM. C. Escher: Waterfall. 18/9/2015 [tsl425 1/29]
M. C. Escher: Waterfall 18/9/2015 [tsl425 1/29] Direct Current Circuit Consider a wire with resistance R = ρl/a connected to a battery. Resistor rule: In the direction of I across a resistor with resistance
More informationGraphing Linear Inequalities
Graphing Linear Inequalities Linear Inequalities in Two Variables: A linear inequality in two variables is an inequality that can be written in the general form Ax + By < C, where A, B, and C are real
More informationPhysics 212 Midterm 2 Form A
1. A wire contains a steady current of 2 A. The charge that passes a cross section in 2 s is: A. 3.2 10-19 C B. 6.4 10-19 C C. 1 C D. 2 C E. 4 C 2. In a Physics 212 lab, Jane measures the current versus
More information3.3. SYSTEMS OF ODES 1. y 0 " 2y" y 0 + 2y = x1. x2 x3. x = y(t) = c 1 e t + c 2 e t + c 3 e 2t. _x = A x + f; x(0) = x 0.
.. SYSTEMS OF ODES. Systems of ODEs MATH 94 FALL 98 PRELIM # 94FA8PQ.tex.. a) Convert the third order dierential equation into a rst oder system _x = A x, with y " y" y + y = x = @ x x x b) The equation
More informationVII Selected Topics. 28 Matrix Operations
VII Selected Topics Matrix Operations Linear Programming Number Theoretic Algorithms Polynomials and the FFT Approximation Algorithms 28 Matrix Operations We focus on how to multiply matrices and solve
More information2.004 Dynamics and Control II Spring 2008
MT OpenCourseWare http://ocwmitedu 200 Dynamics and Control Spring 200 For information about citing these materials or our Terms of Use, visit: http://ocwmitedu/terms Massachusetts nstitute of Technology
More informationMath 2331 Linear Algebra
2.2 The Inverse of a Matrix Math 2331 Linear Algebra 2.2 The Inverse of a Matrix Shang-Huan Chiu Department of Mathematics, University of Houston schiu@math.uh.edu math.uh.edu/ schiu/ Shang-Huan Chiu,
More informationVector Space Basics. 1 Abstract Vector Spaces. 1. (commutativity of vector addition) u + v = v + u. 2. (associativity of vector addition)
Vector Space Basics (Remark: these notes are highly formal and may be a useful reference to some students however I am also posting Ray Heitmann's notes to Canvas for students interested in a direct computational
More informationRC Circuits (32.9) Neil Alberding (SFU Physics) Physics 121: Optics, Electricity & Magnetism Spring / 1
(32.9) We have only been discussing DC circuits so far. However, using a capacitor we can create an RC circuit. In this example, a capacitor is charged but the switch is open, meaning no current flows.
More informationAP Physics C Mechanics Objectives
AP Physics C Mechanics Objectives I. KINEMATICS A. Motion in One Dimension 1. The relationships among position, velocity and acceleration a. Given a graph of position vs. time, identify or sketch a graph
More informationTotal 170. Name. Final Examination M340L-CS
1 10 2 10 3 15 4 5 5 10 6 10 7 20 8 10 9 20 10 25 11 10 12 10 13 15 Total 170 Final Examination Name M340L-CS 1. Use extra paper to determine your solutions then neatly transcribe them onto these sheets.
More informationTHE USE OF ROBUST OBSERVERS IN THE SIMULATION OF GAS SUPPLY NETWORKS
UNIVERSITY OF READING DEPARTMENT OF MATHEMATICS THE USE OF ROBUST OBSERVERS IN THE SIMULATION OF GAS SUPPLY NETWORKS by SIMON M. STRINGER This thesis is submitted forthe degree of Doctorof Philosophy NOVEMBER
More informationGabriel Kron's biography here.
Gabriel Kron, Electric Circuit Model of the Schrödinger Equation, 1945 - Component of :... Page 1 of 12 {This website: Please note: The following article is complete; it has been put into ASCII due to
More informationNotes on Dantzig-Wolfe decomposition and column generation
Notes on Dantzig-Wolfe decomposition and column generation Mette Gamst November 11, 2010 1 Introduction This note introduces an exact solution method for mathematical programming problems. The method is
More informationMEM 255 Introduction to Control Systems: Modeling & analyzing systems
MEM 55 Introduction to Control Systems: Modeling & analyzing systems Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University Outline The Pendulum Micro-machined capacitive accelerometer
More informationANALYTICAL MATHEMATICS FOR APPLICATIONS 2018 LECTURE NOTES 3
ANALYTICAL MATHEMATICS FOR APPLICATIONS 2018 LECTURE NOTES 3 ISSUED 24 FEBRUARY 2018 1 Gaussian elimination Let A be an (m n)-matrix Consider the following row operations on A (1) Swap the positions any
More informationExamples: u = is a vector in 2. is a vector in 5.
3 Vectors and vector equations We'll carefully define vectors, algebraic operations on vectors and geometric interpretations of these operations, in terms of displacements These ideas will eventually give
More informationModeling & Simulation 2018 Lecture 12. Simulations
Modeling & Simulation 2018 Lecture 12. Simulations Claudio Altafini Automatic Control, ISY Linköping University, Sweden Summary of lecture 7-11 1 / 32 Models of complex systems physical interconnections,
More informationLinear Algebra The Inverse of a Matrix
Linear Algebra The Inverse of a Matrix Dr. Bisher M. Iqelan biqelan@iugaza.edu.ps Department of Mathematics The Islamic University of Gaza 2017-2018, Semester 2 Dr. Bisher M. Iqelan (IUG) Sec.2.2: The
More informationLinear graph theory. Basic definitions of linear graphs
Linear graph theory Linear graph theory, a branch of combinatorial mathematics has proved to be a useful tool for the study of large or complex systems. Leonhard Euler wrote perhaps the first paper on
More informationVersion 001 CIRCUITS holland (1290) 1
Version CIRCUITS holland (9) This print-out should have questions Multiple-choice questions may continue on the next column or page find all choices before answering AP M 99 MC points The power dissipated
More information1. Introduction As is well known, the bosonic string can be described by the two-dimensional quantum gravity coupled with D scalar elds, where D denot
RIMS-1161 Proof of the Gauge Independence of the Conformal Anomaly of Bosonic String in the Sense of Kraemmer and Rebhan Mitsuo Abe a; 1 and Noboru Nakanishi b; 2 a Research Institute for Mathematical
More informationSystems of Linear ODEs
P a g e 1 Systems of Linear ODEs Systems of ordinary differential equations can be solved in much the same way as discrete dynamical systems if the differential equations are linear. We will focus here
More informationSymmetric matrices and dot products
Symmetric matrices and dot products Proposition An n n matrix A is symmetric iff, for all x, y in R n, (Ax) y = x (Ay). Proof. If A is symmetric, then (Ax) y = x T A T y = x T Ay = x (Ay). If equality
More information文件中找不不到关系 ID 为 rid3 的图像部件. Review of LINEAR ALGEBRA I. TA: Yujia Xie CSE 6740 Computational Data Analysis. Georgia Institute of Technology
文件中找不不到关系 ID 为 rid3 的图像部件 Review of LINEAR ALGEBRA I TA: Yujia Xie CSE 6740 Computational Data Analysis Georgia Institute of Technology 1. Notations A R $ & : a matrix with m rows and n columns, where
More informationChapter 2. Matrix Arithmetic. Chapter 2
Matrix Arithmetic Matrix Addition and Subtraction Addition and subtraction act element-wise on matrices. In order for the addition/subtraction (A B) to be possible, the two matrices A and B must have the
More informationSolutions of Linear system, vector and matrix equation
Goals: Solutions of Linear system, vector and matrix equation Solutions of linear system. Vectors, vector equation. Matrix equation. Math 112, Week 2 Suggested Textbook Readings: Sections 1.3, 1.4, 1.5
More information