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

Size: px
Start display at page:

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

Transcription

1 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 and Li, etc. For full credit, show all work. Some problems require hand calculations. In those cases, do not use MATLAB except to check your answers. It is OK to talk about the homework beforehand. BUT, once you start writing the answers, MAKE SURE YOU WORK ALONE. The purpose of the Homework is to evaluate you individually, not to evaluate a team. Cheating on the homework will be severely punished. The next page must be signed and turned in at the front of ALL homeworks submitted in this course.

2 EE 533 Nonlinear Control Systems Homework Pledge of Honor On all homeworks in this class - YOU MUST WORK ALONE. Any cheating or collusion will be severely punished. It is very easy to compare your software code and determine if you worked together It does not matter if you change the variable names. Please sign this form and include it as the first page of all of your submitted homeworks Typed Name: Pledge of honor: "On my honor I have neither given nor received aid on this homework. e-signature:

3 EE 533 Homework State Variable Systems, Computer Simulation. Simulate the van der Pol oscillator y" ( y ) y' y 0 using MATLAB for various ICs. Plot y(t) vs. t and also the phase plane plot y'(t) vs. y(t). Use y(0)=0., y'(0)= 0. a. For = b. For= Do MATLAB simulation of the Lorenz Attractor chaotic system. Run for 50 sec. with all initial states equal to 0.4. Plot states versus time, and also make 3-D plot of x, x, x3 using PLOT3(x,x,x3). ( x x ) 3 bx3 xx use = 0, r= 8, b= 8/3. rx x x x 3. Consider the Voltera predator-prey system xxx. x x x 3 Simulate the system using MATLAB for various initial conditions. Take ICs spaced in a uniform mesh in the box x=[-,], x=[-,]. Make one phase plane plot with all the trajectories on it. Plot phase plane on square [-5,5]x[-5,5]. 3

4 Nonlinear Systems and Equilibrium. Consider the Voltera predator-prey system xxx. x x x Find the equilibrium points and their nature. EE 533 Homework. Duffing s equation is interesting in that it exhibits bifurcation, or dependence of stability properties and number of equilibrium points on a parameter. The undamped Duffing equation is 3 x 0 a. Find the equilibrium points. Show that for 0 there is only one e.p. b. For 0 there are 3 eps. Linearize the system and study the nature of these 3 e.p.s c. Simulate the Duffing oscillator and make time plot and phase plane plot. Do for a. b. 0. c. 3. Consider the system y( y ). y x( yx ) Simulate the system using MATLAB for various initial conditions for the two cases: a. Take ICs spaced in a uniform mesh in the box x=[-0,0], x=[-0,0]. Make one phase plane plot with all the trajectories on it. Plot phase plane on square [-5,5]x[-5,5]. b. Take ICs spaced in a uniform mesh in the box x=[-3,3], x=[-3,3]. Make one phase plane plot with all the trajectories on it. Plot phase plane on square [-5,5]x[-5,5]. 4. The system of equations axbxx cx dx exx fx describes the growth of two competing species that prey on each other. The constants are positive parameters. Pick a=c=d=f=, b=e=3. Simulate the system using MATLAB for various initial conditions. Take ICs spaced in a uniform mesh in the box x=[-,], x=[-,]. Make one phase plane plot with all the trajectories on it. Plot phase plane on square [-5,5]x[-5,5]. 4

5 Vector Fields, Flows, First Integrals. Consider the undamped oscillator 0 EE 533 Homework 3 a. Write position-velocity state space form X f( X). b. Plot the trajectories x( t), ( t) vs. time. Use initial conditions of x(0) 0., ( t) 0 f( x, x) c. Plot the vector field f( X) f( x, x) in the phase plane ( x, x ) ( x, ). Plot for points spaced in a uniform mesh in the box x=[-0,0], x=[-0,0]. d. Plot the system trajectories (flows or orbits) in the phase plane. Take ICs spaced in a uniform mesh in the box x=[-0,0], x=[-0,0]. e. Derive the First Integral of Motion F( x, x ) as done in class. Plot the FIM as a 3-D surface over the phase plane on the x=[-0,0], x=[-0,0].. Repeat for the unstable system 0 5

6 EE 533 Homework 4 E.P.s, Lyapunov Stability Analysis. Equilibrium points and linearization System is x( xx ) x ( x x ) a. Find all equilibrium points b. Find Jacobian c. Find the nature of all e.p.s. Use Lyapunov to examine the stability of these systems. Simulate time histories from many uniformly spaced ICs to verify your results. a. b. xx x xx x x ( x ) x 3. Use Lyapunov to show that the system xx x ( x x 3) x x x ( x x 3) is locally asymptotically stable. Find the Region of Asymptotic Stability. Simulate the system from many uniformly spaced ICs. 3. Use Lyapunov to show that the system x x ( x x ) x x ( x x ) is UUB. Simulate the system from many uniformly spaced ICs. 6

7 EE 533 Homework 5 Lyapunov s Method. UUB of system with disturbance. Consider the system on S&L p. 66 with a disturbance d c( x) d 0 Assume that xc( x) ax with a 0 a known positive constant a. Assume that d is unknown but is bounded by d D with D a known positive constant. Prove that the system is UUB and find the bound on x(t). b. Assume that d is unknown but is bounded by d D x with D a known positive constant. Prove that the system is UUB and find the bound on x(t).. UUB Use Lyapunov to show that the system x x x ( x x 3) x x x ( x x 3) is uniformly ultimately bounded UUB. That is, show that the Lyapunov derivative is NEGATIVE OUTSIDE A BOUNDED REGION. Find the radius of the bounded region outside which V <0. Any states outside this region are attracted towards the origin. 3. Use Lyapunov Equation to check the stability of the linear systems a. 0 Ax x 6 5 b. 4 Ax x c. 0 Ax x 4 0 7

8 EE 533 Homework 6 Lyapunov Controls Design, Feedback Linearization. On Slotine & Li p. 7 we used the Lyapunov function 4 V( x) ( x x 0) and we obtained Lyapunov derivative V ( x) (4x x)( x x 0) Plot V( x ) and V ( x) using MATLAB. Pick a region for the domains that reveals a nice plot.. A system is given by xsgn( x) xx u Select Lyapunov function candidate V( x) ( x x) Use Lyapunov to design a controller u(x) to make system SISL. 3. For each of these systems, Use Lyapunov to Design feedback control u(x) to make the system i. SISL, and then ii. AS a. b. c. xx x xu xx u xx xx x u 4. Multi-input Control. Use Lyapunov to design controls u, u to make this system i. SISL, and then ii. AS xx u xx u 3 7 8

9 EE 533 Homework 7 I/O Feedback Linearization. A system is given by sin x 4 x cos x u with output y() t x (t) a. Design a FB linearization controller to make the output follow a desired trajectory yd ( t ) That is, find ut ( ) b. Discuss the internal dynamics. Are they a problem?. NMP System A system is given by 4 y zy u 5 z( y) z z 0 a. Take output y( t ) and find the FB linearization controller ut ( ) to follow the prescribed trajectory yd ( t ). b. Find the internal dynamics. Set yt ( ) 0 to get the zero dynamics. Are the ZD stable? Does the FB linearization controller work? 3. Effect of Output Choice in i/o FB Linearization It is desired to stabilize a system given by x sin x x u x x a. Select the output as y x and use FB lin. design to select the control u(t) to follow the desired trajectory yd ( t ). Check the internal dynamics. Set y=0 to get the zero dynamics. Is the system minimum phase? b. Select the new output y x. Find the FB lin. controller u(t). Does this work? What about the internal dynamics? 9

10 DO NOT WORK BEYOND THIS PAGE The following hwks have not yet been assigned and may change. 0

11 EE 533 Homework 3 Chaos, Phase Plane. A system that exhibits chaos is the logistic function xk ( k ) xk ( xk ) However, chaos only occurs for certain values of sweep through the k values using k k for fixed less than but close to. These two equations form a dynamical system. Perform a MATLAB simulation to reproduce this plot of xk vs k, which was taken for and initial value of 0 3. Interpret the plot with some discussion in terms of bifurcation theory. Plot also k. Show your MATLAB code. It is indeed interesting that the logistic function appears in economic systems and military supply systems. k. Rather than try all values of k, we can. Slotine and Li p. 39 problem.. Simulate and plot phase plane trajectories for various ICs. Do not do the problem requested in the book. 3. Slotine and Li p. 35 Example.7b. Show that this system has an unstable limit cycle.

12 EE 533 Homework 4 Lyapunov Stability Analysis 4. Slotine and Li p. 97 problem Use Lyapunov to show that the system xx x ( x x 3) x x x ( x x 3) is locally asymptotically stable. Find the Region of Asymptotic Stability 6. Use Lyapunov to examine the stability of these systems. Simulate time histories from many uniformly spaced ICs to verify your results. c. d. e. xx x xx x sin x x x sin x x x x ( x ) x 7. Use Lyapunov to show that the system x x ( x x ) x x ( x x ) has a stable limit cycle. Simulate the system from many uniformly spaced ICs.

13 EE 533 Homework 5 Lyapunov s Method. UUB of system with disturbance. Consider the system on S&L p. 66 with a disturbance d c( x) d 0 Assume that xc( x) ax with a 0 a known positive constant c. Assume that d is unknown but is bounded by d D with D a known positive constant. Prove that the system is UUB and find the bound on x(t). d. Assume that d is unknown but is bounded by d D x with D a known positive constant. Prove that the system is UUB and find the bound on x(t).. UUB Use Lyapunov to show that the system x x x ( x x 3) x x x ( x x 3) is uniformly ultimately bounded UUB. That is, show that the Lyapunov derivative is NEGATIVE OUTSIDE A BOUNDED REGION. Find the radius of the bounded region outside which V <0. Any states outside this region are attracted towards the origin. 3. Use Lyapunov Equation to check the stability of the linear systems 0 d. x e. f Ax x 0 Ax x 4 0. Barbalat s Lemma and LaSalle extension a. Use quadratic Lyapunov Function to show this system is locally AS x x( x ) xx Find a ball within which V 0. This region is contained in the region of attraction. b. Use quadratic Lyapunov Function to show this system is locally SISL x x( x ) x Find a region within which V 0. 3

14 c. Use Barbalat s Lemma to verify that V 0. Check uniform continuity of the Lyapunov derivative of the system in part b. d. Use LaSalle s extension to verify that the system in part b is actually AS. Find the equilibrium point. 4

15 EE 533 Homework 6 Lyapunov Controls Design, I/O Feedback Linearization. A system is given by xsgn( x) xx u Select Lyapunov function candidate V( x) ( x x) Use Lyapunov to design a controller u(t) to make system SISL.. A system is given by sin x 4 x cos x u with output y() t x (t) a. Design a FB linearization controller to make the output follow a desired trajectory yd ( t ) That is, find ut ( ) b. Discuss the internal dynamics. Are they a problem? 3. Effect of Output Choice in i/o FB Linearization It is desired to stabilize a system given by x sin x x u x x c. Select the output as y x and use FB lin. design to select the control u(t) to follow the desired trajectory yd ( t ). Check the internal dynamics. Set y=0 to get the zero dynamics. Is the system minimum phase? d. Select the new output y x. Find the FB lin. controller u(t). Does this work? What about the internal dynamics? 5

16 Feedback Linearization, backstepping EE 533 Homework 7. I/O feedback linearization. Slotine and Li Problem 6.3. System is sin x 4 x cos x u y x Do i/o feedback linearization to make output track desired trajectory yd ( t ).. Backstepping. The system is sin x u 4 x cos x x Do backstepping to stabilize this system. Select the desired value x d to yield the first step dynamics of x 0. Compare this to Problem, which uses i/o FB linearization. 3. Backstepping. Slotine and Li Problem 6.. Globally stabilize means backstepping. System is 4 5 y zy u 5 z( y) z z 0 6

17 The following homeworks are out of date. Do not do homeworks until they are assigned. They may change. 7

18 Old hwk Obtain the linear model of the system described by y( y ) y y 0 around the equilibrium point. 9. The system of equations axbxx cx dx exx fx describes the growth of two competing species that prey on each other. The constants are positive parameters and it is assumed that the two states are positive. Determine the linear model of the system around the equilibrium point (0,0). 8

19 EE 533 Homework Nonlinear Systems and Equilibrium 0. Consider the Voltera predator-prey system xxx. x x x a. Find the equilibrium points and their nature. b. Simulate the system using MATLAB for various initial conditions. Take ICs spaced in a uniform mesh in the box x=[-,], x=[-,]. Make one phase plane plot with all the trajectories on it. Plot phase plane on square [-5,5]x[-5,5].. Consider the system y( y ). y x( yx ) Simulate the system using MATLAB for various initial conditions for the two cases: c. Take ICs spaced in a uniform mesh in the box x=[-0,0], x=[-0,0]. Make one phase plane plot with all the trajectories on it. Plot phase plane on square [-5,5]x[-5,5]. d. Take ICs spaced in a uniform mesh in the box x=[-3,3], x=[-3,3]. Make one phase plane plot with all the trajectories on it. Plot phase plane on square [-5,5]x[-5,5].. The system of equations axbxx cx dx exx fx describes the growth of two competing species that prey on each other. The constants are positive parameters. Simulate the system using MATLAB for various initial conditions. Make one phase plane plot with all the trajectories on it. 3. Duffing s equation is interesting in that it exhibits bifurcation, or dependence of stability properties and number of equilibrium points on a parameter. The undamped Duffing equation is 3 x 0 d. Find the equilibrium points. Show that for 0 there is only one e.p. e. For 0 there are 3 e.p.s Linearize the system and study the nature of these 3 e.p.s f. Simulate the Duffing oscillator for. Make time plot and phase plane plot. 9

20 EE 533 Homework 3 Chaos, Phase Plane. A system that exhibits chaos is the logistic function xk ( k ) xk ( xk ) However, chaos only occurs for certain values of sweep through the k values using k k for fixed less than but close to. These two equations form a dynamical system. Perform a MATLAB simulation to reproduce this plot of xk vs k, which was taken for and initial value of 0 3. Interpret the plot with some discussion in terms of bifurcation theory. Plot also k. Show your MATLAB code. It is indeed interesting that the logistic function appears in economic systems and military supply systems. k. Rather than try all values of. For Slotine & Li Example. on P. 0- a. Find equilibrium points b. Linearize the system about each equilibrium point. Find poles in each case. c. Simulate the system to find the Region of Attraction. k, we can 3. Slotine and Li p. 39 problem.. Simulate and plot phase plane trajectories for various ICs. Do not do the problem requested in the book. 4. Slotine and Li p. 35 Example.7b. Show that this system has an unstable limit cycle. 0

21 EE 533 Homework 4 Lyapunov Stability Analysis 4. Slotine and Li p. 97 problem Use Lyapunov to show that the system xx x ( x x 3) x x x ( x x 3) is locally asymptotically stable. Find the Region of Asymptotic Stability 6. Use Lyapunov to examine the stability of these systems. Simulate time histories from many uniformly spaced ICs to verify your results. f. g. h. xx x xx x sin x x x sin x x x x ( x ) x 7. Use Lyapunov to show that the system x x ( x x ) x x ( x x ) has a stable limit cycle. Simulate the system from many uniformly spaced ICs.

22 EE 533 Homework 5 Lyapunov s Method. UUB of system with disturbance. Consider the system on S&L p. 66 with a disturbance d c( x) d 0 Assume that xc( x) ax with a 0 a known positive constant e. Assume that d is unknown but is bounded by d D with D a known positive constant. Prove that the system is UUB and find the bound on x(t). f. Assume that d is unknown but is bounded by d D x with D a known positive constant. Prove that the system is UUB and find the bound on x(t).. UUB Use Lyapunov to show that the system x x x ( x x 3) x x x ( x x 3) is uniformly ultimately bounded UUB. That is, show that the Lyapunov derivative is NEGATIVE OUTSIDE A BOUNDED REGION. Find the radius of the bounded region outside which V <0. Any states outside this region are attracted towards the origin. 3. Lyapunov Theorem for Control Design. A system is given by x xu 3 x x u a. Use Lyapunov Linearization Method to show that the open-loop system with u(t)= 0 is unstable about the origin. b. Select the nonlinear feedback control input u x. Find the closed-loop system. Use a Lyapunov extension to show that the nonlinear closed-loop system is UUB. That is, select the quadratic Lyapunov function and find V along the closed-loop system trajectories. Then show that V is negative outside a region (i.e. if x is large enough). If you cannot solve for x such that V ( x) is negative, then plot V ( x) using MATLAB and draw conclusions about stability. c. Discuss the stability. When is the Lyapunov derivative negative? Can you use a LaSalle Extension to show AS? 4. Use Lyapunov Equation on p. 8 to check the stability of the linear systems 0 g. x h. Ax x

23 Redo Exam and turn in as Hwk 6. Lyapunov Functions- limit cycle, UUB EE 533 Homework 6 Use Lyapunov functions to examine the stability of the following systems. Be clear and show all steps. a. Pick a suitable complicated Lyapunov function to study the stability of the limit cycle for the system x x( x 3x 3) 4 x x ( x 3x 3) b. Use the standard simple quadratic Lyapunov function to show that this system is UUB. Describe the region outside which the Lyapunov derivative is negative x x( x 3x 3) 4 6 x x x ( x 3x 3). Barbalat s Lemma and LaSalle extension e. Use quadratic Lyapunov Function to show this system is locally AS x x( x ) xx Find a ball within which V 0. This region is contained in the region of attraction. f. Use quadratic Lyapunov Function to show this system is locally SISL x x( x ) x Find a region within which V 0. g. Use Barbalat s Lemma to verify that V 0. Check uniform continuity of the Lyapunov derivative of the system in part b. h. Use LaSalle s extension to verify that the system in part b is actually AS. Find the equilibrium point. 3. Equilibrium points and linearization System is x ( x x ) x ( x x ) d. Find all equilibrium points 3

24 e. Find Jacobian f. Find the nature of all e.p.s 4. Lyapunov Equation, AS, SISL a. Use Lyapunov equation to show this system is AS 0 Ax x 8 6 Use Q=I>0. Is the solution P to the Lyapunov equation unique? b. Use Lyapunov equation to show this system is SISL 0 Ax x Use Q 0 0. Now, you have to find ANY positive definite P that solves the Lyapunov equation. Is the solution P unique for this case? 4

25 EE 533 Homework 7 i/o Feedback Linearization. A system is given by sin x 4 x cos x u with output y() t x (t) a. Design a FB linearization controller to make the output follow a desired trajectory yd ( t ) That is, find ut () b. Discuss the internal dynamics. Are they a problem?. NMP System A system is given by 4 y zy u z ( y ) z 5 z 5 0 a. Take output y( t ) and find the FB linearization controller ut ( ) to follow the prescribed trajectory yd ( t ). b. Find the internal dynamics. Set yt ( ) 0 to get the zero dynamics. Are the ZD stable? Does the FB linearization controller work? 3. Effect of Output Choice in i/o FB Linearization It is desired to stabilize a system given by x sin x x u x x e. Select the out as y x and use FB lin. design to select the control u(t) to follow the desired trajectory yd ( t ). Check the internal dynamics. Set y=0 to get the zero dynamics. Is the system minimum phase? f. Select the new output y x. Find the FB lin. controller u(t). Does this work? What about the internal dynamics? 5

26 Feedback Linearization, backstepping EE 533 Homework 8 Slotine and Li. I/O feedback linearization. Slotine and Li Problem 6.3. System is sin x 4 x cos x u y x Do i/o feedback linearization to make output track desired trajectory yd ( t ).. Backstepping. The system is sin x u 4 x cos x x Do backstepping to stabilize this system. Select the desired value x d to yield the first step dynamics of x 0. Compare this to Problem, which uses i/o FB linearization. 3. Backstepping. Slotine and Li Problem 6.. Globally stabilize means backstepping. 4. I/O fb linearization. It is desired to stabilize a system given by x sin x x u x x a. Select the out as y x and use FB lin. design to select the control to make output go to zero. Is the system minimum phase? b. Select the new output y x. Does this work? c. Design a backstepping controller. 6

27 The following homeworks are out of date. Do not do homeworks until they are assigned. They may change. 7

28 Nonlinear Systems and Equilibrium EE 533 Homework Fall 009. Obtain the linear model of the system described by y( y ) y y 0 around the equilibrium point.. The system of equations axbxx cx dx exx fx describes the growth of two competing species that prey on each other. The constants are positive parameters and it is assumed that the two states are positive. Determine the linear model of the system around the equilibrium point. Simulate the system using MATLAB for various initial conditions. Make phase plane plot. 3. Determine the equilibrium points and their nature for the predator-prey system xxx. x x x Simulate the system using MATLAB for various initial conditions. Make phase plane plot. 4. Determine the equilibrium points and their nature for the system y( y ). y x( yx ) 8

29 EE 533 Homework Fall 009 State Variable Systems, Computer Simulation, Linearization. Simulate the van der Pol oscillator y" ( y ) y' y 0 using MATLAB for various ICs. Plot y(t) vs. t and also the phase plane plot y'(t) vs. y(t). Use y(0)=0., y'(0)= 0. c. For = d. For= Do MATLAB simulation of the Lorenz Attractor chaotic system. Run for 50 sec. with all initial states equal to 0.3. Plot states versus time, and also make 3-D plot of x, x, x3 using PLOT3(x,x,x3). ( x x ) 3 bx3 xx use = 0, r= 8, b= 8/3. rx x x x 3 3. A system has transfer function H ( s) s s 4 4s 3 a. Use MATLAB to make a 3-D plot of the magnitude of H(s) b. Use MATLAB to make a 3-D plot of the phase of H(s) c. Use MATLAB to draw magnitude and phase Bode plots 4. Use separation of variables to verify the formula for x(t) in Slotine & Li ex.. on p Duffing s equation is interesting in that it exhibits bifurcation, or dependence of stability properties and number of equilibrium points on a parameter. The undamped Duffing equation is 3 x 0 g. Find the equilibrium points. Show that for 0 there is only one e.p. h. For 0 there are 3 e.p.s Linearize the system and study the nature of these 3 e.p.s i. Simulate the Duffing oscillator for. Make time plot and phase plane plot. 9

30 EE 533 Homework 3 Fall 009 Chaos, Phase Plane. A system that exhibits chaos is the logistic function xk ( k ) xk ( xk ) However, chaos only occurs for certain values of sweep through the k values using k k for fixed less than but close to. These two equations form a dynamical system. Perform a MATLAB simulation to reproduce this plot of xk vs k, which was taken for Interpret the plot with some discussion in terms of bifurcation theory. Plot also. Show your MATLAB code. It is indeed interesting that the logistic function appears in economic systems and military supply systems. k k. Rather than try all values of. For Slotine & Li Example. on P. 0- a. Find equilibrium points b. Linearize the system about each equilibrium point. Find poles in each case. c. Simulate the system to find the Region of Attraction. k, we can 3. Slotine and Li p. 39 problem.. Simulate and plot phase plane trajectories for various ICs. Do not do the problem requested in the book. 4. Slotine and Li p. 39 problem.3 Simulate using MATLAB using various initial conditions. Do not do the problem requested in the book. 30

31 EE 533 Homework 4 Fall 009 Slotine and Li Lyapunov s Method. Slotine and Li p. 97 problem 3... Use Lyapunov Equation on p. 8 to prove asymptotic stability of the system 0 x Use Lyapunov to show that the system xx x ( x x 3) x x x ( x x 3) is locally asymptotically stable. Find the Region of Asymptotic Stability 4. UUB Use Lyapunov to show that the system x x x ( x x 3) x x x ( x x 3) is uniformly ultimately bounded UUB. That is, show that the Lyapunov derivative is NEGATIVE OUTSIDE A BOUNDED REGION. Find the radius of the bounded region outside which V <0. Any states outside this region are attracted towards the origin. 5. Lyapunov Theorem for Control Design. A system is given by x xu 3 x x u a. Use Lyapunov Linearization Method to show that the open-loop system with u(t)= 0 is unstable about the origin. b. Select the nonlinear feedback control input u x. Find the closed-loop system. Use a Lyapunov extension to show that the closed-loop system is UUB. That is, select the quadratic Lyapunov function and find V along the closed-loop system trajectories. Then show that V is negative outside a region (i.e. if x is large enough). c. Discuss the stability. When is the Lyapunov derivative negative? Can you use a LaSalle Extension to show AS? 3

32 EE 533 Homework 5 Fall 009, Slotine and Li Lyapunov. S&L p. 03, Example 4.. a. For the 3 systems given, prove the stability claimed by verifying the 3 conditions given. b. Integrate the state equations to find the solutions x(t) of the three systems.. S&L p. 05, Example 4.3. Integrate the state equation to verify the solution given. 3. S&L p. 55 problem 4.9, parts a and b. 4. Consider the nonlinear dynamics for an m-link robot manipulator, M qq Cq, q q Dq gq, m where q, R. M q accounts for the robot inertia. Cq, q q accounts for centrifugal and Coriolis forces. Dq accounts for viscous damping. g q accounts for gravity forces. In addition, we have the following properties: i. M q is a symmetric positive definite matrix of q. ii. M C is a skew symmetric matrix of q, q. iii. g qw q q, where W q is a positive definite function of q. Show that for D 0, the map from to q is passive lossless. And when D is positive definite, the map from to q is passive dissipative. T Hint: Use the total energy V q M qq Pq as a storage function. Select an appropriate Pq, a positive definite function of q. 3

33 Feedback Linearization, backstepping EE 533 Homework 6 Fall 009, Slotine and Li. I/O feedback linearization. Slotine and Li Problem 6.3. Backstepping. The system is sin x u 4 x cos x x Do backstepping to stabilize this system. Select the desired value x 0 to yield the first step dynamics of x 0. Compare this to Problem 6.3, which uses i/o FB linearization. 3. Backstepping. Slotine and Li Problem 6.. Globally stabilize means backstepping. 4. I/O linearization. It is desired to stabilize a system given by sin x x x x x u d. Select the out as y x and use FB lin. design to select the control. Is the system minimum phase? e. Select the new output y x. Does this work. f. Design a backstepping controller. 5. Input-State Linearization. Slotine and Li Problem 6.7. a. Write f ( x ), gx, ( ) ad f g. b. Is the system input-state linearizable. c. Check the given z ( ) x. Does it work? 33

34 Bifurcations EE 533 Homework 7 Fall 009, Strogatz book x. Plot with MATLABthe following vector fields as 3-D surfaces in the (x,r) plane. Also plot the bifurcation diagrams in the plane. r a. b. c. d. e. f. f( x, r) rx f( x, r) rx f( x, r) rx f( x, r) rx f( x, r) rx f( x, r) rx x

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 5 Homewors Fall 04 Updated: Sunday, October 6, 04 Some homewor assignments refer to the textboos: Slotine and Li, etc. For full credit, show all wor. Some problems require hand calculations. In those

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 5322 Intelligent Control Systems Homewor for Spring 209 Updated: Tuesday, March 26, 209 DO NOT DO HOMEWORK UNTIL IT IS ASSIGNED. THE ASSIGNMENTS MAY CHANGE UNTIL ANNOUNCED For full credit, show all

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

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 537 Homewors Friedland Text Updated: Wednesday November 8 Some homewor assignments refer to Friedland s text For full credit show all wor. Some problems require hand calculations. In those cases do

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 Lyapunov Stability - I Hanz Richter Mechanical Engineering Department Cleveland State University Definition of Stability - Lyapunov Sense Lyapunov

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

Automatic Control 2. Nonlinear systems. Prof. Alberto Bemporad. University of Trento. Academic year

Automatic Control 2. Nonlinear systems. Prof. Alberto Bemporad. University of Trento. Academic year Automatic Control 2 Nonlinear systems Prof. Alberto Bemporad University of Trento Academic year 2010-2011 Prof. Alberto Bemporad (University of Trento) Automatic Control 2 Academic year 2010-2011 1 / 18

More information

Math 266: Phase Plane Portrait

Math 266: Phase Plane Portrait Math 266: Phase Plane Portrait Long Jin Purdue, Spring 2018 Review: Phase line for an autonomous equation For a single autonomous equation y = f (y) we used a phase line to illustrate the equilibrium solutions

More information

Nonlinear System Analysis

Nonlinear System Analysis Nonlinear System Analysis Lyapunov Based Approach Lecture 4 Module 1 Dr. Laxmidhar Behera Department of Electrical Engineering, Indian Institute of Technology, Kanpur. January 4, 2003 Intelligent Control

More information

CDS 101/110a: Lecture 2.1 Dynamic Behavior

CDS 101/110a: Lecture 2.1 Dynamic Behavior CDS 11/11a: Lecture 2.1 Dynamic Behavior Richard M. Murray 6 October 28 Goals: Learn to use phase portraits to visualize behavior of dynamical systems Understand different types of stability for an equilibrium

More information

Practice Problems for Final Exam

Practice Problems for Final Exam Math 1280 Spring 2016 Practice Problems for Final Exam Part 2 (Sections 6.6, 6.7, 6.8, and chapter 7) S o l u t i o n s 1. Show that the given system has a nonlinear center at the origin. ẋ = 9y 5y 5,

More information

CDS 101/110a: Lecture 2.1 Dynamic Behavior

CDS 101/110a: Lecture 2.1 Dynamic Behavior CDS 11/11a: Lecture.1 Dynamic Behavior Richard M. Murray 6 October 8 Goals: Learn to use phase portraits to visualize behavior of dynamical systems Understand different types of stability for an equilibrium

More information

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

EN Nonlinear Control and Planning in Robotics Lecture 3: Stability February 4, 2015 EN530.678 Nonlinear Control and Planning in Robotics Lecture 3: Stability February 4, 2015 Prof: Marin Kobilarov 0.1 Model prerequisites Consider ẋ = f(t, x). We will make the following basic assumptions

More information

LMI Methods in Optimal and Robust Control

LMI Methods in Optimal and Robust Control LMI Methods in Optimal and Robust Control Matthew M. Peet Arizona State University Lecture 15: Nonlinear Systems and Lyapunov Functions Overview Our next goal is to extend LMI s and optimization to nonlinear

More information

Video 8.1 Vijay Kumar. Property of University of Pennsylvania, Vijay Kumar

Video 8.1 Vijay Kumar. Property of University of Pennsylvania, Vijay Kumar Video 8.1 Vijay Kumar 1 Definitions State State equations Equilibrium 2 Stability Stable Unstable Neutrally (Critically) Stable 3 Stability Translate the origin to x e x(t) =0 is stable (Lyapunov stable)

More information

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

Output Feedback and State Feedback. EL2620 Nonlinear Control. Nonlinear Observers. Nonlinear Controllers. ẋ = f(x,u), y = h(x) Output Feedback and State Feedback EL2620 Nonlinear Control Lecture 10 Exact feedback linearization Input-output linearization Lyapunov-based control design methods ẋ = f(x,u) y = h(x) Output feedback:

More information

1 The pendulum equation

1 The pendulum equation Math 270 Honors ODE I Fall, 2008 Class notes # 5 A longer than usual homework assignment is at the end. The pendulum equation We now come to a particularly important example, the equation for an oscillating

More information

EECS 700: Exam # 1 Tuesday, October 21, 2014

EECS 700: Exam # 1 Tuesday, October 21, 2014 EECS 700: Exam # 1 Tuesday, October 21, 2014 Print Name and Signature The rules for this exam are as follows: Write your name on the front page of the exam booklet. Initial each of the remaining pages

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

3 Stability and Lyapunov Functions

3 Stability and Lyapunov Functions CDS140a Nonlinear Systems: Local Theory 02/01/2011 3 Stability and Lyapunov Functions 3.1 Lyapunov Stability Denition: An equilibrium point x 0 of (1) is stable if for all ɛ > 0, there exists a δ > 0 such

More information

Nonlinear Control. Nonlinear Control Lecture # 25 State Feedback Stabilization

Nonlinear Control. Nonlinear Control Lecture # 25 State Feedback Stabilization Nonlinear Control Lecture # 25 State Feedback Stabilization Backstepping η = f a (η)+g a (η)ξ ξ = f b (η,ξ)+g b (η,ξ)u, g b 0, η R n, ξ, u R Stabilize the origin using state feedback View ξ as virtual

More information

Controls Problems for Qualifying Exam - Spring 2014

Controls Problems for Qualifying Exam - Spring 2014 Controls Problems for Qualifying Exam - Spring 2014 Problem 1 Consider the system block diagram given in Figure 1. Find the overall transfer function T(s) = C(s)/R(s). Note that this transfer function

More information

Solutions to homework assignment #7 Math 119B UC Davis, Spring for 1 r 4. Furthermore, the derivative of the logistic map is. L r(x) = r(1 2x).

Solutions to homework assignment #7 Math 119B UC Davis, Spring for 1 r 4. Furthermore, the derivative of the logistic map is. L r(x) = r(1 2x). Solutions to homework assignment #7 Math 9B UC Davis, Spring 0. A fixed point x of an interval map T is called superstable if T (x ) = 0. Find the value of 0 < r 4 for which the logistic map L r has a

More information

EML5311 Lyapunov Stability & Robust Control Design

EML5311 Lyapunov Stability & Robust Control Design EML5311 Lyapunov Stability & Robust Control Design 1 Lyapunov Stability criterion In Robust control design of nonlinear uncertain systems, stability theory plays an important role in engineering systems.

More information

STABILITY. Phase portraits and local stability

STABILITY. Phase portraits and local stability MAS271 Methods for differential equations Dr. R. Jain STABILITY Phase portraits and local stability We are interested in system of ordinary differential equations of the form ẋ = f(x, y), ẏ = g(x, y),

More information

Section 5. Graphing Systems

Section 5. Graphing Systems Section 5. Graphing Systems 5A. The Phase Plane 5A-1. Find the critical points of each of the following non-linear autonomous systems. x = x 2 y 2 x = 1 x + y a) b) y = x xy y = y + 2x 2 5A-2. Write each

More information

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

E209A: Analysis and Control of Nonlinear Systems Problem Set 6 Solutions E9A: Analysis and Control of Nonlinear Systems Problem Set 6 Solutions Michael Vitus Gabe Hoffmann Stanford University Winter 7 Problem 1 The governing equations are: ẋ 1 = x 1 + x 1 x ẋ = x + x 3 Using

More information

ME 680- Spring Representation and Stability Concepts

ME 680- Spring Representation and Stability Concepts ME 680- Spring 014 Representation and Stability Concepts 1 3. Representation and stability concepts 3.1 Continuous time systems: Consider systems of the form x F(x), x n (1) where F : U Vis a mapping U,V

More information

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

Prof. Krstic Nonlinear Systems MAE281A Homework set 1 Linearization & phase portrait 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

More information

PROJECTS on. Multiple Bifurcations of Sample Dynamical Systems

PROJECTS on. Multiple Bifurcations of Sample Dynamical Systems Course on Bifurcation Theory, a.y. 2010/11 PROJECTS on Multiple Bifurcations of Sample Dynamical Systems Students of the Bifurcation Theory course can carry out an individual homework, to be discussed

More information

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

Half of Final Exam Name: Practice Problems October 28, 2014 Math 54. Treibergs Half of Final Exam Name: Practice Problems October 28, 24 Half of the final will be over material since the last midterm exam, such as the practice problems given here. The other half

More information

Question: Total. Points:

Question: Total. Points: MATH 308 May 23, 2011 Final Exam Name: ID: Question: 1 2 3 4 5 6 7 8 9 Total Points: 0 20 20 20 20 20 20 20 20 160 Score: There are 9 problems on 9 pages in this exam (not counting the cover sheet). Make

More information

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

154 Chapter 9 Hints, Answers, and Solutions The particular trajectories are highlighted in the phase portraits below. 54 Chapter 9 Hints, Answers, and Solutions 9. The Phase Plane 9.. 4. The particular trajectories are highlighted in the phase portraits below... 3. 4. 9..5. Shown below is one possibility with x(t) and

More information

Converse Lyapunov theorem and Input-to-State Stability

Converse Lyapunov theorem and Input-to-State Stability Converse Lyapunov theorem and Input-to-State Stability April 6, 2014 1 Converse Lyapunov theorem In the previous lecture, we have discussed few examples of nonlinear control systems and stability concepts

More information

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

MCE/EEC 647/747: Robot Dynamics and Control. Lecture 8: Basic Lyapunov Stability Theory MCE/EEC 647/747: Robot Dynamics and Control Lecture 8: Basic Lyapunov Stability Theory Reading: SHV Appendix Mechanical Engineering Hanz Richter, PhD MCE503 p.1/17 Stability in the sense of Lyapunov A

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

System Control Engineering 0

System Control Engineering 0 System Control Engineering 0 Koichi Hashimoto Graduate School of Information Sciences Text: Nonlinear Control Systems Analysis and Design, Wiley Author: Horacio J. Marquez Web: http://www.ic.is.tohoku.ac.jp/~koichi/system_control/

More information

ADAPTIVE DESIGN OF CONTROLLER AND SYNCHRONIZER FOR LU-XIAO CHAOTIC SYSTEM

ADAPTIVE DESIGN OF CONTROLLER AND SYNCHRONIZER FOR LU-XIAO CHAOTIC SYSTEM ADAPTIVE DESIGN OF CONTROLLER AND SYNCHRONIZER FOR LU-XIAO CHAOTIC SYSTEM WITH UNKNOWN PARAMETERS Sundarapandian Vaidyanathan 1 1 Research and Development Centre, Vel Tech Dr. RR & Dr. SR Technical University

More information

COSMOS: Making Robots and Making Robots Intelligent Lecture 3: Introduction to discrete-time dynamics

COSMOS: Making Robots and Making Robots Intelligent Lecture 3: Introduction to discrete-time dynamics COSMOS: Making Robots and Making Robots Intelligent Lecture 3: Introduction to discrete-time dynamics Jorge Cortés and William B. Dunbar June 3, 25 Abstract In this and the coming lecture, we will introduce

More information

Math Exam 03 Review

Math Exam 03 Review Math 10350 Exam 03 Review 1. The statement: f(x) is increasing on a < x < b. is the same as: 1a. f (x) is on a < x < b. 2. The statement: f (x) is negative on a < x < b. is the same as: 2a. f(x) is on

More information

ADAPTIVE SYNCHRONIZATION FOR RÖSSLER AND CHUA S CIRCUIT SYSTEMS

ADAPTIVE SYNCHRONIZATION FOR RÖSSLER AND CHUA S CIRCUIT SYSTEMS Letters International Journal of Bifurcation and Chaos, Vol. 12, No. 7 (2002) 1579 1597 c World Scientific Publishing Company ADAPTIVE SYNCHRONIZATION FOR RÖSSLER AND CHUA S CIRCUIT SYSTEMS A. S. HEGAZI,H.N.AGIZA

More information

MATH 415, WEEKS 14 & 15: 1 Recurrence Relations / Difference Equations

MATH 415, WEEKS 14 & 15: 1 Recurrence Relations / Difference Equations MATH 415, WEEKS 14 & 15: Recurrence Relations / Difference Equations 1 Recurrence Relations / Difference Equations In many applications, the systems are updated in discrete jumps rather than continuous

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

EE222 - Spring 16 - Lecture 2 Notes 1

EE222 - Spring 16 - Lecture 2 Notes 1 EE222 - Spring 16 - Lecture 2 Notes 1 Murat Arcak January 21 2016 1 Licensed under a Creative Commons Attribution-NonCommercial-ShareAlike 4.0 International License. Essentially Nonlinear Phenomena Continued

More information

28. Pendulum phase portrait Draw the phase portrait for the pendulum (supported by an inextensible rod)

28. Pendulum phase portrait Draw the phase portrait for the pendulum (supported by an inextensible rod) 28. Pendulum phase portrait Draw the phase portrait for the pendulum (supported by an inextensible rod) θ + ω 2 sin θ = 0. Indicate the stable equilibrium points as well as the unstable equilibrium points.

More information

- - - - - - - - - - - - - - - - - - DISCLAIMER - - - - - - - - - - - - - - - - - - General Information: This midterm is a sample midterm. This means: The sample midterm contains problems that are of similar,

More information

EECS C128/ ME C134 Final Wed. Dec. 15, am. Closed book. Two pages of formula sheets. No calculators.

EECS C128/ ME C134 Final Wed. Dec. 15, am. Closed book. Two pages of formula sheets. No calculators. Name: SID: EECS C28/ ME C34 Final Wed. Dec. 5, 2 8- am Closed book. Two pages of formula sheets. No calculators. There are 8 problems worth points total. Problem Points Score 2 2 6 3 4 4 5 6 6 7 8 2 Total

More information

Example Chaotic Maps (that you can analyze)

Example Chaotic Maps (that you can analyze) Example Chaotic Maps (that you can analyze) Reading for this lecture: NDAC, Sections.5-.7. Lecture 7: Natural Computation & Self-Organization, Physics 256A (Winter 24); Jim Crutchfield Monday, January

More information

LECTURE 8: DYNAMICAL SYSTEMS 7

LECTURE 8: DYNAMICAL SYSTEMS 7 15-382 COLLECTIVE INTELLIGENCE S18 LECTURE 8: DYNAMICAL SYSTEMS 7 INSTRUCTOR: GIANNI A. DI CARO GEOMETRIES IN THE PHASE SPACE Damped pendulum One cp in the region between two separatrix Separatrix Basin

More information

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.

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. Entrance Exam, Differential Equations April, 7 (Solve exactly 6 out of the 8 problems). Consider the following initial value problem: { y + y + y cos(x y) =, y() = y. Find all the values y such that the

More information

Adaptive Robust Tracking Control of Robot Manipulators in the Task-space under Uncertainties

Adaptive Robust Tracking Control of Robot Manipulators in the Task-space under Uncertainties Australian Journal of Basic and Applied Sciences, 3(1): 308-322, 2009 ISSN 1991-8178 Adaptive Robust Tracking Control of Robot Manipulators in the Task-space under Uncertainties M.R.Soltanpour, M.M.Fateh

More information

Introduction to Dynamical Systems Basic Concepts of Dynamics

Introduction to Dynamical Systems Basic Concepts of Dynamics Introduction to Dynamical Systems Basic Concepts of Dynamics A dynamical system: Has a notion of state, which contains all the information upon which the dynamical system acts. A simple set of deterministic

More information

MA 262, Fall 2017, Final Version 01(Green)

MA 262, Fall 2017, Final Version 01(Green) INSTRUCTIONS MA 262, Fall 2017, Final Version 01(Green) (1) Switch off your phone upon entering the exam room. (2) Do not open the exam booklet until you are instructed to do so. (3) Before you open the

More information

CHALMERS, GÖTEBORGS UNIVERSITET. EXAM for DYNAMICAL SYSTEMS. COURSE CODES: TIF 155, FIM770GU, PhD

CHALMERS, GÖTEBORGS UNIVERSITET. EXAM for DYNAMICAL SYSTEMS. COURSE CODES: TIF 155, FIM770GU, PhD CHALMERS, GÖTEBORGS UNIVERSITET EXAM for DYNAMICAL SYSTEMS COURSE CODES: TIF 155, FIM770GU, PhD Time: Place: Teachers: Allowed material: Not allowed: August 22, 2018, at 08 30 12 30 Johanneberg Jan Meibohm,

More information

Math 341 Fall 2006 Final Exam December 12th, Name:

Math 341 Fall 2006 Final Exam December 12th, Name: Math 341 Fall 2006 Final Exam December 12th, 2006 Name: You may use a calculator, your note card and something to write with. You must attach your notecard to the exam when you turn it in. You cannot use

More information

Prüfung Regelungstechnik I (Control Systems I) Übersetzungshilfe / Translation aid (English) To be returned at the end of the exam!

Prüfung Regelungstechnik I (Control Systems I) Übersetzungshilfe / Translation aid (English) To be returned at the end of the exam! Prüfung Regelungstechnik I (Control Systems I) Prof. Dr. Lino Guzzella 29. 8. 2 Übersetzungshilfe / Translation aid (English) To be returned at the end of the exam! Do not mark up this translation aid

More information

Differential Equations Handout A

Differential Equations Handout A Differential Equations Handout A Math 1b November 11, 2005 These are problems for the differential equations unit of the course, some topics of which are not covered in the main textbook. Some of these

More information

Dynamical Systems & Lyapunov Stability

Dynamical Systems & Lyapunov Stability Dynamical Systems & Lyapunov Stability Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University Outline Ordinary Differential Equations Existence & uniqueness Continuous dependence

More information

Analysis of Dynamical Systems

Analysis of Dynamical Systems 2 YFX1520 Nonlinear phase portrait Mathematical Modelling and Nonlinear Dynamics Coursework Assignments 1 ϕ (t) 0-1 -2-6 -4-2 0 2 4 6 ϕ(t) Dmitri Kartofelev, PhD 2018 Variant 1 Part 1: Liénard type equation

More information

HOMEWORK 3: Phase portraits for Mathmatical Models, Analysis and Simulation, Fall 2010 Report due Mon Oct 4, Maximum score 7.0 pts.

HOMEWORK 3: Phase portraits for Mathmatical Models, Analysis and Simulation, Fall 2010 Report due Mon Oct 4, Maximum score 7.0 pts. HOMEWORK 3: Phase portraits for Mathmatical Models, Analysis and Simulation, Fall 2010 Report due Mon Oct 4, 2010. Maximum score 7.0 pts. Three problems are to be solved in this homework assignment. The

More information

Lyapunov Design for Controls

Lyapunov Design for Controls F.L. Lewis, NAI Moncrief-O Donnell Chair, UTA Research Institute (UTARI) The University of Texas at Arlington, USA and Foreign Professor, Chongqing University, China Supported by : China Qian Ren Program,

More information

Chaos. Lendert Gelens. KU Leuven - Vrije Universiteit Brussel Nonlinear dynamics course - VUB

Chaos. Lendert Gelens. KU Leuven - Vrije Universiteit Brussel   Nonlinear dynamics course - VUB Chaos Lendert Gelens KU Leuven - Vrije Universiteit Brussel www.gelenslab.org Nonlinear dynamics course - VUB Examples of chaotic systems: the double pendulum? θ 1 θ θ 2 Examples of chaotic systems: the

More information

2D-Volterra-Lotka Modeling For 2 Species

2D-Volterra-Lotka Modeling For 2 Species Majalat Al-Ulum Al-Insaniya wat - Tatbiqiya 2D-Volterra-Lotka Modeling For 2 Species Alhashmi Darah 1 University of Almergeb Department of Mathematics Faculty of Science Zliten Libya. Abstract The purpose

More information

Why are Discrete Maps Sufficient?

Why are Discrete Maps Sufficient? Why are Discrete Maps Sufficient? Why do dynamical systems specialists study maps of the form x n+ 1 = f ( xn), (time is discrete) when much of the world around us evolves continuously, and is thus well

More information

APPM 2460 CHAOTIC DYNAMICS

APPM 2460 CHAOTIC DYNAMICS APPM 2460 CHAOTIC DYNAMICS 1. Introduction Today we ll continue our exploration of dynamical systems, focusing in particular upon systems who exhibit a type of strange behavior known as chaos. We will

More information

Fundamentals of Dynamical Systems / Discrete-Time Models. Dr. Dylan McNamara people.uncw.edu/ mcnamarad

Fundamentals of Dynamical Systems / Discrete-Time Models. Dr. Dylan McNamara people.uncw.edu/ mcnamarad Fundamentals of Dynamical Systems / Discrete-Time Models Dr. Dylan McNamara people.uncw.edu/ mcnamarad Dynamical systems theory Considers how systems autonomously change along time Ranges from Newtonian

More information

Lecture 9 Nonlinear Control Design

Lecture 9 Nonlinear Control Design Lecture 9 Nonlinear Control Design Exact-linearization Lyapunov-based design Lab 2 Adaptive control Sliding modes control Literature: [Khalil, ch.s 13, 14.1,14.2] and [Glad-Ljung,ch.17] Course Outline

More information

Dynamical Systems and Chaos Part I: Theoretical Techniques. Lecture 4: Discrete systems + Chaos. Ilya Potapov Mathematics Department, TUT Room TD325

Dynamical Systems and Chaos Part I: Theoretical Techniques. Lecture 4: Discrete systems + Chaos. Ilya Potapov Mathematics Department, TUT Room TD325 Dynamical Systems and Chaos Part I: Theoretical Techniques Lecture 4: Discrete systems + Chaos Ilya Potapov Mathematics Department, TUT Room TD325 Discrete maps x n+1 = f(x n ) Discrete time steps. x 0

More information

6.2 Brief review of fundamental concepts about chaotic systems

6.2 Brief review of fundamental concepts about chaotic systems 6.2 Brief review of fundamental concepts about chaotic systems Lorenz (1963) introduced a 3-variable model that is a prototypical example of chaos theory. These equations were derived as a simplification

More information

Lecture - 11 Bendixson and Poincare Bendixson Criteria Van der Pol Oscillator

Lecture - 11 Bendixson and Poincare Bendixson Criteria Van der Pol Oscillator Nonlinear Dynamical Systems Prof. Madhu. N. Belur and Prof. Harish. K. Pillai Department of Electrical Engineering Indian Institute of Technology, Bombay Lecture - 11 Bendixson and Poincare Bendixson Criteria

More information

EMA 545 Final Exam - Prof. M. S. Allen Spring 2011

EMA 545 Final Exam - Prof. M. S. Allen Spring 2011 EMA 545 Final Exam - Prof. M. S. Allen Spring 2011 Honor Pledge: On my honor, I pledge that this exam represents my own work, and that I have neither given nor received inappropriate aid in the preparation

More information

Handout 2: Invariant Sets and Stability

Handout 2: Invariant Sets and Stability Engineering Tripos Part IIB Nonlinear Systems and Control Module 4F2 1 Invariant Sets Handout 2: Invariant Sets and Stability Consider again the autonomous dynamical system ẋ = f(x), x() = x (1) with state

More information

ASTATISM IN NONLINEAR CONTROL SYSTEMS WITH APPLICATION TO ROBOTICS

ASTATISM IN NONLINEAR CONTROL SYSTEMS WITH APPLICATION TO ROBOTICS dx dt DIFFERENTIAL EQUATIONS AND CONTROL PROCESSES N 1, 1997 Electronic Journal, reg. N P23275 at 07.03.97 http://www.neva.ru/journal e-mail: diff@osipenko.stu.neva.ru Control problems in nonlinear systems

More information

Solution of Additional Exercises for Chapter 4

Solution of Additional Exercises for Chapter 4 1 1. (1) Try V (x) = 1 (x 1 + x ). Solution of Additional Exercises for Chapter 4 V (x) = x 1 ( x 1 + x ) x = x 1 x + x 1 x In the neighborhood of the origin, the term (x 1 + x ) dominates. Hence, the

More information

A plane autonomous system is a pair of simultaneous first-order differential equations,

A plane autonomous system is a pair of simultaneous first-order differential equations, Chapter 11 Phase-Plane Techniques 11.1 Plane Autonomous Systems A plane autonomous system is a pair of simultaneous first-order differential equations, ẋ = f(x, y), ẏ = g(x, y). This system has an equilibrium

More information

Math 215/255 Final Exam (Dec 2005)

Math 215/255 Final Exam (Dec 2005) Exam (Dec 2005) Last Student #: First name: Signature: Circle your section #: Burggraf=0, Peterson=02, Khadra=03, Burghelea=04, Li=05 I have read and understood the instructions below: Please sign: Instructions:.

More information

APPM 2360: Final Exam 10:30am 1:00pm, May 6, 2015.

APPM 2360: Final Exam 10:30am 1:00pm, May 6, 2015. APPM 23: Final Exam :3am :pm, May, 25. ON THE FRONT OF YOUR BLUEBOOK write: ) your name, 2) your student ID number, 3) lecture section, 4) your instructor s name, and 5) a grading table for eight questions.

More information

Predictability: Does the Flap of a Butterfly s Wings in Brazil set off a Tornado

Predictability: Does the Flap of a Butterfly s Wings in Brazil set off a Tornado Chapter 4 Dynamic Behavior Predictability: Does the Flap of a Butterfly s Wings in Brazil set off a Tornado in Texas? Talk given by Edward Lorenz, December 972 meeting of the American Association for the

More information

Control of Mobile Robots

Control of Mobile Robots Control of Mobile Robots Regulation and trajectory tracking Prof. Luca Bascetta (luca.bascetta@polimi.it) Politecnico di Milano Dipartimento di Elettronica, Informazione e Bioingegneria Organization and

More information

Announcements. Topics: Homework: - sections , 6.1 (extreme values) * Read these sections and study solved examples in your textbook!

Announcements. Topics: Homework: - sections , 6.1 (extreme values) * Read these sections and study solved examples in your textbook! Announcements Topics: - sections 5.2 5.7, 6.1 (extreme values) * Read these sections and study solved examples in your textbook! Homework: - review lecture notes thoroughly - work on practice problems

More information

Edward Lorenz. Professor of Meteorology at the Massachusetts Institute of Technology

Edward Lorenz. Professor of Meteorology at the Massachusetts Institute of Technology The Lorenz system Edward Lorenz Professor of Meteorology at the Massachusetts Institute of Technology In 1963 derived a three dimensional system in efforts to model long range predictions for the weather

More information

Een vlinder in de wiskunde: over chaos en structuur

Een vlinder in de wiskunde: over chaos en structuur Een vlinder in de wiskunde: over chaos en structuur Bernard J. Geurts Enschede, November 10, 2016 Tuin der Lusten (Garden of Earthly Delights) In all chaos there is a cosmos, in all disorder a secret

More information

Nonlinear Dynamics. Moreno Marzolla Dip. di Informatica Scienza e Ingegneria (DISI) Università di Bologna.

Nonlinear Dynamics. Moreno Marzolla Dip. di Informatica Scienza e Ingegneria (DISI) Università di Bologna. Nonlinear Dynamics Moreno Marzolla Dip. di Informatica Scienza e Ingegneria (DISI) Università di Bologna http://www.moreno.marzolla.name/ 2 Introduction: Dynamics of Simple Maps 3 Dynamical systems A dynamical

More information

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

2 Lyapunov Stability. x(0) x 0 < δ x(t) x 0 < ɛ 1 2 Lyapunov Stability Whereas I/O stability is concerned with the effect of inputs on outputs, Lyapunov stability deals with unforced systems: ẋ = f(x, t) (1) where x R n, t R +, and f : R n R + R n.

More information

Synchronization of different chaotic systems and electronic circuit analysis

Synchronization of different chaotic systems and electronic circuit analysis Synchronization of different chaotic systems and electronic circuit analysis J.. Park, T.. Lee,.. Ji,.. Jung, S.M. Lee epartment of lectrical ngineering, eungnam University, Kyongsan, Republic of Korea.

More information

More Details Fixed point of mapping is point that maps into itself, i.e., x n+1 = x n.

More Details Fixed point of mapping is point that maps into itself, i.e., x n+1 = x n. More Details Fixed point of mapping is point that maps into itself, i.e., x n+1 = x n. If there are points which, after many iterations of map then fixed point called an attractor. fixed point, If λ

More information

EG4321/EG7040. Nonlinear Control. Dr. Matt Turner

EG4321/EG7040. Nonlinear Control. Dr. Matt Turner EG4321/EG7040 Nonlinear Control Dr. Matt Turner EG4321/EG7040 [An introduction to] Nonlinear Control Dr. Matt Turner EG4321/EG7040 [An introduction to] Nonlinear [System Analysis] and Control Dr. Matt

More information

Math 122 Test 3. April 15, 2014

Math 122 Test 3. April 15, 2014 SI: Math 1 Test 3 April 15, 014 EF: 1 3 4 5 6 7 8 Total Name Directions: 1. No books, notes or 6 year olds with ear infections. You may use a calculator to do routine arithmetic computations. You may not

More information

ADAPTIVE FEEDBACK LINEARIZING CONTROL OF CHUA S CIRCUIT

ADAPTIVE FEEDBACK LINEARIZING CONTROL OF CHUA S CIRCUIT International Journal of Bifurcation and Chaos, Vol. 12, No. 7 (2002) 1599 1604 c World Scientific Publishing Company ADAPTIVE FEEDBACK LINEARIZING CONTROL OF CHUA S CIRCUIT KEVIN BARONE and SAHJENDRA

More information

Problem set 7 Math 207A, Fall 2011 Solutions

Problem set 7 Math 207A, Fall 2011 Solutions Problem set 7 Math 207A, Fall 2011 s 1. Classify the equilibrium (x, y) = (0, 0) of the system x t = x, y t = y + x 2. Is the equilibrium hyperbolic? Find an equation for the trajectories in (x, y)- phase

More information

Unit Ten Summary Introduction to Dynamical Systems and Chaos

Unit Ten Summary Introduction to Dynamical Systems and Chaos Unit Ten Summary Introduction to Dynamical Systems Dynamical Systems A dynamical system is a system that evolves in time according to a well-defined, unchanging rule. The study of dynamical systems is

More information

GLOBAL CHAOS SYNCHRONIZATION OF UNCERTAIN SPROTT J AND K SYSTEMS BY ADAPTIVE CONTROL

GLOBAL CHAOS SYNCHRONIZATION OF UNCERTAIN SPROTT J AND K SYSTEMS BY ADAPTIVE CONTROL GLOBAL CHAOS SYNCHRONIZATION OF UNCERTAIN SPROTT J AND K SYSTEMS BY ADAPTIVE CONTROL Sundarapandian Vaidyanathan 1 1 Research and Development Centre, Vel Tech Dr. RR & Dr. SR Technical University Avadi,

More information

Stability of Nonlinear Systems An Introduction

Stability of Nonlinear Systems An Introduction Stability of Nonlinear Systems An Introduction Michael Baldea Department of Chemical Engineering The University of Texas at Austin April 3, 2012 The Concept of Stability Consider the generic nonlinear

More information

Problem Value Score Total 100/105

Problem Value Score Total 100/105 RULES This is a closed book, closed notes test. You are, however, allowed one piece of paper (front side only) for notes and definitions, but no sample problems. The top half is the same as from the first

More information

Chaos Suppression in Forced Van Der Pol Oscillator

Chaos Suppression in Forced Van Der Pol Oscillator International Journal of Computer Applications (975 8887) Volume 68 No., April Chaos Suppression in Forced Van Der Pol Oscillator Mchiri Mohamed Syscom laboratory, National School of Engineering of unis

More information

Nonlinear dynamics & chaos BECS

Nonlinear dynamics & chaos BECS Nonlinear dynamics & chaos BECS-114.7151 Phase portraits Focus: nonlinear systems in two dimensions General form of a vector field on the phase plane: Vector notation: Phase portraits Solution x(t) describes

More information

Chapter 13 Internal (Lyapunov) Stability 13.1 Introduction We have already seen some examples of both stable and unstable systems. The objective of th

Chapter 13 Internal (Lyapunov) Stability 13.1 Introduction We have already seen some examples of both stable and unstable systems. The objective of th Lectures on Dynamic Systems and Control Mohammed Dahleh Munther A. Dahleh George Verghese Department of Electrical Engineering and Computer Science Massachuasetts Institute of Technology 1 1 c Chapter

More information

Tracking Control of a Class of Differential Inclusion Systems via Sliding Mode Technique

Tracking Control of a Class of Differential Inclusion Systems via Sliding Mode Technique International Journal of Automation and Computing (3), June 24, 38-32 DOI: 7/s633-4-793-6 Tracking Control of a Class of Differential Inclusion Systems via Sliding Mode Technique Lei-Po Liu Zhu-Mu Fu Xiao-Na

More information

Stability and Robustness Analysis of Nonlinear Systems via Contraction Metrics and SOS Programming

Stability and Robustness Analysis of Nonlinear Systems via Contraction Metrics and SOS Programming arxiv:math/0603313v1 [math.oc 13 Mar 2006 Stability and Robustness Analysis of Nonlinear Systems via Contraction Metrics and SOS Programming Erin M. Aylward 1 Pablo A. Parrilo 1 Jean-Jacques E. Slotine

More information

ADAPTIVE CHAOS SYNCHRONIZATION OF UNCERTAIN HYPERCHAOTIC LORENZ AND HYPERCHAOTIC LÜ SYSTEMS

ADAPTIVE CHAOS SYNCHRONIZATION OF UNCERTAIN HYPERCHAOTIC LORENZ AND HYPERCHAOTIC LÜ SYSTEMS ADAPTIVE CHAOS SYNCHRONIZATION OF UNCERTAIN HYPERCHAOTIC LORENZ AND HYPERCHAOTIC LÜ SYSTEMS Sundarapandian Vaidyanathan 1 1 Research and Development Centre, Vel Tech Dr. RR & Dr. SR Technical University

More information