Pattern Formation in Networks of Nonlinear Systems with Delay Couplings
|
|
- Mark Golden
- 5 years ago
- Views:
Transcription
1 Pattern Formation in Networks of Nonlinear Systems with Delay Couplings Eiichi Uchida Toshiki Oguchi Tokyo Metropolitan University - Minami-osawa, Hachioji-shi, Tokyo Japan ( t.oguchi@tmu.ac.jp) Abstract: This paper considers pattern formation in networks of identical nonlinear systems with time-delays. First of all, applying the harmonic balance methods, we derive the corresponding harmonic balance equations for networks of identical nonlinear systems with delay couplings. Then, solving the equations by reducing to the stability problem of linear retarded systems, we can estimate the oscillation profile such as the frequency, amplitudes and phases of coupled systems. Based on this analysis method, we also develop a design method of networks of nonlinear systems that can achieve a prescribed oscillation profile. The effectiveness of the proposed methods is shown by numerical examples. Keywords: Delay systems, Nonlinear systems, Pattern generation, Complex systems, Harmonic balance. INTRODUCTION In recent years, network systems have attracted attention in applied physics, mathematical biology, social sciences, control theory and interdisciplinary fields. In particular, synchronization and pattern formation of coupled systems has been the subject of intense study. From the standpoint of control engineering, synchronization and pattern formation of coupled systems are important notions to realize decentralized control techniques and biomimetic control approach in increasingly complex applications. For instance, central pattern generators (CPGs), which produce the rhythmic movements such as locomotion, breathing and scratching, are considered to consist of a group of neurons. The study of the CPGs has been carried out not only for understanding of biological behavior but also for producing rhythmic locomotion of multi-legged robots. The behaviors of coupled identical systems and networks have been studied by a large numbers of researchers in broad areas (Golubitsky and Stewart (4); Newman et al. (6) and so on). Golubitsky al. showed that the existence of periodic solutions in symmetric coupled systems and the pattern classification of gaits based on the symmetric network structure of the CPG. Pogromsky et al. (999, ) considered the synchronization problem in symmetric networks of chaotic systems and derived sufficient conditions for partial synchronization in networks and existence for periodic solutions. However these researchesfocus ononlythe symmetric structure of networks. Iwasaki (8) proposed a systematic approach for the analysis and synthesis of CPGs by applying the multivariable harmonic balance method. The proposed approach has no restriction on the network structure. This work was partially supported by Japan Society for the Promotion of Science (JSPS) Grant-in-Aid for Scientific Research (No. 5644). On the other hand, the research interest on synchronization of coupled systems are shifting to the synchronization problem in delayed networks during recent years (Oguchi et al. (8, 9); Neefs et al. ()). Since time-delays caused by signal transmission affect the behavior of coupled system in practical situations, it is therefore important to study the effect of time-delay in network systems. The objective of this paper is to extend the systematic approach based on the harmonic balance method (Iwasaki (8)) for delayed networks. First, we derive the corresponding harmonic balance equations for networks of identicalnonlinearsystemswithdelaycouplings.then, solving the equationsbyreducingtothe stabilityproblemoflinear retarded systems, we can estimate the oscillation profile such as the frequency, amplitudes and phases of coupled systems. Based on this analysis method, we also develop a design method of networks of nonlinear systems that can achieve a prescribed oscillation profile. The effectiveness of the proposed methods is shown by numerical examples.. NONLINEAR NETWORK SYSTEMS. Model Description In this paper, we consider n identical nonlinear oscillators interconnected as follows (Iwasaki (8)). n v i = ψ i (q i ),q i = f(s)u i,u i = µ ir (s)v j () r= whereu i andv i arethe input andtheoutputofoscillatori, ψ i is a static nonlinear function and is defined by ψ i (x)= tanh(x)andµ ir means the transfer function denoting a connection from oscillator r to oscillator i, which is given byµ ir (s)=k ir e τms.hereeachnonlinearoscillatorcanbe consideredas a model of the electricalactivity of a neuron, k ir R means a coupling gain between i and r,andτ m is Copyright by the International Federation of Automatic Control (IFAC) 58
2 a commensurate delay of τ, i.e. τ m = mτ for m =,...,h. f(s) is a linear time-invariant part of oscillator given by f(s)= ω. () s+ω Then the network of N coupled oscillators is summarized as the following equations. ψ v =Ψ(q)=. ψ n, F(s)= f(s)... = f(s)i n f(s) q = F(s)M(s)Ψ(q) f(s)µ (s)ψ (q )+ +f(s)µ n (s)ψ n (q n ) =. f(s)µ n (s)ψ (q )+ +f(s)µ nn (s)ψ n (q n ) where M(s) is the transfer matrix whose (i,j)entryis µ ij (s). For the above system, we deal with the following two problems: (I) Analysis problem : Given a network system with delays,determine whether the coupled systems have oscillatory trajectories, and if so, estimate the oscillation profile (frequency, amplitudes, phase) without actually simulating the differential equations. (II) Synthesis problem : Given a network structure with delays and a desired oscillation profile, determinethecouplingstrengthk ir inthe couplingtransfer function µ ir (s) sothattheresultingnetwork system achieves the given profile.. Analysis of periodic solutions We assume that the system () has a periodic solution. Applying Fourier series expansion, q is described by q(t)= a k sin(ωkt)+b k cos(ωkt) (3) k= where a k,b k R n and ω R. Supposing that ω is sufficiently higher than the band-pass frequency of f(s), q i (t) can be approximated by q i (t) α i sin(ωt+φ i ),i=,...,n where α i and φ i denote the amplitude and phase, respectively. Furthermore, using a describing function κ i, ψ i can be approximated by ψ i (q i ) κ i (α i )q i,κ i (x):= π ψ i (xsinθ)sinθdθ. πx The describing function κ i (x) is a monotonically decreasing function satisfying κ i () = and κ i ( ) =since ψ i (x) = tanh(x). Using these approximations, we obtain the corresponding multivariable harmonic balance equation for the nonlinear oscillator delay network as follows. (M(jω)K(α) /f(jω)i n )q = (4) where j denotes the imaginary unit, I n the n n identity matrix, q i := α i e jφi and K(α) := diag(κ i (α i )). The triplet (ω,α,φ) satisfying the equation(4) is the oscillation profile to be solved. In this case, however, not only are the solution unique in general, there may exist infinite number of solutions due to the existence of time-delays. Among the solutions, we have to choose a triplet (ω,α,π) so that the estimated oscillation is stable. Here, replacing K(α)witha constant matrix K(α):= K( q ), the stability of oscillation can be expected by checking if the following characteristic quasi-polynomial of linear system has a pair of solutions s = ±jω on the imaginary axis and the rest in the open right half plane. ) det (M(s)K(α) /f(s)i n =. (5) If the couplings have no delay, the stability analysis can be reduced to the eigenvalue problem of the constant matrix MK(α)sincethe matrixm isaconstantmatrix.however, if the couplings have delays, then the entries of M(s) are not constant but functions of e sτm. Therefore the stability analysis cannot be done in the same way. Nowwe decompose the transfer function matrix M(s)into a constant matrix part and the rest, i.e. M(s)=M + M m e sτm (6) where M and M m R n n for m =,...,hare constant coupling matrices with delay τ...,τ h. Substituting () and (6) into the quasi-polynomial (5) and setting M c = I n +M K(α), we obtain ( ) det M m e sτm K(α)+M c λi n = (7) where λ := s/ω.therootλsatisfying the equation (7) is identical to the pole of the linear retarded system ẋ(t)=m c x(t)+ M m K(α)x(t τ m ). (8) Therefore, seeking the rightmost characteristic roots of the system (8), we can confirm whether the characteristic quasi-polynomial (5) has at least one pair of solutions on theimaginaryaxisandtherestintheopenrighthalfplane. If there exists a solution λ on the imaginary axis, it is a solution of the modified harmonic balance equation given by ( M m e sτm K(α)+M c λi n )q =. (9) Assuming that there exists a q C n satisfying (9), we can expect that the trajectory of each oscillator q i (t)has a phase φ i and frequency ω, i.e. q i (t) α i sin(ωt+ φ i ), where a triplet (ω,α i,φ i ) is determined by jω/ω = λ, q i = α i e jφi,ω> () The analysis on the oscillation profile can be accomplished by finding q C n satisfying the equation (5), but K(α) depends on q. So, in this paper, extending the algorithm introduced by Iwasaki (8) for delay systems, we propose a heuristic calculation method. Algorithm: Step Set the initial values k =andv k =[,...,]. Step Let K k = K( v k ) and find the rightmost roots λ k which has the maximum imaginary part and satisfies the characteristic quasi polynomial ( ) det M m e λ kω mτ K k +M c λ k I n = 59
3 Step 3 Find the corresponding eigenvector x k to λ k obtained in Step : ( M m e λ kω mτ K k +M c λ k I n )x k = x k = Step 4 Update the eigenvector v k+ = e σ k y k,σ k := R(λ k ),y k := x k x kv k Step 5 If v k+ v k ǫ, then the algorithmterminates, and q := v k+ is the solution of the equation (9). So, we obtain a triplet (ω,α,φ) from (). Otherwise go to Step and iterate. Remark. For computing the rightmost roots of the characteristic quasi polynomial, there are several numerical techniques and some useful softwareslike DDE-BIFTOOL are also developed..3 Numerical Examples Example In this analytical method, we do not assume the symmetric structure of the network. To show the availability of this method for complex network systems, we consider a network of n = 5 oscillators and the network structure is assumed to be given in Fig.. Each vertex denotes a oscillator and the number i inside of each vertex indicates the i-th oscillator. The edge (i,j) represents the existence of coupling between the i-th oscillator and the j-th oscillator. Now we assume that all coupling has identical delay τ and the matrix M(s)in(6) is given by M(s)=M e sτ where M =. For this network system, we estimate the profiles of oscillators for the following two cases: (i) τ =., (ii) τ =.3. Figures and 3 show the comparison between the estimated amplitudes and phases for q i and those obtained by numerical simulations for τ =.andτ =.3, respectively. Here the phaseφ i means the phase difference between the st oscillatorand the i-th oscillator. In addition, the estimated frequency ω for each delay is summarized in Table. Table. Comparison between the estimation and simulation results of the frequency of coupled oscillator τ =. τ =.3 imation.6.77 ulation Fig.. A Small World Network with n =5 Phase φi [deg] Fig.. imated and simulated phase of each oscillator (τ =.) Amplitude αi Fig. 3. imated and simulated amplitude of each oscillator (τ =.) Example Consider n = 5 identical nonlinear systems with delay couplings. The network structure is shown in Fig. 6. Each vertex denotes an oscillator, and the directed edge i,j means that there exists a coupling from the i-th oscillatorto the j-th oscillator.in particular, the red edges mean couplings without delay and others are couplings with delay τ =.3. Now the transfer function matrix of the coupling M(s) isassumedto be givenby where M(s)=M +M e sτ 4 3 5
4 5.8.6 q q 5.4 q 3 q 4 Phase φi [deg] 5 5 q(t)...4 q Fig. 4. imated and simulated phase of each oscillator (τ =.3) Amplitude αi Fig. 5. imated and simulated amplitude of each oscillator (τ =.3) Fig. 6. Network structure in Example M = and M = Thenonzeroentriesofthesematricesaregivenby random numbers. Figure 7 shows the steady state behaviors of coupled oscillators q(t) and the oscillation of each system is stable. Then, applying the foregoing algorithm, we Time Fig. 7. Oscillatory trajectories in Example Table. imated and simulated oscillation profiles in Example q q q 3 q 4 q 5 imation Amp Phase ulation Amp Phase estimate the profile of each oscillator. The comparison between the estimated values and the simulation results is summarized in Table. From Table, we can conclude that the estimation of a triplet (ω,α,φ) can be almost perfectly accomplished. 3. SYNTHESIS OF NETWORKS WITH DELAYS 3. Design method In the previous section, we considered how to estimate the profile of the coupled oscillators under given network structureandcouplings.inthissection,weconsiderhowto design the coupling gain k ij to achieve given specifications of oscillation profile. Fromthe discussionin the foregoingsection,itisnecessary that the coupling gain matrices M,...,M h satisfy the equation (5) for a given triplet (ω,α,φ). Furthermore, the roots of eq. (7) or equivalently the poles of the linear retarded system (8) have to be on the imaginary axis and in the left half-plane of the complex plane only. Based on the above idea, we obtain the following results on the design of coupling gains. Theorem. Let ω R, α R n,φ R n and mτ R for m =,...,h be given, where ω>, α i > for i =,...,nand τ<π/ω.furthermore,assumethat φ is chosen such that there exists at least one φ i satisfying φ i φ j kπ for any j i and k Z. Then the delaycoupled oscillator is expected to have oscillatory behavior if there exist M i R n n for i =,...,h and a positive symmetric matrix P satisfying where [R(Aq) I(Aq)]=RΩ () B T P +PB< () 5
5 A := B := M m e jωτm K(α)+M K(α) I M m K(α)+M K(α) I n R := [c s] R n,c i := α i cosφ i,s i := α i sinφ i [ ] ω/ω Ω:=. ω/ω Proof. For given ω, α i, φ i and τ m, the coupling gain matrices M i for i =,...,h must satisfy the following harmonic balance equation: ( M m e jωτm K(α)+M c )q = λq (3) where λ = jω/ω and M c = I n + M K(α). Now if we suppose that the characteristic quasi-polynomial has the characteristic roots at ±jω, the following equation equivalently holds: ( M m e jωτm K(α)+M c )q = j ω q (4) ω Payingattentionto that q i = α i e jφi = α i cosφ i +jα i sinφ i and rewriting the above equation, we obtain eq. (). On the( other hand, if the system (8) with τ =, i.e. ẋ(t)= M c + ) h A m x(t),isasymptoticallystable,the corresponding matrix pencil Λ is regular. Then Theorem 3 in Appendix A shows a necessary and sufficient condition for λi n M c M m e λ kω mτ K(α) = to have at least one nonzero root on the imaginary axis and the corresponding delay values. From this result, if we choose the minimum value in the delay values, the characteristic quasi polynomial has at least one pair of roots ±jω on the imaginary axis and the rest on the left half plane. Therefore it is necessary that τ< π ω and there exists a positive definite matrix P satisfying B T P +PB<. 3. Examples Consider network systems with the network structure given in Figure 6 again. Example 3 Firstly, we consider a case in which all couplings have the same length of time-delay τ.thenwe design the coupling gain matrix M such that the coupled oscillators meet the specification given in Table 3. Table 3. Specifications and simulation results (Example 3) ω α i φ φ φ 3 φ 4 φ 5 Spec Case ) Case ) It should be noted that the structures of M i have been already fixed due to the network structure. Solving the LMI q(t) Time Fig. 8. Behaviors of coupled oscillators (τ =.) () with respect to M subject to the specific structure and the equation (), we can obtain the following matrix M for τ =.and.3, respectively. ) Case τ =.: M = ) Case τ =.3: M = The comparison between the given specification and the simulation results with M mentioned above is summarized in Table 3. Figure 8 shows the behavior of each oscillator in case of τ =. From these results, we see that the design specifications are almost fulfilled. Example 4 Consider a case in which the network has both delay free couplings and time-delay couplings with delay τ =.. As in Example, the red edges indicate delay free couplings and others time-delay couplings. Based on the specifications given in Table 4, the coupling gain matrices M and M are obtained by applying Theorem as follows..6 M = M = The comparison between the given specification and the design result is summarized in Table 4 and the simulation result is shown in Fig. 9. Although the phases are short by about 5 degrees in whole, the amplitudes almost coincide with the specification. q q q 3 q 4 q 5 5
6 q(t) Table 4. Specifications and simulation results (Example 4) ω α α α 3 α 4 α 5 Spec φ φ φ 3 φ 4 φ 5 Spec Time Fig. 9. Behaviors of coupled oscillators (Example 4) 4. CONCLUSIONS In this paper, weconsideredpatternformationinnetworks of nonlinear systems with delays. The approach employed in this paper is based on the harmonic balance method and the obtained results are extensions of Iwasaki (8) for delayed networks. Although these methods contain approximations and the proposed algorithms are not guaranteed to converge, the validity of the proposed methods were supported by numerical examples. REFERENCES Chen, J., Gu, G., and Nett, C.N. (995). A new method for computing delay margins for stability of linear delay systems. Systems & Control Letters, 6, 7 7. Golubitsky, M. and Stewart, I. (4). The Symmetry Perspective - From Equilibrium to Chaos in Phase Space and Physical Space. Birkhäuser. Iwasaki, T. (8). Multivariable harmonic balance for central pattern generators. Automatica, 44(), Neefs, P., Steur, E., and Nijmeijer, H. (). Network complexity and synchronous behavior an experimental approach. Int. J. of Neural Systems,(3), Newman, M., Barabási, A.L., and Watts, D.J. (6). The Structure and Dynamics of Networks. Princeton University Press. Niculescu, S.I., Fu, P., and Chen, J. (998). Stability switches and reversalsof linear systems with commensurate delays: A matrix pencil characterization. IMA Journal of Mathematical Control & Information, 5, Oguchi, T., Nijmeijer, H., and Yamamoto, T. (8). Synchronization in networks of chaotic systems with time-delay coupling. Chaos, 8(3), 378. Oguchi, T., Yamamoto, T., and Nijmeijer, H. (9). Synchronization of bidirectionally coupled nonlinear systems with time-varying delay. In J.J. Loiseau et al.(ed.), q q q 3 q 4 q 5 Topics in Time Delay Systems: Analysis, Algorithms, and Control, volume 388 of LNCIS, Springerverlag. Pogromsky, A., Glad, T., and Nijmeijer, H. (999). On diffusion driven oscillations in coupled dynamical systems. International Journal of Bifurcation and Chaos, 9(4), Pogromsky, A., Santoboni, G., and Nijmeijer, H. (). Partial synchronization: from symmetry towards stability. Physica D, 7(), Appendix A. STABILITY FOR LINEAR RETARDED SYSTEMS Regarding to the characteristic roots on the imaginary, the following result holds regarding to the stability of linearretardedsystems (Chenet al.(995),niculescuet al. (998)). Consider the following linear delay system: m ẋ(t)=a x(t)+ A i x(t iτ) i= (A.) The corresponding characteristic function is given by ( m p(λ;τ)=det λi n A A i e λiτ) (A.) i= For the system (A.), the matrix pencil is defined by Λ(z):=zM +N (A.3) where M,N R (mn ) (mn ) are given by I n I n M :=... I n B m L n I n N :=... I n B m B m+ B m+ B m and B k for k =,...,m and B i for i =,...,m are defined as B k = I n A T k,b i = A i I n,b = A A T. The operator and denote the Kronecker product and sum. Then the following theorem holds. Theorem 3. (Niculescu et al. (998)) Assume that the matrix pencil Λ is regular. Then the quasi polynomial p(λ;τ)=has a crossing root on the imaginary axis for some positive delay value τ if and only if there exists a complex number z σ(λ) C(,) such that ( m ) jω σ A + A i zo i where σ(a) denotes a point spectrum of A.Furthermore, for some z satisfying the condition (ii) above, the set of delays corresponding to the induced crossing is given by: T (z )= i= { z ω + πl ω > :jω ( A + m } A i ),l Z i= 53
Invariance properties in the root sensitivity of time-delay systems with double imaginary roots
Invariance properties in the root sensitivity of time-delay systems with double imaginary roots Elias Jarlebring, Wim Michiels Department of Computer Science, KU Leuven, Celestijnenlaan A, 31 Heverlee,
More informationReliably computing all characteristic roots of delay differential equations in a given right half plane using a spectral method
Reliably computing all characteristic roots of delay differential equations in a given right half plane using a spectral method Zhen Wu Wim Michiels Report TW 596, May 2011 Katholieke Universiteit Leuven
More informationPolynomial two-parameter eigenvalue problems and matrix pencil methods for stability of delay-differential equations
Polynomial two-parameter eigenvalue problems and matrix pencil methods for stability of delay-differential equations Elias Jarlebring a, Michiel E. Hochstenbach b,1, a Technische Universität Braunschweig,
More informationPolynomial two-parameter eigenvalue problems and matrix pencil methods for stability of delay-differential equations
Polynomial two-parameter eigenvalue problems and matrix pencil methods for stability of delay-differential equations Elias Jarlebring a, Michiel E. Hochstenbach b,1, a Technische Universität Braunschweig,
More informationResearch Article Mathematical Model and Cluster Synchronization for a Complex Dynamical Network with Two Types of Chaotic Oscillators
Applied Mathematics Volume 212, Article ID 936, 12 pages doi:1.11/212/936 Research Article Mathematical Model and Cluster Synchronization for a Complex Dynamical Network with Two Types of Chaotic Oscillators
More informationZero controllability in discrete-time structured systems
1 Zero controllability in discrete-time structured systems Jacob van der Woude arxiv:173.8394v1 [math.oc] 24 Mar 217 Abstract In this paper we consider complex dynamical networks modeled by means of state
More informationRobust Entrainment to Natural Oscillations of Asymmetric Systems arising from Animal Locomotion
Joint 48th IEEE Conference on Decision and Control and 28th Chinese Control Conference Shanghai, P.R. China, December 16-18, 2009 WeCIn4.1 Robust Entrainment to Natural Oscillations of Asymmetric Systems
More informationLecture 5: Linear Systems. Transfer functions. Frequency Domain Analysis. Basic Control Design.
ISS0031 Modeling and Identification Lecture 5: Linear Systems. Transfer functions. Frequency Domain Analysis. Basic Control Design. Aleksei Tepljakov, Ph.D. September 30, 2015 Linear Dynamic Systems Definition
More informationComputing partial synchronization manifolds of delay-coupled systems
ENOC 2017 June 25-30 2017 Budapest Hungary Computing partial synchronization manifolds of delay-coupled systems L Su W Michiels E Steur and H Nijmeijer Department of Computer Science KU Leuven Institute
More informationLecture Note #6 (Chap.10)
System Modeling and Identification Lecture Note #6 (Chap.) CBE 7 Korea University Prof. Dae Ryoo Yang Chap. Model Approximation Model approximation Simplification, approximation and order reduction of
More informationSTABILITY ANALYSIS OF DELAY DIFFERENTIAL EQUATIONS WITH TWO DISCRETE DELAYS
CANADIAN APPLIED MATHEMATICS QUARTERLY Volume 20, Number 4, Winter 2012 STABILITY ANALYSIS OF DELAY DIFFERENTIAL EQUATIONS WITH TWO DISCRETE DELAYS XIHUI LIN AND HAO WANG ABSTRACT. We use an algebraic
More informationController Design for Delay-Independent Stability of Multiple Time-Delay Systems via Déscartes s Rule of Signs
Controller Design for Delay-Independent Stability of Multiple Time-Delay Systems via Déscartes s Rule of Signs Ismail Ilker Delice Rifat Sipahi,1 Department of Mechanical and Industrial Engineering, Northeastern
More informationLinearization of Differential Equation Models
Linearization of Differential Equation Models 1 Motivation We cannot solve most nonlinear models, so we often instead try to get an overall feel for the way the model behaves: we sometimes talk about looking
More informationSTUDY OF SYNCHRONIZED MOTIONS IN A ONE-DIMENSIONAL ARRAY OF COUPLED CHAOTIC CIRCUITS
International Journal of Bifurcation and Chaos, Vol 9, No 11 (1999) 19 4 c World Scientific Publishing Company STUDY OF SYNCHRONIZED MOTIONS IN A ONE-DIMENSIONAL ARRAY OF COUPLED CHAOTIC CIRCUITS ZBIGNIEW
More informationTHE EIGENVALUE PROBLEM
THE EIGENVALUE PROBLEM Let A be an n n square matrix. If there is a number λ and a column vector v 0 for which Av = λv then we say λ is an eigenvalue of A and v is an associated eigenvector. Note that
More informationADAPTIVE 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 informationConsensus Protocols for Networks of Dynamic Agents
Consensus Protocols for Networks of Dynamic Agents Reza Olfati Saber Richard M. Murray Control and Dynamical Systems California Institute of Technology Pasadena, CA 91125 e-mail: {olfati,murray}@cds.caltech.edu
More informationDesign of Oscillator Networks for Generating Signal with Prescribed Statistical Property
Journal of Physics: Conference Series PAPER OPEN ACCESS Design of Oscillator Networks for Generating Signal with Prescribed Statistical Property To cite this article: Tatsuo Yanagita 2017 J. Phys.: Conf.
More informationStability of Hybrid Control Systems Based on Time-State Control Forms
Stability of Hybrid Control Systems Based on Time-State Control Forms Yoshikatsu HOSHI, Mitsuji SAMPEI, Shigeki NAKAURA Department of Mechanical and Control Engineering Tokyo Institute of Technology 2
More informationMath Ordinary Differential Equations
Math 411 - Ordinary Differential Equations Review Notes - 1 1 - Basic Theory A first order ordinary differential equation has the form x = f(t, x) (11) Here x = dx/dt Given an initial data x(t 0 ) = x
More informationSynchronization of a General Delayed Complex Dynamical Network via Adaptive Feedback
Synchronization of a General Delayed Complex Dynamical Network via Adaptive Feedback Qunjiao Zhang and Junan Lu College of Mathematics and Statistics State Key Laboratory of Software Engineering Wuhan
More information21 Linear State-Space Representations
ME 132, Spring 25, UC Berkeley, A Packard 187 21 Linear State-Space Representations First, let s describe the most general type of dynamic system that we will consider/encounter in this class Systems may
More informationSolution via Laplace transform and matrix exponential
EE263 Autumn 2015 S. Boyd and S. Lall Solution via Laplace transform and matrix exponential Laplace transform solving ẋ = Ax via Laplace transform state transition matrix matrix exponential qualitative
More informationSynchronization and Bifurcation Analysis in Coupled Networks of Discrete-Time Systems
Commun. Theor. Phys. (Beijing, China) 48 (2007) pp. 871 876 c International Academic Publishers Vol. 48, No. 5, November 15, 2007 Synchronization and Bifurcation Analysis in Coupled Networks of Discrete-Time
More informationECE 388 Automatic Control
Lead Compensator and PID Control Associate Prof. Dr. of Mechatronics Engineeering Çankaya University Compulsory Course in Electronic and Communication Engineering Credits (2/2/3) Course Webpage: http://ece388.cankaya.edu.tr
More informationANALYSIS AND CONTROLLING OF HOPF BIFURCATION FOR CHAOTIC VAN DER POL-DUFFING SYSTEM. China
Mathematical and Computational Applications, Vol. 9, No., pp. 84-9, 4 ANALYSIS AND CONTROLLING OF HOPF BIFURCATION FOR CHAOTIC VAN DER POL-DUFFING SYSTEM Ping Cai,, Jia-Shi Tang, Zhen-Bo Li College of
More informationAdaptive synchronization of chaotic neural networks with time delays via delayed feedback control
2017 º 12 È 31 4 ½ Dec. 2017 Communication on Applied Mathematics and Computation Vol.31 No.4 DOI 10.3969/j.issn.1006-6330.2017.04.002 Adaptive synchronization of chaotic neural networks with time delays
More informationLecture Note #6 (Chap.10)
System Modeling and Identification Lecture Note #6 (Chap.) CBE 7 Korea University Prof. Dae Ryook Yang Chap. Model Approximation Model approximation Simplification, approximation and order reduction of
More informationDECENTRALIZED CONTROL DESIGN USING LMI MODEL REDUCTION
Journal of ELECTRICAL ENGINEERING, VOL. 58, NO. 6, 2007, 307 312 DECENTRALIZED CONTROL DESIGN USING LMI MODEL REDUCTION Szabolcs Dorák Danica Rosinová Decentralized control design approach based on partial
More informationPhase Synchronization
Phase Synchronization Lecture by: Zhibin Guo Notes by: Xiang Fan May 10, 2016 1 Introduction For any mode or fluctuation, we always have where S(x, t) is phase. If a mode amplitude satisfies ϕ k = ϕ k
More informationTime Response Analysis (Part II)
Time Response Analysis (Part II). A critically damped, continuous-time, second order system, when sampled, will have (in Z domain) (a) A simple pole (b) Double pole on real axis (c) Double pole on imaginary
More informationDependence of Delay Margin on Network Topology: Single Delay Case
Dependence of Delay Margin on Network Topology: Single Delay Case Wei Qiao Rifat Sipahi,1 Department of Mechanical and Industrial Engineering, Northeastern University, Boston, MA 02115 USA. e-mail: qiao.w@husky.neu.edu,
More informationSPEED-GRADIENT-BASED CONTROL OF POWER NETWORK: CASE STUDY
CYBERNETICS AND PHYSICS, VOL. 5, NO. 3, 2016, 85 90 SPEED-GRADIENT-BASED CONTROL OF POWER NETWORK: CASE STUDY Igor Furtat Control of Complex Systems ITMO University Russia cainenash@mail.ru Nikita Tergoev
More informationStability and Stabilization of Systems with Time Delay. Limitations and Opportunities
University of New Mexico UNM Digital Repository Electrical & Computer Engineering Technical Reports Engineering Publications 10-26-2010 Stability and Stabilization of Systems with Time Delay. Limitations
More informationAutomatic Control Systems theory overview (discrete time systems)
Automatic Control Systems theory overview (discrete time systems) Prof. Luca Bascetta (luca.bascetta@polimi.it) Politecnico di Milano Dipartimento di Elettronica, Informazione e Bioingegneria Motivations
More informationDESYNCHRONIZATION TRANSITIONS IN RINGS OF COUPLED CHAOTIC OSCILLATORS
Letters International Journal of Bifurcation and Chaos, Vol. 8, No. 8 (1998) 1733 1738 c World Scientific Publishing Company DESYNCHRONIZATION TRANSITIONS IN RINGS OF COUPLED CHAOTIC OSCILLATORS I. P.
More informationChaotic motion. Phys 420/580 Lecture 10
Chaotic motion Phys 420/580 Lecture 10 Finite-difference equations Finite difference equation approximates a differential equation as an iterative map (x n+1,v n+1 )=M[(x n,v n )] Evolution from time t
More informationModelling biological oscillations
Modelling biological oscillations Shan He School for Computational Science University of Birmingham Module 06-23836: Computational Modelling with MATLAB Outline Outline of Topics Van der Pol equation Van
More informationDelay-Dependent Stability Criteria for Linear Systems with Multiple Time Delays
Delay-Dependent Stability Criteria for Linear Systems with Multiple Time Delays Yong He, Min Wu, Jin-Hua She Abstract This paper deals with the problem of the delay-dependent stability of linear systems
More informationLecture 4: Analysis of MIMO Systems
Lecture 4: Analysis of MIMO Systems Norms The concept of norm will be extremely useful for evaluating signals and systems quantitatively during this course In the following, we will present vector norms
More informationUnderstanding interactions between networks controlling distinct behaviours: Escape and swimming in larval zebrafish
Understanding interactions between networks controlling distinct behaviours: Escape and swimming in larval zebrafish P. Dwight Kuo 1 and Chris Eliasmith 1,2 20th September 2003 Systems Design Engineering
More informationPeriod-One Rotating Solutions of Horizontally Excited Pendulum Based on Iterative Harmonic Balance
Advances in Pure Mathematics, 015, 5, 413-47 Published Online June 015 in Scies. http://www.scirp.org/journal/apm http://dx.doi.org/10.436/apm.015.58041 Period-One otating Solutions of Horizontally Excited
More informationSynchronization Transitions in Complex Networks
Synchronization Transitions in Complex Networks Y. Moreno 1,2,3 1 Institute for Biocomputation and Physics of Complex Systems (BIFI) University of Zaragoza, Zaragoza 50018, Spain 2 Department of Theoretical
More informationDynamics on a General Stage Structured n Parallel Food Chains
Memorial University of Newfoundland Dynamics on a General Stage Structured n Parallel Food Chains Isam Al-Darabsah and Yuan Yuan November 4, 2013 Outline: Propose a general model with n parallel food chains
More informationLesson 4: Non-fading Memory Nonlinearities
Lesson 4: Non-fading Memory Nonlinearities Nonlinear Signal Processing SS 2017 Christian Knoll Signal Processing and Speech Communication Laboratory Graz University of Technology June 22, 2017 NLSP SS
More informationarxiv: v1 [math.oc] 22 Jan 2008
arxiv:0801.3390v1 [math.oc] 22 Jan 2008 LQR-based coupling gain for synchronization of linear systems S. Emre Tuna tuna@eee.metu.edu.tr May 28, 2018 Abstract Synchronization control of coupled continuous-time
More informationControl Systems I. Lecture 5: Transfer Functions. Readings: Emilio Frazzoli. Institute for Dynamic Systems and Control D-MAVT ETH Zürich
Control Systems I Lecture 5: Transfer Functions Readings: Emilio Frazzoli Institute for Dynamic Systems and Control D-MAVT ETH Zürich October 20, 2017 E. Frazzoli (ETH) Lecture 5: Control Systems I 20/10/2017
More informationAnalysis of Neural Networks with Chaotic Dynamics
Chaos, Solitonr & Fructals Vol. 3, No. 2, pp. 133-139, 1993 Printed in Great Britain @60-0779/93$6.00 + 40 0 1993 Pergamon Press Ltd Analysis of Neural Networks with Chaotic Dynamics FRANCOIS CHAPEAU-BLONDEAU
More informationPatterned Linear Systems: Rings, Chains, and Trees
Patterned Linear Systems: Rings Chains and Trees Sarah C Hamilton and Mireille E Broucke Abstract In a first paper we studied system theoretic properties of patterned systems and solved classical control
More informationStability and hybrid synchronization of a time-delay financial hyperchaotic system
ISSN 76-7659 England UK Journal of Information and Computing Science Vol. No. 5 pp. 89-98 Stability and hybrid synchronization of a time-delay financial hyperchaotic system Lingling Zhang Guoliang Cai
More informationOscillatory Motion. Simple pendulum: linear Hooke s Law restoring force for small angular deviations. Oscillatory solution
Oscillatory Motion Simple pendulum: linear Hooke s Law restoring force for small angular deviations d 2 θ dt 2 = g l θ θ l Oscillatory solution θ(t) =θ 0 sin(ωt + φ) F with characteristic angular frequency
More informationOscillatory Motion. Simple pendulum: linear Hooke s Law restoring force for small angular deviations. small angle approximation. Oscillatory solution
Oscillatory Motion Simple pendulum: linear Hooke s Law restoring force for small angular deviations d 2 θ dt 2 = g l θ small angle approximation θ l Oscillatory solution θ(t) =θ 0 sin(ωt + φ) F with characteristic
More informationAnalysis of a lumped model of neocortex to study epileptiform ac
of a lumped model of neocortex to study epileptiform activity Sid Visser Hil Meijer Stephan van Gils March 21, 2012 What is epilepsy? Pathology Neurological disorder, affecting 1% of world population Characterized
More informationCalculus and Differential Equations II
MATH 250 B Second order autonomous linear systems We are mostly interested with 2 2 first order autonomous systems of the form { x = a x + b y y = c x + d y where x and y are functions of t and a, b, c,
More informationResearch Article Numerical Stability Test of Neutral Delay Differential Equations
Mathematical Problems in Engineering Volume 2008, Article ID 698043, 10 pages doi:10.1155/2008/698043 Research Article Numerical Stability Test of Neutral Delay Differential Equations Z. H. Wang Institute
More informationSynchronization of Limit Cycle Oscillators by Telegraph Noise. arxiv: v1 [cond-mat.stat-mech] 5 Aug 2014
Synchronization of Limit Cycle Oscillators by Telegraph Noise Denis S. Goldobin arxiv:148.135v1 [cond-mat.stat-mech] 5 Aug 214 Department of Physics, University of Potsdam, Postfach 61553, D-14415 Potsdam,
More informationSmall-Signal Stability Analysis of Neutral Delay Differential Equations
1 Small-Signal Stability Analysis of Neutral Delay Differential Equations Muyang Liu, Ioannis Dassios, and Federico Milano, Fellow, IEEE Abstract This paper focuses on the small-signal stability analysis
More informationResearch Article Stabilization Analysis and Synthesis of Discrete-Time Descriptor Markov Jump Systems with Partially Unknown Transition Probabilities
Research Journal of Applied Sciences, Engineering and Technology 7(4): 728-734, 214 DOI:1.1926/rjaset.7.39 ISSN: 24-7459; e-issn: 24-7467 214 Maxwell Scientific Publication Corp. Submitted: February 25,
More informationConstant Amplitude and Zero Autocorrelation Sequences and Single Pixel Camera Imaging
Constant Amplitude and Zero Autocorrelation Sequences and Single Pixel Camera Imaging Mark Magsino mmagsino@math.umd.edu Norbert Wiener Center for Harmonic Analysis and Applications Department of Mathematics
More informationNeural Excitability in a Subcritical Hopf Oscillator with a Nonlinear Feedback
Neural Excitability in a Subcritical Hopf Oscillator with a Nonlinear Feedback Gautam C Sethia and Abhijit Sen Institute for Plasma Research, Bhat, Gandhinagar 382 428, INDIA Motivation Neural Excitability
More informationarxiv: v1 [nlin.ao] 8 Dec 2017
Nonlinearity 31(1):R1-R23 (2018) Identical synchronization of nonidentical oscillators: when only birds of different feathers flock together arxiv:1712.03245v1 [nlin.ao] 8 Dec 2017 Yuanzhao Zhang 1 and
More information2nd Symposium on System, Structure and Control, Oaxaca, 2004
263 2nd Symposium on System, Structure and Control, Oaxaca, 2004 A PROJECTIVE ALGORITHM FOR STATIC OUTPUT FEEDBACK STABILIZATION Kaiyang Yang, Robert Orsi and John B. Moore Department of Systems Engineering,
More informationAutonomous Helicopter Landing A Nonlinear Output Regulation Perspective
Autonomous Helicopter Landing A Nonlinear Output Regulation Perspective Andrea Serrani Department of Electrical and Computer Engineering Collaborative Center for Control Sciences The Ohio State University
More informationPhysics 106b: Lecture 7 25 January, 2018
Physics 106b: Lecture 7 25 January, 2018 Hamiltonian Chaos: Introduction Integrable Systems We start with systems that do not exhibit chaos, but instead have simple periodic motion (like the SHO) with
More informationOn the convergence of normal forms for analytic control systems
Chapter 1 On the convergence of normal forms for analytic control systems Wei Kang Department of Mathematics, Naval Postgraduate School Monterey, CA 93943 wkang@nps.navy.mil Arthur J. Krener Department
More informationAn efficient algorithm to compute the real perturbation values of a matrix
An efficient algorithm to compute the real perturbation values of a matrix Simon Lam and Edward J. Davison 1 Abstract In this paper, an efficient algorithm is presented for solving the nonlinear 1-D optimization
More informationThe Tuning of Robust Controllers for Stable Systems in the CD-Algebra: The Case of Sinusoidal and Polynomial Signals
The Tuning of Robust Controllers for Stable Systems in the CD-Algebra: The Case of Sinusoidal and Polynomial Signals Timo Hämäläinen Seppo Pohjolainen Tampere University of Technology Department of Mathematics
More informationNonsmooth systems: synchronization, sliding and other open problems
John Hogan Bristol Centre for Applied Nonlinear Mathematics, University of Bristol, England Nonsmooth systems: synchronization, sliding and other open problems 2 Nonsmooth Systems 3 What is a nonsmooth
More informationOn the simultaneous diagonal stability of a pair of positive linear systems
On the simultaneous diagonal stability of a pair of positive linear systems Oliver Mason Hamilton Institute NUI Maynooth Ireland Robert Shorten Hamilton Institute NUI Maynooth Ireland Abstract In this
More informationEE482: Digital Signal Processing Applications
Professor Brendan Morris, SEB 3216, brendan.morris@unlv.edu EE482: Digital Signal Processing Applications Spring 2014 TTh 14:30-15:45 CBC C222 Lecture 05 IIR Design 14/03/04 http://www.ee.unlv.edu/~b1morris/ee482/
More informationStabilization of Higher Periodic Orbits of Discrete-time Chaotic Systems
ISSN 749-3889 (print), 749-3897 (online) International Journal of Nonlinear Science Vol.4(27) No.2,pp.8-26 Stabilization of Higher Periodic Orbits of Discrete-time Chaotic Systems Guoliang Cai, Weihuai
More informationV.B. LARIN 1. Keywords: unilateral quadratic matrix equation, matrix sign function, parameters updating.
TWMS Jour. Pure Appl. Math., V.3, N.2, 2012, pp.202-209 THE UNILATERAL QUADRATIC MATRIX EQUATION AND PROBLEM OF UPDATING OF PARAMETERS OF MODEL V.B. LARIN 1 Abstract. The algorithm of construction of solutions
More informationStability Theory for Nonnegative and Compartmental Dynamical Systems with Time Delay
1 Stability Theory for Nonnegative and Compartmental Dynamical Systems with Time Delay Wassim M. Haddad and VijaySekhar Chellaboina School of Aerospace Engineering, Georgia Institute of Technology, Atlanta,
More informationEigenstructure Assignment for Helicopter Hover Control
Proceedings of the 17th World Congress The International Federation of Automatic Control Eigenstructure Assignment for Helicopter Hover Control Andrew Pomfret Stuart Griffin Tim Clarke Department of Electronics,
More informationHOPF BIFURCATION WITH S 3 -SYMMETRY
PORTUGALIAE MATHEMATICA Vol. 63 Fasc. 006 Nova Série HOPF BIFURCATION WITH S 3 -SYMMETRY Ana Paula S. Dias and Rui C. Paiva Recommended by José Basto-Gonçalves Abstract: The aim of this paper is to study
More informationHomework Set 5 Solutions
MATH 672-010 Vector Spaces Prof. D. A. Edwards Due: Oct. 7, 2014 Homework Set 5 Solutions 1. Let S, T L(V, V finite-dimensional. (a (5 points Prove that ST and T S have the same eigenvalues. Solution.
More informationMaster Thesis. Equivariant Pyragas control. Isabelle Schneider Freie Universität Berlin Fachbereich Mathematik und Informatik Berlin
Master Thesis Equivariant Pyragas control Isabelle Schneider Freie Universität Berlin Fachbereich Mathematik und Informatik 14195 Berlin February 3, 2014 Abstract In this thesis we show how the concept
More informationDigital Control Engineering Analysis and Design
Digital Control Engineering Analysis and Design M. Sami Fadali Antonio Visioli AMSTERDAM BOSTON HEIDELBERG LONDON NEW YORK OXFORD PARIS SAN DIEGO SAN FRANCISCO SINGAPORE SYDNEY TOKYO Academic Press is
More informationA MATRIX INEQUALITY BASED DESIGN METHOD FOR CONSENSUS PROBLEMS IN MULTI AGENT SYSTEMS
Int. J. Appl. Math. Comput. Sci., 2009, Vol. 19, No. 4, 639 646 DOI: 10.2478/v10006-009-0051-1 A MATRIX INEQUALITY BASED DESIGN METHOD FOR CONSENSUS PROBLEMS IN MULTI AGENT SYSTEMS GUISHENG ZHAI, SHOHEI
More informationStability crossing curves of SISO systems controlled by delayed output feedback
Accepted to Dynamics of Continuous, Discrete and Impulsive Systems http:monotone.uwaterloo.ca/ journal Stability crossing curves of SISO systems controlled by delayed output feedback Constantin-Irinel
More informationChaotic motion. Phys 750 Lecture 9
Chaotic motion Phys 750 Lecture 9 Finite-difference equations Finite difference equation approximates a differential equation as an iterative map (x n+1,v n+1 )=M[(x n,v n )] Evolution from time t =0to
More informationObservations on the Stability Properties of Cooperative Systems
1 Observations on the Stability Properties of Cooperative Systems Oliver Mason and Mark Verwoerd Abstract We extend two fundamental properties of positive linear time-invariant (LTI) systems to homogeneous
More informationStability and Hopf Bifurcation for a Discrete Disease Spreading Model in Complex Networks
International Journal of Difference Equations ISSN 0973-5321, Volume 4, Number 1, pp. 155 163 (2009) http://campus.mst.edu/ijde Stability and Hopf Bifurcation for a Discrete Disease Spreading Model in
More informationAssignment # 4. K-WTA neural networks are relatively simple. N neurons, with all the connection strengths equal to -1 (with no self-connections).
Farrukh Jabeen COM 58 Mathematical Methods In AI Assignment # 4 This problem asks you to analyze a K-WTA network (K winner take all). First read the K-WTA reference article, then answer the questions below.
More informationCPG CONTROL OF A TENSEGRITY MORPHING STRUCTURE FOR BIOMIMETIC APPLICATIONS
CPG CONTROL OF A TENSEGRITY MORPHING STRUCTURE FOR BIOMIMETIC APPLICATIONS T. K. BLISS, T. IWASAKI, and H. BART-SMITH University of Virginia Department of Mechanical and Aerospace Engineering 122 Engineer
More informationSection 3: Complex numbers
Essentially: Section 3: Complex numbers C (set of complex numbers) up to different notation: the same as R 2 (euclidean plane), (i) Write the real 1 instead of the first elementary unit vector e 1 = (1,
More informationMath 331 Homework Assignment Chapter 7 Page 1 of 9
Math Homework Assignment Chapter 7 Page of 9 Instructions: Please make sure to demonstrate every step in your calculations. Return your answers including this homework sheet back to the instructor as a
More informationSimple Chaotic Oscillator: From Mathematical Model to Practical Experiment
6 J. PERŽELA, Z. KOLKA, S. HANUS, SIMPLE CHAOIC OSCILLAOR: FROM MAHEMAICAL MODEL Simple Chaotic Oscillator: From Mathematical Model to Practical Experiment Jiří PERŽELA, Zdeněk KOLKA, Stanislav HANUS Dept.
More informationStabilization 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 informationarxiv: v1 [nlin.ao] 3 May 2015
UNI-DIRECTIONAL SYNCHRONIZATION AND FREQUENCY DIFFERENCE LOCKING INDUCED BY A HETEROCLINIC RATCHET arxiv:155.418v1 [nlin.ao] 3 May 215 Abstract A system of four coupled oscillators that exhibits unusual
More informationDiscrete Time Signals and Systems Time-frequency Analysis. Gloria Menegaz
Discrete Time Signals and Systems Time-frequency Analysis Gloria Menegaz Time-frequency Analysis Fourier transform (1D and 2D) Reference textbook: Discrete time signal processing, A.W. Oppenheim and R.W.
More information8 Ecosystem stability
8 Ecosystem stability References: May [47], Strogatz [48]. In these lectures we consider models of populations, with an emphasis on the conditions for stability and instability. 8.1 Dynamics of a single
More informationNETWORK control is of utmost importance in many
Exploiting nodes symmetries to control synchronization and consensus patterns in multiagent systems Davide Fiore, Student Member, IEEE, Giovanni Russo, Member, IEEE, and Mario di Bernardo, Fellow, IEEE
More informationStability and bifurcation of a simple neural network with multiple time delays.
Fields Institute Communications Volume, 999 Stability and bifurcation of a simple neural network with multiple time delays. Sue Ann Campbell Department of Applied Mathematics University of Waterloo Waterloo
More informationBoundedness of the Lorenz Stenflo system
Manuscript Click here to download Manuscript: Boundedness of the Lorenz-Stenflo system.tex Boundedness of the Lorenz Stenflo system Zeraoulia Elhadj 1, J. C. Sprott 1 Department of Mathematics, University
More informationGlocal Control for Hierarchical Dynamical Systems
Lecture Series, TU Munich October 22, 29 & November 5, 2013 Glocal Control for Hierarchical Dynamical Systems Theoretical Foundations with Applications in Energy Networks Shinji HARA The University of
More informationConsensus of Multi-Agent Systems with
Consensus of Multi-Agent Systems with 1 General Linear and Lipschitz Nonlinear Dynamics Using Distributed Adaptive Protocols arxiv:1109.3799v1 [cs.sy] 17 Sep 2011 Zhongkui Li, Wei Ren, Member, IEEE, Xiangdong
More informationCONTROLLING IN BETWEEN THE LORENZ AND THE CHEN SYSTEMS
International Journal of Bifurcation and Chaos, Vol. 12, No. 6 (22) 1417 1422 c World Scientific Publishing Company CONTROLLING IN BETWEEN THE LORENZ AND THE CHEN SYSTEMS JINHU LÜ Institute of Systems
More informationA SYSTEMATIC APPROACH TO GENERATING n-scroll ATTRACTORS
International Journal of Bifurcation and Chaos, Vol. 12, No. 12 (22) 297 2915 c World Scientific Publishing Company A SYSTEMATIC APPROACH TO ENERATIN n-scroll ATTRACTORS UO-QUN ZHON, KIM-FUN MAN and UANRON
More informationOn Positive Real Lemma for Non-minimal Realization Systems
Proceedings of the 17th World Congress The International Federation of Automatic Control On Positive Real Lemma for Non-minimal Realization Systems Sadaaki Kunimatsu Kim Sang-Hoon Takao Fujii Mitsuaki
More information