Sample Solutions of Assignment 9 for MAT3270B

Similar documents
Sample Solutions of Assignment 10 for MAT3270B

Math 266: Phase Plane Portrait

Solutions to the Calculus and Linear Algebra problems on the Comprehensive Examination of January 28, 2011

Lecture 38. Almost Linear Systems

Math Ordinary Differential Equations

Nonlinear Autonomous Systems of Differential

Math 3301 Homework Set Points ( ) ( ) I ll leave it to you to verify that the eigenvalues and eigenvectors for this matrix are, ( ) ( ) ( ) ( )

1. Diagonalize the matrix A if possible, that is, find an invertible matrix P and a diagonal

Stability of Nonlinear Systems An Introduction

6 Linear Equation. 6.1 Equation with constant coefficients

FINAL EXAM MATH303 Theory of Ordinary Differential Equations. Spring dx dt = x + 3y dy dt = x y.

Old Math 330 Exams. David M. McClendon. Department of Mathematics Ferris State University

ENGI Duffing s Equation Page 4.65

ODE, part 2. Dynamical systems, differential equations

Section 5.4 (Systems of Linear Differential Equation); 9.5 Eigenvalues and Eigenvectors, cont d

Def. (a, b) is a critical point of the autonomous system. 1 Proper node (stable or unstable) 2 Improper node (stable or unstable)

Autonomous Systems and Stability

4 Second-Order Systems

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

Changing coordinates to adapt to a map of constant rank

STABILITY. Phase portraits and local stability

Chapter 6 Nonlinear Systems and Phenomena. Friday, November 2, 12

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

2.10 Saddles, Nodes, Foci and Centers

Module 6 : Solving Ordinary Differential Equations - Initial Value Problems (ODE-IVPs) Section 3 : Analytical Solutions of Linear ODE-IVPs

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

Preliminary Exam 2016 Solutions to Morning Exam

Stability of Dynamical systems

ENGI Linear Approximation (2) Page Linear Approximation to a System of Non-Linear ODEs (2)

MATH 5720: Unconstrained Optimization Hung Phan, UMass Lowell September 13, 2018

Problem set 7 Math 207A, Fall 2011 Solutions

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

Calculus for the Life Sciences II Assignment 6 solutions. f(x, y) = 3π 3 cos 2x + 2 sin 3y

Math 322. Spring 2015 Review Problems for Midterm 2

Problem set 6 Math 207A, Fall 2011 Solutions. 1. A two-dimensional gradient system has the form

Math 273 (51) - Final

Section 5. Graphing Systems

Sample Solutions of Assignment 3 for MAT3270B: 2.8,2.3,2.5,2.7

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

MATH 1700 FINAL SPRING MOON

Solutions to Math 53 First Exam April 20, 2010

Solutions to the Multivariable Calculus and Linear Algebra problems on the Comprehensive Examination of January 31, 2014

Hello everyone, Best, Josh

EC Control Engineering Quiz II IIT Madras

Properties of Linear Transformations from R n to R m

MATH 251 Final Examination December 19, 2012 FORM A. Name: Student Number: Section:

Department of Mathematics IIT Guwahati

MATH 4B Differential Equations, Fall 2016 Final Exam Study Guide

Solution: In standard form (i.e. y + P (t)y = Q(t)) we have y t y = cos(t)

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

Designing Information Devices and Systems II Spring 2018 J. Roychowdhury and M. Maharbiz Discussion 6B

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.

Math 1270 Honors ODE I Fall, 2008 Class notes # 14. x 0 = F (x; y) y 0 = G (x; y) u 0 = au + bv = cu + dv

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

Ordinary Differential Equations II

B. Differential Equations A differential equation is an equation of the form

Exercise 8.1 We have. the function is differentiable, with. f (x 0, y 0 )(u, v) = (2ax 0 + 2by 0 )u + (2bx 0 + 2cy 0 )v.

ECEEN 5448 Fall 2011 Homework #5 Solutions

ENGI 9420 Lecture Notes 4 - Stability Analysis Page Stability Analysis for Non-linear Ordinary Differential Equations

THE SEPARATRIX FOR A SECOND ORDER ORDINARY DIFFERENTIAL EQUATION OR A 2 2 SYSTEM OF FIRST ORDER ODE WHICH ALLOWS A PHASE PLANE QUANTITATIVE ANALYSIS

Lyapunov Stability Theory

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

Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games

Physics: spring-mass system, planet motion, pendulum. Biology: ecology problem, neural conduction, epidemics

n n. ( t) ( ) = = Ay ( ) a y

1. < 0: the eigenvalues are real and have opposite signs; the fixed point is a saddle point

Systems of Algebraic Equations and Systems of Differential Equations

Understand the existence and uniqueness theorems and what they tell you about solutions to initial value problems.

Practice Problems for Final Exam

Dynamical Systems Solutions to Exercises

6.3. Nonlinear Systems of Equations

Linear ODEs. Existence of solutions to linear IVPs. Resolvent matrix. Autonomous linear systems

Math 312 Lecture Notes Linear Two-dimensional Systems of Differential Equations

Note: Each problem is worth 14 points except numbers 5 and 6 which are 15 points. = 3 2

1 The relation between a second order linear ode and a system of two rst order linear odes

Lecture Notes for Math 251: ODE and PDE. Lecture 27: 7.8 Repeated Eigenvalues

Solutions to Dynamical Systems 2010 exam. Each question is worth 25 marks.

Dynamical Systems & Lyapunov Stability

CDS 101 Precourse Phase Plane Analysis and Stability

1 The pendulum equation

REVIEW NOTES FOR MATH 266

CDS 101/110a: Lecture 2.1 Dynamic Behavior

Math 312 Lecture Notes Linearization

Input-to-state stability and interconnected Systems

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

Some recent results on controllability of coupled parabolic systems: Towards a Kalman condition

Probability Density (1)

LECTURE 8: DYNAMICAL SYSTEMS 7

NAME: MA Sample Final Exam. Record all your answers on the answer sheet provided. The answer sheet is the only thing that will be graded.

Do not write below here. Question Score Question Score Question Score

MATH 215/255 Solutions to Additional Practice Problems April dy dt

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

MATH 251 Final Examination August 10, 2011 FORM A. Name: Student Number: Section:

B5.6 Nonlinear Systems

Math 4B Notes. Written by Victoria Kala SH 6432u Office Hours: T 12:45 1:45pm Last updated 7/24/2016

1 Existence of Travelling Wave Fronts for a Reaction-Diffusion Equation with Quadratic- Type Kinetics

Nonlinear Control Lecture 2:Phase Plane Analysis

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

Section 9.3 Phase Plane Portraits (for Planar Systems)

MAT 22B - Lecture Notes

Transcription:

Sample Solutions of Assignment 9 for MAT370B. For the following ODEs, find the eigenvalues and eigenvectors, and classify the critical point 0,0 type and determine whether it is stable, asymptotically stable, or unstable. a. b. c. d. 3 5 3 5 0 Answer: a. deta λi = 3 λ λ = λ λ + = 0 λ =, λ = The eigenvectors are X = and X = 0,0 is an unstable critical point of the system.

b. deta λi = λ 5 3 λ = λ + λ + = 0 λ = + i, λ = i The eigenvectors are X = 5 i and X = 5 i 0,0 is an asymptotically stable critical point of the system. c. deta λi = The eigenvectors are λ 5 0 λ = λ = 0 X = λ = 0 and 0 X = t + 0 5 0,0 is an unstable critical point of the system. d. deta λi = The eigenvectors are λ λ = λ + = 0 λ = i, λ = i X = i

3 and X = i 0,0 is a stable critical point of the system.. Transform the following problems into the form λ µ µ λ a. b. 3 4 3 Answer: a. deta λi = 3 λ 4 λ = λ λ + = 0 λ = 7 + i, λ = 7 i The problem can be transformed to b. deta λi = 7 7 λ 3 λ = λ λ + = 0

4 λ = 7 + i, λ = 7 i The problem can be transformed to 7 7 3. Determine the critical points for each of the following systems a. d = 3y, = y y b. d = y, = y 4 y 3y 4 c. d = y, = µ y, µ > 0 d. d = y + y, = + y y Answer: a. The critical points of the system are 0,0, 0, and 3,. 5 5 b. The critical points of the system are 0,0, 0, and,0. c. The critical points of the system is 0,0. d. The critical points of the system is 0,0. 4. Find out an equation for the form H, y = c satisfies for the trajectories: a. d = + y +, = y y b. d = y 3 4y, = y + 6y

5 Answer: a. We get d = y y + y + The ODE can be rewritten as y yd + + y + = 0 M, y y So this ODE is eact for M,y y that = N, y = ψ = M, y = y y Integrating above equation, we obtain = N,y. Thus there is a ψ, y such ψ, y = y y + fy Setting ψ y = N gives then we get Hence, f y = y fy = y + c y y y = c b. We get d = y + 6y y 3 4y The ODE can be rewritten as y 6yd + y 3 4y = 0 M, y y N, y = 4y 6 = 4y 6

6 So this ODE is eact for M,y y that = N,y. Thus there is a ψ, y such ψ = M, y = y 6y Integrating above equation, we obtain ψ, y = y 3 y + fy Setting ψ y = N gives then we get Hence, f y = 4y fy = y + c y 3 y y = c 5. Show that the trajectories of the nonlinear undamped pendulum equation are given by where = θ, y = dθ. d θ + g l sin θ = 0 g y cos + l = c Answer: Let then = θ, y = dθ d = dθ = y = d θ = g l sin θ = g l sin

7 We get Hence, d = g sin l y g l cos + y = c. 6. Prove that for the system dx = fx there is at most one trajectory passing through a given point X 0. Answer: Let φ t, φ t be two trajectories passing through X 0, i.e. φ t = X 0, φ t = X 0, then Set dφ = fφ, φ t = X 0 dφ = fφ, φ t = X 0 φ = φ t + t t, then By uniqueness, d φ = dφ f φ, φ t = X 0 φ t = φ t = φ t + t t Hence, φ and φ are the same trajectories. 7. Prove that if a trajectory starts at a noncritical point of the system dx = fx

8 the it can not reach a critical point in a finite length of time. Answer: Assume the contrary. we suppose that 0, y 0 is a critical point of the system, and the solution = φt, y = ψt satisfies φa = 0, ψa = y 0. On the other hand, = 0, y = y 0 is a solution of the given system satisfying the initial condition = 0, y = y 0 at t = a. By the uniqueness, φt = 0, ψt = y 0, t 0, hence, φ0 = 0, ψt = y 0. This contradict the assumption of the problem. 8. Assuming the trajectory corresponding to solution = φt, y = ψt, < t < +, of an autonomous system is closed, show that the solution must be periodic. Answer: Suppose that φt = φt = X 0, t t. By the same proof as in Problem 6, we get φt + t t = φt, t Let T = t t, the φt + T = φt, t Hence φt is periodic with period T. 9. Write the following spring-mass system as a system of equations by introducing = u, u = du m d u + cdu + ku = 0 Find out the critical point and analyze the stability of the critical point.

9 Answer: Let then = u, y = du d = du = y = d u = c m y k m 0,0 is a critical point of the system. deta λi = λ c λ m k m = mλ + cλ + k = 0 Case i c 4mk > 0 and k 0 In this case, there are two distinct real eigenvalues, say r, r. Since m, c, k > 0, then r, r < 0 0,0 is asymptotically stable critical point of the system. Case ii c 4mk = 0 and k 0 In this case, there are repeated eigenvalues, say r. Since m, c, k > 0, then r < 0 0,0 is asymptotically stable critical point of the system. Case iii c 4mk < 0 and k 0 In this case, there are two comple eigenvalues, say r = λ + µi, r = λ µi where λ = c. Since m, c > 0, 0,0 is asymptotically stable m critical point of the system if c 0. Similarly, 0,0 is stable critical point of the system if c = 0.

0 Case iv k = 0 In this case, r = 0, r = c. 0,0 is stable critical point of the m system. 0. Show that the system is almost linear and 0,0 is a stable critical point of the system d = y, = y y Answer: d = 0 0 Let = r cos θ, y = r sin θ then y y + y lim y r = lim r cos θ sin θ = 0 lim y = 0 r det A =, the system is almost linear. deta λi = λ 0 0 λ λ = = λ + = 0 0,0 is stable critical point of the system.. Show that the system is almost linear and 0,0 is an asymptotically stable critical point of the system d = 3 + y, = y3

Answer: Let = r cos θ, y = r sin θ then lim 3 + y r = lim r cos 3 θ + r cos θ sin θ = 0 lim y3 r Hence, the system is almost linear. = limr sin 3 θ = 0 As for the type of the critical point, we will give the solution on Problem 5 of Assignment 0.. Determine all real critical points of the following system of equations a. d = + y, = + y b. d = y, = y Answer: a The critical points of the system are 0,0 and,-. b The critical points of the system are, and,-. 3. Consider the following problem. Determine the eigenvalues and critical points, classify the type of the critical point and determine whether it is stable, asymptotically stable, or unstable. ε = ε Here ε is a real number.

Answer: Since det A = ε + > 0, 0,0 is the only critical point of the system. deta λi = Case A: ε = 0 ε λ ε λ = λ ελ + ε + = 0 λ = i, λ = i 0,0 is a stable center. Case B: ε > 0 λ = ε + i, λ = ε i 0,0 is an unstable spiral point. Case C: ε < 0 λ = ε + i, λ = ε i 0,0 is an asymptotical stable spiral point. 4. Consider the following Lienard equation where g0 = 0. d + cd + g = 0 a. Write it as a system of two first order equations by introducing y = d. b. Show that 0,0 is a critical point and that the system is almost linear in the neighborhood of 0,0.

c. Show that if c0 > 0, g 0 > 0 then the critical point is asymptotically stable, and that if c0 < 0 or g 0 < 0, then the critical point is unstable. 3 Answer: a.let then y = d = d From the original equation, we get Let = cy g. X = y We get X = 0 0 c X + 0 g Let b. Since g0 = 0, then 0,0 is the critical point of the system. where η 0, 0 where η / 0, 0 The system can written as X 0 = g 0 c0 c = c0 + η g = g 0 + η X + 0 η y η

4 Take G T = 0, η y η Then G 0,, y 0, 0. Hence the system is almost linear in the neighborhood of 0,0. c. deta λi = Hence Case A: c0 > 0, g 0 > 0 λ g 0 c0 λ λ = c0 + c 0 4g 0 λ = c0 c 0 4g 0 = λ + c0λ + g 0 = 0 If = c 0 4g > 0, then λ < λ < 0 and 0,0 is an asymptotical stable critical point. If = c 0 4g < 0, then λ = α + βi, λ = α βi where α = c0 < 0 and 0,0 is an asymptotical stable critical point. If = c 0 4g = 0, then λ = λ = c0 < 0 and 0,0 is an asymptotical stable critical point. Case B: c0 < 0 If = c 0 4g critical point. > 0, then λ > λ > 0 and 0,0 is an unstable If = c 0 4g < 0, then λ = α + βi, λ = α βi where α = c0 > 0 and 0,0 is an unstable critical point.

If = c 0 4g = 0, then λ = λ = c0 > 0 and 0,0 is unstable critical point. 5 Case C: g 0 < 0 In this case, λ > 0 > λ > 0 and 0,0 is an unstable critical point.