Hopf Bifurcation Analysis for the Comprehensive National Strength Model with Time Delay Xiao-hong Wang 1, Yan-hui Zhai 2

Similar documents
Chapter 7 Stead St y- ate Errors

CHAPTER-2. S.No Name of the Sub-Title No. 2.5 Use of Modified Heffron Phillip's model in Multi- Machine Systems 32

Hypernormal Form at Quintic of a Class of Four-dimensional Vector Fields with Linear Part has Two Pairs of Pure Imaginary Eigenvalues

Conventional Hot-Wire Anemometer

Generalized Half Linear Canonical Transform And Its Properties

Consider a system of 2 simultaneous first order linear equations

CIVL 8/ D Boundary Value Problems - Triangular Elements (T6) 1/8

Supplementary Figure 1. Experiment and simulation with finite qudit. anharmonicity. (a), Experimental data taken after a 60 ns three-tone pulse.

Engineering Circuit Analysis 8th Edition Chapter Nine Exercise Solutions

Instructors Solution for Assignment 3 Chapter 3: Time Domain Analysis of LTIC Systems

Summary: Solving a Homogeneous System of Two Linear First Order Equations in Two Unknowns

Lecture 4 : Backpropagation Algorithm. Prof. Seul Jung ( Intelligent Systems and Emotional Engineering Laboratory) Chungnam National University

Problem 1: Consider the following stationary data generation process for a random variable y t. e t ~ N(0,1) i.i.d.

Chapter 13 Laplace Transform Analysis

9. Simple Rules for Monetary Policy

An N-Component Series Repairable System with Repairman Doing Other Work and Priority in Repair

Fractal diffusion retrospective problems

The Variance-Covariance Matrix

NONLOCAL BOUNDARY VALUE PROBLEM FOR SECOND ORDER ANTI-PERIODIC NONLINEAR IMPULSIVE q k INTEGRODIFFERENCE EQUATION

Advanced Queueing Theory. M/G/1 Queueing Systems

The Mathematics of Harmonic Oscillators

Wave Phenomena Physics 15c

Final Exam : Solutions

Euler-Maruyama Approximation for Mean-Reverting Regime Switching CEV Process

Improved Exponential Estimator for Population Variance Using Two Auxiliary Variables

Chapter 9 Transient Response

14.02 Principles of Macroeconomics Fall 2005 Quiz 3 Solutions

Valuation and Analysis of Basket Credit Linked Notes with Issuer Default Risk

KNOWLEDGE AND TECHNOLOGY TRANSFER: NOVEL RESULTS WHILE TEACHING ADVANCED COURSES

Cmd> data<-matread("hwprobs.dat","exmpl8.2") This is the data shown in Table 8.5, the copper and diet factors.

Chapter 5 The Laplace Transform. x(t) input y(t) output Dynamic System

10.5 Linear Viscoelasticity and the Laplace Transform

SIMEON BALL AND AART BLOKHUIS

Valuing Credit Derivatives Using Gaussian Quadrature : A Stochastic Volatility Framework

Valuation of a Basket Loan Credit Default Swap

Implementation of the Extended Conjugate Gradient Method for the Two- Dimensional Energized Wave Equation

Poisson process Markov process

Improved Exponential Estimator for Population Variance Using Two Auxiliary Variables

1 Finite Automata and Regular Expressions

4.1 The Uniform Distribution Def n: A c.r.v. X has a continuous uniform distribution on [a, b] when its pdf is = 1 a x b

Transfer function and the Laplace transformation

1. Inverse Matrix 4[(3 7) (02)] 1[(0 7) (3 2)] Recall that the inverse of A is equal to:

THE STOCHASTIC SEASONAL BEHAVIOR OF THE NATURAL GAS PRICE

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

DEPARTMENT OF ELECTRICAL &ELECTRONICS ENGINEERING SIGNALS AND SYSTEMS. Assoc. Prof. Dr. Burak Kelleci. Spring 2018

CONTINUOUS TIME DYNAMIC PROGRAMMING

Enhanced DBCC for high-speed permanent magnet synchronous motor drives

Response of LTI Systems to Complex Exponentials

UNIT #5 EXPONENTIAL AND LOGARITHMIC FUNCTIONS

Frequency Response. Response of an LTI System to Eigenfunction

Cost Effective Multi-Period Spraying for Routing in Delay Tolerant Networks

Chapter 8 Theories of Systems

EE243 Advanced Electromagnetic Theory Lec # 10: Poynting s Theorem, Time- Harmonic EM Fields

(,,, ) (,,, ). In addition, there are three other consumers, -2, -1, and 0. Consumer -2 has the utility function

Homework: Introduction to Motion

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

Ma/CS 6a Class 15: Flows and Bipartite Graphs

An Indian Journal FULL PAPER. Trade Science Inc. A stage-structured model of a single-species with density-dependent and birth pulses ABSTRACT

EXERCISE - 01 CHECK YOUR GRASP

Midterm exam 2, April 7, 2009 (solutions)

FL/VAL ~RA1::1. Professor INTERVI of. Professor It Fr recru. sor Social,, first of all, was. Sys SDC? Yes, as a. was a. assumee.

Lecture 1: Numerical Integration The Trapezoidal and Simpson s Rule

t=0 t>0: + vr - i dvc Continuation

A universal saturation controller design for mobile robots

Response of MDOF systems

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

AQUIFER DRAWDOWN AND VARIABLE-STAGE STREAM DEPLETION INDUCED BY A NEARBY PUMPING WELL

Lecture 4: Laplace Transforms

Control Systems. Transient and Steady State Response.

I M P O R T A N T S A F E T Y I N S T R U C T I O N S W h e n u s i n g t h i s e l e c t r o n i c d e v i c e, b a s i c p r e c a u t i o n s s h o

On the Derivatives of Bessel and Modified Bessel Functions with Respect to the Order and the Argument

innovations shocks white noise

Improved Ratio Estimators for Population Mean Based on Median Using Linear Combination of Population Mean and Median of an Auxiliary Variable

Control of a Lake Network Invasion: Bioeconomics Approach

t + t sin t t cos t sin t. t cos t sin t dt t 2 = exp 2 log t log(t cos t sin t) = Multiplying by this factor and then integrating, we conclude that

Some Transcendental Elements in Positive Characteristic

Fourier Series: main points

Gauge Theories. Elementary Particle Physics Strong Interaction Fenomenology. Diego Bettoni Academic year

Introduction to Laplace Transforms October 25, 2017

Chapter 12 Introduction To The Laplace Transform

V A. V-A ansatz for fundamental fermions

AR(1) Process. The first-order autoregressive process, AR(1) is. where e t is WN(0, σ 2 )

Pricing Arithmetic Average Reset Options With Control Variates

State Observer Design

A. Inventory model. Why are we interested in it? What do we really study in such cases.

Chapter 5 Transient Analysis

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

Economics 302 (Sec. 001) Intermediate Macroeconomic Theory and Policy (Spring 2011) 3/28/2012. UW Madison

ELEN E4830 Digital Image Processing

Stability and Bifurcation in a Neural Network Model with Two Delays

Modified-Zero-Fourth-Order-Cumulant-Expansion Approximation Method for Investigation of Turbulence with Reacting and Mixing Chemical Elements

2. The Laplace Transform

Relation between Fourier Series and Transform

A Simple Representation of the Weighted Non-Central Chi-Square Distribution

Analysis of Laser-Driven Particle Acceleration from Planar Transparent Boundaries *

Wave Equation (2 Week)

English Made Easy: Foundation Book 1 Notes for parents

Graphene nanoplatelets induced heterogeneous bimodal structural magnesium matrix composites with enhanced mechanical properties

T ansmission Line Theory Smith Chart

MEM 355 Performance Enhancement of Dynamical Systems A First Control Problem - Cruise Control

Transcription:

opf Bfurcaon Analy for h Comprhnv Naonal Srnh ol wh m Dlay Xao-hon an Yan-hu Zha School of Scnc ann Polychnc Unvry ann Abrac h papr manly mof an furhr vlop h comprhnv naonal rnh mol By mofyn h bac comprhnv naonal rnh mol can mor accuraly llura h ocy phnomna wh m lay Fr w rarch h ynamc of h mof wh m lay By mployn h normal form hory an cnr manfol mho w oban om abl rul on h u h concluon confrm ha a opf bfurcaon occur u o h xnc of ably wch whn h lay var Fnally om numrcal mulaon ar vn o llura h ffcvn of our rul Kywor opf bfurcaon Sably Comprhnv naonal rnh mol Cnr manfol Normal form I E ESABLISEN OF E ODEL Ornary ffrnal uaon u o uy h comprhnv naonal rnh of a counry ovrall uaon whch raually vlop n rcn yar h arcl brn n m lay n h comprhnv naonal rnh mol B h comprhnv naonal rnh mol mol ar al analy x x har powr rnh n whch a compo ncaor h lvl of h maral cvlaon rourc conomy mlary cnc an chnoloy c of a counry A x br h maral cvlaon mor proprou y of powr rnh n whch h prual cvlaon of a counry y an for ocal vl con-makn rror ucaon falur offcal corrupon popl al c A h m of powr ha obacl o h ocal vlopmn y man of powr upror w con-makn naonal ualy c ha ha a promon ffc on ocal vlopmn m conan So n h lraur [] comprhnv naonal rnh hown by h follown ffrnal uaon: x x x y y y m x x A known o all om chan of h of powr can b rflc on h har powr afr a cran pro of m h lraur [] ha ablh h follown lay mol of h comprhnv naonal rnh x x x y y y m x x pov an h ohr paramr ar h am a of hl mol of ha conr h m lay n h ral worl har powr rnh no only rly on har powr rnh of of powr n h pa bu alo rly on of powr In orr o mak h mol mor accura w mofy h ym of o h follown form: x x x ky y y y m x x 3 Pa 74

k h oranaon of h papr a follow: rarn a bfurcaon paramr w uy h ably of h ulbrum pon of h ym 3 an opf bfurcaon of h ulbrum pnn on hn ba on h nw normal form of h ffrnal-albrac ym nrouc by Chn al [3] an h normal form approach hory an cnr manfol hory nrouc by aar al [4] w rv h formula for rmnn h propr of opf bfurcaon of h ym n h hr con Numrcal mulaon am a ufyn h horcal analy wll b rpor n Scon 4 Fnally h papr n wh a cuon II SABILIY AND LOCAL OPF BIFURCAION ANALYSIS h ably an opf concluon for h ym of can b oban rcly from h lraur [56] From ym 3 w can ha hr x ulbrum E m m m Y E X Accorn o h praccal nfcanc of h mol hr w only cu h problm of h hopf bfurcaon an ably for h ol pov ulbrum pon E S u x X u y Y h ym 3 bcom u p u u ku u u u m p u u m hr p An h lnaraon of ym a E u p u ku u u u m p u h characrc uaon of ym a follow c f 3 c p m p f p m p k o uy h ably of h ulbrum pon E an bfurcaon ym w only n o cu h rbuon of h roo for characrc E 3 If h all roo of uaon 3 hav nav ral par h ulbrum pon E ay If h uaon ha a roo ha conan pov ral par h ulbrum pon E no abl In orr o uy h rbuon of h roo for E 3 Conrn n fr characrc uaon 3 c 4 hr p m p c Pa 75

Obvouly all roo of uaon of 3 ha nav ral par f af So h ym n ulbrum E locally aympocally abl for Now w nva h local ably aroun h pov ulbrum pon for h E whn ym 3 an h xnc of opf bfurcaon occurrn a h ulbrum pon Lmma For h ym 3 E ha a par of purly manary roo for whn af Proof If a oluon of h characrc uaon whn an only whn m c co n h paraon of h ral an manary par yl f f co c n 5 whch la o v c f v f 6 v aum ha h coffcn af h follown conon If h conon c f 4 f c f c f 4 f or roo From 6 w oban A a rul whn af hn E 6 ha pov roo hrfor E hav purly manary f c c f 4 f f arcco h characrc uaon hav a par of purly manary roo h proof compl Lmma h ranvraly conon R or af Proof By ffrnan boh of E 3 wh rar o an olvn hav hn c R c R Pa 76

can c f by conon co cn c f or hu R h proof compl Lmma 3 For E 3 f all of h roo hav nav ral par Pov ulbrum aympocally [ abl an h pov ulbrum prouc opf bfurcaon n III E DIRECION OF E OPF BIFURCAION AND E SABILIY OF PERIODIC SOLUIONS In w oban h conon of opf bfurcaon In h con w cu h rcon of opf bfurcaon an h ably of h bfurcan proc oluon ba on h normal form nrouc by Chn al [3] an h cnr manfol hory nrouc by aar al [4] In h follown par w aum ha h ym 3 unro opf bfurcaon a h ulbrum E for an w l h corrponn purly manary roo of h characrc uaon a h ulbrum E S R clarly h opf bfurcaon valu of ym 3 S u u for convnnc w connu o u u a u Dffrnal Euaon FDE ym n C C[ ] R hn h ym 3 uvaln o h follown Funconal u L F u 3 u u u L : C R F : RC R ar vn rpcvly by L B B F L 为 C[ ] R p B m p k B By h R rprnaon horm hr x a marx funcon who componn ar boun varaon funcon :[ ] R uch ha In fac w ak L C 3 Pa 77

B B 33 a Dla funcon For C [ ] R w fn opraor A an R a follow [ A 34 F hn h ym can b wrn a h follown form: [ R 35 u A u R u 36 Sn C [] R h aon opraor A an a blnar nnr prouc vn by A of A fn a ] 37 38 From h cuon n Scon w know ha ar nvalu of A hu hy ar alo nvalu of w calcula h nvcor of A blonn o an nvcor By h fnon of A w hav A an In aon of A Nx A blonn o h nvalu B B A 39 L E39 bcom p k m p p k By h fnon of A w hav A an hn Pa 78

Pa 79 A B B 3 L D E3 bcom p m k p p ] [ D ] [ B D ] [ k D o ] [ k D Nx w uy h pcfc paramr of h rcon an for bfurca proc oluon Un h am noaon a n aar al [4] w fr compu h coorna o crb h cnr manfol C a Dfn } R{ u u 3 On h cnr manfol C w hav 3 hr In fac an ar local coorna for cnr manfol C n h rcon of an No ha ral f u ral conr only ral oluon For h oluon of E36 n cnr manfol C w hav R A * R A R A F 33

L from 33 an 35 w hav By 3 w hav hn 34 35 F F V V F F 36 F u u u 37 u u D[ [ ] Comparn h coffcn wh 36 follow ha D D D D[ ] By E35 36 an fnon 3 of w hav Pa 8

A R u A R F F F [ A 38 L A hr From h fnon 3 of w hav 39 3 By 38 39 an 3 w hav A A 3 For [] by 34 38 an 39 w R[ F ] Comparn h coffcn wh 39 w v ha 3 33 I follow from 3 ha A 34 Solvn for an w oban Pa 8

3 E E 35 By 34 an 3 w 36 From 37 an 38 w Accorn o E35 36 an 37 w hn E E E3 3 E E E 37 I E E E E 38 3 p E E p k m p m p Accorn o h abov Propoon an [7] w can compu h follown paramr: C R{ C} R{ } R{ C } 3 Im{ C} Im{ } IV NUERICAL SIULAION Numrcal mulaon how from abl o unabl complx ranformaon proc h ym 3 v a concr xampl o how h ynamc bhavor of comprhnv naonal rnh mol ak 3 k 8 m 9 n h ym 3 So w can oban Pa 8

ulbrum E 893 hrouh calculan of h ahmacal ofwar w oban 6 88 By lmma an abov concluon w can oban ha ulbrum E abl f 7 Fur By conra h ulbrum E unabl f Fur B whn 886 h proc oluon occur from h ulbrum E Fur 3-4 owvr our analy nca ha h ynamc of h comprhnv naonal rnh y mol wh m lay can b much mor complca han w may y hav xpc I ll nrn an nprn o rarch 3 5 8 6 5 4 5 x x 5 5 5 y 8 4 6 8 FIGURE E EQUILIBRIU E FIGURE E EQUILIBRIU E OF SYSE 3 IS SABLE I 7 OF SYSE 3 IS SABLE I 5 5 5 5 x FIGURE 3 SYSE 3 PRODUCE E PERIODIC SOLUIONS I 886 Pa 83

Inrnaonal Journal x of Ennrn Rarch & Scnc IJOER ISSN: [395-699] [Vol- Iu-3 arch- 6] y 5 5 5 5 4 6 8 4 5 4 6 8 4 FIGURE 4 IE SEQUENCE DIAGRA OF FIGURE 3 V CONCLUSION h papr apply lay n h comprhnv naonal rnh mol whch how rch ynamc bhavor Dffrn from prvou u w a h nflunc of m lay fback n h ym 3 Dynamc bhavor of comprhnv naonal rnh mol wh m lay analy by un h mho of uanav hn lay $\au$ acro a r of crcal valu nonlnar ynamc ym nra h opf bfurcaon In aon by mployn h normal form hory an cnr manfol mho w oban om abl rul on h u u normav hory an cnr of popular horm oban h calculan mho of h rcon of opf bfurcaon an ably of proc oluon Fnally h abov horcal analy vrf by numrcal mulaon h ynamc bhavor of h comprhnv naonal rnh mol rch any apc no mnn y o b Rfrnc CONFLIC OF INERESS h auhor clar ha hr no conflc of nr rarn h publcaon of h papr ACKNOLEDGEENS h auhor ar raful o h rfr for hr hlpful commn an conrucv uon AUOR CONRIBUIONS Concv: Xaohon an Drawn raphc: Xaohon an Calcula: Xaohon an of: Yanhu Zha ro h papr: Xaohon an REFERENCES [] S an Dffrnal uaon mol an chao Journal of Chna Scnc an chnoloy Unvry pp 3 34 999 [] X Xn opf Bfurcaon Analy n Comprhnv Naonal Powr ol wh wo Dlay[J]Journal of Yl Normal UnvryNaural Scnc Eo vol no3 pp 4 7 [3] BS Chn XX Lao YQ Lu Normal form an bfurcaon for h ffrnc-albrac ym Aca ah Appl Sn 3 49 443 n Chn [4] B aar D Kaarnoff Y an hory an Applcaon of opf Bfurcaon Cambr Unvry Pr Cambr 98 [5] L an Z Lu opf bfurcaon of prc Rylh uaon wh m lay Journal of Jln Unvry vol47 no3 9 [6] GD Zhan LL Zhu BS Chn opf bfurcaon an ably for a ffrnal albrac bolocal conomc ym Appl ah Compu 7 33 338 [7] Y an an Y Z Sably an opf bfurcaon of ffrnal uaon mol of prc wh m lay hlh of Scncpapr Onln vol4 no [8] Y Z Y B Y Xon an XN opf bfurcaon analy for h mof Raylh prc mol wh m lay Abrac an Appl Analy vol3arcl ID 9497 6 pa 3 [9] D D aar N D Kaarnoff an Y an hory an Applcaon of opf Bfurcaon Cambr Unvry Pr Cambr UK 98 [] J a YQ Chn LX Luopf bfurcaon an chao of fnancal ym on conon of pcfc combnaon of paramr Journal of Sym Scnc an Complxy 8 Vol no pp5-59 4 [] J al hory of Funconal Dffrnal Euaon Sprnr Nw York NY USA n on 977 Pa 84