Claude Elysée Lobry Université de Nice, Faculté des Sciences, parc Valrose, NICE, France.

Similar documents
D.S.G. POLLOCK: TOPICS IN TIME-SERIES ANALYSIS STATISTICAL FOURIER ANALYSIS

Chapter 3.1: Polynomial Functions

Multi-objective Programming Approach for. Fuzzy Linear Programming Problems

Intermediate Division Solutions

Ch. 1 Introduction to Estimation 1/15

Quantum Mechanics for Scientists and Engineers. David Miller

Study of Energy Eigenvalues of Three Dimensional. Quantum Wires with Variable Cross Section

MATH Midterm Examination Victor Matveev October 26, 2016

are specified , are linearly independent Otherwise, they are linearly dependent, and one is expressed by a linear combination of the others

ENGI 4421 Central Limit Theorem Page Central Limit Theorem [Navidi, section 4.11; Devore sections ]

Fourier Series & Fourier Transforms

5.1 Two-Step Conditional Density Estimator

Every gas consists of a large number of small particles called molecules moving with very high velocities in all possible directions.

Fourier Method for Solving Transportation. Problems with Mixed Constraints

5.80 Small-Molecule Spectroscopy and Dynamics

A New Method for Finding an Optimal Solution. of Fully Interval Integer Transportation Problems

Author. Introduction. Author. o Asmir Tobudic. ISE 599 Computational Modeling of Expressive Performance

K [f(t)] 2 [ (st) /2 K A GENERALIZED MEIJER TRANSFORMATION. Ku(z) ()x) t -)-I e. K(z) r( + ) () (t 2 I) -1/2 e -zt dt, G. L. N. RAO L.

BIO752: Advanced Methods in Biostatistics, II TERM 2, 2010 T. A. Louis. BIO 752: MIDTERM EXAMINATION: ANSWERS 30 November 2010

, the random variable. and a sample size over the y-values 0:1:10.

The Excel FFT Function v1.1 P. T. Debevec February 12, The discrete Fourier transform may be used to identify periodic structures in time ht.

Axial Temperature Distribution in W-Tailored Optical Fibers

A Study on Estimation of Lifetime Distribution with Covariates Under Misspecification

Markov processes and the Kolmogorov equations

MATHEMATICS 9740/01 Paper 1 14 Sep hours

Solutions. Definitions pertaining to solutions

Distributed Trajectory Generation for Cooperative Multi-Arm Robots via Virtual Force Interactions

ENGI 4421 Central Limit Theorem Page Central Limit Theorem [Navidi, section 4.11; Devore sections ]

IJISET - International Journal of Innovative Science, Engineering & Technology, Vol. 2 Issue 12, December

Mean residual life of coherent systems consisting of multiple types of dependent components

Lecture 18. MSMPR Crystallization Model

UNIVERSITY OF TECHNOLOGY. Department of Mathematics PROBABILITY THEORY, STATISTICS AND OPERATIONS RESEARCH GROUP. Memorandum COSOR 76-10

Comparative analysis of bayesian control chart estimation and conventional multivariate control chart

Solutions to Midterm II. of the following equation consistent with the boundary condition stated u. y u x y

ON FREE RING EXTENSIONS OF DEGREE N

On the affine nonlinearity in circuit theory

x 2 x 3 x b 0, then a, b, c log x 1 log z log x log y 1 logb log a dy 4. dx As tangent is perpendicular to the x axis, slope

The Molecular Diffusion of Heat and Mass from Two Spheres

Active redundancy allocation in systems. R. Romera; J. Valdés; R. Zequeira*

A Hartree-Fock Calculation of the Water Molecule

The Complexity of Translation Membership for Macro Tree Transducers

Physical Chemistry Laboratory I CHEM 445 Experiment 2 Partial Molar Volume (Revised, 01/13/03)

Preliminary Test Single Stage Shrinkage Estimator for the Scale Parameter of Gamma Distribution

An S-type upper bound for the largest singular value of nonnegative rectangular tensors

Electrostatics. . where,.(1.1) Maxwell Eqn. Total Charge. Two point charges r 12 distance apart in space

Grade 3 Mathematics Course Syllabus Prince George s County Public Schools

Design and Implementation of Cosine Transforms Employing a CORDIC Processor

Study in Cylindrical Coordinates of the Heat Transfer Through a Tow Material-Thermal Impedance

Function representation of a noncommutative uniform algebra

Matching a Distribution by Matching Quantiles Estimation

Fluctuating Groundwater Flow in Homogeneous Finite Porous Domain

AP Statistics Notes Unit Eight: Introduction to Inference

Copyright 1978, by the author(s). All rights reserved.

Review for cumulative test

Lecture 21: Signal Subspaces and Sparsity

Portfolio Performance Evaluation in a Modified Mean-Variance-Skewness Framework with Negative Data

Christensen, Mads Græsbøll; Vera-Candeas, Pedro; Somasundaram, Samuel D.; Jakobsson, Andreas

Review of Important Concepts

Chapter 4. Problem Solutions

The generalized marginal rate of substitution

Identical Particles. We would like to move from the quantum theory of hydrogen to that for the rest of the periodic table

Declarative approach to cyclic steady state space refinement: periodic process scheduling

CHAPTER 24: INFERENCE IN REGRESSION. Chapter 24: Make inferences about the population from which the sample data came.

[1 & α(t & T 1. ' ρ 1

General Chemistry 1 (CHEM1141) Shawnee State University Fall 2016

Unifying the Derivations for. the Akaike and Corrected Akaike. Information Criteria. from Statistics & Probability Letters,

On the structure of space-time and matter as obtained from the Planck scale by period doubling in three and four dimensions

Hº = -690 kj/mol for ionization of n-propylene Hº = -757 kj/mol for ionization of isopropylene

A Simplified Nonlinear Generalized Maxwell Model for Predicting the Time Dependent Behavior of Viscoelastic Materials

Modelling of Clock Behaviour. Don Percival. Applied Physics Laboratory University of Washington Seattle, Washington, USA

RMO Sample Paper 1 Solutions :

Result on the Convergence Behavior of Solutions of Certain System of Third-Order Nonlinear Differential Equations

Frequency-Domain Study of Lock Range of Injection-Locked Non- Harmonic Oscillators

ALE 26. Equilibria for Cell Reactions. What happens to the cell potential as the reaction proceeds over time?

ABOUT CHAOS AND SENSITIVITY IN TOPOLOGICAL DYNAMICS

WEST VIRGINIA UNIVERSITY

Recovery of Third Order Tensors via Convex Optimization

ROOT FINDING FOR NONLINEAR EQUATIONS

MODIFIED LEAKY DELAYED LMS ALGORITHM FOR IMPERFECT ESTIMATE SYSTEM DELAY

Dynamic Response of Second Order Mechanical Systems with Viscous Dissipation forces

EECE 301 Signals & Systems

Directional Duality Theory

On natural cubic splines, with an application to numerical integration formulae Schurer, F.

The generation of successive approximation methods for Markov decision processes by using stopping times

STRUCTURES IN MIKE 21. Flow over sluice gates A-1

Homework 2. Show that if h is a bounded sesquilinear form on the Hilbert spaces X and Y, then h has the representation

E o and the equilibrium constant, K

International Journal of Solids and Structures

MATHEMATICS SYLLABUS SECONDARY 5th YEAR

Complex-Analytic Theory of the -Function*

This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and

MULTIDIMENSIONAL EXTENICS THEORY

Computational modeling techniques

(b) y(t) is not periodic although sin t and 4 cos 2πt are independently periodic.

Abstract: The asympttically ptimal hypthesis testig prblem with the geeral surces as the ull ad alterative hyptheses is studied uder expetial-type err

Journal of Applied and Computational Mechanics, Vol. 3, No. 1, (2017), DOI: /jacm

ACTIVE FILTERS EXPERIMENT 2 (EXPERIMENTAL)

ESWW-2. Israeli semi-underground great plastic scintillation multidirectional muon telescope (ISRAMUTE) for space weather monitoring and forecasting

Cold mirror based on High-Low-High refractive index dielectric materials

THE ASYMPTOTIC PERFORMANCE OF THE LOG LIKELIHOOD RATIO STATISTIC FOR THE MIXTURE MODEL AND RELATED RESULTS

Transcription:

CHAOS AND CELLULAR AUTOMATA Claude Elysée Lbry Uiversité de Nice, Faculté des Scieces, parc Valrse, 06000 NICE, Frace. Keywrds: Chas, bifurcati, cellularautmata, cmputersimulatis, dyamical system, ifectius disease. Ctets 1. Chas 1.1. Itrducti 1.2. Oe Dimesial Discrete Chatic Systems 1. 3. Tw Dimesial Discrete Chatic System 1.4. The Cascade f Bifurcati 1.5. Ctiuus Chatic Systems 1.6. The Stry f Chas i Dyamical Systems 2. Cellular autmata 2.1. Itrducti 2.2. The S.E.I.R. Mdel fr Epidemics 2.3. Aalysis f the I.R.N autmat 2.3.1. Defiitis 2.3.2. Sme simulatis 2.4. A partial Thery fr the (I, R, N, G) Autmat 2.5. The Need fr a Thery whe e uses Cellular Autmata i Mdelig. 3. Cclusi Glssary Bibligraphy Bigraphical Sketch Summary I this chapter, we are ccered with Chas thery ad Cellular Autmata thery. These tw theries have i cmm the fact that they were very ppular sme time ag ad the fact that they are itimately cected t the develpmet f digital cmputig. The chapter is divided i tw parts the first e develped t chas, the secd e devted t cellular autmata. We d t try t make a survey f all aspects f these theries. We just ccetrate few aspects which seem t be relevat t mdelig ad give cautiary advises relative t these theries i the prcess f mdelig. 1. Chas 1.1. Itrducti Cmpare the tw sigals i Figure 1a:

Figure 1a: Tw sigals They lk very differet; the e belw is icely peridic sice the e abve has apparetly regularity. It lks radm. But it is t! Bth sigals are btaied by iterati f the very simple scheme: (1) U U ( U ) + 1 =λ 1 the first e with λ= 3.8 the secd e with λ = 3.5. This is a defiitely determiistic prcess! The first sigal is called chatic because it is differet frm a pure radm sigal. T make the differece betwee chas ad pure radm we cmpare the chatic sigal abve with the sigal btaied by plttig successively the result f a call t the radm fucti f a cmputer. Oe btais the fucti shw i Figure 1b: Figure 1b: A radm fucti This sigal lks similar t the chatic e but if we plt, t the value f the sigal agaist time, but the value at time + 1 agaist the value at time we btai the fucti i Figure 2: Figure 2: Chatic ad radm sigals This illustrates the essetial differece betwee chatic sigals ad radm sigals. The

ideal mathematical ccept f radmess requires that tw successive ccurrece f U must be idepedet, which meas that the prbability f appearace f ay value f 1 U + is t affected by the kwledge f U. The visual csequece f this is the apparetly uifrm repartiti the square f the ccurreces f the pits f crdiates ( U+ 1, U). Actually defiiti f perfect radmess is mre striget ad requires that all fiite sequeces are mutually idepedet, but we d t g further i this directi. Cmpletely differet is the chatic sigal. Sice the ext U + 1 is related t U by the fuctial relati: ( ) U + 1 = 3.8 U 1 U it is t surprisig at all that all the pit f crdiates ( U U ) graph f the mappig: ( ) U 3.8 U 1 U + 1, are ctaied i the Ntice that the radm prcess f ur cmputer is t able t realize a true mathematical pure radm prcess, but ly a apprximati f it by sme determiistic device. 1.2. Oe Dimesial Discrete Chatic Systems Csider a mappig f frm sme buded subset K f dyamical system defied by: + 1 = ( ) x f x ad dete by ( ) R it itself ad csider the x x the trajectry issued frm x (this meas the sequece issued frm x ad defied by iducti with the abve frmula. Defiiti : Sesitivity t iitial cditis Oe says that the system pssesses the prperty f sesitivity t iitial cditis sme subset A f K if there exist sme psitive umber M such that fr every x i A ad fr every (small) psitive ε e ca fid sme iitial cditi x ad a iteger such that: ( ) ( ) x x x x M We illustrate this by the fllwig experimets with the system (1) fr the value 0.8 f

the parameter. We have pltted the differece U ( U 0) U ( U 0 ) differig successively frm 6 9 12 10, 10, 10 i Figure 3. fr iitial cditis Figure 3: The fucti U ( U 0) U ( U 0 ) fr differet iitial cditis We see that fr few iterates (15,25,35) the tw trajectries are impssible t distiguish (up t the precisi f the pixel f the scree) ad the, abruptly, the differece becmes visible at upredictable. The differece grws expetially which meas that t add sme cstat umber t the umber f iterates where the trajectries are udistiguishable e has t divide by sme cstat factr the distace f the tw iitial cditis; i ur example e have t divide by 1000 t btai te mre iteratis. Defiiti : We say that system (2) is chatic if it has a trajectry which is dese a subset A which is a attractr, if it has the prperty f sesitivity t iitial cditi subset A ad, mrever, peridic rbits are dese i A. Because f the sesibility t iitial cditi, what is the meaig f ur umerical cmputatis is t clear. All what we kw is that ur cmputer simulati has thig t d with the actual trajectry! Frtuately there is a deep mathematical result which state (i a precise mathematical way) that every cmputer simulati is actually clse t the true trajectry issued frm a iitial cditi which is t the e we used i the cmputati but is clse t it. By the way, i a chatic system, it is impssible t predict the trajectry issued frm a iitial cditi but what we cmpute is a typical pssible utcme f the system. 1. 3. Tw Dimesial Discrete Chatic System S far we have csidered a e dimesial chatic system. A example f a famus chatic system i tw dimesis is the Hé mappig, prpsed i 1976. Csider the mappig: ( x, y) ( y + 1 ax 2, bx) ad the iteratis:

x = y + 1 ax + 1 y + 1 = bx Figure 4: Tw successive elargemets fr the iterates f the Hé mappig 7 fr the values a = 1, 4 ad b = 0.3. Figure 4 shws 10 iterates f this mappig startig frm sme iitial cditi. Each elargemet is a cpy f the previus picture. This is ather feature f chatic systems. The fractal ature f each trajectry. - - - Bibligraphy TO ACCESS ALL THE 19 PAGES OF THIS CHAPTER, Visit: http://www.elss.et/elss-sampleallchapter.aspx Alluche J-P, ad C. Reder (1984), Oscillatis spati-temprelles egedrées par u autmate cellulaire. Discrete Applied Mathematics 8 : 215-254. Cartwright M.L. ad J.E. Littlewd. (1945) O liear equatis f the secd rder; the equati 2 y k 1 y y + y = bλk cs λ t+α, k large. J; Ld Math. Sc. 20: 180-189 ( ) ( ) Checier A. (1989). Systèmes Dyamiques Différetiables. Ecycpedia Uiveralis. Emetrut G-B ad L. Edelstei-Keshet (1993), Cellular Autmata Apprach t Bilgical Mdelig. Jural f Theretical Bilgy 160 : 97-133. Feigebaum M. (1978), Quatitative uiversality fr a class f liear trasfrmatis. J. Stat. Physics. 19: 25-52. Lbry C. a H. Elmzi (2000). Cmbiatrial Prperties f sme Cellular Autmata Related t the Msaïc Cycle Ccept. Acta Bitheretica 48: 219-242. Lrez E. (1963). J. Atms. Sci. 20: 130-141.

May R. (1976). Simple mathematical mdels with very cmplicated dyamics, Nature 261, 459-467 Picaré H. (Euvres cmplètes. Gauthier- Vilars 1951 Smale S. (1976). Differetial Dyamical Systems. Bull. A.M.S. Wieer N., ad A. Rseblueth (1946). The mathematical frmulati f the prblem f cducti f impulses i a etwrk f cected excitable elemets. Archives f the Istitute f Cardilgy f Mexic 16 : 202-2065. Wlfram S. (1994) Cellular Autmata ad Cmplexity. Addis-Wesley Bigraphical Sketch Claude Lbry, was br i 1943 ad gt a thèse d état at the uiversity f Greble, i ptimal ctrl, i 1972. He was the appited as prfessr at the Uiversity f Brdeaux ad jied the Uiversity f Nice i 1981. His mai iterest is mathematics applied t autmatic ctrl ad atural systems where he became a specialist f ctrllability. Durig the eighties he used the tls f Nstadard Aalysis i sme prblems f sigular perturbatis f differetial systems.frm 1990 t w he has bee very active i creatig iterdiscipliary teams r etwrks with the Cetre Natial de la Recherche Scietifique, the Istitut Natial de Recherche Agrmique, the Istitut Natial de Recherche e Ifrmatique et Autmatique. Besides his scietific activities he has a strg ivlvemet i the prmti f mathematical research i develpig cutries, especially Africa.