Prof. Krstic Nonlinear Systems MAE281A Homework set 1 Linearization & phase portrait

Similar documents
Solution of Additional Exercises for Chapter 4

Nonlinear Systems and Control Lecture # 12 Converse Lyapunov Functions & Time Varying Systems. p. 1/1

Lecture 4. Chapter 4: Lyapunov Stability. Eugenio Schuster. Mechanical Engineering and Mechanics Lehigh University.

LMI Methods in Optimal and Robust Control

EN Nonlinear Control and Planning in Robotics Lecture 3: Stability February 4, 2015

Introduction to Nonlinear Control Lecture # 3 Time-Varying and Perturbed Systems

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

Complex Dynamic Systems: Qualitative vs Quantitative analysis

Nonlinear Systems February 12, Late submissions will not be accepted. Collaboration not allowed. ẋ = cx + y 2m x cos 2 x

Nonlinear Control. Nonlinear Control Lecture # 8 Time Varying and Perturbed Systems

MCE693/793: Analysis and Control of Nonlinear Systems

Part II. Dynamical Systems. Year

Modeling and Analysis of Dynamic Systems

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

1 Lyapunov theory of stability

Nonlinear systems. Lyapunov stability theory. G. Ferrari Trecate

CDS Solutions to the Midterm Exam

Entrance Exam, Differential Equations April, (Solve exactly 6 out of the 8 problems) y + 2y + y cos(x 2 y) = 0, y(0) = 2, y (0) = 4.

3 Stability and Lyapunov Functions

Applied Math Qualifying Exam 11 October Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems.

UNIVERSITY OF MANITOBA

MCE/EEC 647/747: Robot Dynamics and Control. Lecture 8: Basic Lyapunov Stability Theory

Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games

Dynamical Systems & Lyapunov Stability

MCE693/793: Analysis and Control of Nonlinear Systems

Converse Lyapunov theorem and Input-to-State Stability

Lecture 8. Chapter 5: Input-Output Stability Chapter 6: Passivity Chapter 14: Passivity-Based Control. Eugenio Schuster.

1.7. Stability and attractors. Consider the autonomous differential equation. (7.1) ẋ = f(x),

STABILITY. Phase portraits and local stability

Problem set 7 Math 207A, Fall 2011 Solutions

Nonlinear Systems Theory

Nonlinear Control. Nonlinear Control Lecture # 2 Stability of Equilibrium Points

ODE Homework 1. Due Wed. 19 August 2009; At the beginning of the class

Chapter III. Stability of Linear Systems

Half of Final Exam Name: Practice Problems October 28, 2014

Georgia Institute of Technology Nonlinear Controls Theory Primer ME 6402

EE363 homework 7 solutions

E209A: Analysis and Control of Nonlinear Systems Problem Set 6 Solutions

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

Nonlinear Control. Nonlinear Control Lecture # 3 Stability of Equilibrium Points

154 Chapter 9 Hints, Answers, and Solutions The particular trajectories are highlighted in the phase portraits below.

Stabilization of persistently excited linear systems

Nonlinear Control Systems

High-Gain Observers in Nonlinear Feedback Control. Lecture # 3 Regulation

L2 gains and system approximation quality 1

Lyapunov Stability Theory

Nonlinear Control Lecture 5: Stability Analysis II

EE222 - Spring 16 - Lecture 2 Notes 1

System Control Engineering 0

CDS 101/110a: Lecture 2.1 Dynamic Behavior

Stability in the sense of Lyapunov

TTK4150 Nonlinear Control Systems Solution 6 Part 2

Pattern generation, topology, and non-holonomic systems

2 Lyapunov Stability. x(0) x 0 < δ x(t) x 0 < ɛ

Output Feedback and State Feedback. EL2620 Nonlinear Control. Nonlinear Observers. Nonlinear Controllers. ẋ = f(x,u), y = h(x)

Problem set 3: Solutions Math 207A, Fall where a, b are positive constants, and determine their linearized stability.

Now I switch to nonlinear systems. In this chapter the main object of study will be

Chapter 4: First-order differential equations. Similarity and Transport Phenomena in Fluid Dynamics Christophe Ancey

Nonlinear Control Theory. Lecture 9

Linear Differential Equations. Problems

EG4321/EG7040. Nonlinear Control. Dr. Matt Turner

CDS 101/110a: Lecture 2.1 Dynamic Behavior

Introduction to Nonlinear Control Lecture # 4 Passivity

2.10 Saddles, Nodes, Foci and Centers

Nonlinear Control Systems

Solution of Linear State-space Systems

BIBO STABILITY AND ASYMPTOTIC STABILITY

Nonlinear Control. Nonlinear Control Lecture # 24 State Feedback Stabilization

Nonlinear Autonomous Systems of Differential

Lecture 9. Systems of Two First Order Linear ODEs

Nonlinear Control. Nonlinear Control Lecture # 8 Time Varying and Perturbed Systems

The first order quasi-linear PDEs

High-Gain Observers in Nonlinear Feedback Control. Lecture # 2 Separation Principle

B5.6 Nonlinear Systems

Lecture 10: Singular Perturbations and Averaging 1

Applied Dynamical Systems

7.1 Linear Systems Stability Consider the Continuous-Time (CT) Linear Time-Invariant (LTI) system

ME Fall 2001, Fall 2002, Spring I/O Stability. Preliminaries: Vector and function norms

Nonlinear Control. Nonlinear Control Lecture # 6 Passivity and Input-Output Stability

+ i. cos(t) + 2 sin(t) + c 2.

Lecture 7. Please note. Additional tutorial. Please note that there is no lecture on Tuesday, 15 November 2011.

CDS 101 Precourse Phase Plane Analysis and Stability

ASTATISM IN NONLINEAR CONTROL SYSTEMS WITH APPLICATION TO ROBOTICS

Poincaré Map, Floquet Theory, and Stability of Periodic Orbits

MATH 220: MIDTERM OCTOBER 29, 2015

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

Chapter 3 Second Order Linear Equations

Feedback Linearization

Analysis of Dynamical Systems

7 Planar systems of linear ODE

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

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

On the Stability of the Best Reply Map for Noncooperative Differential Games

11 Chaos in Continuous Dynamical Systems.

Asymptotic Stability by Linearization

High-Gain Observers in Nonlinear Feedback Control

Nonlinear Dynamical Systems Lecture - 01

Practice Problems for Final Exam

EE363 homework 8 solutions

PROJECTS on. Multiple Bifurcations of Sample Dynamical Systems

Transcription:

Prof. Krstic Nonlinear Systems MAE28A Homework set Linearization & phase portrait. For each of the following systems, find all equilibrium points and determine the type of each isolated equilibrium. Use Matlab to compute the eigenvalues. () ẋ = x 2 (2) ẋ = x +x 2 ẋ 2 = x + 6 x3 x 2 ẋ 2 = 0.x 2x 2 x 2 0.x 3 (3) ẋ = x +x 2 (+x ) (4) ẋ = x 3 +x 2 ẋ 2 = x (+x ) ẋ 2 = x x 3 2 2. The phase portrait of the following four systems are shown in Figure : parts (a), (b), (c), and (d), respectively. Mark the arrowheads and discuss the qualitative behavior of each system. () ẋ = x 2 (2) ẋ = x 2 ẋ 2 = x x 2 ( x 2 +0.x 4 ) ẋ 2 = x +x 2 3tan (x +x 2 ) (3) ẋ = x 2 (4) ẋ = x 2 ẋ 2 = (0.5x +x 3 ) ẋ 2 = x 2 ψ(x x 2 ) Where ψ(y) = y 3 +0.5y if y and ψ(y) = 2y 0.5 if y >. Figure : Exercise 2

Prof. Krstic Nonlinear Systems MAE28A Homework set 2 Gronwall s inequality & comparison principle. Let y(t) be a nonnegative scalar function that satisfies the inequality y(t) k e α(t t 0) + t t 0 e α(t τ)[ k 2 y(τ)+k 3 ] dτ () where k,k 2, and k 3 are nonnegative constants and α is a positive constant that satisfies α > k 2. Using the Gronwall-Bellman inequality, show that y(t) k e (α k 2)(t t 0 ) + k 3 α k 2 [ e (α k 2 )(t t 0 ) ] (2) Hint: Take z(t) = y(t)e α(t t 0) and find the inequality satisfied by z. This exercise shows how to apply Gronwall lemma to the case: y(t) < λ(t)+ρ(t) µ(τ)y(τ)dτ. 2. Using the comparison principle, show that if ν,l, and l 2 are functions that satisfy and if c > 0, then where denotes the L norm defined as ν cν +l (t)ν +l 2 (t), ν(0) 0 (3) ν(t) ( ν(0)e ct + l 2 ) e l, (4) Using Gronwall s lemma show that f = 0 f(t) dt. (5) ν(t) ( )( ν(0)e ct + l 2 + l e l ). (6) Which of the two bounds is less conservative? 3. Consider the system: ẋ = cx+y 2m xcos 2 (x) (7) ẏ = y 3. (8) Using either Gronwall s inequality or the comparison principle, show that a) x(t) is bounded for all t 0 whenever c = 0 and m >. b) x(t) 0 as t whenever c > 0 and m =.

4. Consider the system ẋ = x+yxsin(x) (9) ẏ = y +zysin(y) (0) ż = z. () Using Gronwall s lemma (twice), show that x(t) x 0 e y 0 e z 0 e t, t 0. (2) 2

Prof. Krstic Nonlinear Systems MAE28A Homework set 3 Sensitivity equations. Derive the sensitivity equations for the system ẋ = tan (αx ) x x 2 () ẋ 2 = bx 2 cx 2 (2) as the parameters a,b, and c vary from their nominal values a 0 =,b 0 = 0, and c 0 =. 2. Calculate exactly (in closed form) the sensitivity function at λ 0 = for the system ẋ = λx 3. (3) What is the approximation x(t,λ) x(t,λ 0 )+S(t)(λ λ 0 ) for λ = 7/2?

Prof. Krstic Nonlinear Systems MAE28A Homework set 4 Lyapunov stability - autonomous systems. Prove global stability of the origin of the system ẋ = x 2 () ẋ 2 = x. (2) +x 2 2 2. Prove global asymptotatic stability of the origin of the system ẋ = x 3 2 (3) ẋ 2 = x x 2. (4) 3. Prove global asymptotatic stability of the origin of the system ẋ = x 2 (2x 2 +x 2 2)x (5) ẋ 2 = x 2(2x 2 +x 2 2)x 2. (6) Is the origin locally exponentially stable and why or why not? 4. Consider the system ẋ = x +x x 2 (7) ẋ 2 = x2. (8) +x 2 Show that the equilibrium x = x 2 = 0 is globally stable and that Hint: Seek a Lyapunov function in the form lim x (t) = 0 (9) t V(x,x 2 ) = φ(x )+x 2 2, (0) where the function φ( ) is to be found. Make sure that your Lyapunov function V(, ) is positive definite and radially unbounded.

Prof. Krstic Nonlinear Systems MAE28A Homework set 5 Chetaev s theorem. With Chetaev s theorem, show that the equilibrium at the origin of the following three systems is unstable: a) ẋ = x 3 +xy 3 () ẏ = y +x 2 (2) b) ξ = η +ξ 3 +3ξη 2 (3) η = ξ +η 3 +3ηξ 2 (4) c) ẋ = x x+xy y (5) ẏ = y + x y. (6)

Prof. Krstic Nonlinear Systems MAE28A Homework set 6 LaSalle s theorem. Using the Lyapunov function candidate V = 4 x4 + 2 y2 + 4 z4 () study stability of the origin of the system ẋ = y (2) ẏ = x 3 y 3 z 3 (3) ż = z +y. (4) 2. Consider the system ẋ = x+yx+zcos(x) (5) ẏ = x 2 (6) ż = xcos(x) (7) a) Determine all the equilibria of the system. b) Show that the equilibrium x = y = z = 0 is globally stable. c) Show that x(t) 0 as t. d) Show that z(t) 0 as t. 3. Which of the state variables of the following system are guaranteed to converge to zero from any initial condition? ẋ = x 2 +x x 3 (8) ẋ 2 = x x 2 +x 2 x 3 (9) ẋ 3 = x 2 x 2 2. (0)

Prof. Krstic Nonlinear Systems MAE28A Homework set 7 Perturbation theory. Consider the system ẋ = x 2 () ẋ 2 = x x 2 +ǫx 3. (2) a) Find an O(ǫ) approximation. b) Find an O(ǫ 2 ) approximation. c) Investigate the validity of the approximation on the infinite interval. 2. Repeat exercise for the system ẋ = x +x 2 (+x )+ǫ(+x ) 2 (3) ẋ 2 = x (x +). (4)

Prof. Krstic Nonlinear Systems MAE28A Homework set 8 Averaging theory. Using averaging theory, analyze the following system: ẋ = ǫ [ x+ 2(y +sin(t)) 2] () ẏ = ǫz (2) [ ż = ǫ z sin(t) [ 2 x+(y +sin(t))2]]. (3) 2. Analyze the following system using the method of averaging for large ω: 3. Consider the second-order system ẋ = ( x 2 sin(ωt) 2 ) x x 3 (4) ẋ 2 = x 2 + ( x 2 2sin(ωt) 2x 3 cos(ωt) ) cos(ωt) (5) ẋ 3 = 2x 2 sin(x 3 )+ ( 4x 2 sin(ωt)+x 3 ) sin(ωt) (6) ẋ = sin(ωt)y (7) ẋ 2 = cos(ωt)y 2 (8) y = [ x +sin(ωt) ][ x 2 +cos(ωt) x sin(ωt) ] (9) y 2 = [ x 2 +cos(ωt) ][ x +sin(ωt) x 2 cos(ωt) ]. (0) Show that for sufficiently large ω there exists an exponentially stable periodic orbit in an O(/ω) neighborhood of the origin x = x 2 = 0. Hint: the following functions have a zero mean over the interval [0,2π] : sin(τ),cos(τ), sin(τ)cos(τ),sin 3 (τ),cos 3 (τ),sin 2 (τ)cos(τ), and sin(τ)cos 2 (τ). 4. Consider Rayleigh s equation m d2 u dt 2 +ku = λ [ α where m,k,λ, and α are positive constants. ( du dt ) 2 ] du dt () a) Using the dimensionless variables y = u u,τ = t t, and ǫ = λ λ, where (u ) 2 αk = m 3, t = m k, and λ = km, show that the equation can be normalized to where ẏ denotes the derivative of y with respect to τ. ÿ +y = ǫ(ẏ 3ẏ3 ) (2) b) Apply the averaging method to show that the normalized Rayleigh equation has a stable limit cycle. Estimate the location of the limit cycle in the plane (y,ẏ).

Prof. Krstic Nonlinear Systems MAE28A Homework set 9 Singular perturbation. Using singular perturbation theory, study local exponential stability of the origin of the system Is the origin globally exponentially stable? 2. Consider the following control system: ẋ = y () ẏ = z (2) ǫż = z +sin(x)+y. (3) ẋ = A x+a 2 z +B u (4) ǫż = A 2 x+a 22 z. (5) Assume that the matrix A 22 is Hurwitz and that there exists a matrix/vector K (of appropriate dimensions) such that A A 2 A 22A 2 +B K (6) is also Hurwitz. Prove that the partial-state feedback law u = Kx (7) exponentially stabilizes the equilibrium (x, z) = (0, 0) for sufficiently small ǫ. 3. Find the exact slow manifold of the singularly perturbed system ẋ = xz 3 (8) ǫż = z x 4/3 + 4 3 ǫx6/3. (9) Hint: try substitute the quasi-steady state into the (exact) manifold condition. 4. How many slow manifolds does the following system have? Which of these manifolds will attract trajectories of the system? ẋ = xz (0) ǫż = (z sin 2 (x))(z e ax )(z 2e 2ax ), a > 0. ()

Prof. Krstic Nonlinear Systems MAE28A Homework set 0 Input-to-state stability. Show that the following system is ISS and determine its gain function: ẋ = x 3 +xu. () 2. Show that the following system is ISS and determine its gain function: ẋ = x+u 3. (2) 3. Show that the following system is ISS and guess its gain function: ẋ = x 3 +xy (3) ẏ = y +u 3. (4) 4. Consider the system Show that this system is ISS using the Lyapunov function ẋ = x+y 3 (5) x ẏ = y +x 2 +z2 (6) ż = z +u. (7) V = +x 2 + 4 y4 + 2 z8. (8) 5. Show that the following system is ISS ẋ = x+x /3 y +p 2 (9) ẏ = y x 4/3 +p 3 (0) ṗ = p+u. ()

Prof. Krstic Nonlinear Systems MAE28A Homework set Center manifold theorem. For each of the following systems, investigate the stability of the origin using the center manifold theorem. a) ẋ = x 2 2 () ẋ 2 = x 2 +x 2 +x x 2 (2) b) ẋ = ax 2 x 2 2, a 0 (3) ẋ 2 = x 2 +x 2 +x x 2 (4) 2. Consider the system ẋ = x x 2 +ax 3 +bx x 2 2 (5) ẋ 2 = x 2 +cx 2 +dx 2 x 2 (6) Investigate the stability of the origin using the center manifold theorem for each of the following cases i) a+c > 0 ii) a+c < 0 iii) a+c = 0 and cd+bc 2 < 0 iv) a+c = 0 and cd+bc 2 > 0 3. Using the center manifold theorem, determine whether the origin of the following system is asymptotically stable: ẋ = x 2 +x x 3 (7) ẋ 2 = x +x 2 x 3 (8) ẋ 3 = x 3 (x 2 +x 2 2)+x 2 3. (9) 4. Using the center manifold theorem, determine whether the origin of the following system is asymptotically stable: ẏ = yz +2y 3 (0) ż = z 2y 2 4y 4 2y 2 z. ()

Prof. Krstic Nonlinear Systems MAE28A Homework set 2 Describing function. For each of odd nonlinearity ψ(y) on the following list, verify the given expression of the describing function Ψ(a). i) ψ(y) = y 5 ; Ψ(a) = 5a 4 /8 ii) ψ(y) = y 3 y ; Ψ(a) = 32a 3 /5π iii) ψ(y) : Figure ; Ψ(a) = k +4A/aπ ψ(y) A Slope=k y Figure 0.7 A Figure : ψ(y) Figure 2: Feedback connection 2. Consider the feedback connection of Figure 2 with G(s) = (s+) 6 and ψ(y) = sgn(y). () Using the describing function method, investigate the possibility of existence of periodic solutions and the possible frequency and amplitude of oscillation. 3. Repeat exercise 2 with G(s) = s s 2 s+ and ψ(y) = y 5. (2) 4. Consider the feedback system with a linear block s(s+)(s+2) andanonlinearitysgn(y)+ y y (note that it is an odd nonlinearity, so the describing functions method applies, and note that the first term was already studied in class). First, find the describing function for the nonlinearity. Then, determine if the feedback system is likely to have any periodic solutions.