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

Size: px
Start display at page:

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

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)

2 @ 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)

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

4 ( 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 = 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 = [ ], D = 0

5 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\ ~ 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

6 @! 3),-:S- Ot4'\d. g, O-{re ~ Gee C.OCA k'll\o ':.~ 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

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

8 ;:. '=- +c~ \1:-"-) +-G v-tct-) ~ -=- ~ (~l t:-) + '"Vci-) -r- -=-.,.c-'y' -\- c;~ b;: - j ~~, \ ] ~"+-G"\HLi:-) =F ~~+ Gv-r - l 0)(..?"-_;>t...K\)~ 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~) \ ~ ~=-~(+-)

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

Review: control, feedback, etc. Today s topic: state-space models of systems; linearization Plan of the Lecture Review: control, feedback, etc Today s topic: state-space models of systems; linearization Goal: a general framework that encompasses all examples of interest Once we have mastered

More information

A review of stability and dynamical behaviors of differential equations:

A review of stability and dynamical behaviors of differential equations: A review of stability and dynamical behaviors of differential equations: scalar ODE: u t = f(u), system of ODEs: u t = f(u, v), v t = g(u, v), reaction-diffusion equation: u t = D u + f(u), x Ω, with boundary

More information

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

Time Response of Dynamic Systems! Multi-Dimensional Trajectories Position, velocity, and acceleration are vectors Time Response of Dynamic Systems Robert Stengel Robotics and Intelligent Systems MAE 345, Princeton University, 217 Multi-dimensional trajectories Numerical integration Linear and nonlinear systems Linearization

More information

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

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 B B, ch, 9 B h f h c h l f f Bll lh, c chcl l l k l, h, h ch f h f ll ll l f h lh h c ll k f Bll lh, c ck ll ch,9, c lh lh B B c x l l f f f 9 B l f c l f f l 9 f B c, h f c x ch 9 B B h f f fc Bll c f

More information

Lecture Note 1: Background

Lecture Note 1: Background ECE5463: Introduction to Robotics Lecture Note 1: Background Prof. Wei Zhang Department of Electrical and Computer Engineering Ohio State University Columbus, Ohio, USA Spring 2018 Lecture 1 (ECE5463 Sp18)

More information

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

Chapter #4 EEE8086-EEE8115. Robust and Adaptive Control Systems Chapter #4 Robust and Adaptive Control Systems Nonlinear Dynamics.... Linear Combination.... Equilibrium points... 3 3. Linearisation... 5 4. Limit cycles... 3 5. Bifurcations... 4 6. Stability... 6 7.

More information

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

June 2011 PURDUE UNIVERSITY Study Guide for the Credit Exam in (MA 262) Linear Algebra and Differential Equations June 20 PURDUE UNIVERSITY Study Guide for the Credit Exam in (MA 262) Linear Algebra and Differential Equations The topics covered in this exam can be found in An introduction to differential equations

More information

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

Jsee x dx = In Isec x + tanxl + C Jcsc x dx = - In I cscx + cotxl + C MAC 2312 Final Exam Review Instructions: The Final Exam will consist of 15 questions plus a bonus problem. All questions will be multiple choice, which will be graded partly on whether or not you circle

More information

Linearization problem. The simplest example

Linearization problem. The simplest example Linear Systems Lecture 3 1 problem Consider a non-linear time-invariant system of the form ( ẋ(t f x(t u(t y(t g ( x(t u(t (1 such that x R n u R m y R p and Slide 1 A: f(xu f(xu g(xu and g(xu exist and

More information

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

25. Chain Rule. Now, f is a function of t only. Expand by multiplication: 25. Chain Rule The Chain Rule is present in all differentiation. If z = f(x, y) represents a two-variable function, then it is plausible to consider the cases when x and y may be functions of other variable(s).

More information

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

Proc. of the 23rd Intl. Conf. on Parallel Processing, St. Charles, Illinois, August 1994, vol. 3, pp. 227{ Hanan Samet P. 23 Il. C. Plll P, S. Cl, Ill, 1994, vl. 3,. 227{234 1 DT-PRE SPTI JOI GORITHMS Ek G. Hl y Gy Dv B C W, D.C. 20233 H S C S D C R I v C S Uvy Myl Cll Pk, Myl 20742 { E -lll l j l R-, l,. T l l (.., B

More information

Laplace Transforms Chapter 3

Laplace Transforms Chapter 3 Laplace Transforms Important analytical method for solving linear ordinary differential equations. - Application to nonlinear ODEs? Must linearize first. Laplace transforms play a key role in important

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

Topic # Feedback Control

Topic # Feedback Control Topic #7 16.31 Feedback Control State-Space Systems What are state-space models? Why should we use them? How are they related to the transfer functions used in classical control design and how do we develop

More information

Lecture 1 From Continuous-Time to Discrete-Time

Lecture 1 From Continuous-Time to Discrete-Time Lecture From Continuous-Time to Discrete-Time Outline. Continuous and Discrete-Time Signals and Systems................. What is a signal?................................2 What is a system?.............................

More information

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

A MODEL FOR THE POPULATION DYNAMICS OF INTERNAL PARASITES. Keewhan Choi ABSTRACT .. A MODEL FOR THE POPULATION DYNAMICS OF INTERNAL PARASITES BU-176-M Keewhan Choi ABSTRACT We are interested in the loss of red cells in a parasitized sheep in relation to the contamination of posture

More information

PH.D. PRELIMINARY EXAMINATION MATHEMATICS

PH.D. PRELIMINARY EXAMINATION MATHEMATICS UNIVERSITY OF CALIFORNIA, BERKELEY Dept. of Civil and Environmental Engineering FALL SEMESTER 2014 Structural Engineering, Mechanics and Materials NAME PH.D. PRELIMINARY EXAMINATION MATHEMATICS Problem

More information

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

Topic # /31 Feedback Control Systems. Analysis of Nonlinear Systems Lyapunov Stability Analysis Topic # 16.30/31 Feedback Control Systems Analysis of Nonlinear Systems Lyapunov Stability Analysis Fall 010 16.30/31 Lyapunov Stability Analysis Very general method to prove (or disprove) stability of

More information

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

Problem Set 1. This week. Please read all of Chapter 1 in the Strauss text. Math 425, Spring 2015 Jerry L. Kazdan Problem Set 1 Due: Thurs. Jan. 22 in class. [Late papers will be accepted until 1:00 PM Friday.] This is rust remover. It is essentially Homework Set 0 with a few

More information

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

ACM/CMS 107 Linear Analysis & Applications Fall 2016 Assignment 4: Linear ODEs and Control Theory Due: 5th December 2016 ACM/CMS 17 Linear Analysis & Applications Fall 216 Assignment 4: Linear ODEs and Control Theory Due: 5th December 216 Introduction Systems of ordinary differential equations (ODEs) can be used to describe

More information

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 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. C PROPERTIES OF MATRICES 697 to whether the permutation i 1 i 2 i N is even or odd, respectively Note that I =1 Thus, for a 2 2 matrix, the determinant takes the form A = a 11 a 12 = a a 21 a 11 a 22 a

More information

11.2 Basic First-order System Methods

11.2 Basic First-order System Methods 112 Basic First-order System Methods 797 112 Basic First-order System Methods Solving 2 2 Systems It is shown here that any constant linear system u a b = A u, A = c d can be solved by one of the following

More information

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

Chapter 3. f(x) == [. k 1 =? - == k Chapter 3 3.1 (1) The term Ixl is not continuously differentiable at x == 0, but it is globally Lipschitz. The term x2 is continuously differentiable, but its partial derivative is not globally bounded.

More information

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

Chapter 2: Time-Domain Representations of Linear Time-Invariant Systems. Chih-Wei Liu Chapter : Time-Domain Representations of Linear Time-Invariant Systems Chih-Wei Liu Outline Characteristics of Systems Described by Differential and Difference Equations Block Diagram Representations State-Variable

More information

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

PHYS 410/555 Computational Physics Solution of Non Linear Equations (a.k.a. Root Finding) (Reference Numerical Recipes, 9.0, 9.1, 9. PHYS 410/555 Computational Physics Solution of Non Linear Equations (a.k.a. Root Finding) (Reference Numerical Recipes, 9.0, 9.1, 9.4) We will consider two cases 1. f(x) = 0 1-dimensional 2. f(x) = 0 d-dimensional

More information

Chapter 1 Fundamental Concepts

Chapter 1 Fundamental Concepts Chapter 1 Fundamental Concepts Signals A signal is a pattern of variation of a physical quantity as a function of time, space, distance, position, temperature, pressure, etc. These quantities are usually

More information

Linear Algebra. Session 8

Linear Algebra. Session 8 Linear Algebra. Session 8 Dr. Marco A Roque Sol 08/01/2017 Abstract Linear Algebra Range and kernel Let V, W be vector spaces and L : V W, be a linear mapping. Definition. The range (or image of L is the

More information

Lecture 8 - IIR Filters (II)

Lecture 8 - IIR Filters (II) Lecture 8 - IIR Filters (II) James Barnes (James.Barnes@colostate.edu) Spring 2009 Colorado State University Dept of Electrical and Computer Engineering ECE423 1 / 27 Lecture 8 Outline Introduction Digital

More information

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

T h e C S E T I P r o j e c t T h e P r o j e c t T H E P R O J E C T T A B L E O F C O N T E N T S A r t i c l e P a g e C o m p r e h e n s i v e A s s es s m e n t o f t h e U F O / E T I P h e n o m e n o n M a y 1 9 9 1 1 E T

More information

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

The Derivative. Appendix B. B.1 The Derivative of f. Mappings from IR to IR Appendix B The Derivative B.1 The Derivative of f In this chapter, we give a short summary of the derivative. Specifically, we want to compare/contrast how the derivative appears for functions whose domain

More information

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

176 5 t h Fl oo r. 337 P o ly me r Ma te ri al s A g la di ou s F. L. 462 E l ec tr on ic D ev el op me nt A i ng er A.W.S. 371 C. A. M. A l ex an de r 236 A d mi ni st ra ti on R. H. (M rs ) A n dr ew s P. V. 326 O p ti ca l Tr an sm is si on A p ps

More information

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

Linear Algebra Review (Course Notes for Math 308H - Spring 2016) Linear Algebra Review (Course Notes for Math 308H - Spring 2016) Dr. Michael S. Pilant February 12, 2016 1 Background: We begin with one of the most fundamental notions in R 2, distance. Letting (x 1,

More information

Module 02 CPS Background: Linear Systems Preliminaries

Module 02 CPS Background: Linear Systems Preliminaries Module 02 CPS Background: Linear Systems Preliminaries Ahmad F. Taha EE 5243: Introduction to Cyber-Physical Systems Email: ahmad.taha@utsa.edu Webpage: http://engineering.utsa.edu/ taha/index.html August

More information

Control Systems

Control Systems 6.5 Control Systems Last Time: Introdction Motivation Corse Overview Project Math. Descriptions of Systems ~ Review Classification of Systems Linear Systems LTI Systems The notion of state and state variables

More information

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

f (r) (a) r! (x a) r, r=0 Part 3.3 Differentiation v1 2018 Taylor Polynomials Definition 3.3.1 Taylor 1715 and Maclaurin 1742) If a is a fixed number, and f is a function whose first n derivatives exist at a then the Taylor polynomial

More information

Nonlinear Control Lecture 9: Feedback Linearization

Nonlinear Control Lecture 9: Feedback Linearization Nonlinear Control Lecture 9: Feedback Linearization Farzaneh Abdollahi Department of Electrical Engineering Amirkabir University of Technology Fall 2011 Farzaneh Abdollahi Nonlinear Control Lecture 9 1/75

More information

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

P a g e 5 1 of R e p o r t P B 4 / 0 9 P a g e 5 1 of R e p o r t P B 4 / 0 9 J A R T a l s o c o n c l u d e d t h a t a l t h o u g h t h e i n t e n t o f N e l s o n s r e h a b i l i t a t i o n p l a n i s t o e n h a n c e c o n n e

More information

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

fur \ \,,^N/ D7,,)d.s) 7. The champion and Runner up of the previous year shall be allowed to play directly in final Zone. OUL O GR SODRY DUTO, ODS,RT,SMTUR,USWR.l ntuctin f cnuct f Kbi ( y/gil)tunent f 2L-Lg t. 2.. 4.. 6. Mtche hll be lye e K ule f ene f tie t tie Dutin f ech tch hll be - +0 (Rece)+ = M The ticint f ech Te

More information

Modeling nonlinear systems using multiple piecewise linear equations

Modeling nonlinear systems using multiple piecewise linear equations Nonlinear Analysis: Modelling and Control, 2010, Vol. 15, No. 4, 451 458 Modeling nonlinear systems using multiple piecewise linear equations G.K. Lowe, M.A. Zohdy Department of Electrical and Computer

More information

ECE557 Systems Control

ECE557 Systems Control ECE557 Systems Control Bruce Francis Course notes, Version.0, September 008 Preface This is the second Engineering Science course on control. It assumes ECE56 as a prerequisite. If you didn t take ECE56,

More information

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

Multiple Choice Answers. MA 114 Calculus II Spring 2013 Final Exam 1 May Question MA 114 Calculus II Spring 2013 Final Exam 1 May 2013 Name: Section: Last 4 digits of student ID #: This exam has six multiple choice questions (six points each) and five free response questions with points

More information

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

Classification of Phase Portraits at Equilibria for u (t) = f( u(t)) Classification of Phase Portraits at Equilibria for u t = f ut Transfer of Local Linearized Phase Portrait Transfer of Local Linearized Stability How to Classify Linear Equilibria Justification of the

More information

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

Math 212-Lecture 8. The chain rule with one independent variable Math 212-Lecture 8 137: The multivariable chain rule The chain rule with one independent variable w = f(x, y) If the particle is moving along a curve x = x(t), y = y(t), then the values that the particle

More information

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

Section 9.7 and 9.10: Taylor Polynomials and Approximations/Taylor and Maclaurin Series Section 9.7 and 9.10: Taylor Polynomials and Approximations/Taylor and Maclaurin Series Power Series for Functions We can create a Power Series (or polynomial series) that can approximate a function around

More information

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

Review of Strain Energy Methods and Introduction to Stiffness Matrix Methods of Structural Analysis uke University epartment of Civil and Environmental Engineering CEE 42L. Matrix Structural Analysis Henri P. Gavin Fall, 22 Review of Strain Energy Methods and Introduction to Stiffness Matrix Methods

More information

Lecture 2: Controllability of nonlinear systems

Lecture 2: Controllability of nonlinear systems DISC Systems and Control Theory of Nonlinear Systems 1 Lecture 2: Controllability of nonlinear systems Nonlinear Dynamical Control Systems, Chapter 3 See www.math.rug.nl/ arjan (under teaching) for info

More information

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

EE 637 Final April 30, Spring Each problem is worth 20 points for a total score of 100 points EE 637 Final April 30, Spring 2018 Name: Instructions: This is a 120 minute exam containing five problems. Each problem is worth 20 points for a total score of 100 points You may only use your brain and

More information

Georgia Institute of Technology Nonlinear Controls Theory Primer ME 6402

Georgia Institute of Technology Nonlinear Controls Theory Primer ME 6402 Georgia Institute of Technology Nonlinear Controls Theory Primer ME 640 Ajeya Karajgikar April 6, 011 Definition Stability (Lyapunov): The equilibrium state x = 0 is said to be stable if, for any R > 0,

More information

ECE 314 Signals and Systems Fall 2012

ECE 314 Signals and Systems Fall 2012 ECE 31 ignals and ystems Fall 01 olutions to Homework 5 Problem.51 Determine the impulse response of the system described by y(n) = x(n) + ax(n k). Replace x by δ to obtain the impulse response: h(n) =

More information

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

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 OH BOY! O h Boy!, was or igin a lly cr eat ed in F r en ch an d was a m a jor s u cc ess on t h e Fr en ch st a ge f or young au di enc es. It h a s b een s een by ap pr ox i ma t ely 175,000 sp ect at

More information

Section 5.8. Taylor Series

Section 5.8. Taylor Series Difference Equations to Differential Equations Section 5.8 Taylor Series In this section we will put together much of the work of Sections 5.-5.7 in the context of a discussion of Taylor series. We begin

More information

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

1. Type your solutions. This homework is mainly a programming assignment. THE UNIVERSITY OF TEXAS AT SAN ANTONIO EE 5243 INTRODUCTION TO CYBER-PHYSICAL SYSTEMS H O M E W O R K S # 6 + 7 Ahmad F. Taha October 22, 2015 READ Homework Instructions: 1. Type your solutions. This homework

More information

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

DO NOT DO HOMEWORK UNTIL IT IS ASSIGNED. THE ASSIGNMENTS MAY CHANGE UNTIL ANNOUNCED. EE 533 Homeworks Spring 07 Updated: Saturday, April 08, 07 DO NOT DO HOMEWORK UNTIL IT IS ASSIGNED. THE ASSIGNMENTS MAY CHANGE UNTIL ANNOUNCED. Some homework assignments refer to the textbooks: Slotine

More information

ECE7850 Lecture 7. Discrete Time Optimal Control and Dynamic Programming

ECE7850 Lecture 7. Discrete Time Optimal Control and Dynamic Programming ECE7850 Lecture 7 Discrete Time Optimal Control and Dynamic Programming Discrete Time Optimal control Problems Short Introduction to Dynamic Programming Connection to Stabilization Problems 1 DT nonlinear

More information

Course Summary Math 211

Course Summary Math 211 Course Summary Math 211 table of contents I. Functions of several variables. II. R n. III. Derivatives. IV. Taylor s Theorem. V. Differential Geometry. VI. Applications. 1. Best affine approximations.

More information

Translational and Rotational Dynamics!

Translational and Rotational Dynamics! Translational and Rotational Dynamics Robert Stengel Robotics and Intelligent Systems MAE 345, Princeton University, 217 Copyright 217 by Robert Stengel. All rights reserved. For educational use only.

More information

Least-squares data fitting

Least-squares data fitting EE263 Autumn 2015 S. Boyd and S. Lall Least-squares data fitting 1 Least-squares data fitting we are given: functions f 1,..., f n : S R, called regressors or basis functions data or measurements (s i,

More information

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

Control Systems! Copyright 2017 by Robert Stengel. All rights reserved. For educational use only. Control Systems Robert Stengel Robotics and Intelligent Systems MAE 345, Princeton University, 2017 Analog vs. digital systems Continuous- and Discretetime Dynamic Models Frequency Response Transfer Functions

More information

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

1. Linearity of a Function A function f(x) is defined linear if. f(αx 1 + βx 2 ) = αf(x 1 ) + βf(x 2 ) 1. Linearity of a Function A function f(x) is defined linear if f(αx 1 + βx 2 ) αf(x 1 ) + βf(x 2 ) where α and β are scalars. Example of a linear function: f(x) 2x A nonlinear function: What about f(x)

More information

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

5. Hand in the entire exam booklet and your computer score sheet. WINTER 2016 MATH*2130 Final Exam Last name: (PRINT) First name: Student #: Instructor: M. R. Garvie 19 April, 2016 INSTRUCTIONS: 1. This is a closed book examination, but a calculator is allowed. The test

More information

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.

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. Algebra II Midterm Exam Review Solve: R view of Algebra 1 4. 215x + 31 = 16 /5xt3 1:: & 3 :: 2.1 3" ::. -J5 /,:::-=r\ or L":: -j) 2. II-2xl = 15 / j - ;).'1 -:.115 00( 1-).)(":.-15 - X-: 1"1 by:-t3-8 5

More information

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

Suppose that we have a specific single stage dynamic system governed by the following equation: Dynamic Optimisation Discrete Dynamic Systems A single stage example Suppose that we have a specific single stage dynamic system governed by the following equation: x 1 = ax 0 + bu 0, x 0 = x i (1) where

More information

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:

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: 1. (A). (C). (B). (B) 5. (A) 6. (B) 7. (B) 8. (A) 9. (A) 1. (D) 11. (C) 1. (B) 1. (A) 1. (D) 15. (B) 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 5 XII-IC SPARK Date: 9.1.17

More information

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

AM 205: lecture 14. Last time: Boundary value problems Today: Numerical solution of PDEs AM 205: lecture 14 Last time: Boundary value problems Today: Numerical solution of PDEs ODE BVPs A more general approach is to formulate a coupled system of equations for the BVP based on a finite difference

More information

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

Laplace Transforms. Chapter 3. Pierre Simon Laplace Born: 23 March 1749 in Beaumont-en-Auge, Normandy, France Died: 5 March 1827 in Paris, France Pierre Simon Laplace Born: 23 March 1749 in Beaumont-en-Auge, Normandy, France Died: 5 March 1827 in Paris, France Laplace Transforms Dr. M. A. A. Shoukat Choudhury 1 Laplace Transforms Important analytical

More information

Math 106: Review for Exam II - SOLUTIONS

Math 106: Review for Exam II - SOLUTIONS Math 6: Review for Exam II - SOLUTIONS INTEGRATION TIPS Substitution: usually let u a function that s inside another function, especially if du (possibly off by a multiplying constant) is also present

More information

Math 216 Final Exam 24 April, 2017

Math 216 Final Exam 24 April, 2017 Math 216 Final Exam 24 April, 2017 This sample exam is provided to serve as one component of your studying for this exam in this course. Please note that it is not guaranteed to cover the material that

More information

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

g(t) = f(x 1 (t),..., x n (t)). Reading: [Simon] p. 313-333, 833-836. 0.1 The Chain Rule Partial derivatives describe how a function changes in directions parallel to the coordinate axes. Now we shall demonstrate how the partial derivatives

More information

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

Trade Patterns, Production networks, and Trade and employment in the Asia-US region Trade Patterns, Production networks, and Trade and employment in the Asia-U region atoshi Inomata Institute of Developing Economies ETRO Development of cross-national production linkages, 1985-2005 1985

More information

Ritz Method. Introductory Course on Multiphysics Modelling

Ritz Method. Introductory Course on Multiphysics Modelling Ritz Method Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research of the Polish Academy of Sciences Warsaw Poland

More information

Taylor series - Solutions

Taylor series - Solutions Taylor series - Solutions. f(x) sin(x) sin(0) + x cos(0) + x x ( sin(0)) +!! ( cos(0)) + + x4 x5 (sin(0)) + 4! 5! 0 + x + 0 x x! + x5 5! x! + 0 + x5 (cos(0)) + x6 6! ( sin(0)) + x 7 7! + x9 9! 5! + 0 +

More information

Lecture two. January 17, 2019

Lecture two. January 17, 2019 Lecture two January 17, 2019 We will learn how to solve rst-order linear equations in this lecture. Example 1. 1) Find all solutions satisfy the equation u x (x, y) = 0. 2) Find the solution if we know

More information

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

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 Math 6 Fall 4 Exam. October 3, 4. The following questions have to do with the integral (a) Evaluate dx. Use integration by parts (x 3 dx = ) ( dx = ) x3 x dx = x x () dx = x + x x dx = x + x 3 dx dx =

More information

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

Jim Lambers MAT 285 Spring Semester Practice Exam 2 Solution. y(t) = 5 2 e t 1 2 e 3t. . Solve the initial value problem which factors into Jim Lambers MAT 85 Spring Semester 06-7 Practice Exam Solution y + 4y + 3y = 0, y(0) =, y (0) =. λ + 4λ + 3 = 0, (λ + )(λ + 3) = 0. Therefore, the roots

More information

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

Perspective. ECE 3640 Lecture 11 State-Space Analysis. To learn about state-space analysis for continuous and discrete-time. Objective: systems ECE 3640 Lecture State-Space Analysis Objective: systems To learn about state-space analysis for continuous and discrete-time Perspective Transfer functions provide only an input/output perspective of

More information

Partial differential equations (ACM30220)

Partial differential equations (ACM30220) (ACM3. A pot on a stove has a handle of length that can be modelled as a rod with diffusion constant D. The equation for the temperature in the rod is u t Du xx < x

More information

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

A L A BA M A L A W R E V IE W A L A BA M A L A W R E V IE W Volume 52 Fall 2000 Number 1 B E F O R E D I S A B I L I T Y C I V I L R I G HT S : C I V I L W A R P E N S I O N S A N D TH E P O L I T I C S O F D I S A B I L I T Y I N

More information

Math review. Math review

Math review. Math review Math review 1 Math review 3 1 series approximations 3 Taylor s Theorem 3 Binomial approximation 3 sin(x), for x in radians and x close to zero 4 cos(x), for x in radians and x close to zero 5 2 some geometry

More information

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.

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. ECE 35 Spring - Final Exam 9 May ECE 35 Signals and Systems Spring Final Exam - Solutions Three 8 ½ x sheets of notes, and a calculator are allowed during the exam Write all answers neatly and show your

More information

Differential Equation Types. Moritz Diehl

Differential Equation Types. Moritz Diehl Differential Equation Types Moritz Diehl Overview Ordinary Differential Equations (ODE) Differential Algebraic Equations (DAE) Partial Differential Equations (PDE) Delay Differential Equations (DDE) Ordinary

More information

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

Quasi-linear first order equations. Consider the nonlinear transport equation Quasi-linear first order equations Consider the nonlinear transport equation u t + c(u)u x = 0, u(x, 0) = f (x) < x < Quasi-linear first order equations Consider the nonlinear transport equation u t +

More information

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

ECEN 420 LINEAR CONTROL SYSTEMS. Lecture 6 Mathematical Representation of Physical Systems II 1/67 1/67 ECEN 420 LINEAR CONTROL SYSTEMS Lecture 6 Mathematical Representation of Physical Systems II State Variable Models for Dynamic Systems u 1 u 2 u ṙ. Internal Variables x 1, x 2 x n y 1 y 2. y m Figure

More information

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

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 ~ SOME NEW FINITE ELEMENT METHODS FOR THE ANALYSIS OF SHOCK AND ACCELERATION WAVES IN NONLINEAR MATERIALS J. T. ODEN The University of Teas at Austin L. C. WELLFORD, Jr. University of Southern California

More information

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

Differential Equations 2280 Sample Midterm Exam 3 with Solutions Exam Date: 24 April 2015 at 12:50pm Differential Equations 228 Sample Midterm Exam 3 with Solutions Exam Date: 24 April 25 at 2:5pm Instructions: This in-class exam is 5 minutes. No calculators, notes, tables or books. No answer check is

More information

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

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 fl ) 2 ;;:i c.l l) ( # =S >' 5 ^'R 1? l.y i.i.9 9 P * v ,>f { e e v? 'il v * 5,?: S 'V i: :i (g Y 1,Y iv cg G J :< >,c Z^ /^ c..l Cl i l 1 3 11 5 (' \ h 9 J :'i g > _ ^ j, \= f{ '{l #l * C\? 0l = 5,

More information

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

MATH4406 (Control Theory) Unit 6: The Linear Quadratic Regulator (LQR) and Model Predictive Control (MPC) Prepared by Yoni Nazarathy, Artem MATH4406 (Control Theory) Unit 6: The Linear Quadratic Regulator (LQR) and Model Predictive Control (MPC) Prepared by Yoni Nazarathy, Artem Pulemotov, September 12, 2012 Unit Outline Goal 1: Outline linear

More information

Intro. Computer Control Systems: F8

Intro. Computer Control Systems: F8 Intro. Computer Control Systems: F8 Properties of state-space descriptions and feedback Dave Zachariah Dept. Information Technology, Div. Systems and Control 1 / 22 dave.zachariah@it.uu.se F7: Quiz! 2

More information

MAHALAKSHMI ENGINEERING COLLEGE-TRICHY

MAHALAKSHMI ENGINEERING COLLEGE-TRICHY DIGITAL SIGNAL PROCESSING UNIT-I PART-A DEPT. / SEM.: CSE/VII. Define a causal system? AUC APR 09 The causal system generates the output depending upon present and past inputs only. A causal system is

More information

MATH 23 Exam 2 Review Solutions

MATH 23 Exam 2 Review Solutions MATH 23 Exam 2 Review Solutions Problem 1. Use the method of reduction of order to find a second solution of the given differential equation x 2 y (x 0.1875)y = 0, x > 0, y 1 (x) = x 1/4 e 2 x Solution

More information

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

Exact and high order discretization schemes. Wishart processes and their affine extensions for Wishart processes and their affine extensions CERMICS, Ecole des Ponts Paris Tech - PRES Université Paris Est Modeling and Managing Financial Risks -2011 Plan 1 Motivation and notations 2 Splitting

More information

Lecture Discrete dynamic systems

Lecture Discrete dynamic systems Chapter 3 Low-level io Lecture 3.4 Discrete dynamic systems Lecture 3.4 Discrete dynamic systems Suppose that we wish to implement an embedded computer system that behaves analogously to a continuous linear

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

Control for Coordination of Linear Systems

Control for Coordination of Linear Systems Control for Coordination of Linear Systems André C.M. Ran Department of Mathematics, Vrije Universiteit De Boelelaan 181a, 181 HV Amsterdam, The Netherlands Jan H. van Schuppen CWI, P.O. Box 9479, 19 GB

More information

Lecture 2: Review of Prerequisites. Table of contents

Lecture 2: Review of Prerequisites. Table of contents Math 348 Fall 217 Lecture 2: Review of Prerequisites Disclaimer. As we have a textbook, this lecture note is for guidance and supplement only. It should not be relied on when preparing for exams. In this

More information

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

Optimal Control. Quadratic Functions. Single variable quadratic function: Multi-variable quadratic function: Optimal Control Control design based on pole-placement has non unique solutions Best locations for eigenvalues are sometimes difficult to determine Linear Quadratic LQ) Optimal control minimizes a quadratic

More information

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

ON THE AVERAGE NUMBER OF REAL ROOTS OF A RANDOM ALGEBRAIC EQUATION ON THE AVERAGE NUMBER OF REAL ROOTS OF A RANDOM ALGEBRAIC EQUATION M. KAC 1. Introduction. Consider the algebraic equation (1) Xo + X x x + X 2 x 2 + + In-i^" 1 = 0, where the X's are independent random

More information

Advanced Quantum Mechanics

Advanced Quantum Mechanics Advanced Quantum Mechanics Rajdeep Sensarma! sensarma@theory.tifr.res.in Lecture #22 Path Integrals and QM Recap of Last Class Statistical Mechanics and path integrals in imaginary time Imaginary time

More information

Solution of Linear State-space Systems

Solution of Linear State-space Systems Solution of Linear State-space Systems Homogeneous (u=0) LTV systems first Theorem (Peano-Baker series) The unique solution to x(t) = (t, )x 0 where The matrix function is given by is called the state

More information

Solution: APPM 1350 Final Exam Spring 2014

Solution: APPM 1350 Final Exam Spring 2014 APPM 135 Final Exam Spring 214 1. (a) (5 pts. each) Find the following derivatives, f (x), for the f given: (a) f(x) = x 2 sin 1 (x 2 ) (b) f(x) = 1 1 + x 2 (c) f(x) = x ln x (d) f(x) = x x d (b) (15 pts)

More information