---_._-_ _._._-... _ (56-r1S:_~0.~~\<;).e~~~b::;.~ 11<3<3 G

Similar documents
Review: control, feedback, etc. Today s topic: state-space models of systems; linearization

A review of stability and dynamical behaviors of differential equations:

cc) c, >o,, c, =c7 $"(x\ <o - S ~ Z

Time Response of Dynamic Systems! Multi-Dimensional Trajectories Position, velocity, and acceleration are vectors

TRAN S F O R M E R S TRA N SMI S S I O N. SECTION AB Issue 2, March, March,1958, by American Telephone and Telegraph Company

Lecture Note 1: Background

Chapter #4 EEE8086-EEE8115. Robust and Adaptive Control Systems

June 2011 PURDUE UNIVERSITY Study Guide for the Credit Exam in (MA 262) Linear Algebra and Differential Equations

Jsee x dx = In Isec x + tanxl + C Jcsc x dx = - In I cscx + cotxl + C

Linearization problem. The simplest example

25. Chain Rule. Now, f is a function of t only. Expand by multiplication:

Proc. of the 23rd Intl. Conf. on Parallel Processing, St. Charles, Illinois, August 1994, vol. 3, pp. 227{ Hanan Samet

Laplace Transforms Chapter 3

Chapter 1 Fundamental Concepts

Topic # Feedback Control

Lecture 1 From Continuous-Time to Discrete-Time

A MODEL FOR THE POPULATION DYNAMICS OF INTERNAL PARASITES. Keewhan Choi ABSTRACT

PH.D. PRELIMINARY EXAMINATION MATHEMATICS

Topic # /31 Feedback Control Systems. Analysis of Nonlinear Systems Lyapunov Stability Analysis

Problem Set 1. This week. Please read all of Chapter 1 in the Strauss text.

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

11 a 12 a 21 a 11 a 22 a 12 a 21. (C.11) A = The determinant of a product of two matrices is given by AB = A B 1 1 = (C.13) and similarly.

11.2 Basic First-order System Methods

Chapter 3. f(x) == [. k 1 =? - == k

Chapter 2: Time-Domain Representations of Linear Time-Invariant Systems. Chih-Wei Liu

PHYS 410/555 Computational Physics Solution of Non Linear Equations (a.k.a. Root Finding) (Reference Numerical Recipes, 9.0, 9.1, 9.

Chapter 1 Fundamental Concepts

Linear Algebra. Session 8

Lecture 8 - IIR Filters (II)

T h e C S E T I P r o j e c t

The Derivative. Appendix B. B.1 The Derivative of f. Mappings from IR to IR

176 5 t h Fl oo r. 337 P o ly me r Ma te ri al s

Linear Algebra Review (Course Notes for Math 308H - Spring 2016)

Module 02 CPS Background: Linear Systems Preliminaries

Control Systems

f (r) (a) r! (x a) r, r=0

Nonlinear Control Lecture 9: Feedback Linearization

P a g e 5 1 of R e p o r t P B 4 / 0 9

fur \ \,,^N/ D7,,)d.s) 7. The champion and Runner up of the previous year shall be allowed to play directly in final Zone.

Modeling nonlinear systems using multiple piecewise linear equations

ECE557 Systems Control

Multiple Choice Answers. MA 114 Calculus II Spring 2013 Final Exam 1 May Question

Classification of Phase Portraits at Equilibria for u (t) = f( u(t))

Math 212-Lecture 8. The chain rule with one independent variable

Section 9.7 and 9.10: Taylor Polynomials and Approximations/Taylor and Maclaurin Series

Review of Strain Energy Methods and Introduction to Stiffness Matrix Methods of Structural Analysis

Lecture 2: Controllability of nonlinear systems

EE 637 Final April 30, Spring Each problem is worth 20 points for a total score of 100 points

Georgia Institute of Technology Nonlinear Controls Theory Primer ME 6402

ECE 314 Signals and Systems Fall 2012

OH BOY! Story. N a r r a t iv e a n d o bj e c t s th ea t e r Fo r a l l a g e s, fr o m th e a ge of 9

Section 5.8. Taylor Series

1. Type your solutions. This homework is mainly a programming assignment.

DO NOT DO HOMEWORK UNTIL IT IS ASSIGNED. THE ASSIGNMENTS MAY CHANGE UNTIL ANNOUNCED.

ECE7850 Lecture 7. Discrete Time Optimal Control and Dynamic Programming

Course Summary Math 211

Translational and Rotational Dynamics!

Least-squares data fitting

Control Systems! Copyright 2017 by Robert Stengel. All rights reserved. For educational use only.

1. Linearity of a Function A function f(x) is defined linear if. f(αx 1 + βx 2 ) = αf(x 1 ) + βf(x 2 )

5. Hand in the entire exam booklet and your computer score sheet.

l(- oo)-~j [-I <. )( L6\ \ -J ~ ~ ~~~ ~~L{ ,~:::-=r\ or L":: -j) {fevylemr.eor k, ("p J~ -4" e S ' e,~ :; ij or J iv I 0"'& ~~ a. 11 qa.

Suppose that we have a specific single stage dynamic system governed by the following equation:

NARAYANA I I T / P M T A C A D E M Y. C o m m o n P r a c t i c e T e s t 05 XII-IC SPARK Date:

AM 205: lecture 14. Last time: Boundary value problems Today: Numerical solution of PDEs

Laplace Transforms. Chapter 3. Pierre Simon Laplace Born: 23 March 1749 in Beaumont-en-Auge, Normandy, France Died: 5 March 1827 in Paris, France

Math 106: Review for Exam II - SOLUTIONS

Math 216 Final Exam 24 April, 2017

g(t) = f(x 1 (t),..., x n (t)).

Trade Patterns, Production networks, and Trade and employment in the Asia-US region

Ritz Method. Introductory Course on Multiphysics Modelling

Taylor series - Solutions

Lecture two. January 17, 2019

Math 106 Fall 2014 Exam 2.1 October 31, ln(x) x 3 dx = 1. 2 x 2 ln(x) + = 1 2 x 2 ln(x) + 1. = 1 2 x 2 ln(x) 1 4 x 2 + C

Jim Lambers MAT 285 Spring Semester Practice Exam 2 Solution. y(t) = 5 2 e t 1 2 e 3t.

Perspective. ECE 3640 Lecture 11 State-Space Analysis. To learn about state-space analysis for continuous and discrete-time. Objective: systems

Partial differential equations (ACM30220)

A L A BA M A L A W R E V IE W

Math review. Math review

ECE 350 Signals and Systems Spring 2011 Final Exam - Solutions. Three 8 ½ x 11 sheets of notes, and a calculator are allowed during the exam.

Differential Equation Types. Moritz Diehl

Quasi-linear first order equations. Consider the nonlinear transport equation

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

SOME NEW FINITE ELEMENT METHODS FOR THE ANALYSIS OF SHOCK AND ACCELERATION WAVES IN NONLINEAR MATERIALS. J. T. ODEN The University of Texas at Austin

Differential Equations 2280 Sample Midterm Exam 3 with Solutions Exam Date: 24 April 2015 at 12:50pm

i.ea IE !e e sv?f 'il i+x3p \r= v * 5,?: S i- co i, ==:= SOrq) Xgs'iY # oo .9 9 PE * v E=S s->'d =ar4lq 5,n =.9 '{nl a':1 t F #l *r C\ t-e

MATH4406 (Control Theory) Unit 6: The Linear Quadratic Regulator (LQR) and Model Predictive Control (MPC) Prepared by Yoni Nazarathy, Artem

Intro. Computer Control Systems: F8

MAHALAKSHMI ENGINEERING COLLEGE-TRICHY

MATH 23 Exam 2 Review Solutions

Exact and high order discretization schemes. Wishart processes and their affine extensions

Lecture Discrete dynamic systems

10 Transfer Matrix Models

Control for Coordination of Linear Systems

Lecture 2: Review of Prerequisites. Table of contents

Optimal Control. Quadratic Functions. Single variable quadratic function: Multi-variable quadratic function:

ON THE AVERAGE NUMBER OF REAL ROOTS OF A RANDOM ALGEBRAIC EQUATION

Advanced Quantum Mechanics

Solution of Linear State-space Systems

Solution: APPM 1350 Final Exam Spring 2014

Transcription:

~-~-- ---_._-_.------_._._-... _--- -.-.. -_. -- --- ----,---------- _......1-:( ~~,A ~~AT1~ H.Qr~!j.'-~ f'-\ ;. A-P- S':1S'"[E-r.-t.>_ L~'lt1I\. (56-r1S:_~0.~~\<;).e~~~b::;.~ 11<3<3 G Consider now the general nonlinear dynamic control system in matrix form represented by d dtx(t) = F(x(t), u(t)) 0.66) where x(t), u(t), and F are, respectively, the n-dimensional state space vector, the r-dimensional control vector, and the n-dimensional vector function. Assume that the nominal (operating) system trajectory Xn (t) is known and that the nominal control that keeps the system on the nominal trajectory is given by un(t). Using the same logic as for the scalar case, we can assume that the actual system dynamics in the immediate proximity of the system nominal trajectories can be approximated by the first terms of the Taylor series. That is, starting with x(t) = xn(t) +.6.x(t), u(t) = un(t) +.6.u(t) (1.67) and (1.68) we expand equation (1.66) as follows d d dtxn + dt.6.x = F(xn +.6.x, Un +.6.u) = F(xn' Un) + (~F)..6.x + (~F)..6.u + high-order terms ux I Xn{t) uu "",,,(t) unit) I un(') (1.69) High-order terms contain at least quadratic quantities of.6.x and.6.u. Since.6.x and.6.u are small their squares are even smaller, and hence the high-order terms can be neglected. Using (1.67) and neglecting high-order terms, an approximation is obtained d (OF) (OF) dt.6.x(t)= -;r.6.x(t) + ~.6.u(t) (1.70) ox IXn(t) uu Ix"{') u.. (t) unit)

@ INTRODUCTION 27 Partial derivatives in (1.70) represent the Jacobian matrices given by (OF) ox IXn(t) un(t) _Anxn_ - - 8F, ~ 8F, 8x, 8X2 8xn li 8F2 8x, ~ 8Fi 8xj 8Fn 8Fn 8Fn 8x, 8x, 8xn I Xn(t) Un (t) (l.71a) 8F, 8F, &F, &u, e,;:; 8ur &F2 li &u, BUr BFi (~~) IXn(t) Bu] un(t) =B nxt = (1.71b) BFn BFn BFn I Bu, Xn(t) 8U2 BUr Un (t) Note that the Jacobian matrices have to be evaluated at the nominal points, i.e. at xn(t) and un(t). With this notation, the linearized system (1.70) has the form d dt 6.x(t) = A6.x(t) +B6.u(t), 6.x(to) = x(to) - xn(to) (1.72) The output of a nonlinear system, in general, satisfies a nonlinear algebraic equation, that is y(t) = 9(x(t), u(t)) (1. 73) This equation can be also linearized by expanding its right-hand side into a Taylor series about nominal points Xn ( t) and Un(t). This leads to + high-order term,? (1.74) Note that Yn cancels term 9(xn,yn). By neglecting high-order terms in (1.74), the linearized part of the output equation is given by 6.y(t) = Cax(t) + nau(t) (1.75)

28 INTRODUCTION where the Jacobian matrices C and D satisfy ~ ~ f!flj.. ax, ax2 ax" ag2 ag2 ax, ax" cpxn _ ag, (09) axj - ox Ix,,(t) (1.76a) u,,(t) ag p ag p ag p I x,,(t) ax, ax2 ax" u,,(t) (09) ~ ag~ au, ~ au, DPxr _ ~ - OU auj IX"(') un(t).q... OUr ag2 OUr!2.r.!2.r.!2.r. I x,,(t) au, au2 OUr u" (t) (1.76b) Example 1.2: Let a nonlinear system be represented by dxi. at = Xl SlllX2 + X2 U dx2 2 dt =:: Xle- x2 + u y = 2XIX2 + x~ Assume that the values for the system nominal trajectories and control are known and given by Xl n, X2n, and Un. The linearized state space equation of the above nonlinear system is obtained as [ L\~l(t)] = [ sinx 2n XlncoSX2n+un] [L\XI(t)] + [X2n]L\u(t) L\x2(t) e-x2n -xlne-x2n L\X2(t) 2un L\y(t) = [2X2n 2Xln + 2X2n 1[~:~~~~] + O~u(t) Having obtained the solution of this linearized system under the given control input 6.u( t), the corresponding approximation'of the nonlinear system trajectories is <>

( INTRODUCTION 29 Example 1.3: Consider the mathematical model of a single-link robotic manipulator with a flexible joint (Spong and Vidyasagar, 1989) IBI + mgl sinbi +k(bi - ( 2 ) = 0 JB2 - k(bi - ( 2 ) = U where BI, ()2 are angular positions, I, J are moments of inertia, m and l are, respectively, the link's mass and length, and k is the link's spring constant. Introducing the change of variables as the manipulator's state space nonlinear model equivalent to (1.66) is given by Xl = Xz. mgl. k ( ) X2 = ---smxi - - Xl - X3 I I X3 = X4. k ( ) 1 X4 = J Xl - X3 + JU Take the nominal points as (Xl n, Xzn, X3n, X4n, un), then the matrices A and B defined in (1.71) are given by _ k+m~l o A= [ I k J COSX' n 1 o o o In Spong (1995) the following numerical values are used for system parameters: mgl = 5, I = J = 1, k = 0.08. Assuming that the output variable is equal to the link's angular position, that is the matrices C and D, defined in (1.76), are given by C = [1 0 0 0], D = 0

5-1 ct"s.-rz...e.r '7:pvCU-5 ~. <1 'Pr-oblevtA S+~ Jv\.l- Hcn')-e(.;I)i)~ PLa'\'1t- ~li"~ }-,-tec.s~e.--\'1/le.lf\ -t- 'ModeQs> (JoY\t-c'Yl'U--0U"') t-c5vv1. e ;;;.. ::::- f (,t... \:t; j-+- 6 'V-f? Z = ( C~ 1'=-) +- --I ') ~. C5 y-,ut- ')'\ece-s<:;~ t::h~ ~ ~'\dl r-y ~_ G~u-,.Y5CGU<..>(., vfh~.---.;h 1i..=c.><... =: BO""Oss~-) sc"od,,;g's~ ~C~~ --=P~' e) ex:.) $0 y:. C:-t-) U f \ 1, ~"_L"!"'<;~ ",:H~Gle '()C\L)-e = NOlA 6dt.5-S '5~ '?~stt.c ~~S.~ [t-) 11 1\ ~12... - f C,.c.)L-.J J lj-~..i) + '""Vfl.L-1 :;z:l.2- - Z (:~1i.. \ ),L ') +- "'"~ ~ ~ '':J-I}I-'; --=:. () -t=~ ~bl ~{'J -=- &~ A (;L-L) ~1 ""12-) = D q "'IJ1..vt) -=P~ AL~L) k{ )(.o::r =-~o~'v) C~ t ~:; =-.ilm.(.;>~' t.-1

@! 3),-:S- Ot4'\d. g, O-{re ~ Gee C.OCA k'll\o 1...9-. ':.~ elf:_~_o\ ::r-c..c!'~le, "--7 ~~2~~ f>ross---.1i)/.e (@ 1-IH{;-A~2.A=nc-.t'-f Or Dl'S~4'F-'1lf"'>'(E Ho.t{(...(NC~ S'"fSCE-Ks, (j... A-9601 A- H6}-1\ I H: Ae-- 'TP-A--~~_e::-X '~'-~~--~~,--':::::"-q ') ')(~)J-. =-- -+ c -" "J,L..~? Jl.~) + ~ JL- ~ / ~JL.N6M= ~()<.:~r' :;12--4-) ~pot-'l=-~~ 6i-?Gi/'\E' cf; '-H G~~... ~ r C'il:.V~ die;> U/\ '-fj~ l;2:::) ~_ te.;;a.uz- e 'J-f'rettlsecJ M6M b 'I<-\L = ~ -,zjl. S -;z J.L == X-)L - ~ ( -,i..jl.t-\o M 1.t.l-) ~'Us, ('2-) ~"- --~JL8?<l<--J- -t ~ L6 ~,ej...a,'l~~ m.9c:{~,~ o~, C~').t

i-i f-\ E-A-P-t"ZA-'1l6H or '(l-te DUTrO') ~(:) All on (?> ) f':i ll -= =~ ('>YL. I l.i.. )-';-~J? ~Wl 'Z~ '=- e, 't'-'io1-1_ 0 (+VI--"! ) <2)L- \/".. ~jl! (-J... R ( )ZjJ. I.\L) = ~ I C)'C:CV-\!.LJ -l- :dr::) ( """- +- ~I )4t.=-~OM 4 l:> c~) -K)L. 8 ~ st. 4-~ '" ~J -.. ~

;:. '=- +c~ \1:-"-) +-G v-tct-) ~ -=- ~ (~l t:-) + '"Vci-) -r- -=-.,.c-'y' -\- c;~ b;: - j ~~, \ ] ~"+-G"\HLi:-) =F ~~+ Gv-r - l 0)(..?"-_;>t...K\)~ Mso @:~ J-l'N~j2A'TlC\N ~OlJT~~~S,0T\r--t~_ TAA-~7 ~ 'Y'\01Jli\0Y"\a;..( ~~ <-b re-rlac..j?~ ~ tee e'stvl~~_ ~c.~.. ~ ;-e~t- 'LS ~ 5G\-'Y\e, :d ~ (~\ I.L") ''0-'", '~,t~~c) a~,,- J,,~"=-#+),~~&I~) \ ~ ~=-~(+-)