Phase plane method is an important graphical methods to deal with problems related to a second-order autonomous system.

Similar documents
Variational Equation or Continuous Dependence on Initial Condition or Trajectory Sensitivity & Floquet Theory & Poincaré Map

1973 AP Calculus BC: Section I

Numerical Simulation for the 2-D Heat Equation with Derivative Boundary Conditions

MAT3700. Tutorial Letter 201/2/2016. Mathematics III (Engineering) Semester 2. Department of Mathematical sciences MAT3700/201/2/2016

ENGG 1203 Tutorial. Difference Equations. Find the Pole(s) Finding Equations and Poles

Continous system: differential equations

Anouncements. Conjugate Gradients. Steepest Descent. Outline. Steepest Descent. Steepest Descent

Chapter 3 Linear Equations of Higher Order (Page # 144)

Partial Fraction Expansion

Part B: Transform Methods. Professor E. Ambikairajah UNSW, Australia

Boyce/DiPrima/Meade 11 th ed, Ch 4.1: Higher Order Linear ODEs: General Theory

Fourier Series: main points

(A) 1 (B) 1 + (sin 1) (C) 1 (sin 1) (D) (sin 1) 1 (C) and g be the inverse of f. Then the value of g'(0) is. (C) a. dx (a > 0) is

k of the incident wave) will be greater t is too small to satisfy the required kinematics boundary condition, (19)

( A) ( B) ( C) ( D) ( E)

Note 6 Frequency Response

EEE 303: Signals and Linear Systems

Supplementary Information

Quality Monitoring Calibration Assuring Standardization Among Monitors

ECSE Partial fraction expansion (m<n) 3 types of poles Simple Real poles Real Equal poles

Phys Nov. 3, 2017 Today s Topics. Continue Chapter 2: Electromagnetic Theory, Photons, and Light Reading for Next Time

An Asymptotic Expansion for the Non-Central Chi-square Distribution. By Jinan Hamzah Farhood Department of Mathematics College of Education

Lecture 14. Time Harmonic Fields

1. Mathematical tools which make your life much simpler 1.1. Useful approximation formula using a natural logarithm

Linear Systems Analysis in the Time Domain

GUIDE FOR SUPERVISORS 1. This event runs most efficiently with two to four extra volunteers to help proctor students and grade the student

ODEs II, Supplement to Lectures 6 & 7: The Jordan Normal Form: Solving Autonomous, Homogeneous Linear Systems. April 2, 2003

A New Skew Linear Interpolation Characteristic Difference Method for Sobolev Equation

Response of LTI Systems to Complex Exponentials

Boyce/DiPrima 9 th ed, Ch 7.6: Complex Eigenvalues

Study of Tyre Damping Ratio and In-Plane Time Domain Simulation with Modal Parameter Tyre Model (MPTM)

Mixing time with Coupling

Poisson Arrival Process

Poisson Arrival Process

Chapter 21: Connecting with a Network

Ideal Amplifier/Attenuator. Memoryless. where k is some real constant. Integrator. System with memory

Math-303 Chapter 7 Linear systems of ODE November 16, Chapter 7. Systems of 1 st Order Linear Differential Equations.

Linear Time Invariant Systems

E F. and H v. or A r and F r are dual of each other.

and integrated over all, the result is f ( 0) ] //Fourier transform ] //inverse Fourier transform

Mathematical Preliminaries for Transforms, Subbands, and Wavelets

Math 2414 Homework Set 7 Solutions 10 Points

MODERN CONTROL SYSTEMS

Math 6710, Fall 2016 Final Exam Solutions

David Randall. ( )e ikx. k = u x,t. u( x,t)e ikx dx L. x L /2. Recall that the proof of (1) and (2) involves use of the orthogonality condition.

Update. Continuous-Time Frequency Analysis. Frequency Analysis. Frequency Synthesis. prism. white light. spectrum. Coverage:

ECEN620: Network Theory Broadband Circuit Design Fall 2014

Advanced Engineering Mathematics, K.A. Stroud, Dexter J. Booth Engineering Mathematics, H.K. Dass Higher Engineering Mathematics, Dr. B.S.

A TRANSIENT HEAT CONDUCTION PROBLEM OF SEMI-INFINITE SOLID CIRCULAR CYLINDER AND ITS THERMAL DEFLECTION BY QUASI-STATIC APPROACH

( ) ( ) ( ) 2011 HSC Mathematics Solutions ( 6) ( ) ( ) ( ) π π. αβ = = 2. α β αβ. Question 1. (iii) 1 1 β + (a) (4 sig. fig.

Department of Mathematics. Birla Institute of Technology, Mesra, Ranchi MA 2201(Advanced Engg. Mathematics) Session: Tutorial Sheet No.

State-Space Model. In general, the dynamic equations of a lumped-parameter continuous system may be represented by

Today s topic 2 = Setting up the Hydrogen Atom problem. Schematic of Hydrogen Atom

Modeling of the CML FD noise-to-jitter conversion as an LPTV process

2 Dirac delta function, modeling of impulse processes. 3 Sine integral function. Exponential integral function

One of the common descriptions of curvilinear motion uses path variables, which are measurements made along the tangent t and normal n to the path of

Week 8 Lecture 3: Problems 49, 50 Fourier analysis Courseware pp (don t look at French very confusing look in the Courseware instead)

Mathematical Statistics

The Development of Suitable and Well-founded Numerical Methods to Solve Systems of Integro- Differential Equations,

Calculus Limits. Limit of a function.. 1. One-Sided Limits...1. Infinite limits 2. Vertical Asymptotes...3. Calculating Limits Using the Limit Laws.

Web-appendix 1: macro to calculate the range of ( ρ, for which R is positive definite

λiv Av = 0 or ( λi Av ) = 0. In order for a vector v to be an eigenvector, it must be in the kernel of λi

DEVELOPMENTOF FERRITE BASED ABSORBING PAINTS FOR REDUCTION OF EMI /EMC PROBLEMS IN MICROWAVE RANGE

ECE351: Signals and Systems I. Thinh Nguyen

Dielectric Waveguide 1

Outline. Review Homework Problem. Review Homework Problem II. Review Dimensionless Problem. Review Convection Problem

A Dash of Maxwell s. A Maxwell s Equations Primer. Chapter V Radiation from a Small Wire Element

FIXED FUZZY POINT THEOREMS IN FUZZY METRIC SPACE

EXERCISE - 01 CHECK YOUR GRASP

F.Y. Diploma : Sem. II [CE/CR/CS] Applied Mathematics

Lecture 15: Three-tank Mixing and Lead Poisoning

ECE-314 Fall 2012 Review Questions

SHINGLETON FOREST AREA Stand Level Information Compartment: 44 Entry Year: 2009

Suggested Solution for Pure Mathematics 2011 By Y.K. Ng (last update: 8/4/2011) Paper I. (b) (c)

Comparison between Fourier and Corrected Fourier Series Methods

SHAPE DESIGN SENSITIVITY ANALYSIS OF CONTACT PROBLEM WITH FRICTION

ECE 570 Session 7 IC 752-E Computer Aided Engineering for Integrated Circuits. Transient analysis. Discuss time marching methods used in SPICE

TAKA KUSANO. laculty of Science Hrosh tlnlersty 1982) (n-l) + + Pn(t)x 0, (n-l) + + Pn(t)Y f(t,y), XR R are continuous functions.

Numerical Method for Ordinary Differential Equation

BE.430 Tutorial: Linear Operator Theory and Eigenfunction Expansion

EXISTENCE THEORY OF RANDOM DIFFERENTIAL EQUATIONS D. S. Palimkar

ON H-TRICHOTOMY IN BANACH SPACES

N! AND THE GAMMA FUNCTION

A TAUBERIAN THEOREM FOR THE WEIGHTED MEAN METHOD OF SUMMABILITY

CHAPTER 5 CIRCULAR MOTION

Ch. 6 Free Electron Fermi Gas

1 Notes on Little s Law (l = λw)

Physics 111. Lecture 38 (Walker: ) Phase Change Latent Heat. May 6, The Three Basic Phases of Matter. Solid Liquid Gas

Lecture 2: Bayesian inference - Discrete probability models

Dynamics of Bloch Electrons 1

ELEC9721: Digital Signal Processing Theory and Applications

1985 AP Calculus BC: Section I

The Central Limit Theorems for Sums of Powers of Function of Independent Random Variables

K3 p K2 p Kp 0 p 2 p 3 p

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

Boyce/DiPrima 9 th ed, Ch 7.8: Repeated Eigenvalues

STK4080/9080 Survival and event history analysis

COS 522: Complexity Theory : Boaz Barak Handout 10: Parallel Repetition Lemma

Actuarial Society of India

(Reference: sections in Silberberg 5 th ed.)

Transcription:

NCTU Dpam of Elcical ad Compu Egiig Sio Cous <Dyamic Sysm alysis ad Simulaio By Pof. Yo-Pig Ch 8. Phas Pla Mhod. Phas Pla B. Naual Mods of TI Sysms C. Phas Poais Phas pla mhod is a impoa gaphical mhods o dal wih poblms lad o a scod-od auoomous sysm.. Phas Pla Fis, l s ioduc h basic cocps of h phas pla by a scod-od auoomous sysm which is mahmaically pssd as blow: & ( ) f ( ( ) ( ) ), ( ), & ( ) f( ( ), ( ) ), ( ) (8-) (8-) wihou ay ipus. I is kow ha h sysm sas () ad () a uiquly dmid by iiial alus ( ) ad ( ) codiios ( ) ad ( ). Claly, diff iiial will sul i diff sas () ad () fo. If w adop h alus of h wo sa aiabls as wo as, h a coodia is fomd i Figu 8-, calld h phas pla. W ca daw all h pois ( (), ()) fo o obai a ajcoy wih spc o spcifid iiial codiios ( ) ad ( ). Figu 8- shows a ampl of h ajcoy ( (), ()) fo <<5, wh a aow is usd o show h dicio of h ajcoy fom o 5 sc. ( (), ()) ( (3), (3)) 3 ( (4), (4)) 4 ( (), ()) 5 ( (5), (5)) ( (), ()) Figu 8-8-

NCTU Dpam of Elcical ad Compu Egiig Sio Cous <Dyamic Sysm alysis ad Simulaio By Pof. Yo-Pig Ch Sic ach ajcoy is uiquly dmid by a s of iiial codiios, ha mas o iscio iss amog all h ajcois, cp a a sigula poi (, ) saisfyig & ( ) ad ( ) ( ) & fo all, i.., f (8-3),, ( ) f (8-4),, I fac, (, ) is a quilibium poi ad o cssay o b uiqu. s of ajcois obaid fom diff iiial codiios is calld a phas poai, as show i Figu 8-, wh P is a quilibium poi o a sigula poi. P Figu 8- B. Naual Mods of TI Sysms Th aual mods flcs h iisic popis of sysms which a o cid by ay al ipus. I oh wods, h aual mods a h sysm bhaio causd oly by is iiial codiios. s mploy h simpls d od dyamic sysm as a ampl o plai how o dmi h aual mods of a TI sysm. Cosid h followig d od TI sysm wihou ay ipu: o i a mai fom ( ) a ( ) a ( ) &, () (8-5) ( ) a ( ) a ( ) &, () (8-6) &, (8-7) 8-

wh a a, a a ( ) ( ) NCTU Dpam of Elcical ad Compu Egiig Sio Cous <Dyamic Sysm alysis ad Simulaio By Pof. Yo-Pig Ch ad. Claly, h oigi (,) i h phas pla is a quilibium poi sic h uh of lads o &. Fom (-8), h sa co ca b sold as (8-8) wih dfid as I (8-9) a fom simila o. I is kow ha ay squa mai ca b asfomd io a mai i diagoal fom o oda fom, pssd as (8-) wh is iibl. I is asy o show ha Th a h ypical foms of, dod as,,, (8-) c, o (8-) c wh h disgoal ms a h igalus of ad ca b sold fom is chaacisic quaio a a I (8-3) a a Tha mas hy a h soluios of ( a a ) a a a a (8-4) ad boh of hm may b al o compl i cojuga ad may b disic o pad. ccodigly, h a h cass of h igalus: Cas- disic al igalus, Cas- pad al igalus, Cas-3 compl igalus i cojuga, α jω ad α jω. c I is kow ha fo a igalu i h iss a las o igco i such ha (8-5) i i i c 8-3

NCTU Dpam of Elcical ad Compu Egiig Sio Cous <Dyamic Sysm alysis ad Simulaio By Pof. Yo-Pig Ch 8-4 Fo Cas-, basd o (8-5) w ha (8-6) which suls i [ ] [ ] 443 443 443 (8-7) No ha is h fis ypical fom i (8-) ad is iibl ad calld h igcomai sic i coais all h igcos of. Claly, (8-7) ca b aagd as, a fom gi i (8-). Bsids, w ha (8-8) Hc, h sa co i (8-8) is obaid as I I (8-9) o i h followig fom ( ) α α (8-) ( ) α α (8-) wh α, α, α ad α dpd o () ad (). Th fucios ad a calld h aual mods of his sysm. Fo Cas- of pad igalu, w may obai o o wo igcos cospodigly. s assum oly o igco is foud, h w ha (8-)

NCTU Dpam of Elcical ad Compu Egiig Sio Cous <Dyamic Sysm alysis ad Simulaio By Pof. Yo-Pig Ch 8-5 I is kow ha w ca ca a w co which saisfis (8-3) ad mos impoaly ad a idpd. s a sul, ww aai [ ] [ ] 443 443 443 (8-4) No ha is h scod ypical fom i (8-) ad is iibl sic ad a idpd. Claly, (8-4) ca b aagd as, sam as (8-). Bsids, w ha (8-5) Hc, h sa co i (8-8) is obaid as 3 I I (8-6) o i h followig fom ( ) β γ, (8-7) ( ) β γ, (8-8) wh γ, γ, β ad β dpd o () ad (). Th fucios ad a calld h aual mods of his sysm. Fo Cas-3 of compl igalus i cojuga, w will skip h diaio ad jus show h suld sa aiabls which a ( ) si B cos ω ω α α (8-9) ( ) si B cos ω ω α α (8-3)

NCTU Dpam of Elcical ad Compu Egiig Sio Cous <Dyamic Sysm alysis ad Simulaio By Pof. Yo-Pig Ch wh,, B ad B dpd o () ad (). Th fucios α siω α cosω a calld h aual mods of his sysm. No ha if α, i.., h igalus a pu imagiay, h h aual mods a cosω ad siω. I is cla ha all h aual mods dpd o h igalus of h sysm mai. If a igalu is locad o h lf half compl pla, h h aual mod cospodig his igalu will cog o as. O h oh had, if a igalu is locad o h igh half compl pla, h h lad aual mod will icas o as. s fo h igalu o h imagiay ais, is aual mod will oscilla. ad C. Phas Poais Now, l s us som ampls o daw phas poais, a s of ajcois, ad discuss hi lad popis. cually, diff aual mods will sul i diff ajcois ad diff phas poais as wll. Cosid h scod-od TI sysm i caoical fom, which is dscibd as blow: ( ) ( ) &, () (8-3) ( ) a ( ) a ( ) &, () (8-3) Th sysm mai is a a ad h chaacisic quaio is I a a (8-33) a a whos oos a a ± a 4a, (8-34) If a a, h ad a al ad disic. Wihou loss of galiy, 4 assum, h h a h cass lisd as blow: I. < <, i.., a a, a < ad a <. 4 II., i.., a a, a ad a <. 4 8-6

III., i.., a a ad a. 4 NCTU Dpam of Elcical ad Compu Egiig Sio Cous <Dyamic Sysm alysis ad Simulaio By Pof. Yo-Pig Ch Thi phas poais a dpicd i Figu 8-3(a)(b)(c), spcily. Fom (8-), h fis sa aiabl ( ) is ( ) α α (8-35) ad fom (8-3), h scod sa aiabl is ( ) ( ) α α & (8-36) Fo som iiial codiios () ad (), h coffici α i (8-35) may b qual o, h w ha ( ) α (8-37) which implis ( ) α ( ) (8-38) ( ) ( ) (8-39) i.., h lad ajcoy is h saigh li show i Figu 8-3(a)(b)(c). Similaly, if iiial codiios sul i α, h ( ) α (8-4) which implis ( ) α ( ) (8-4) ( ) ( ) (8-4) i.., h lad ajcoy is h saigh li show i Figu 8-3(a)(b)(c). I Cas-I wih < <, is phas poai is dpicd i Figu 8-3(a). s, du o h fac ha h phas poai mus saisfy ad, all h ajcois cp i ( ) α (8-43) ( ) α ( ) (8-44) i.., ( ) ( ). I oh wods, all h ajcois cp i h phas poai will appoach as. Claly, i is a sabl sysm sic h sysm sa will cog o h quilibium poi (,) as. 8-7

NCTU Dpam of Elcical ad Compu Egiig Sio Cous <Dyamic Sysm alysis ad Simulaio By Pof. Yo-Pig Ch (a) (b) (c) Figu 8-3 I Cas-II ad Cas III, hi phas poais a dpicd i Figu 8-3(b) ad Figu 8-3(c). s, du o h fac ha ajcois cp i h phas poai mus saisfy ( ) α ad, all h (8-45) ( ) α ( ) (8-46) i.., ( ) ( ). I oh wods, all h ajcois cp i h phas poai will appoach as. Claly, i is a usabl sysm sic h sysm sa will mo o (, ) o (, ) as, cp i Cas-II. How, Cas-II is sill ad as usabl bcaus ay mly small disubac will di h sysm o la ad mo o (, ) o (, ). (a) Figu 8-4 (b) If a a, h fom (8-34) w ha 4 lisd as blow: I. <, i.., a a ad a <. 4 a. Th a wo cass 8-8

., i.., a a ad a. 4 NCTU Dpam of Elcical ad Compu Egiig Sio Cous <Dyamic Sysm alysis ad Simulaio By Pof. Yo-Pig Ch Boh phas poais a dpicd i Figu 8-4(a)(b), spcily. I fac, hy a spcial cass of Cas-I ad Cas-II wh ad a mgd ogh. Fom h phas poais, i is asy o coclud ha Cas-I is sabl ad Cas- is usabl. (a) (b) (c) Figu 8-5 If a a, h ad a compl i cojuga. Th a wo cass lisd as blow: 4 < I., ±jω wh ω, ω o ω<. II., α±jω wh α o α<. Thi phas poais a dpicd i Figu-5(a)(b)(c) ad Figu-6(a)(b). Fo Cas-I, fom (8-9) wih α, h fis sa aiabl is ( ) cosω B siω C cos( ω θ ) (8-47) ad fom (8-3) h scod sa aiabl is ( ) ( ) Cω si( ω θ ) & (8-48) wh C is dmid by h iiial codiios () ad (). Claly, i is a oscillaoy sysm. Fom (8-47) ad (8-48), w ha ( ) (8-49) ω ( ) C whos phas poais a llipss o cicls, show i Figu 8-5(a)(b)(c) spcily fo ω, ω ad ω<. No ha h aow i hs figus is poiig fom lf o igh i h gio () sic & ( ) ( ), i.., ( ) is icasd as icass, ad fom igh o lf i h gio ()< sic & ( ) ( ), i.., ( ) < is dcasd as icass. Whil passig hough h ( ) ais, i.., ( ) ( ) h aow gos dow o up ically. &, 8-9

NCTU Dpam of Elcical ad Compu Egiig Sio Cous <Dyamic Sysm alysis ad Simulaio By Pof. Yo-Pig Ch (a) Figu 8-6 (b) Th abo dscipio of h aow dicio is also suiabl fo Cas-I o Cas-. Fo Cas-II, fom (3), h fis sa aiabl is α α ( ) ( cosω B siω) C cos( ω θ ) (8-5) ad fom (8-3) h scod sa aiabl is α α ( ) ( ) C ω si( ω θ ) Cα cos( ω θ ) α C cos( ω θ ) & (8-5) Claly, i is a dig oscillaoy sysm fo α as show i Figu 8-6(a) ad a auad oscillaoy sysm fo α< as show i Figu 8-6(b). Ulik a TI sysm, h phas poai of a olia sysm is of mo complicad sic i may ha mulipl quilibium pois o limi cycls. How, h local bhaio of a coiuous olia sysm aoud a quilibium poi ca b dmid by liaizaio which suls i a appoimd TI sysm a o h quilibium poi. If h appoimad TI sysm is sabl h h quilibium poi is sabl, ohwis i is usabl. Now, cosid a olia sysm ad dmi is sabiliy by h us of phas pla mhod. h sa quaios b & & ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) (8-5) If ( ) ( ), h i is asy o sol h sa aiabls as ( ) cos( θ ) ad ( ) si( θ ) ( ) ( ). I oh wods, if h sysm is iiially locad o h cicl, h h sysm sa will kp i a cicula ajcoy. Whil h iiial codiio is a lil diad fom h cicl, fo ampl ( ) ( ) ε, 8-

NCTU Dpam of Elcical ad Compu Egiig Sio Cous <Dyamic Sysm alysis ad Simulaio By Pof. Yo-Pig Ch is bhaio is h god by & & ( ) ε( ) ( ) ( ) ( ) ε ( ) (8-53) ε whos sysm mai is. Cospodigly, is chaacisic quaio is ε ε (8-54) ε ± 4ε 4 ad h igalus a, ε ± j. Claly, fo ε h sysm is i h gio ousid h cicl ad bcoms usabl, i.., i mos away h cicl o h ifiiy. Fo ε< h sysm is i h gio isid h cicl ad bcoms sabl, i.., i mos away h cicl o h oigi. Figu 8-7 shows h phas poai icludig a limi cycl ( ) ( ). ( ) ( ) Figu 8-7 quaios N, l s cosid a aiabl sucu sysm, dscibd by h followig sa & & ( ) ( ) ( ) (. 5. 5sig( ( ) ( ) )) ( ) (8-55) I is obious ha h a wo sucus of h sysm, o fo ( ) ( ) ad h oh fo ( ) ( ). No ha h sysm is udfid fo ( ) ( ) < ha ( ) ( ) is god by. I cas, i.., i h fis ad hid quadas of h phas pla, h sysm 8-

& & ( ) ( ) ( ) 4 ( ) which pfoms as a llips. Whil ( ) ( ) NCTU Dpam of Elcical ad Compu Egiig Sio Cous <Dyamic Sysm alysis ad Simulaio By Pof. Yo-Pig Ch < quadas of h phas pla, h sysm is god by & & ( ) ( ) ( ) ( ) (8-56), i.., i h scod ad fouh (8-57) which is a cicl. Th suld phas poai is plod by h followig MTB simulaio: Ca m-fil: ss.m fucio dss(,) dzos(,); d()(); d()-(.5.5*sg(()*())*(); % ky i h followig isucios [,]od45(@ss,[:.:],[ ]) plo((:,),(:,)); labl( ); ylabl( ); gid 5-5 - -5 5 I ca b cocludd ha h aiabl sucu sysm is sabl, hough boh of is sub-sysms a usabl. Poblms P.8- Cosid h followig sysm: ( ) ( ) &, ()5 ( ) ( ) ( ) &, () 4 8-

NCTU Dpam of Elcical ad Compu Egiig Sio Cous <Dyamic Sysm alysis ad Simulaio By Pof. Yo-Pig Ch Daw h ajcoy o h phas pla fom o i MTB. 8-3