NLControl - package for modelling, analysis and synthesis of nonlinear control systems

Size: px
Start display at page:

Download "NLControl - package for modelling, analysis and synthesis of nonlinear control systems"

Transcription

1 NLControl - package for modelling, analysis and synthesis of nonlinear control systems Maris Tõnso Institute of Cybernetics Tallinn University of Technology Estonia maris@cc.ioc.ee May, 2013

2 Contents 1 Difference/differential field Sequence {H k } Problem 1: Linearization 2 Non-commutative polynomials Polynomial system Problem 2: Reduction Problem 3: Realization 2 / 15

3 The main idea Allows to address a wide range of nonlinear control problems 1 Equations One-: The nonlinear system is globally linearised and described by its infinitesimal representation. 2 Subspaces of differential one- are computed. 3 Intrinsic necessary and sufficient conditions for the existence of the solution are formulated in terms of the one-. 4 One- Eqations: Integration of one-. Integrability restrictions. 5 Final solution is found. 3 / 15

4 [E. Aranda-Bricaire, Ü. Kotta, C. Moog, 1996] x(t + 1) = f (x(t), u(t)), x(t) R n, u(t) R m K - the field of meromorphic functions in a finite number of variables from the infinite set {x(0), u j (t), j = 1,..., m, t 0} Forward-shift operator δ : K K δx i (t) = x i (t + 1) = f i ( ), i = 1,..., n δu j (t) = u j (t + 1), j = 1,..., m δ is injective (K, δ) is a difference field 4 / 15

5 Define a difference vector space E := span K {dϕ ϕ K}. The elements of E are called one-. δ : K K induces δ : E E by a i dϕ i (δa i )d(δϕ i ), a i, ϕ i K. i i δ and d are commutative The relative degree r of a one-form ω E is defined to be the least integer such that δ r ω / span K {dx(0)}. If such an integer does not exist, we set r =. 5 / 15

6 Sequence of subspaces {H k } Define a sequence of subspaces {H k } of E: H 1 = span K {dx(0)} H k+1 = {ω H k δω H k }, k 1. Each H k contains the one- whose relative degree r k {H k } is a decreasing sequence. H 1 H k H k +1 = H k +2 = =: H 6 / 15

7 Problem 1: Static state feedback linearization x(t + 1) = f (x(t), u(t)), x(t) R n, u(t) R m Theorem: The system is linearizable by static stat feedback iff All subspaces H k are integrable H = {0}. To check integrability: Frobenius Theorem: The subspace span{ω 1,..., ω k } is integrable iff dω l ω 1 ω k = 0 for all l = 1,..., k. 7 / 15

8 Non-commutative polynomial ring Noncommutative polynomials act on one- as operators allow to repersent computations in a more compact way. Introduce the left polynomial ring K[, δ]. A polynomial from K[, δ] has the form a = n a i n i, i=0 a i K polynomial indeterminate. 8 / 15

9 Non-commutative polynomial ring (continuation) Multiplication in the ring K[, δ] is not commutative since a a, a K K[, δ]. Multiplication is defined by the commutation rule: Example: a = δ(a), a K u(t) = u(t + 1) u(t) = u(t) The ring K[, δ] is an Ore ring. In the following we interprete as δ. 9 / 15

10 Polynomial system i/o equation y(t + n) = φ(y(t),..., y(t + n 1), u(t),..., u(t + s)) Differentiate the i/o equations Use definitions dy(t +j) = j dy(t), du(t +r) = r du(t) p( )dy(t) = q( )du(t) n 1 p = n φ y(t + j) j q = j=0 s r=0 φ u(t + r) r 10 / 15

11 Problem II - Reduction Problem setting y(t + n) = φ(y(t),..., y(t + n 1), u(t),..., u(t + s)) Def: A function ϕ r K is said to be an autonomous variable for a system, if there exist an integer µ 1 and a non-zero function F such that F (ϕ r, δϕ r,..., δ µ ϕ r ) = 0. The system is said to be irreducible if there does not exist any non-zero autonomous variable for it. The system is reducible if it there exists a lower order system which is transfer equivalent with the original system. 11 / 15

12 Problem II - Reduction Solution p( )dy(t) q( )du(t) = 0 Theorem: Nonlinear system is irreducible, iff p and q are relatively left prime. Let g = LeftGCD(p, q). Then p( )dy(t) q( )du(t) = g(s)[ p( )dy(t) q( )du(t)] Denote dψ := [ p( )dy(t) q( )du(t)]. Then ψ = 0 is the reduced equation. 12 / 15

13 Problem 3 - Realization Problem setting Input-output equations y(t + n) = φ(y(t),..., y(t + n 1), u(t),..., u(t + s)) State equations x(t + 1) = f (x(t), u(t)) y(t) = h(x(t)) We are looking for minimal and observable realization For arbitrary i/o equation the state space form does not necessarily exist. 13 / 15

14 Problem III - Realization Solution p( )dy(t) q( )du(t) = 0 For r = r k k + + r 1 + r 0 Define cut-and-shift opeartor δc 1 (r) := δ 1 (r k ) k δ 1 (r 2 ) + δ 1 (r 1 ) [ ] dy ω l := δc l [p, q] for l = 1,..., n du Theorem: The i/o equation is realizable if the subspace H = {ω 1,..., ω n } is integrable. s of the state coordinates dx i can be found as integrable linear combinations of ω 1,..., ω n. 14 / 15

15 More Problems Addressed by algebraic approach and implemented in NLControl Tranforming the state equations into observer form Model matching problems Flatness Input-output linearization Disturbance decoupling. 15 / 15

Transfer functions of discrete-time nonlinear control systems

Transfer functions of discrete-time nonlinear control systems Proc. Estonian Acad. Sci. Phys. Math., 2007, 56, 4, 322 335 a Transfer functions of discrete-time nonlinear control systems Miroslav Halás a and Ülle Kotta b Institute of Control and Industrial Informatics,

More information

Disturbance decoupling by measurement feedback

Disturbance decoupling by measurement feedback Preprints of the 19th Word Congress The Internationa Federation of Automatic Contro Disturbance decouping by measurement feedback Arvo Kadmäe, Üe Kotta Institute of Cybernetics at TUT, Akadeemia tee 21,

More information

WEBMATHEMATICA BASED TOOLS FOR NONLINEAR CONTROL SYSTEMS

WEBMATHEMATICA BASED TOOLS FOR NONLINEAR CONTROL SYSTEMS WEBMATHEMATICA BASED TOOLS FOR NONLINEAR CONTROL SYSTEMS Heli Rennik, Maris Tõnso and Ülle Kotta Institute of Cybernetics, Tallinn University of Technology, Akadeemia tee 21, Tallinn, 12818, Estonia heli.rennik@mail.ee,

More information

Chapter 6. Differentially Flat Systems

Chapter 6. Differentially Flat Systems Contents CAS, Mines-ParisTech 2008 Contents Contents 1, Linear Case Introductory Example: Linear Motor with Appended Mass General Solution (Linear Case) Contents Contents 1, Linear Case Introductory Example:

More information

Chap. 3. Controlled Systems, Controllability

Chap. 3. Controlled Systems, Controllability Chap. 3. Controlled Systems, Controllability 1. Controllability of Linear Systems 1.1. Kalman s Criterion Consider the linear system ẋ = Ax + Bu where x R n : state vector and u R m : input vector. A :

More information

Solutions. Problems of Chapter 1

Solutions. Problems of Chapter 1 Solutions Problems of Chapter 1 1.1 A real square matrix A IR n n is invertible if and only if its determinant is non zero. The determinant of A is a polynomial in the entries a ij of A, whose set of zeros

More information

Vector Space Concepts

Vector Space Concepts Vector Space Concepts ECE 174 Introduction to Linear & Nonlinear Optimization Ken Kreutz-Delgado ECE Department, UC San Diego Ken Kreutz-Delgado (UC San Diego) ECE 174 Fall 2016 1 / 25 Vector Space Theory

More information

Spin(10,1)-metrics with a parallel null spinor and maximal holonomy

Spin(10,1)-metrics with a parallel null spinor and maximal holonomy Spin(10,1)-metrics with a parallel null spinor and maximal holonomy 0. Introduction. The purpose of this addendum to the earlier notes on spinors is to outline the construction of Lorentzian metrics in

More information

COMPLEX MULTIPLICATION: LECTURE 13

COMPLEX MULTIPLICATION: LECTURE 13 COMPLEX MULTIPLICATION: LECTURE 13 Example 0.1. If we let C = P 1, then k(c) = k(t) = k(c (q) ) and the φ (t) = t q, thus the extension k(c)/φ (k(c)) is of the form k(t 1/q )/k(t) which as you may recall

More information

WebMathematica-based tools for discrete-time nonlinear control systems

WebMathematica-based tools for discrete-time nonlinear control systems Proceedings of the Estonian Academy of Sciences, 2009, 58, 4, 224 240 doi: 10.3176/proc.2009.4.04 Available online at www.eap.ee/proceedings WebMathematica-based tools for discrete-time nonlinear control

More information

WELL-FORMED DYNAMICS UNDER QUASI-STATIC STATE FEEDBACK

WELL-FORMED DYNAMICS UNDER QUASI-STATIC STATE FEEDBACK GEOMETRY IN NONLINEAR CONTROL AND DIFFERENTIAL INCLUSIONS BANACH CENTER PUBLICATIONS, VOLUME 32 INSTITUTE OF MATHEMATICS POLISH ACADEMY OF SCIENCES WARSZAWA 1995 WELL-FORMED DYNAMICS UNDER QUASI-STATIC

More information

Chap. 1. Some Differential Geometric Tools

Chap. 1. Some Differential Geometric Tools Chap. 1. Some Differential Geometric Tools 1. Manifold, Diffeomorphism 1.1. The Implicit Function Theorem ϕ : U R n R n p (0 p < n), of class C k (k 1) x 0 U such that ϕ(x 0 ) = 0 rank Dϕ(x) = n p x U

More information

TANGENT VECTORS. THREE OR FOUR DEFINITIONS.

TANGENT VECTORS. THREE OR FOUR DEFINITIONS. TANGENT VECTORS. THREE OR FOUR DEFINITIONS. RMONT We define and try to understand the tangent space of a manifold Q at a point q, as well as vector fields on a manifold. The tangent space at q Q is a real

More information

Math 121 Homework 4: Notes on Selected Problems

Math 121 Homework 4: Notes on Selected Problems Math 121 Homework 4: Notes on Selected Problems 11.2.9. If W is a subspace of the vector space V stable under the linear transformation (i.e., (W ) W ), show that induces linear transformations W on W

More information

Di erential Algebraic Geometry, Part I

Di erential Algebraic Geometry, Part I Di erential Algebraic Geometry, Part I Phyllis Joan Cassidy City College of CUNY Fall 2007 Phyllis Joan Cassidy (Institute) Di erential Algebraic Geometry, Part I Fall 2007 1 / 46 Abstract Di erential

More information

Qualifying Exams I, 2014 Spring

Qualifying Exams I, 2014 Spring Qualifying Exams I, 2014 Spring 1. (Algebra) Let k = F q be a finite field with q elements. Count the number of monic irreducible polynomials of degree 12 over k. 2. (Algebraic Geometry) (a) Show that

More information

INTRODUCTION TO ALGEBRAIC GEOMETRY

INTRODUCTION TO ALGEBRAIC GEOMETRY INTRODUCTION TO ALGEBRAIC GEOMETRY WEI-PING LI 1 Preliminary of Calculus on Manifolds 11 Tangent Vectors What are tangent vectors we encounter in Calculus? (1) Given a parametrised curve α(t) = ( x(t),

More information

MCE693/793: Analysis and Control of Nonlinear Systems

MCE693/793: Analysis and Control of Nonlinear Systems MCE693/793: Analysis and Control of Nonlinear Systems Input-Output and Input-State Linearization Zero Dynamics of Nonlinear Systems Hanz Richter Mechanical Engineering Department Cleveland State University

More information

Chap 4. State-Space Solutions and

Chap 4. State-Space Solutions and Chap 4. State-Space Solutions and Realizations Outlines 1. Introduction 2. Solution of LTI State Equation 3. Equivalent State Equations 4. Realizations 5. Solution of Linear Time-Varying (LTV) Equations

More information

Stabilization and Passivity-Based Control

Stabilization and Passivity-Based Control DISC Systems and Control Theory of Nonlinear Systems, 2010 1 Stabilization and Passivity-Based Control Lecture 8 Nonlinear Dynamical Control Systems, Chapter 10, plus handout from R. Sepulchre, Constructive

More information

B 1 = {B(x, r) x = (x 1, x 2 ) H, 0 < r < x 2 }. (a) Show that B = B 1 B 2 is a basis for a topology on X.

B 1 = {B(x, r) x = (x 1, x 2 ) H, 0 < r < x 2 }. (a) Show that B = B 1 B 2 is a basis for a topology on X. Math 6342/7350: Topology and Geometry Sample Preliminary Exam Questions 1. For each of the following topological spaces X i, determine whether X i and X i X i are homeomorphic. (a) X 1 = [0, 1] (b) X 2

More information

Algebra Homework, Edition 2 9 September 2010

Algebra Homework, Edition 2 9 September 2010 Algebra Homework, Edition 2 9 September 2010 Problem 6. (1) Let I and J be ideals of a commutative ring R with I + J = R. Prove that IJ = I J. (2) Let I, J, and K be ideals of a principal ideal domain.

More information

Topic # Feedback Control

Topic # Feedback Control Topic #11 16.31 Feedback Control State-Space Systems State-space model features Observability Controllability Minimal Realizations Copyright 21 by Jonathan How. 1 Fall 21 16.31 11 1 State-Space Model Features

More information

Robust Control of Linear Quantum Systems

Robust Control of Linear Quantum Systems B. Ross Barmish Workshop 2009 1 Robust Control of Linear Quantum Systems Ian R. Petersen School of ITEE, University of New South Wales @ the Australian Defence Force Academy, Based on joint work with Matthew

More information

Krull Dimension and Going-Down in Fixed Rings

Krull Dimension and Going-Down in Fixed Rings David Dobbs Jay Shapiro April 19, 2006 Basics R will always be a commutative ring and G a group of (ring) automorphisms of R. We let R G denote the fixed ring, that is, Thus R G is a subring of R R G =

More information

Computable Differential Fields

Computable Differential Fields Computable Differential Fields Russell Miller Queens College & CUNY Graduate Center New York, NY Mathematisches Forschungsinstitut Oberwolfach 10 February 2012 (Joint work with Alexey Ovchinnikov & Dmitry

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

ORDINARY DIFFERENTIAL EQUATIONS

ORDINARY DIFFERENTIAL EQUATIONS ORDINARY DIFFERENTIAL EQUATIONS GABRIEL NAGY Mathematics Department, Michigan State University, East Lansing, MI, 48824. JANUARY 3, 25 Summary. This is an introduction to ordinary differential equations.

More information

Lecture 5: Linear Systems. Transfer functions. Frequency Domain Analysis. Basic Control Design.

Lecture 5: Linear Systems. Transfer functions. Frequency Domain Analysis. Basic Control Design. ISS0031 Modeling and Identification Lecture 5: Linear Systems. Transfer functions. Frequency Domain Analysis. Basic Control Design. Aleksei Tepljakov, Ph.D. September 30, 2015 Linear Dynamic Systems Definition

More information

Exam in Systems Engineering/Process Control

Exam in Systems Engineering/Process Control Department of AUTOMATIC CONTROL Exam in Systems Engineering/Process Control 7-6- Points and grading All answers must include a clear motivation. Answers may be given in English or Swedish. The total number

More information

ẋ n = f n (x 1,...,x n,u 1,...,u m ) (5) y 1 = g 1 (x 1,...,x n,u 1,...,u m ) (6) y p = g p (x 1,...,x n,u 1,...,u m ) (7)

ẋ n = f n (x 1,...,x n,u 1,...,u m ) (5) y 1 = g 1 (x 1,...,x n,u 1,...,u m ) (6) y p = g p (x 1,...,x n,u 1,...,u m ) (7) EEE582 Topical Outline A.A. Rodriguez Fall 2007 GWC 352, 965-3712 The following represents a detailed topical outline of the course. It attempts to highlight most of the key concepts to be covered and

More information

Serge Ballif January 18, 2008

Serge Ballif January 18, 2008 ballif@math.psu.edu The Pennsylvania State University January 18, 2008 Outline Rings Division Rings Noncommutative Rings s Roots of Rings Definition A ring R is a set toger with two binary operations +

More information

Linear Algebra. Paul Yiu. Department of Mathematics Florida Atlantic University. Fall 2011

Linear Algebra. Paul Yiu. Department of Mathematics Florida Atlantic University. Fall 2011 Linear Algebra Paul Yiu Department of Mathematics Florida Atlantic University Fall 2011 Linear Algebra Paul Yiu Department of Mathematics Florida Atlantic University Fall 2011 1A: Vector spaces Fields

More information

We are IntechOpen, the world s leading publisher of Open Access books Built by scientists, for scientists. International authors and editors

We are IntechOpen, the world s leading publisher of Open Access books Built by scientists, for scientists. International authors and editors We are IntechOpen, the world s leading publisher of Open Access books Built by scientists, for scientists 3,500 08,000.7 M Open access books available International authors and editors Downloads Our authors

More information

NONCOMMUTATIVE POLYNOMIAL EQUATIONS. Edward S. Letzter. Introduction

NONCOMMUTATIVE POLYNOMIAL EQUATIONS. Edward S. Letzter. Introduction NONCOMMUTATIVE POLYNOMIAL EQUATIONS Edward S Letzter Introduction My aim in these notes is twofold: First, to briefly review some linear algebra Second, to provide you with some new tools and techniques

More information

2. Time-Domain Analysis of Continuous- Time Signals and Systems

2. Time-Domain Analysis of Continuous- Time Signals and Systems 2. Time-Domain Analysis of Continuous- Time Signals and Systems 2.1. Continuous-Time Impulse Function (1.4.2) 2.2. Convolution Integral (2.2) 2.3. Continuous-Time Impulse Response (2.2) 2.4. Classification

More information

NOTES ON FINITE FIELDS

NOTES ON FINITE FIELDS NOTES ON FINITE FIELDS AARON LANDESMAN CONTENTS 1. Introduction to finite fields 2 2. Definition and constructions of fields 3 2.1. The definition of a field 3 2.2. Constructing field extensions by adjoining

More information

Linear System Theory

Linear System Theory Linear System Theory Wonhee Kim Lecture 3 Mar. 21, 2017 1 / 38 Overview Recap Nonlinear systems: existence and uniqueness of a solution of differential equations Preliminaries Fields and Vector Spaces

More information

1. Find the solution of the following uncontrolled linear system. 2 α 1 1

1. Find the solution of the following uncontrolled linear system. 2 α 1 1 Appendix B Revision Problems 1. Find the solution of the following uncontrolled linear system 0 1 1 ẋ = x, x(0) =. 2 3 1 Class test, August 1998 2. Given the linear system described by 2 α 1 1 ẋ = x +

More information

A Lie-Group Approach for Nonlinear Dynamic Systems Described by Implicit Ordinary Differential Equations

A Lie-Group Approach for Nonlinear Dynamic Systems Described by Implicit Ordinary Differential Equations A Lie-Group Approach for Nonlinear Dynamic Systems Described by Implicit Ordinary Differential Equations Kurt Schlacher, Andreas Kugi and Kurt Zehetleitner kurt.schlacher@jku.at kurt.zehetleitner@jku.at,

More information

Section III.6. Factorization in Polynomial Rings

Section III.6. Factorization in Polynomial Rings III.6. Factorization in Polynomial Rings 1 Section III.6. Factorization in Polynomial Rings Note. We push several of the results in Section III.3 (such as divisibility, irreducibility, and unique factorization)

More information

Signals & Systems. Lecture 5 Continuous-Time Fourier Transform. Alp Ertürk

Signals & Systems. Lecture 5 Continuous-Time Fourier Transform. Alp Ertürk Signals & Systems Lecture 5 Continuous-Time Fourier Transform Alp Ertürk alp.erturk@kocaeli.edu.tr Fourier Series Representation of Continuous-Time Periodic Signals Synthesis equation: x t = a k e jkω

More information

NATIONAL BOARD FOR HIGHER MATHEMATICS. Research Awards Screening Test. February 25, Time Allowed: 90 Minutes Maximum Marks: 40

NATIONAL BOARD FOR HIGHER MATHEMATICS. Research Awards Screening Test. February 25, Time Allowed: 90 Minutes Maximum Marks: 40 NATIONAL BOARD FOR HIGHER MATHEMATICS Research Awards Screening Test February 25, 2006 Time Allowed: 90 Minutes Maximum Marks: 40 Please read, carefully, the instructions on the following page before you

More information

Let X be a topological space. We want it to look locally like C. So we make the following definition.

Let X be a topological space. We want it to look locally like C. So we make the following definition. February 17, 2010 1 Riemann surfaces 1.1 Definitions and examples Let X be a topological space. We want it to look locally like C. So we make the following definition. Definition 1. A complex chart on

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

(3) Let Y be a totally bounded subset of a metric space X. Then the closure Y of Y

(3) Let Y be a totally bounded subset of a metric space X. Then the closure Y of Y () Consider A = { q Q : q 2 2} as a subset of the metric space (Q, d), where d(x, y) = x y. Then A is A) closed but not open in Q B) open but not closed in Q C) neither open nor closed in Q D) both open

More information

Problem 4 (Wed Jan 29) Let G be a finite abelian group. Prove that the following are equivalent

Problem 4 (Wed Jan 29) Let G be a finite abelian group. Prove that the following are equivalent Last revised: May 16, 2014 A.Miller M542 www.math.wisc.edu/ miller/ Problem 1 (Fri Jan 24) (a) Find an integer x such that x = 6 mod 10 and x = 15 mod 21 and 0 x 210. (b) Find the smallest positive integer

More information

Group Theory. Problem Set 3, Solution

Group Theory. Problem Set 3, Solution Problem Set 3, Solution Symmetry group of H From O b φ n = m T mnbφ m and O a φ q = p T pqaφ p we have O b O a φ q = O b T pq aφ p = T pq a T mp bφ m p p m =! T mp bt pq a φ m = O ba φ q m p = m T mq baφ

More information

University of Ottawa

University of Ottawa University of Ottawa Department of Mathematics and Statistics MAT3143: Ring Theory Professor: Hadi Salmasian Final Exam April 21, 2015 Surname First Name Instructions: (a) You have 3 hours to complete

More information

Homogeneous Linear Systems of Differential Equations with Constant Coefficients

Homogeneous Linear Systems of Differential Equations with Constant Coefficients Objective: Solve Homogeneous Linear Systems of Differential Equations with Constant Coefficients dx a x + a 2 x 2 + + a n x n, dx 2 a 2x + a 22 x 2 + + a 2n x n,. dx n = a n x + a n2 x 2 + + a nn x n.

More information

MATH32062 Notes. 1 Affine algebraic varieties. 1.1 Definition of affine algebraic varieties

MATH32062 Notes. 1 Affine algebraic varieties. 1.1 Definition of affine algebraic varieties MATH32062 Notes 1 Affine algebraic varieties 1.1 Definition of affine algebraic varieties We want to define an algebraic variety as the solution set of a collection of polynomial equations, or equivalently,

More information

Presentation of Normal Bases

Presentation of Normal Bases Presentation of Normal Bases Mohamadou Sall mohamadou1.sall@ucad.edu.sn University Cheikh Anta Diop, Dakar (Senegal) Pole of Research in Mathematics and their Applications in Information Security (PRMAIS)

More information

Modern Control Systems

Modern Control Systems Modern Control Systems Matthew M. Peet Illinois Institute of Technology Lecture 18: Linear Causal Time-Invariant Operators Operators L 2 and ˆL 2 space Because L 2 (, ) and ˆL 2 are isomorphic, so are

More information

Abstract Vector Spaces

Abstract Vector Spaces CHAPTER 1 Abstract Vector Spaces 1.1 Vector Spaces Let K be a field, i.e. a number system where you can add, subtract, multiply and divide. In this course we will take K to be R, C or Q. Definition 1.1.

More information

Prime and irreducible elements of the ring of integers modulo n

Prime and irreducible elements of the ring of integers modulo n Prime and irreducible elements of the ring of integers modulo n M. H. Jafari and A. R. Madadi Department of Pure Mathematics, Faculty of Mathematical Sciences University of Tabriz, Tabriz, Iran Abstract

More information

A Little Beyond: Linear Algebra

A Little Beyond: Linear Algebra A Little Beyond: Linear Algebra Akshay Tiwary March 6, 2016 Any suggestions, questions and remarks are welcome! 1 A little extra Linear Algebra 1. Show that any set of non-zero polynomials in [x], no two

More information

Lecture 7.5: Euclidean domains and algebraic integers

Lecture 7.5: Euclidean domains and algebraic integers Lecture 7.5: Euclidean domains and algebraic integers Matthew Macauley Department of Mathematical Sciences Clemson University http://www.math.clemson.edu/~macaule/ Math 4120, Modern Algebra M. Macauley

More information

Power series solutions for 2nd order linear ODE s (not necessarily with constant coefficients) a n z n. n=0

Power series solutions for 2nd order linear ODE s (not necessarily with constant coefficients) a n z n. n=0 Lecture 22 Power series solutions for 2nd order linear ODE s (not necessarily with constant coefficients) Recall a few facts about power series: a n z n This series in z is centered at z 0. Here z can

More information

10 Transfer Matrix Models

10 Transfer Matrix Models MIT EECS 6.241 (FALL 26) LECTURE NOTES BY A. MEGRETSKI 1 Transfer Matrix Models So far, transfer matrices were introduced for finite order state space LTI models, in which case they serve as an important

More information

ECE504: Lecture 10. D. Richard Brown III. Worcester Polytechnic Institute. 11-Nov-2008

ECE504: Lecture 10. D. Richard Brown III. Worcester Polytechnic Institute. 11-Nov-2008 ECE504: Lecture 10 D. Richard Brown III Worcester Polytechnic Institute 11-Nov-2008 Worcester Polytechnic Institute D. Richard Brown III 11-Nov-2008 1 / 25 Lecture 10 Major Topics We are finishing up Part

More information

Rings. EE 387, Notes 7, Handout #10

Rings. EE 387, Notes 7, Handout #10 Rings EE 387, Notes 7, Handout #10 Definition: A ring is a set R with binary operations, + and, that satisfy the following axioms: 1. (R, +) is a commutative group (five axioms) 2. Associative law for

More information

1. If 1, ω, ω 2, -----, ω 9 are the 10 th roots of unity, then (1 + ω) (1 + ω 2 ) (1 + ω 9 ) is A) 1 B) 1 C) 10 D) 0

1. If 1, ω, ω 2, -----, ω 9 are the 10 th roots of unity, then (1 + ω) (1 + ω 2 ) (1 + ω 9 ) is A) 1 B) 1 C) 10 D) 0 4 INUTES. If, ω, ω, -----, ω 9 are the th roots of unity, then ( + ω) ( + ω ) ----- ( + ω 9 ) is B) D) 5. i If - i = a + ib, then a =, b = B) a =, b = a =, b = D) a =, b= 3. Find the integral values for

More information

Qualifying Examination HARVARD UNIVERSITY Department of Mathematics Tuesday, January 19, 2016 (Day 1)

Qualifying Examination HARVARD UNIVERSITY Department of Mathematics Tuesday, January 19, 2016 (Day 1) Qualifying Examination HARVARD UNIVERSITY Department of Mathematics Tuesday, January 19, 2016 (Day 1) PROBLEM 1 (DG) Let S denote the surface in R 3 where the coordinates (x, y, z) obey x 2 + y 2 = 1 +

More information

Complex Dynamic Systems: Qualitative vs Quantitative analysis

Complex Dynamic Systems: Qualitative vs Quantitative analysis Complex Dynamic Systems: Qualitative vs Quantitative analysis Complex Dynamic Systems Chiara Mocenni Department of Information Engineering and Mathematics University of Siena (mocenni@diism.unisi.it) Dynamic

More information

ON FLATNESS OF NONLINEAR IMPLICIT SYSTEMS

ON FLATNESS OF NONLINEAR IMPLICIT SYSTEMS ON FLATNESS OF NONLINEAR IMPLICIT SYSTEMS Paulo Sergio Pereira da Silva, Simone Batista Escola Politécnica da USP Av. Luciano Gualerto trav. 03, 158 05508-900 Cidade Universitária São Paulo SP BRAZIL Escola

More information

1. Group Theory Permutations.

1. Group Theory Permutations. 1.1. Permutations. 1. Group Theory Problem 1.1. Let G be a subgroup of S n of index 2. Show that G = A n. Problem 1.2. Find two elements of S 7 that have the same order but are not conjugate. Let π S 7

More information

Algebra Review 2. 1 Fields. A field is an extension of the concept of a group.

Algebra Review 2. 1 Fields. A field is an extension of the concept of a group. Algebra Review 2 1 Fields A field is an extension of the concept of a group. Definition 1. A field (F, +,, 0 F, 1 F ) is a set F together with two binary operations (+, ) on F such that the following conditions

More information

Control Systems. Internal Stability - LTI systems. L. Lanari

Control Systems. Internal Stability - LTI systems. L. Lanari Control Systems Internal Stability - LTI systems L. Lanari outline LTI systems: definitions conditions South stability criterion equilibrium points Nonlinear systems: equilibrium points examples stable

More information

Honors Algebra 4, MATH 371 Winter 2010 Assignment 4 Due Wednesday, February 17 at 08:35

Honors Algebra 4, MATH 371 Winter 2010 Assignment 4 Due Wednesday, February 17 at 08:35 Honors Algebra 4, MATH 371 Winter 2010 Assignment 4 Due Wednesday, February 17 at 08:35 1. Let R be a commutative ring with 1 0. (a) Prove that the nilradical of R is equal to the intersection of the prime

More information

Nov : Lecture 22: Differential Operators, Harmonic Oscillators

Nov : Lecture 22: Differential Operators, Harmonic Oscillators 14 Nov. 3 005: Lecture : Differential Operators, Harmonic Oscillators Reading: Kreyszig Sections:.4 (pp:81 83),.5 (pp:83 89),.8 (pp:101 03) Differential Operators The idea of a function as something that

More information

Discussion on: Measurable signal decoupling with dynamic feedforward compensation and unknown-input observation for systems with direct feedthrough

Discussion on: Measurable signal decoupling with dynamic feedforward compensation and unknown-input observation for systems with direct feedthrough Discussion on: Measurable signal decoupling with dynamic feedforward compensation and unknown-input observation for systems with direct feedthrough H.L. Trentelman 1 The geometric approach In the last

More information

Page Points Possible Points. Total 200

Page Points Possible Points. Total 200 Instructions: 1. The point value of each exercise occurs adjacent to the problem. 2. No books or notes or calculators are allowed. Page Points Possible Points 2 20 3 20 4 18 5 18 6 24 7 18 8 24 9 20 10

More information

3.4. ZEROS OF POLYNOMIAL FUNCTIONS

3.4. ZEROS OF POLYNOMIAL FUNCTIONS 3.4. ZEROS OF POLYNOMIAL FUNCTIONS What You Should Learn Use the Fundamental Theorem of Algebra to determine the number of zeros of polynomial functions. Find rational zeros of polynomial functions. Find

More information

ASSIGNMENT - 1, DEC M.Sc. (FINAL) SECOND YEAR DEGREE MATHEMATICS. Maximum : 20 MARKS Answer ALL questions. is also a topology on X.

ASSIGNMENT - 1, DEC M.Sc. (FINAL) SECOND YEAR DEGREE MATHEMATICS. Maximum : 20 MARKS Answer ALL questions. is also a topology on X. (DM 21) ASSIGNMENT - 1, DEC-2013. PAPER - I : TOPOLOGY AND FUNCTIONAL ANALYSIS Maimum : 20 MARKS 1. (a) Prove that every separable metric space is second countable. Define a topological space. If T 1 and

More information

COPRIME FACTORIZATIONS AND WELL-POSED LINEAR SYSTEMS

COPRIME FACTORIZATIONS AND WELL-POSED LINEAR SYSTEMS SIAM J. CONTROL OPTIM. c 1998 Society for Industrial and Applied Mathematics Vol. 36, No. 4, pp. 1268 1292, July 1998 009 COPRIME FACTORIZATIONS AND WELL-POSED LINEAR SYSTEMS OLOF J. STAFFANS Abstract.

More information

ON TWO DEFINITIONS OF OBSERVATION SPACES 1

ON TWO DEFINITIONS OF OBSERVATION SPACES 1 ON TWO DEFINITIONS OF OBSERVATION SPACES Yuan Wang Eduardo D. Sontag Department of Mathematics Rutgers University, New Brunswick, NJ 08903 E-mail: sycon@fermat.rutgers.edu ABSTRACT This paper establishes

More information

Flatness conditions for control systems

Flatness conditions for control systems The Diffiety Institute Preprint Series Preprint DIPS 01 January 24, 2002 Flatness conditions for control systems by V N CHETVERIKOV Available via INTERNET: http://diffietyacru/ The Diffiety Institute Polevaya

More information

5 Set Operations, Functions, and Counting

5 Set Operations, Functions, and Counting 5 Set Operations, Functions, and Counting Let N denote the positive integers, N 0 := N {0} be the non-negative integers and Z = N 0 ( N) the positive and negative integers including 0, Q the rational numbers,

More information

A differential geometric setting for dynamic equivalence and dynamic linearization J.-B. Pomet 3

A differential geometric setting for dynamic equivalence and dynamic linearization J.-B. Pomet 3 A differential geometric setting for dynamic equivalence and dynamic linearization J-B Pomet 3 Infinitesimal Brunovsý form for nonlinear systems, with applications to dynamic linearization E Aranda-Bricaire

More information

ORDINARY DIFFERENTIAL EQUATIONS

ORDINARY DIFFERENTIAL EQUATIONS ORDINARY DIFFERENTIAL EQUATIONS GABRIEL NAGY Mathematics Department, Michigan State University, East Lansing, MI, 4884 NOVEMBER 9, 7 Summary This is an introduction to ordinary differential equations We

More information

On the Universal Enveloping Algebra: Including the Poincaré-Birkhoff-Witt Theorem

On the Universal Enveloping Algebra: Including the Poincaré-Birkhoff-Witt Theorem On the Universal Enveloping Algebra: Including the Poincaré-Birkhoff-Witt Theorem Tessa B. McMullen Ethan A. Smith December 2013 1 Contents 1 Universal Enveloping Algebra 4 1.1 Construction of the Universal

More information

Special Conformal Invariance

Special Conformal Invariance Chapter 6 Special Conformal Invariance Conformal transformation on the d-dimensional flat space-time manifold M is an invertible mapping of the space-time coordinate x x x the metric tensor invariant up

More information

EC Control Engineering Quiz II IIT Madras

EC Control Engineering Quiz II IIT Madras EC34 - Control Engineering Quiz II IIT Madras Linear algebra Find the eigenvalues and eigenvectors of A, A, A and A + 4I Find the eigenvalues and eigenvectors of the following matrices: (a) cos θ sin θ

More information

Monotone Control System. Brad C. Yu SEACS, National ICT Australia And RSISE, The Australian National University June, 2005

Monotone Control System. Brad C. Yu SEACS, National ICT Australia And RSISE, The Australian National University June, 2005 Brad C. Yu SEACS, National ICT Australia And RSISE, The Australian National University June, 005 Foreword The aim of this presentation is to give a (primitive) overview of monotone systems and monotone

More information

Feedback Linearization Lectures delivered at IIT-Kanpur, TEQIP program, September 2016.

Feedback Linearization Lectures delivered at IIT-Kanpur, TEQIP program, September 2016. Feedback Linearization Lectures delivered at IIT-Kanpur, TEQIP program, September 216 Ravi N Banavar banavar@iitbacin September 24, 216 These notes are based on my readings o the two books Nonlinear Control

More information

arxiv: v1 [nlin.si] 27 Aug 2012

arxiv: v1 [nlin.si] 27 Aug 2012 Characteristic Lie rings and symmetries of differential Painlevé I and Painlevé III equations arxiv:208.530v [nlin.si] 27 Aug 202 O. S. Kostrigina Ufa State Aviation Technical University, 2, K. Marx str.,

More information

Formal power series rings, inverse limits, and I-adic completions of rings

Formal power series rings, inverse limits, and I-adic completions of rings Formal power series rings, inverse limits, and I-adic completions of rings Formal semigroup rings and formal power series rings We next want to explore the notion of a (formal) power series ring in finitely

More information

Math 581 Problem Set 3 Solutions

Math 581 Problem Set 3 Solutions Math 581 Problem Set 3 Solutions 1. Prove that complex conjugation is a isomorphism from C to C. Proof: First we prove that it is a homomorphism. Define : C C by (z) = z. Note that (1) = 1. The other properties

More information

ON IDENTIFIABILITY OF NONLINEAR ODE MODELS AND APPLICATIONS IN VIRAL DYNAMICS

ON IDENTIFIABILITY OF NONLINEAR ODE MODELS AND APPLICATIONS IN VIRAL DYNAMICS ON IDENTIFIABILITY OF NONLINEAR ODE MODELS AND APPLICATIONS IN VIRAL DYNAMICS HONGYU MIAO, XIAOHUA XIA, ALAN S. PERELSON, AND HULIN WU Abstract. Ordinary differential equations (ODE) are a powerful tool

More information

ALGEBRA QUALIFYING EXAM PROBLEMS LINEAR ALGEBRA

ALGEBRA QUALIFYING EXAM PROBLEMS LINEAR ALGEBRA ALGEBRA QUALIFYING EXAM PROBLEMS LINEAR ALGEBRA Kent State University Department of Mathematical Sciences Compiled and Maintained by Donald L. White Version: August 29, 2017 CONTENTS LINEAR ALGEBRA AND

More information

The flat trefoil and other oddities

The flat trefoil and other oddities The flat trefoil and other oddities Joel Langer Case Western Reserve University ICERM June, 2015 Plane curves with compact polyhedral geometry Thm Assume: C CP 2 is an irreducible, real algebraic curve

More information

Definition 1. A set V is a vector space over the scalar field F {R, C} iff. there are two operations defined on V, called vector addition

Definition 1. A set V is a vector space over the scalar field F {R, C} iff. there are two operations defined on V, called vector addition 6 Vector Spaces with Inned Product Basis and Dimension Section Objective(s): Vector Spaces and Subspaces Linear (In)dependence Basis and Dimension Inner Product 6 Vector Spaces and Subspaces Definition

More information

MATH 54 - TOPOLOGY SUMMER 2015 FINAL EXAMINATION. Problem 1

MATH 54 - TOPOLOGY SUMMER 2015 FINAL EXAMINATION. Problem 1 MATH 54 - TOPOLOGY SUMMER 2015 FINAL EXAMINATION ELEMENTS OF SOLUTION Problem 1 1. Let X be a Hausdorff space and K 1, K 2 disjoint compact subsets of X. Prove that there exist disjoint open sets U 1 and

More information

The Generalized Laplace Transform: Applications to Adaptive Control*

The Generalized Laplace Transform: Applications to Adaptive Control* The Transform: Applications to Adaptive * J.M. Davis 1, I.A. Gravagne 2, B.J. Jackson 1, R.J. Marks II 2, A.A. Ramos 1 1 Department of Mathematics 2 Department of Electrical Engineering Baylor University

More information

Quantum Mechanics for Mathematicians: Energy, Momentum, and the Quantum Free Particle

Quantum Mechanics for Mathematicians: Energy, Momentum, and the Quantum Free Particle Quantum Mechanics for Mathematicians: Energy, Momentum, and the Quantum Free Particle Peter Woit Department of Mathematics, Columbia University woit@math.columbia.edu November 28, 2012 We ll now turn to

More information

The goal of this chapter is to study linear systems of ordinary differential equations: dt,..., dx ) T

The goal of this chapter is to study linear systems of ordinary differential equations: dt,..., dx ) T 1 1 Linear Systems The goal of this chapter is to study linear systems of ordinary differential equations: ẋ = Ax, x(0) = x 0, (1) where x R n, A is an n n matrix and ẋ = dx ( dt = dx1 dt,..., dx ) T n.

More information

. Consider the linear system dx= =! = " a b # x y! : (a) For what values of a and b do solutions oscillate (i.e., do both x(t) and y(t) pass through z

. Consider the linear system dx= =! =  a b # x y! : (a) For what values of a and b do solutions oscillate (i.e., do both x(t) and y(t) pass through z Preliminary Exam { 1999 Morning Part Instructions: No calculators or crib sheets are allowed. Do as many problems as you can. Justify your answers as much as you can but very briey. 1. For positive real

More information

TEST CODE: PMB SYLLABUS

TEST CODE: PMB SYLLABUS TEST CODE: PMB SYLLABUS Convergence and divergence of sequence and series; Cauchy sequence and completeness; Bolzano-Weierstrass theorem; continuity, uniform continuity, differentiability; directional

More information

ALGEBRA HW 3 CLAY SHONKWILER

ALGEBRA HW 3 CLAY SHONKWILER ALGEBRA HW 3 CLAY SHONKWILER (a): Show that R[x] is a flat R-module. 1 Proof. Consider the set A = {1, x, x 2,...}. Then certainly A generates R[x] as an R-module. Suppose there is some finite linear combination

More information

Cheng Soon Ong & Christian Walder. Canberra February June 2018

Cheng Soon Ong & Christian Walder. Canberra February June 2018 Cheng Soon Ong & Christian Walder Research Group and College of Engineering and Computer Science Canberra February June 2018 (Many figures from C. M. Bishop, "Pattern Recognition and ") 1of 89 Part II

More information