Power Spectral Analysis of Elementary Cellular Automata
|
|
- Gavin Weaver
- 6 years ago
- Views:
Transcription
1 Power Spectral Analysis o Elementary Cellular Automata Shigeru Ninagawa Division o Inormation and Computer Science, Kanazawa Institute o Technology, 7- Ohgigaoka, Nonoichi, Ishikawa , Japan Spectral analysis o elementary cellular automata is perormed. A power spectrum is calculated rom the evolution o 88 independent rules starting rom random initial conigurations. As a result, it is ound that rule 0 exhibits / noise during the longest time steps. Rule 0 has proved to be capable o supporting universal computation. These results suggest that there is a relationship between computational universality and / noise in cellular automata.. Introduction Cellular automata (CAs) are spatially and temporally discrete dynamical systems with large degrees o reedom. Since spectral analysis is one useul method or investigating the behavior o dynamical systems [], it is reasonable to apply it to studying the behavior o CAs. In this paper we deal with elementary CAs (ECAs), namely one-dimensional and twostate, three-neighbor CAs. ECAs have been investigated in detail or their simplicity ([2], appendix in [3]). Although the ECA rule space is not large, there is a wide variety o rules which exhibit regular, chaotic, or complex behavior. Spectral analysis was used or the investigation o the spatial structures produced by ECAs in connection with regular languages in [4], whereas we ocus on the temporal behavior o ECAs. 2. Spectral analysis o elementary cellular automata Let x i (t) be the value o site i at time step t in an ECA. The value o each site is speciied as 0 or. The site value evolves by iterating the mapping x i (t ) F(x i (t), x i (t), x i (t)). () Here F is an arbitrary unction speciying the ECA rule. The ECA rule is determined by a binary sequence with length 2 3 8, F(,,), F(,,0),..., F(0,0,0). (2) Electronic mail address: ninagawa@inor.kanazawa-it.ac.jp. ; 2008 Complex Systems Publications, Inc.
2 400 S. Ninagawa Thereore the total number o possible distinct ECA rules is and each rule is abbreviated by the decimal representation o the binary sequence (2) as used in [2]. A parameter Λ is deined as the density o ones in the binary sequence (2) [5]. Out o the 256 ECA rules, 88 o them remain independent (appendix in [6]). With a ew exceptions, we use as representative rules those with smaller Λ values and the smaller decimal rule numbers adopted in [6]. The discrete Fourier transorm o the time series o state x i (t) or t 0,,..., T isgivenby T ˆx i ( ) x T i (t)exp i 2Πt T t 0. (3) It is natural to deine the power spectrum o a CA as N S( ) ˆx N i ( ) 2, (4) i where N means array size and the summation is taken over all cells in the array. The power S( ) at requency intuitively means the strength o the periodic vibration with period T/ in the evolution o observation length T. Another deinition o power spectrum is given by calculating spatial density beore the Fourier transorm. It might be useul to investigate the behavior o CAs with an oscillating density in time, although we do not ocus on them in this research. Since we are interested in the dynamics that a CA rule brings, rather than the behavior rom a particular initial coniguration which is elaborately designed, throughout this paper we use random initial conigurations in which each site takes state 0 or randomly with independent equal probabilities. Generally, the evolution o one-dimensional CAs tends to depend more heavily on an initial coniguration than that o two-dimensional CAs. So we use some dierent initial conigurations to grasp the typical behavior o a CA. In addition, an array size larger than 400 seems required to avoid singular behavior. We adopt periodic boundary conditions where each end o the array is connected like a ring. Spectral analysis revealed that 88 independent ECAs can be classiied into several categories according to the shape o their power spectra. Theresultsareasollows. 2. Category : Extremely low power density The power spectra in this category are characterized by an extremely low power density at almost all requencies. Figure shows the spacetime patterns (let) and the power spectra (right) o rule 32 (top) and rule
3 Power Spectral Analysis o Elementary Cellular Automata 40 Rule #32 S() Rule #08 S() Figure. Space-time patterns (let) and power spectra (right) o rule 32 (top) and rule 08 (bottom). The space-time patterns consist o 00 cells or 00 time steps. The power spectra are calculated rom the evolution o 700 cells or 024 time steps. 08 (bottom). Space-time patterns show conigurations obtained at successive time steps on successive horizontal lines in which black squares represent sites with value, white squares value 0. The power spectrum is calculated rom the evolution o 700 cells or 024 time steps. Only hal o the components o the spectrum are shown since the other hal are redundant. The y-axis is plotted on a logarithmic scale. The power density in the power spectrum o rule 32 is extremely low at all requencies except 0. This is ascribed to the initial transient behavior which vanishes ater the irst ew time steps. Other rules with similar power spectra are rules 0, 4, 8, 2, 3, 26, 40, 44, 72, 76, 77, 78, 04, 28, 32, 36, 40, 60, 64, 68, 72, 200, and 232. These rules are considered to be in class I, or a part o class II, with evolutions that lead to stable structures according to Wolram s classiication scheme in [7]. The power spectrum o rule 08 has a peak at 52 which is caused by periodic structures with period 2. Rules like this are a part o class II with evolutions that lead to periodic structures, namely rules, 5, 9, 23, 28, 29, 33, 37, 50, 5, 33, 56, and 78.
4 402 S. Ninagawa 2.2 Category 2: Broad-band noise The power spectra in this category are characterized by broad-band noise. In general, broad-band noise in a power spectrum means that the original time series is nonperiodic. This category is divided into two subcategories. A typical example o the irst subcategory 2-A is shown in the top o Figure 2. The power spectrum o rule 6 has a higher power density at all requencies compared with those in category. The space-time pattern o rule 6 yields uniormly shiting structures which move one cell let every time step. So the evolution with array size N becomes periodic with period N under periodic boundary conditions. Thereore the power spectrum with observation length T N has a undamental peak component at T/N and several harmonic peak components, while the power spectrum with T N is similar to white noise, which means the time series is random. But its randomness is mainly included in an initial coniguration rather than generated by a rule. Subcategory 2-A consists o rules 2, 3, 6, 7, 9, 0,, 4, 5, 24, 25, 26, 27, 34, 35, 38, 4, 42, 43, 46, 56, 57, 58, 74, 30, 34, 38, 42, 52, 54, 62, 70, and 84. Rule #6 S() Rule 90 S() 0.0 # Figure 2. Space-time patterns (let) and power spectra (right) o rule 6 (top) and rule 90 (bottom).
5 Power Spectral Analysis o Elementary Cellular Automata 403 The power spectra in subcategory 2-B are characterized by white noise. The bottom o Figure 2 shows the space-time pattern (let) and the power spectrum (right) o rule 90. The power spectrum in subcategory 2-B has an almost equal power density at all requencies except or rules 8 and 46, which have a peak with period 2. This result implies that evolution in subcategory 2-B is virtually orderless, although it is causally determined by a rule. Randomness in the evolution o subcategory 2-B owes much to the rule itsel rather than initial conigurations. Subcategory 2-B consists o rules 8, 22, 30, 45, 60, 90, 05, 06, 29, 46, 50, and 6 which coincide with class III [6]. By using an initial coniguration containing a single site with value, the dierence between subcategories 2-A and 2-B becomes clear. Subcategory 2-A resembles category in the shape o its power spectrum with the exceptions o rules 26 and 54, which exhibit broad-band noise. The power spectrum in subcategory 2-B, however, remains broad-band noise with the exception o rule 06, having an extremely low power density. Rules 26, 5, and 06 are considered as intermediates between these two subcategories. 2.3 Category 3: Power law spectrum The power spectrum in this category is characterized by power law at low requencies. Figure 3 shows the space-time patterns (let) and the power spectra (right) o rules 54 (top), 62 (middle), and 0 (bottom). The x and y axes in the power spectra are plotted on a logarithmic scale and the broken line represents the least square itting o the power spectrum rom to 0 by ln S( ) Α Βln,withΒ 0.5 or rule 54, Β.2 or rule 62, and Β.5 or rule 0. Rules 54 and 0 are considered to be in class IV, which is supposed to be capable o supporting universal computation [7] while there is a guess that the evolution o rule 62 is too simple to be universal [8]. 2.4 Exceptional rules There are two exceptional rules which do not belong to any o the three categories. Rule 73 is called locally chaotic because the array is divided into some independent domains by stable walls and chaotic patterns are generated in each domain [6]. Each domain has an individual period which contributes the peak in the power spectrum. Figure 4 shows the space-time pattern (let) and the power spectrum (right) o rule 73. The power spectrum o rule 204 has zero power density at all requencies except or 0 because the evolution o rule 204 retains the initial coniguration.
6 404 S. Ninagawa Rule 54 #54 ^(-0.5) 0.0 S() Rule 62 #62 ^(-.2) 0.0 S() Rule ^(-.5) S() # Figure 3. Space-time patterns (let) and the power spectra (right) o rules 54 (top), 62 (middle), and 0 (bottom). The x and y axes in the power spectra are plotted on a logarithmic scale. The broken line in the power spectra represents the least square itting o the power spectrum rom to 0 by ln S( ) Α Βln, Β 0.5 or rule 54, Β.2 or rule 62, and Β.5 or rule 0.
7 Power Spectral Analysis o Elementary Cellular Automata 405 Figure 4. Space-time pattern (let) and power spectrum (right) o rule / noise in elementary cellular automata A luctuation with a power spectrum S( ) that is inversely proportional to requency is called / noise [9]. / noise can be observed in a wide variety o phenomena such as the voltage o vacuum tubes, the rate o traic low, and the loudness o music. However, its origin is not well understood. The most controversial problem in / noise is whether / noise lasts orever or not. Generally speaking, i there is inite correlation time Τ in a luctuation, the power spectrum with an observation length o T > Τ has an almost equal power density at requencies smaller than /Τ. Likewise i there is inite correlation time Τ in the evolution o a CA, the power spectrum with an observation length o T > Τ has an almost equal power density at requencies smaller than T/Τ. So the power spectrum becomes close to a lat line at low requencies as the observation length T becomes longer than the correlation time Τ. Moreover we can guess that the correlation time Τ depends on array size N because the average transient time steps T ave o rule 0, which is supposed to be relevant to Τ, increases algebraically with array size N, T ave N Α,withΑ.08 [0]. Thereore we investigate the value o the exponent Β o the power spectra S( ) Β or various observation lengths and array sizes to compare the behaviors o rules 0, 54, and 62. Since one-dimensional CAs in general have larger variations in the value o Β with initial conigurations than two-dimensional CAs, we calculate the value o Β averaged over some dierent random initial conigurations. Figure 5 shows the average Β o the exponent estimated by the least square itting o the power spectrum rom to 0 by ln S( ) Α Βln( ) over 400 random initial conigurations as a unction o array size N and observation lengtht or rules 0 and 54 (top), and rules 0 and 62 (bottom). The solid line represents Β o rule 0 and the broken line rules 54 (top) and 62 (bottom). Β o rules 0, 54,
8 406 S. Ninagawa <β> #0 # N T <β> #0 # N T Figure 5. Average Β o the exponent estimated by the least square itting o the power spectrum rom to 0 by ln S( ) Α Βln( ) or 400 random initial conigurations as a unction o array size N and observation length T or rules 0 and 54 (top), and or rules 0 and 62 (bottom). The solid line represents Β o rule 0 and the broken line rules 54 (top) and 62 (bottom). Β o rules 0, 54, and 62 at N 200, T 8000 equal.2, 0.4, and 0.6 respectively. and62atn 200, T 6000 equal.2, 0.4, and 0.6 respectively. On the whole, Β o rules 0, 54, and 62 increases as the array size N decreases or the observation length T increases. The reason Β o rule 0 at N 300 is close to zero is that singular behavior tends to occur requently in one-dimensional CAs with small array sizes. When the array size is larger than 300, Β o rule 0 is smaller than those o rules 54 and 62. This result means that the evolution o rule 0 exhibits / noise during the longest observation length o the ECAs.
9 Power Spectral Analysis o Elementary Cellular Automata Origin o / noise in rule 0 In this section we compare two examples o the evolution o rule 0 rom distinct random initial conigurations to study the origin o / noise in the evolution o rule 0. Figure 6 shows a set o space-time patterns o the evolution o rule 0 rom a random initial coniguration o 200 cells or 000 time steps. The irst 200 time steps are shown on the let and the last 200 time steps are on the right. The evolution at irst exhibits transient behavior but becomes periodic at t 306 with period 750. The space-time pattern on the right shows that several gliders are moving monotonously within a periodic background. Figure 7 shows another set o space-time patterns rom dierent random initial conigurations. The evolution keeps transient behavior at t 999 although it eventually becomes periodic at t 246 with period 400. The space-time pattern on the right shows that several gliders interact complexly within a periodic background. Figure 8 shows the power spectra calculated rom the space-time patterns shown in Figures 6 and 7. The let one is calculated rom the evolution shown in Figure 6 and the right one rom Figure 7. The broken line represents the least square itting o the power spectrum rom to 0 by ln S( ) Α Βln, with Β 0.5 (let) and Β.5 (right). The let one is too lat to be considered as type /. Comparing these two examples o the evolution suggests that the dierence in the exponent o power spectra at low requencies depends on the duration o transient behavior. t 0 99 t Figure 6. Space-time patterns o rule 0 rom a random initial coniguration. Each picture is 200 cells wide. The irst 200 steps o evolution are shown on the let and the last 200 steps on the right. The evolution becomes periodic at t 306 with period 750. Several gliders are moving monotonously in the last 200 time step evolution.
10 408 S. Ninagawa t 0 99 t Figure 7. Space-time patterns o rule 0 rom a dierent random initial coniguration than the one used in Figure 6. Several gliders interact complexly in the last 200 steps o evolution. The evolution eventually becomes periodic at t 246 with period 400. Figure 8. Power spectra o rule 0 rom two random initial conigurations. The let one is calculated rom the evolution in Figure 6 and the right one rom Figure 7. The broken line represents the least square itting o the spectrum rom to 0 by ln S( ) Α Βln,withΒ 0.5 (let) and Β.5 (right). Intermittency in chaotic dynamical systems is one o the main mechanisms o / noise []. A system with a particular parameter value exhibits periodic behavior which is disrupted occasionally and irregularly by a burst. This burst persists or a inite duration, it stops and a new periodic behavior starts. Intermittent chaos occurs when the transition rom periodic to chaotic behavior takes place as the parameter is varied. During the transient behavior o rule 0 there are alternating periodic phases in which gliders are shiting in a periodic background without collisions and a burst that is caused by colliding gliders. It seems likely that the recurrence o the periodic phase and the burst in the transient behavior generates intermittency and causes / noise in rule 0.
11 Power Spectral Analysis o Elementary Cellular Automata Conclusion In this research we perormed spectral analysis on the evolution o 88 independent elementary cellular automata (ECAs) starting rom random initial conigurations and concluded that rule 0 exhibits / noise during the longest time steps. Rule 0 has proved to be capable o supporting universal computation [8, 2]. These results suggest that there is a relationship between computational universality and / noise in CAs. But there remains the problem o how to decide array size and observation length to detect / noise in CAs. The average Β o the exponent o the power spectrum o rules 0, 54, and 62 in Figure 5 at N 00, T 6000 equal 0.4, 0.5, and 0.04 respectively. These power spectra are not regarded as type /. These examples make it clear that / noise in ECAs is observed only in a speciic range o array size N and observation length T. Since the guess is that / noise in CAs is caused by transient behavior, as mentioned in section 4, we have to estimate experimentally a proper range o array size and observation length in which the evolution remains transient to observe a type / power spectrum. Another example which possesses both the exhibition o / noise and the capability o supporting universal computation is the Game o Lie (Lie) which is a two-dimensional and two-state, nine-neighbor outer totalistic CA [3]. It is supposed that a universal computer can be constructed on the array o Lie by considering a glider as a pulse in a digital circuit. In addition, the evolution rom random initial conigurations has a type / spectrum [4]. Another important property o Lie is sel-organized criticality where the distribution D(T) otime steps required or the lattice to return to stability ollowing a random single-site perturbation is characterized by D(T) T.6 [5]. Since / noise and sel-organized criticality seem to be interrelated phenomena, rule 0 might have the property o sel-organized criticality as well. It must be noted that the power law in the power spectrum o rule 0 holds or one decade in the range o requencies ( 0) as shown in Figure 8 (right), whereas it holds or three decades ( 000) in Lie [4]. The periodic background which is observed in rule 0, but not in Lie, accounts or this dierence. The three peaks o the power spectrum o rule 0 at 43, 286, and 429 in Figure 8 (right) are the undamental, the second, and the third harmonic components caused by the periodic background with period seven respectively. The interaction between the periodic background and other structures broadens the peaks and contributes to the power in high requencies. The hypothesis o the edge o chaos has evoked considerable controversy [6]. This hypothesis says the ability to perorm universal computation in a system arises near a transition rom regular to chaotic
12 40 S. Ninagawa behavior such as class IV. So ar various statistical quantities, such as entropy and dierence pattern spreading rate, have been proposed to detect class IV quantitatively [7]. The exponent o the power spectrum at low requencies might be another useul index. Acknowledgments I would like to thank the reeree or instructive comments on the original manuscript. Reerences [] J. Crutchield, D. Farmer, N. Packard, R. Shaw, G. Jones and R. J. Donnelly, Power Spectral Analysis o a Dynamical System, Physics Letters, 76A (980) 4. [2] S. Wolram, Statistical Mechanics o Cellular Automata, Reviews o Modern Physics, 55 (983) [3] S. Wolram (editor), Theory and Applications o Cellular Automata (World Scientiic, Singapore, 986). [4] W. Li, Power Spectra o Regular Languages and Cellular Automata, Complex Systems, (987) [5] C. G. Langton, Studying Artiicial Lie with Cellular Automata, Physica D, 22 (986) [6] W. Li and N. Packard, The Structure o the Elementary Cellular Automata Rule Space, Complex Systems, 4 (990) [7] S. Wolram, Universality and Complexity in Cellular Automata, Physica D, 0 (984) 35. [8] S. Wolram, A New Kind o Science (Wolram Media, Champaign, IL, 2002). [9] M. S. Keshner, / Noise, Proceedings o the IEEE, 70 (982) [0] W. Li and M. G. Nordahl, Transient Behavior o Cellular Automaton Rule 0, Physics Letters A, 66 (992) [] Y. Pomeau and P. Manneville, Intermittent Transition to Turbulence in Dissipative Dynamical Systems, Communications in Mathematical Physics, 74 (980) [2] M. Cook, Universality in Elementary Cellular Automata, Complex Systems, 5 (2004) 40. [3] E. R. Berlekamp, J. H. Conway, and R. K. Guy, Winning Ways or Your Mathematical Plays, Volume 4 (A. K. Peters, Wellesley, second edition, 2004).
13 Power Spectral Analysis o Elementary Cellular Automata 4 [4] S. Ninagawa, M. Yoneda, and S. Hirose, / Fluctuation in the Game o Lie, Physica D, 8 (998) [5] P. Bak, K. Chen, and M. Creutz, Sel-organized Criticality in the Game o Lie, Nature, 342 (989) [6] C. Langton, Computation at the Edge o Chaos: Phase Transitions and Emergent Computation, Physica D, 42 (990) [7] W. Li, N. H. Packard, and C. G. Langton, Transition Phenomena in Cellular Automata Rule Space, Physica D, 45 (990)
Extension of cellular automata by introducing an algorithm of recursive estimation of neighbors
Extension of cellular automata by introducing an algorithm of recursive estimation of neighbors Yoshihiko Kayama BAIKA Women s University, Japan (Tel: 81-72-643-6221, Fax: 81-72-643-8473) kayama@baika.ac.jp
More informationENSC327 Communications Systems 2: Fourier Representations. School of Engineering Science Simon Fraser University
ENSC37 Communications Systems : Fourier Representations School o Engineering Science Simon Fraser University Outline Chap..5: Signal Classiications Fourier Transorm Dirac Delta Function Unit Impulse Fourier
More informationAPPENDIX 1 ERROR ESTIMATION
1 APPENDIX 1 ERROR ESTIMATION Measurements are always subject to some uncertainties no matter how modern and expensive equipment is used or how careully the measurements are perormed These uncertainties
More informationLeast-Squares Spectral Analysis Theory Summary
Least-Squares Spectral Analysis Theory Summary Reerence: Mtamakaya, J. D. (2012). Assessment o Atmospheric Pressure Loading on the International GNSS REPRO1 Solutions Periodic Signatures. Ph.D. dissertation,
More informationAdditional exercises in Stationary Stochastic Processes
Mathematical Statistics, Centre or Mathematical Sciences Lund University Additional exercises 8 * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
More informationThe achievable limits of operational modal analysis. * Siu-Kui Au 1)
The achievable limits o operational modal analysis * Siu-Kui Au 1) 1) Center or Engineering Dynamics and Institute or Risk and Uncertainty, University o Liverpool, Liverpool L69 3GH, United Kingdom 1)
More informationEstimation and detection of a periodic signal
Estimation and detection o a periodic signal Daniel Aronsson, Erik Björnemo, Mathias Johansson Signals and Systems Group, Uppsala University, Sweden, e-mail: Daniel.Aronsson,Erik.Bjornemo,Mathias.Johansson}@Angstrom.uu.se
More information2. ETA EVALUATIONS USING WEBER FUNCTIONS. Introduction
. ETA EVALUATIONS USING WEBER FUNCTIONS Introduction So ar we have seen some o the methods or providing eta evaluations that appear in the literature and we have seen some o the interesting properties
More informationCategories and Natural Transformations
Categories and Natural Transormations Ethan Jerzak 17 August 2007 1 Introduction The motivation or studying Category Theory is to ormalise the underlying similarities between a broad range o mathematical
More information8. INTRODUCTION TO STATISTICAL THERMODYNAMICS
n * D n d Fluid z z z FIGURE 8-1. A SYSTEM IS IN EQUILIBRIUM EVEN IF THERE ARE VARIATIONS IN THE NUMBER OF MOLECULES IN A SMALL VOLUME, SO LONG AS THE PROPERTIES ARE UNIFORM ON A MACROSCOPIC SCALE 8. INTRODUCTION
More informationFluctuationlessness Theorem and its Application to Boundary Value Problems of ODEs
Fluctuationlessness Theorem and its Application to Boundary Value Problems o ODEs NEJLA ALTAY İstanbul Technical University Inormatics Institute Maslak, 34469, İstanbul TÜRKİYE TURKEY) nejla@be.itu.edu.tr
More informationTLT-5200/5206 COMMUNICATION THEORY, Exercise 3, Fall TLT-5200/5206 COMMUNICATION THEORY, Exercise 3, Fall Problem 1.
TLT-5/56 COMMUNICATION THEORY, Exercise 3, Fall Problem. The "random walk" was modelled as a random sequence [ n] where W[i] are binary i.i.d. random variables with P[W[i] = s] = p (orward step with probability
More informationAn Improved Expression for a Classical Type of Explicit Approximation of the Colebrook White Equation with Only One Internal Iteration
International Journal o Hydraulic Engineering 06, 5(): 9-3 DOI: 0.593/j.ijhe.06050.03 An Improved Expression or a Classical Type o Explicit Approximation o the Colebrook White Equation with Only One Internal
More informationCellular automata are idealized models of complex systems Large network of simple components Limited communication among components No central
Cellular automata are idealized models of complex systems Large network of simple components Limited communication among components No central control Complex dynamics from simple rules Capability of information
More informationIn many diverse fields physical data is collected or analysed as Fourier components.
1. Fourier Methods In many diverse ields physical data is collected or analysed as Fourier components. In this section we briely discuss the mathematics o Fourier series and Fourier transorms. 1. Fourier
More information1/f noise from the nonlinear transformations of the variables
Modern Physics Letters B Vol. 9, No. 34 (015) 15503 (6 pages) c World Scientiic Publishing Company DOI: 10.114/S0179849155031 1/ noise rom the nonlinear transormations o the variables Bronislovas Kaulakys,
More informationThe concept of limit
Roberto s Notes on Dierential Calculus Chapter 1: Limits and continuity Section 1 The concept o limit What you need to know already: All basic concepts about unctions. What you can learn here: What limits
More informationCurve Sketching. The process of curve sketching can be performed in the following steps:
Curve Sketching So ar you have learned how to ind st and nd derivatives o unctions and use these derivatives to determine where a unction is:. Increasing/decreasing. Relative extrema 3. Concavity 4. Points
More informationRoberto s Notes on Differential Calculus Chapter 8: Graphical analysis Section 1. Extreme points
Roberto s Notes on Dierential Calculus Chapter 8: Graphical analysis Section 1 Extreme points What you need to know already: How to solve basic algebraic and trigonometric equations. All basic techniques
More informationControlling chaos in random Boolean networks
EUROPHYSICS LETTERS 20 March 1997 Europhys. Lett., 37 (9), pp. 597-602 (1997) Controlling chaos in random Boolean networks B. Luque and R. V. Solé Complex Systems Research Group, Departament de Fisica
More informationNumerical Solution of Ordinary Differential Equations in Fluctuationlessness Theorem Perspective
Numerical Solution o Ordinary Dierential Equations in Fluctuationlessness Theorem Perspective NEJLA ALTAY Bahçeşehir University Faculty o Arts and Sciences Beşiktaş, İstanbul TÜRKİYE TURKEY METİN DEMİRALP
More information2 Frequency-Domain Analysis
2 requency-domain Analysis Electrical engineers live in the two worlds, so to speak, o time and requency. requency-domain analysis is an extremely valuable tool to the communications engineer, more so
More informationPhysics 5153 Classical Mechanics. Solution by Quadrature-1
October 14, 003 11:47:49 1 Introduction Physics 5153 Classical Mechanics Solution by Quadrature In the previous lectures, we have reduced the number o eective degrees o reedom that are needed to solve
More informationSpatial Vector Algebra
A Short Course on The Easy Way to do Rigid Body Dynamics Roy Featherstone Dept. Inormation Engineering, RSISE The Australian National University Spatial vector algebra is a concise vector notation or describing
More informationSpecial types of Riemann sums
Roberto s Notes on Subject Chapter 4: Deinite integrals and the FTC Section 3 Special types o Riemann sums What you need to know already: What a Riemann sum is. What you can learn here: The key types o
More informationPart I: Thin Converging Lens
Laboratory 1 PHY431 Fall 011 Part I: Thin Converging Lens This eperiment is a classic eercise in geometric optics. The goal is to measure the radius o curvature and ocal length o a single converging lens
More informationFigure 3.1 Effect on frequency spectrum of increasing period T 0. Consider the amplitude spectrum of a periodic waveform as shown in Figure 3.2.
3. Fourier ransorm From Fourier Series to Fourier ransorm [, 2] In communication systems, we oten deal with non-periodic signals. An extension o the time-requency relationship to a non-periodic signal
More informationA Fourier Transform Model in Excel #1
A Fourier Transorm Model in Ecel # -This is a tutorial about the implementation o a Fourier transorm in Ecel. This irst part goes over adjustments in the general Fourier transorm ormula to be applicable
More informationCellular Automata. Jason Frank Mathematical Institute
Cellular Automata Jason Frank Mathematical Institute WISM484 Introduction to Complex Systems, Utrecht University, 2015 Cellular Automata Game of Life: Simulator: http://www.bitstorm.org/gameoflife/ Hawking:
More informationApplication of Wavelet Transform Modulus Maxima in Raman Distributed Temperature Sensors
PHOTONIC SENSORS / Vol. 4, No. 2, 2014: 142 146 Application o Wavelet Transorm Modulus Maxima in Raman Distributed Temperature Sensors Zongliang WANG, Jun CHANG *, Sasa ZHANG, Sha LUO, Cuanwu JIA, Boning
More informationSIO 211B, Rudnick. We start with a definition of the Fourier transform! ĝ f of a time series! ( )
SIO B, Rudnick! XVIII.Wavelets The goal o a wavelet transorm is a description o a time series that is both requency and time selective. The wavelet transorm can be contrasted with the well-known and very
More informationOn Elementary and Algebraic Cellular Automata
Chapter On Elementary and Algebraic Cellular Automata Yuriy Gulak Center for Structures in Extreme Environments, Mechanical and Aerospace Engineering, Rutgers University, New Jersey ygulak@jove.rutgers.edu
More informationCHAPTER 1: INTRODUCTION. 1.1 Inverse Theory: What It Is and What It Does
Geosciences 567: CHAPTER (RR/GZ) CHAPTER : INTRODUCTION Inverse Theory: What It Is and What It Does Inverse theory, at least as I choose to deine it, is the ine art o estimating model parameters rom data
More informationSignals & Linear Systems Analysis Chapter 2&3, Part II
Signals & Linear Systems Analysis Chapter &3, Part II Dr. Yun Q. Shi Dept o Electrical & Computer Engr. New Jersey Institute o echnology Email: shi@njit.edu et used or the course:
More informationII. Cellular Automata 8/27/03 1
II. Cellular Automata 8/27/03 1 Cellular Automata (CAs) Invented by von Neumann in 1940s to study reproduction He succeeded in constructing a self-reproducing CA Have been used as: massively parallel computer
More informationSemideterministic Finite Automata in Operational Research
Applied Mathematical Sciences, Vol. 0, 206, no. 6, 747-759 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/0.2988/ams.206.62 Semideterministic Finite Automata in Operational Research V. N. Dumachev and
More informationPhiladelphia University Faculty of Engineering Communication and Electronics Engineering
Module: Electronics II Module Number: 6503 Philadelphia University Faculty o Engineering Communication and Electronics Engineering Ampliier Circuits-II BJT and FET Frequency Response Characteristics: -
More informationScattered Data Approximation of Noisy Data via Iterated Moving Least Squares
Scattered Data Approximation o Noisy Data via Iterated Moving Least Squares Gregory E. Fasshauer and Jack G. Zhang Abstract. In this paper we ocus on two methods or multivariate approximation problems
More informationRATIONAL FUNCTIONS. Finding Asymptotes..347 The Domain Finding Intercepts Graphing Rational Functions
RATIONAL FUNCTIONS Finding Asymptotes..347 The Domain....350 Finding Intercepts.....35 Graphing Rational Functions... 35 345 Objectives The ollowing is a list o objectives or this section o the workbook.
More informationIntrinsic Small-Signal Equivalent Circuit of GaAs MESFET s
Intrinsic Small-Signal Equivalent Circuit o GaAs MESFET s M KAMECHE *, M FEHAM M MELIANI, N BENAHMED, S DALI * National Centre o Space Techniques, Algeria Telecom Laboratory, University o Tlemcen, Algeria
More informationModelling with cellular automata
Modelling with cellular automata Shan He School for Computational Science University of Birmingham Module 06-23836: Computational Modelling with MATLAB Outline Outline of Topics Concepts about cellular
More informationTwo-step self-tuning phase-shifting interferometry
Two-step sel-tuning phase-shiting intererometry J. Vargas, 1,* J. Antonio Quiroga, T. Belenguer, 1 M. Servín, 3 J. C. Estrada 3 1 Laboratorio de Instrumentación Espacial, Instituto Nacional de Técnica
More information(Anti-)Stable Points and the Dynamics of Extended Systems
(Anti-)Stable Points and the Dynamics of Extended Systems P.-M. Binder SFI WORKING PAPER: 1994-02-009 SFI Working Papers contain accounts of scientific work of the author(s) and do not necessarily represent
More informationCHAPTER 8 ANALYSIS OF AVERAGE SQUARED DIFFERENCE SURFACES
CAPTER 8 ANALYSS O AVERAGE SQUARED DERENCE SURACES n Chapters 5, 6, and 7, the Spectral it algorithm was used to estimate both scatterer size and total attenuation rom the backscattered waveorms by minimizing
More informationA cellular-automata model of flow in ant trails: non-monotonic variation of speed with density
INSTITUTE OF PHYSICSPUBLISHING JOURNAL OFPHYSICSA: MATHEMATICAL AND GENERAL J. Phys. A: Math. Gen. 35 (22) L573 L577 PII: S35-447(2)52965-6 LETTER TO THE EDITOR A cellular-automata model o low in ant trails:
More informationImprovement of Sparse Computation Application in Power System Short Circuit Study
Volume 44, Number 1, 2003 3 Improvement o Sparse Computation Application in Power System Short Circuit Study A. MEGA *, M. BELKACEMI * and J.M. KAUFFMANN ** * Research Laboratory LEB, L2ES Department o
More informationAnalysis Scheme in the Ensemble Kalman Filter
JUNE 1998 BURGERS ET AL. 1719 Analysis Scheme in the Ensemble Kalman Filter GERRIT BURGERS Royal Netherlands Meteorological Institute, De Bilt, the Netherlands PETER JAN VAN LEEUWEN Institute or Marine
More informationDecibels, Filters, and Bode Plots
23 db Decibels, Filters, and Bode Plots 23. LOGAITHMS The use o logarithms in industry is so extensive that a clear understanding o their purpose and use is an absolute necessity. At irst exposure, logarithms
More informationNETWORK REPRESENTATION OF THE GAME OF LIFE
JAISCR, 2011, Vol.1, No.3, pp. 233 240 NETWORK REPRESENTATION OF THE GAME OF LIFE Yoshihiko Kayama and Yasumasa Imamura Department of Media and Information, BAIKA Women s University, 2-19-5, Shukuno-sho,
More informationSimpler Functions for Decompositions
Simpler Functions or Decompositions Bernd Steinbach Freiberg University o Mining and Technology, Institute o Computer Science, D-09596 Freiberg, Germany Abstract. This paper deals with the synthesis o
More informationAH 2700A. Attenuator Pair Ratio for C vs Frequency. Option-E 50 Hz-20 khz Ultra-precision Capacitance/Loss Bridge
0 E ttenuator Pair Ratio or vs requency NEEN-ERLN 700 Option-E 0-0 k Ultra-precision apacitance/loss ridge ttenuator Ratio Pair Uncertainty o in ppm or ll Usable Pairs o Taps 0 0 0. 0. 0. 07/08/0 E E E
More informationPhysics 2B Chapter 17 Notes - First Law of Thermo Spring 2018
Internal Energy o a Gas Work Done by a Gas Special Processes The First Law o Thermodynamics p Diagrams The First Law o Thermodynamics is all about the energy o a gas: how much energy does the gas possess,
More informationDynamic Modeling and Control of Shunt Active Power Filter
Dynamic Modeling and Control o Shunt Active Power Filter Utkal Ranjan Muduli Electrical Engineering Indian Institute o Technology Gandhinagar Gujarat, India-82424 Email: utkalranjan muduli@iitgn.ac.in
More informationObjectives. By the time the student is finished with this section of the workbook, he/she should be able
FUNCTIONS Quadratic Functions......8 Absolute Value Functions.....48 Translations o Functions..57 Radical Functions...61 Eponential Functions...7 Logarithmic Functions......8 Cubic Functions......91 Piece-Wise
More informationMEASUREMENT UNCERTAINTIES
MEASUREMENT UNCERTAINTIES What distinguishes science rom philosophy is that it is grounded in experimental observations. These observations are most striking when they take the orm o a quantitative measurement.
More informationIntroduction to Analog And Digital Communications
Introduction to Analog And Digital Communications Second Edition Simon Haykin, Michael Moher Chapter Fourier Representation o Signals and Systems.1 The Fourier Transorm. Properties o the Fourier Transorm.3
More informationDefinition: Let f(x) be a function of one variable with continuous derivatives of all orders at a the point x 0, then the series.
2.4 Local properties o unctions o several variables In this section we will learn how to address three kinds o problems which are o great importance in the ield o applied mathematics: how to obtain the
More informationAn Alternative Poincaré Section for Steady-State Responses and Bifurcations of a Duffing-Van der Pol Oscillator
An Alternative Poincaré Section or Steady-State Responses and Biurcations o a Duing-Van der Pol Oscillator Jang-Der Jeng, Yuan Kang *, Yeon-Pun Chang Department o Mechanical Engineering, National United
More informationLecture 25: Heat and The 1st Law of Thermodynamics Prof. WAN, Xin
General Physics I Lecture 5: Heat and he 1st Law o hermodynamics Pro. WAN, Xin xinwan@zju.edu.cn http://zimp.zju.edu.cn/~xinwan/ Latent Heat in Phase Changes Latent Heat he latent heat o vaporization or
More informationA Systematic Approach to Frequency Compensation of the Voltage Loop in Boost PFC Pre- regulators.
A Systematic Approach to Frequency Compensation o the Voltage Loop in oost PFC Pre- regulators. Claudio Adragna, STMicroelectronics, Italy Abstract Venable s -actor method is a systematic procedure that
More informationGaussian Process Regression Models for Predicting Stock Trends
Gaussian Process Regression Models or Predicting Stock Trends M. Todd Farrell Andrew Correa December 5, 7 Introduction Historical stock price data is a massive amount o time-series data with little-to-no
More informationReview D: Potential Energy and the Conservation of Mechanical Energy
MSSCHUSETTS INSTITUTE OF TECHNOLOGY Department o Physics 8. Spring 4 Review D: Potential Energy and the Conservation o Mechanical Energy D.1 Conservative and Non-conservative Force... D.1.1 Introduction...
More informationFORMAL LANGUAGES, AUTOMATA AND COMPUTABILITY. FLAC (15-453) Spring l. Blum TIME COMPLEXITY AND POLYNOMIAL TIME;
15-453 TIME COMPLEXITY AND POLYNOMIAL TIME; FORMAL LANGUAGES, AUTOMATA AND COMPUTABILITY NON DETERMINISTIC TURING MACHINES AND NP THURSDAY Mar 20 COMPLEXITY THEORY Studies what can and can t be computed
More informationarxiv: v1 [nlin.cg] 23 Sep 2010
Complex networks derived from cellular automata Yoshihiko Kayama Department of Media and Information, BAIKA Women s University, 2-9-5, Shukuno-sho, Ibaraki-city, Osaka-pref., Japan arxiv:009.4509v [nlin.cg]
More informationChaos-low periodic orbits transition in a synchronous switched circuit
Chaos-low periodic orbits transition in a synchronous switched circuit N. JABLI 1, H. KHAMMARI 2 and M.F. MIMOUNI 1 Electrical Department (1) National Engineering School o Monastir Rue Ibn aljazzar, 5019
More informationDecibels, Filters, and Bode Plots
Decibels, Filters, and Bode Plots 2 Objectives Develop conidence in the use o logarithms and decibels in the description o power and voltage levels. Become amiliar with the requency response o high- and
More informationComplexity Classes in the Two-dimensional Life Cellular Automata Subspace
Complexity Classes in the Two-dimensional Life Cellular Automata Subspace Michael Magnier Claude Lattaud Laboratoire d Intelligence Artificielle de Paris V, Université René Descartes, 45 rue des Saints
More informationOn High-Rate Cryptographic Compression Functions
On High-Rate Cryptographic Compression Functions Richard Ostertág and Martin Stanek Department o Computer Science Faculty o Mathematics, Physics and Inormatics Comenius University Mlynská dolina, 842 48
More informationStochastic Model for Adaptation Using Basin Hopping Dynamics
Stochastic Model for Adaptation Using Basin Hopping Dynamics Peter Davis NTT Communication Science Laboratories 2-4 Hikaridai, Keihanna Science City, Kyoto, Japan 619-0237 davis@cslab.kecl.ntt.co.jp Abstract
More informationAnderson impurity in a semiconductor
PHYSICAL REVIEW B VOLUME 54, NUMBER 12 Anderson impurity in a semiconductor 15 SEPTEMBER 1996-II Clare C. Yu and M. Guerrero * Department o Physics and Astronomy, University o Caliornia, Irvine, Caliornia
More informationarxiv: v1 [cond-mat.stat-mech] 6 Mar 2008
CD2dBS-v2 Convergence dynamics of 2-dimensional isotropic and anisotropic Bak-Sneppen models Burhan Bakar and Ugur Tirnakli Department of Physics, Faculty of Science, Ege University, 35100 Izmir, Turkey
More informationNONLINEAR CONTROL OF POWER NETWORK MODELS USING FEEDBACK LINEARIZATION
NONLINEAR CONTROL OF POWER NETWORK MODELS USING FEEDBACK LINEARIZATION Steven Ball Science Applications International Corporation Columbia, MD email: sball@nmtedu Steve Schaer Department o Mathematics
More informationAnalysis of self-induced vibrations in a pushing V-belt CVT
4CVT-32 Analysis o sel-induced vibrations in a pushing V-belt CVT Copyright 24 SAE International Wolram Lebrecht Institute o Applied Mechanics, Technical University o Munich Friedrich Peier, Heinz Ulbrich
More informationDerivation, f-derivation and generalized derivation of KUS-algebras
PURE MATHEMATICS RESEARCH ARTICLE Derivation, -derivation and generalized derivation o KUS-algebras Chiranjibe Jana 1 *, Tapan Senapati 2 and Madhumangal Pal 1 Received: 08 February 2015 Accepted: 10 June
More informationCOMPARISON OF THERMAL CHARACTERISTICS BETWEEN THE PLATE-FIN AND PIN-FIN HEAT SINKS IN NATURAL CONVECTION
HEFAT014 10 th International Conerence on Heat Transer, Fluid Mechanics and Thermodynamics 14 6 July 014 Orlando, Florida COMPARISON OF THERMA CHARACTERISTICS BETWEEN THE PATE-FIN AND PIN-FIN HEAT SINKS
More informationA Rough Sets Based Classifier for Induction Motors Fault Diagnosis
A ough Sets Based Classiier or Induction Motors Fault Diagnosis E. L. BONALDI, L. E. BOGES DA SILVA, G. LAMBET-TOES AND L. E. L. OLIVEIA Electronic Engineering Department Universidade Federal de Itajubá
More informationAnalog Computing Technique
Analog Computing Technique by obert Paz Chapter Programming Principles and Techniques. Analog Computers and Simulation An analog computer can be used to solve various types o problems. It solves them in
More informationINPUT GROUND MOTION SELECTION FOR XIAOWAN HIGH ARCH DAM
3 th World Conerence on Earthquake Engineering Vancouver, B.C., Canada August -6, 24 Paper No. 2633 INPUT GROUND MOTION LECTION FOR XIAOWAN HIGH ARCH DAM CHEN HOUQUN, LI MIN 2, ZHANG BAIYAN 3 SUMMARY In
More informationUNIVERSITY OF CALIFORNIA - SANTA CRUZ DEPARTMENT OF PHYSICS PHYS 112. Homework #4. Benjamin Stahl. February 2, 2015
UIERSIY OF CALIFORIA - SAA CRUZ DEPARME OF PHYSICS PHYS Homework #4 Benjamin Stahl February, 05 PROBLEM It is given that the heat absorbed by a mole o ideal gas in a uasi-static process in which both its
More information( x) f = where P and Q are polynomials.
9.8 Graphing Rational Functions Lets begin with a deinition. Deinition: Rational Function A rational unction is a unction o the orm ( ) ( ) ( ) P where P and Q are polynomials. Q An eample o a simple rational
More informationThe Deutsch-Jozsa Problem: De-quantization and entanglement
The Deutsch-Jozsa Problem: De-quantization and entanglement Alastair A. Abbott Department o Computer Science University o Auckland, New Zealand May 31, 009 Abstract The Deustch-Jozsa problem is one o the
More informationMixed Signal IC Design Notes set 6: Mathematics of Electrical Noise
ECE45C /8C notes, M. odwell, copyrighted 007 Mied Signal IC Design Notes set 6: Mathematics o Electrical Noise Mark odwell University o Caliornia, Santa Barbara rodwell@ece.ucsb.edu 805-893-344, 805-893-36
More information2.6 Two-dimensional continuous interpolation 3: Kriging - introduction to geostatistics. References - geostatistics. References geostatistics (cntd.
.6 Two-dimensional continuous interpolation 3: Kriging - introduction to geostatistics Spline interpolation was originally developed or image processing. In GIS, it is mainly used in visualization o spatial
More informationAnalysis of Friction-Induced Vibration Leading to Eek Noise in a Dry Friction Clutch. Abstract
The 22 International Congress and Exposition on Noise Control Engineering Dearborn, MI, USA. August 19-21, 22 Analysis o Friction-Induced Vibration Leading to Eek Noise in a Dry Friction Clutch P. Wickramarachi
More informationReview of Prerequisite Skills for Unit # 2 (Derivatives) U2L2: Sec.2.1 The Derivative Function
UL1: Review o Prerequisite Skills or Unit # (Derivatives) Working with the properties o exponents Simpliying radical expressions Finding the slopes o parallel and perpendicular lines Simpliying rational
More informationClicker questions. Clicker question 2. Clicker Question 1. Clicker question 2. Clicker question 1. the answers are in the lower right corner
licker questions the answers are in the lower right corner question wave on a string goes rom a thin string to a thick string. What picture best represents the wave some time ater hitting the boundary?
More informationYURI LEVIN AND ADI BEN-ISRAEL
Pp. 1447-1457 in Progress in Analysis, Vol. Heinrich G W Begehr. Robert P Gilbert and Man Wah Wong, Editors, World Scientiic, Singapore, 003, ISBN 981-38-967-9 AN INVERSE-FREE DIRECTIONAL NEWTON METHOD
More informationProducts and Convolutions of Gaussian Probability Density Functions
Tina Memo No. 003-003 Internal Report Products and Convolutions o Gaussian Probability Density Functions P.A. Bromiley Last updated / 9 / 03 Imaging Science and Biomedical Engineering Division, Medical
More informationVirtually everything you and I know about the Cosmos has been discovered via electromagnetic
Gravitational Waves Gravitational wave astronomy Virtually everything you and I know about the Cosmos has been discovered via electromagnetic observations. Some inormation has recently been gleaned rom
More informationRELIABILITY OF BURIED PIPELINES WITH CORROSION DEFECTS UNDER VARYING BOUNDARY CONDITIONS
REIABIITY OF BURIE PIPEIES WITH CORROSIO EFECTS UER VARYIG BOUARY COITIOS Ouk-Sub ee 1 and ong-hyeok Kim 1. School o Mechanical Engineering, InHa University #53, Yonghyun-ong, am-ku, Incheon, 40-751, Korea
More informationBINARY MORPHOLOGY AND CELLULAR AUTOMATA
BINARY MORPHOLOGY AND CELLULAR AUTOMATA I can't leave this subject without mentioning cellular automata (CAs). Conway's "Game of Life" is an example of a cellular automaton (CA). In each generation (or
More informationGrover Algorithm Applied to Four Qubits System
www.ccsenet.org/cis Computer and Inormation Science Vol., No. ; May Grover Algorithm Applied to Four Qubits System Z. Sakhi (Corresponding author) Laboratory o Inormation Technology and Modelisation, and
More informationMODULE 6 LECTURE NOTES 1 REVIEW OF PROBABILITY THEORY. Most water resources decision problems face the risk of uncertainty mainly because of the
MODULE 6 LECTURE NOTES REVIEW OF PROBABILITY THEORY INTRODUCTION Most water resources decision problems ace the risk o uncertainty mainly because o the randomness o the variables that inluence the perormance
More informationEx x xf xdx. Ex+ a = x+ a f x dx= xf x dx+ a f xdx= xˆ. E H x H x H x f x dx ˆ ( ) ( ) ( ) μ is actually the first moment of the random ( )
Fall 03 Analysis o Eperimental Measurements B Eisenstein/rev S Errede The Epectation Value o a Random Variable: The epectation value E[ ] o a random variable is the mean value o, ie ˆ (aa μ ) For discrete
More informationCellular Automata CS 591 Complex Adaptive Systems Spring Professor: Melanie Moses 2/02/09
Cellular Automata CS 591 Complex Adaptive Systems Spring 2009 Professor: Melanie Moses 2/02/09 Introduction to Cellular Automata (CA) Invented by John von Neumann (circa~1950). A cellular automata consists
More informationbiologically-inspired computing lecture 12 Informatics luis rocha 2015 INDIANA UNIVERSITY biologically Inspired computing
lecture 12 -inspired Sections I485/H400 course outlook Assignments: 35% Students will complete 4/5 assignments based on algorithms presented in class Lab meets in I1 (West) 109 on Lab Wednesdays Lab 0
More informationAggregate Growth: R =αn 1/ d f
Aggregate Growth: Mass-ractal aggregates are partly described by the mass-ractal dimension, d, that deines the relationship between size and mass, R =αn 1/ d where α is the lacunarity constant, R is the
More informationNumerical Methods - Lecture 2. Numerical Methods. Lecture 2. Analysis of errors in numerical methods
Numerical Methods - Lecture 1 Numerical Methods Lecture. Analysis o errors in numerical methods Numerical Methods - Lecture Why represent numbers in loating point ormat? Eample 1. How a number 56.78 can
More informationGENERALIZED ABSTRACT NONSENSE: CATEGORY THEORY AND ADJUNCTIONS
GENERALIZED ABSTRACT NONSENSE: CATEGORY THEORY AND ADJUNCTIONS CHRIS HENDERSON Abstract. This paper will move through the basics o category theory, eventually deining natural transormations and adjunctions
More information