ANALYSIS OF THE CHAOS DYNAMICS IN (X n,x n+1) PLANE

Similar documents
Numerical Method for Ordinary Differential Equation

Comparison between Fourier and Corrected Fourier Series Methods

Paper 3A3 The Equations of Fluid Flow and Their Numerical Solution Handout 1

Calculus BC 2015 Scoring Guidelines

2 f(x) dx = 1, 0. 2f(x 1) dx d) 1 4t t6 t. t 2 dt i)

Review Exercises for Chapter 9

BEST LINEAR FORECASTS VS. BEST POSSIBLE FORECASTS

EE757 Numerical Techniques in Electromagnetics Lecture 8

N! AND THE GAMMA FUNCTION

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

The analysis of the method on the one variable function s limit Ke Wu

Approximating Solutions for Ginzburg Landau Equation by HPM and ADM

Mean Square Convergent Finite Difference Scheme for Stochastic Parabolic PDEs

F D D D D F. smoothed value of the data including Y t the most recent data.

Some Properties of Semi-E-Convex Function and Semi-E-Convex Programming*

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

1. Solve by the method of undetermined coefficients and by the method of variation of parameters. (4)

SUMMATION OF INFINITE SERIES REVISITED

Linear System Theory

Time Dependent Queuing

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

If boundary values are necessary, they are called mixed initial-boundary value problems. Again, the simplest prototypes of these IV problems are:

Four equations describe the dynamic solution to RBC model. Consumption-leisure efficiency condition. Consumption-investment efficiency condition

Fresnel Dragging Explained

Energy Density / Energy Flux / Total Energy in 1D. Key Mathematics: density, flux, and the continuity equation.

C(p, ) 13 N. Nuclear reactions generate energy create new isotopes and elements. Notation for stellar rates: p 12

BE.430 Tutorial: Linear Operator Theory and Eigenfunction Expansion

NEWTON METHOD FOR DETERMINING THE OPTIMAL REPLENISHMENT POLICY FOR EPQ MODEL WITH PRESENT VALUE

ECE-314 Fall 2012 Review Questions

Electrical Engineering Department Network Lab.

Review Answers for E&CE 700T02

A Complex Neural Network Algorithm for Computing the Largest Real Part Eigenvalue and the corresponding Eigenvector of a Real Matrix

Big O Notation for Time Complexity of Algorithms

CS623: Introduction to Computing with Neural Nets (lecture-10) Pushpak Bhattacharyya Computer Science and Engineering Department IIT Bombay

S n. = n. Sum of first n terms of an A. P is

Prakash Chandra Rautaray 1, Ellipse 2

Available online at J. Math. Comput. Sci. 4 (2014), No. 4, ISSN:

Fourier transform. Continuous-time Fourier transform (CTFT) ω ω

INTEGER INTERVAL VALUE OF NEWTON DIVIDED DIFFERENCE AND FORWARD AND BACKWARD INTERPOLATION FORMULA

Exploring and Simulating Chaotic Advection: A Difference Equations Approach

INVESTMENT PROJECT EFFICIENCY EVALUATION

FIXED FUZZY POINT THEOREMS IN FUZZY METRIC SPACE

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

Notes 03 largely plagiarized by %khc

BAYESIAN ESTIMATION METHOD FOR PARAMETER OF EPIDEMIC SIR REED-FROST MODEL. Puji Kurniawan M

Pure Math 30: Explained!

Homotopy Analysis Method for Solving Fractional Sturm-Liouville Problems

Extremal graph theory II: K t and K t,t

REDUCED DIFFERENTIAL TRANSFORM METHOD FOR GENERALIZED KDV EQUATIONS. Yıldıray Keskin and Galip Oturanç

Clock Skew and Signal Representation

Problems and Solutions for Section 3.2 (3.15 through 3.25)

Inference of the Second Order Autoregressive. Model with Unit Roots

The Solution of the One Species Lotka-Volterra Equation Using Variational Iteration Method ABSTRACT INTRODUCTION

STK4080/9080 Survival and event history analysis

Fermat Numbers in Multinomial Coefficients

Class 36. Thin-film interference. Thin Film Interference. Thin Film Interference. Thin-film interference

Lateral torsional buckling of rectangular beams using variational iteration method

FRACTIONAL SYNCHRONIZATION OF CHAOTIC SYSTEMS WITH DIFFERENT ORDERS

Electromagnetic Waves: Outline. Electromagnetic wave propagation in Particle-In-Cell codes. 1D discrete propagation equation in vacuum

The Uniqueness Theorem for Inverse Nodal Problems with a Chemical Potential

Inverse Heat Conduction Problem in a Semi-Infinite Circular Plate and its Thermal Deflection by Quasi-Static Approach

A Note on Random k-sat for Moderately Growing k

OLS bias for econometric models with errors-in-variables. The Lucas-critique Supplementary note to Lecture 17

An interesting result about subset sums. Nitu Kitchloo. Lior Pachter. November 27, Abstract

A Novel Approach for Solving Burger s Equation

Comparisons Between RV, ARV and WRV

Optimization of Rotating Machines Vibrations Limits by the Spring - Mass System Analysis

MCR3U FINAL EXAM REVIEW (JANUARY 2015)

Academic Forum Cauchy Confers with Weierstrass. Lloyd Edgar S. Moyo, Ph.D. Associate Professor of Mathematics

Supplement for SADAGRAD: Strongly Adaptive Stochastic Gradient Methods"

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

EXISTENCE THEORY OF RANDOM DIFFERENTIAL EQUATIONS D. S. Palimkar

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

Effects of Forces Applied in the Middle Plane on Bending of Medium-Thickness Band

International journal of Engineering Research-Online A Peer Reviewed International Journal Articles available online

A TAUBERIAN THEOREM FOR THE WEIGHTED MEAN METHOD OF SUMMABILITY

CLOSED FORM EVALUATION OF RESTRICTED SUMS CONTAINING SQUARES OF FIBONOMIAL COEFFICIENTS

MODIFIED ADOMIAN DECOMPOSITION METHOD FOR SOLVING RICCATI DIFFERENTIAL EQUATIONS

Samuel Sindayigaya 1, Nyongesa L. Kennedy 2, Adu A.M. Wasike 3

A Study On (H, 1)(E, q) Product Summability Of Fourier Series And Its Conjugate Series

Lecture 15: Three-tank Mixing and Lead Poisoning

Conditional Probability and Conditional Expectation

CSE 241 Algorithms and Data Structures 10/14/2015. Skip Lists

Four equations describe the dynamic solution to RBC model. Consumption-leisure efficiency condition. Consumption-investment efficiency condition

The Connection between the Basel Problem and a Special Integral

14.02 Principles of Macroeconomics Fall 2005

Numerical Solution of Parabolic Volterra Integro-Differential Equations via Backward-Euler Scheme

Using Linnik's Identity to Approximate the Prime Counting Function with the Logarithmic Integral

K3 p K2 p Kp 0 p 2 p 3 p

LINEAR APPROXIMATION OF THE BASELINE RBC MODEL SEPTEMBER 17, 2013

The Moment Approximation of the First Passage Time for the Birth Death Diffusion Process with Immigraton to a Moving Linear Barrier

Application of the Adomian Decomposition Method (ADM) and the SOME BLAISE ABBO (SBA) method to solving the diffusion-reaction equations

6/10/2014. Definition. Time series Data. Time series Graph. Components of time series. Time series Seasonal. Time series Trend

An Efficient Method to Reduce the Numerical Dispersion in the HIE-FDTD Scheme

Variational Iteration and Homotopy Perturbation Method for Solving a Three-Species Food Chain Model

th m m m m central moment : E[( X X) ] ( X X) ( x X) f ( x)

ECE 340 Lecture 19 : Steady State Carrier Injection Class Outline:

A Bayesian Approach for Detecting Outliers in ARMA Time Series

Moment Generating Function

Section 8 Convolution and Deconvolution

Transcription:

ANALYSIS OF THE CHAOS DYNAMICS IN (X,X ) PLANE Soegiao Soelisioo, The Houw Liog Badug Isiue of Techolog (ITB) Idoesia soegiao@sude.fi.ib.ac.id Absrac I he las decade, sudies of chaoic ssem are more ofe used for classical chaoic ssem ha for quaum chaoic ssem. There are ma was of observig he chaoic ssem such as aalzig he frequec wih Fourier rasform or aalze iiial codiio disace wih Liapuov epoe. This paper eplais damic chaoic process b observig rajecor of damic ssem i (, ) pla. Lab.Fisika Theor da Kompuasi, Fisika UNAIR, Phsics Deparme, ITB, Idoesia Lecurer ad Professor of Phsics, Phsics Deparme, ITB, Idoesia

INTRODUCTION This paper ries o eplai some resuls of researches i aswerig he quesio of wha is chaos?. Chaos as characerisic of a ssem is defied as he ssem ha is suscepible o iiial codiio so ha he resul is oo much upredicable. This research ries o defie he cause of chaoic ssem. To aswer he quesio, he ivesigaors cosider (X, X ) plae. (, ) plae is o help udersadig he behaviors of ssem. (See picure ). X (X, X ) PLANE IS TOOL TO INVESTIGATE THE ITERATION OF SYSTEM DYNAMICS For eplaaio coveiece, he damic graph from ha ieraig process is used. Cosider he ieraio process for he logisic mappig: X = AX ( X + () for X =, ad A = 3 Afer imes ieraio he values are obaied: Value,,48 3,7488 4,5649568 5,73759897 6,58649 7,7349947 8,59677 9,7536384,59763444,74548,694736 3,78838 4,675539 5,755539 6,6649 7,7395999 8,635445 9,734349,6689 Wih ha above umber, i is difficul o udersad behavior ha ssem. Eplaaio wih ) Picure. visualizaio of damic ieraio i (X,X ) plae To clarif he coecig lie is draw bewee do ad add lie X =X ad lie of he fucio F() = 3 ( - ) Wih ha picure he damical of he ieraio is kow. Afer his he curve F() is called ieraed fucio curve ad lie X =X c is called reversal lie. GUIDANCE DYNAMICAL SYSTEM IN (X,X ) PLANE The research ries o udersad he behavior of chaoic ssem i (, ) plae ad discussio i his paper is o heoreical basis o eperimeal for ODE ad PDE ha is wellkow havig chaoic ssem The resul will esimae ha: The ssem damics will be guided o move oward he poi or he regio. Afer he ssem damics is moved io his poi or regio he he ssem damics cao move o aoher place. X

The poi or he regios ha arac he damic ssem are called aracor. The Ssem will move damicall ad coiuousl sice here s o poi or regio o arac. I his case, he ssem damics is called he chaoic ssem. The ssem damics will be guided o repel ad escape o ifii. This codiio is o good because o iformaio ca be obaied. The ssem damics will be guided o move i circular moio like plaear moio, graig elecro, ec. I reali, deails of ssem damics guidace are: curve cross he reversal lie, ad saisf > > 4. The damic ssem will be guided o repel he crossig poi if he ieraig fucio curve is crossig he reversal lie, ad saisf <. The ssem damics will be guided o move dow if he ieraig fucio curve is locaed uder reversal lie Tha rule of guidace ssem damics will be used o aalsis he logisic equaio (equaio ).. The ssem damics will be guided o move up if ieraig fucio curve is locaed upper reversal lie. < A < Case 3. The damic ssem will be guide o oward crossig poi if ieraed fucio Because ieraig fucio curve is uder he reversal lie, so he ssem damics will be guided o move dow oward he crossig poi (aracor) ha is (,).

< A < 3 Case: Because ieraig fucio curve is locaed above he reversal lie so ha he ssem damics will be guided o move up oward he crossig poi ha is called aracor. curve ha is locaed above reversal lie so ha ssem damics will be guided o move up ad secod ieraig fucio curve ha is crossig reversal lie ha ad < Tha ssem damics will be guided o repel crossig poi. 3.5699... < A < 4 Case: 3< A < 3.449... Case: I he of case he curve likes above, here are wo rule of guidace: firs ieraig fucio curve ha is locaed above reversal lie. The ssem damics will be guided o move up ad secod ieraig fucio curve ha is crossig reversal lie ad saisf > > Tha ssem damics will be guided o move forward o crossig poi. Tha appeared soluio is - /A 3.449... < A < 3.5699... Case: I he case of above curve, here s o aracor o arac he ssem damics so ha ssem damics will keep o movig Usig guidace rule, ssem damics ca defie aracor ad chaos, ha is: Aracor is he poi or regio ha arac he ssem damics, if he ssem damics is rapped io his poi or regio, he ssem damics will ever move awhere else. Chaos is a codiio ha here s o aracor o arac he ssem damics, ssem will keep movig. Wih he above defiiio, equaio of chaoic ssem ca be cosruced ha is equaio wih o aracor. Eample : I he case of above curve, here are also wo rule of guidace: firs ieraig fucio f ( ) = ( )... for... <.4 (,5 )... if.. else

Solvig ODE usig Numerical Aalsis i Chaoic Ssem Time series he damical above is I aure, phsical pheomea should be formulaed usig differeial equaios. I is ofe difficul o solve differeial equaio usig aalical mehod. Aoher mehods are usig umerical compuaios ha use ieraios. As i is udersood ha ieraio ca ield chaoic soluio. Oe eample below ca be solved b usig aalical mehod ha geerae eac soluio. Usig some umerical mehod for he same problem produce chaos. Eample : f ( ) = (.+. )... for... >.3 (,5 )... if.. else Eample : Cosider he equaio d( ) = a. b. d () wih a ad b as cosa he eac soluio ca be obaied a. e a + b. e a. = a. (3) The above ODE ca be re-wrie b usig Euler umerical soluio as : h = a. b. (4) Time series he damical above is Equaio (4) also ca be re-wrie as = α. + β. + (5) wih : From he above eplaaio, i ca be show ha ieraio ca produce chaos. Therefore, oe should be careful whe usig ieraio process. α = a. h + β = b. h (6) if α =β = A he equaio (5) ca be wrie :

= A. ( ) + (7) Time series equaio (7) is he same as equaio (), ad i produces chaoic damics. alhough i is kow o have eac soluio (equaio 3). This umerical soluio is o iflueced b h value, because h is ol used for α ad β. Suppose i is chose a = 3 ad b = + a, wih h =., Usig he umerical Euler mehod Profile ((),()).75 YEPC.5 Y_.43..4.6.8 solusi umerik solusi aaliik Usig Ruge - Kua mehod Profile ((),()) fmrk( ) fe( ) YE k.5 YRK k.5.5,, XE k.5,, XRKk Time series.998 Time series VROXVLSHUVDPDDQ.8.75 YE Y_.5 YRK Y_.6 6.48 3..4.6.8 solusi umerik solusi aaliik.4.4..4.6.8 solusi umerik solusi aaliik Usig predicor correcor mehod Profile ((),()) fepc( ) YEPC k.5.5,, XEPC k Wih ha codiio, ew problem arises, i.e. if oe differeial equaio ca geerae eac soluio ad some umerical mehod for he same problems produce chaos, he i is assumed ha he eac soluio is correc. Wha abou differeial equaios ha have o eac soluio ad he umerical mehod o solve hem produce chaos? Ca he observed ssem be called chaoic ssem? Oe of he eample is Lorez equaio.

Coclusio From above sud, i ca be eplaied ha ieraio mehod coribue o chaoic soluio Some umerical mehod coai problems: ieraio compuaio, error samplig frequec, roudig he resul. To overcome hose ieraio, error samplig ad roudig he resul problems, oe mus provide o ieraig mehod. I is kow ha eural ework, fuzz logic, cellular auomaa, ec are o-ieraig mehod such as Euler, Predicor-Correcor, ec. Neural ework, fuzz logic ad cellular auomaa are kow as compee uiversal approimaor as well as adapive mehod. Equaios from uiversal approimaio ca be used o solve umerical problem of ODE ad PDE. B miimizig he hree umerical problems oe ca be sure ha compuig umerical process do o cause chaoic ssem. Refereces [] DeVries, Paul L.,993, A Firs Course i Compuaioal Phsics. Joh Wile & Sos, Ic [] Fu, Li Mi, 994, Neural Nework i Compuer Ielligece. McGraw-Hill [3]. H Nagashima & Y Baba,999, Iroducio o Chaos, IOP Publishig Ld. [4] Lagaris, e al, 997, Arificial Neural Nework for solvig ordiar ad parial differeial equaio. ArXiv:phsics / 9753