Final State Sensitivity: an Obstruction to Predictability
Total Page:16
File Type:pdf, Size:1020Kb
Volume 99A, number 9 PHYSICS LETTERS 26 December 1983 FINAL STATE SENSITIVITY: AN OBSTRUCTION TO PREDICTABILITY Celso GREBOGI 1 Steven W. McDONALD 1 , Edward OTT 1'2 and James A. YORKE 3 University of Maryland, College Park, MD 20742, USA Received 21 October 1983 It is shown that nonlinear systems with multiple attractors commonly require very accurate initial conditions for the reliable prediction of final states. A scaling exponent for the final-state-uncertain phase space volume dependence on un- certainty in initial conditions is defined and related to the fractai dimension of basin boundaries. Typical nonlinear dynamical systems may have more is the "basin bounctary". It the initial conaitions are than one possible time-asymptotic final state. In such uncertain by an amount e, then (cf. fig. 1), for those cases the final state that is eventually reached depends initial conditions within e of the boundary, we cannot on the initial state of the system. In this letter we say a priori to which attractor the orbit eventually wish to consider the extent to which uncertainty in tends. For example, in fig. 1, 1 and 2 represent two initial conditions leads to uncertainty in the final state. initial conditions with an uncertainty e. The orbit To orient the discussion, consider the simple two- generated by initial condition 1 is attracted to attrac- dimensional phase space diagram schematically de- tor B. Initial condition 2, however, is uncertain in the picted in fig. 1. There are two possible final states sense that the orbit generated by 2 may be attracted ("attractors") denoted A and B. Initial conditions on either to A or to B. In particular, consider the fraction one side of the boundary, Z, eventually asymptote to of the uncertain phase space volume within the rec- B, while those on the other side of Z eventually go to tangle shown and denote this fraction f. For the case A. The region to the left (right) of Z is the "basin of shown in fig. 1, we clearly have f ~ e. It is one of the attraction" for attractor A (or B, respectively) and main points of this letter that, from the point of view of prediction, much worse scalings of f with e fre- quently occur in nonlinear dynamics. In particular, if 1 Laboratory for Plasma and Fusion Energy Studies, and De- partment of Physics and Astronomy. f ~e ~ , (1) 2 Also, Department of Electrical Engineering. 3 Institute for Physical Science and Technology, and Depart- with c~ < 1, we shall say that there is final state sensi- ment of Mathematics. tivity. In fact, a substantially less than one is, we be- lieve, fairly common. In such a case, a substantial im- provement in the initial condition uncertainty, e, yields only a relatively small decrease in the uncertainty of the final state as measured by f. While ¢x is equal to one for simple basin boundaries, such as that depicted in fig. 1, highly convoluted boundaries with noninteger (fractal) dimension also occur. We use here the capacity definition of dimen- sion [ 1 ], Fig. 1. A region of phase space divided by the basin boundary Y. into basins of attraction for the two attractors A and B. 1 In N(8) and 2 represent two initial conditions with uncertainty e. d = lim (2) 6--,0 ln(1/8) ' 0.031-9163/83/0000-0000/$ 03.00 © 1983 North-Holland 415 Volume 99A, number 9 PHYSICS LETTERS 26 December 1983 where, if the dimensionality of the relevant phase space is D, then N(6) is the minimum number of D ( dimensional cubes of side 6 needed to cover the basin boundary. In general, since the basin boundary di- c vides the phase space, its dimension d must satisfy d ~> D-1. It can be proven ,1 [2] that the following c relation between the index c~ and the basin boundary dimension holds -( c~ = O - d. (3) -( For a simple boundary, such as that depicted in fig. 1, we have d = D-l, and eq. (3) then gives c~ = 1, as ex- O rr/2 Jr pected. For a fractal basin boundary, d > D-l, and Fig. 2. Basins of attraction for eqs. (4). Arrows denote the eq. (3) gives a < 1 (i.e., there is final state sensitivity). two attractors. To see heuristically how (3) comes about from (2), set the cube edge 6 equal to the initial condition un- certainty e. The volume of the uncertain region of these parameters we find numerically that the only phase space will be of the order of the total volume attractors are the fixed points (0, -0.3) and (Tr, 0.3). of all the N(e) D-dimensional cubes of side e needed Fig. 2 shows a computer generated picture of the basins to cover the basin boundary. Since the volume of one of attraction for the two fixed point attractors. This of these D-dimensional cubes is e D, the uncertain vol- figure is constructed using a 256 X 256 grid of initial ume of phase space is of the order of eDN(e). Noting conditions. Each initial condition is iterated until it is that N(e) ~ e -d satisfies (2), we use this to estimate close to one of the two attractors (100 iterates of the the uncertain phase space volume as eDN(e) ~ e D-d, map is always sufficient to accomplish this). If an orbit which is the result predicted by eqs. (1) and (3). goes to the attractor at 0 = 0, a black dot is plotted at We now illustrate the above with a concrete exam- the corresponding initial condition. If the orbit goes ple. Subsequently, we argue that the demonstrated to the other attractor, no dot is plotted. Thus the phenomena occur commonly in nonlinear dynamics black and blank regions are essentially pictures of the (e.g., in the Lorenz system [3] and in experiments on basins of attractions for the two attractors to the ac- fluids such as those of Berg6 and Dubois [4]). curacy of the grid used and of the computer plotter. We consider the two-dimensional map (Due to the symmetry of the map we have only shown 0 ~< 0 ~ rr in fig. 2). Fine scale structure in the basins On+ 1 = On + a sin 20 n - b sin 40 n - x n sin On , (4a) of attraction is evident. This is a consequence of the Cantor set nature of the basin boundary * 2. In fact, Xn+l = -Jo cos On , (4b) magnifications of the basin boundary show that, as where 0 and 0 + 27r are identified as equivalent. This we examine it on a smaller and smaller scale, it con- map has two fixed points (0, x) = (0, -J0), and (Tr, J0), tinues to have structure. which are attracting for I1 + 2a - 4bl < 1. Numerical experiments on eqs. (4) with different sets of parame- ters have been performed. Here, as an example, we re- ~:2 The Cantor set structure of the basin boundary for eqs. (4) is due to the presence of a horseshoe in the dynamics port results forJ 0 = 0.3, a = 1.32, and b = 0.90. For [5]. In fact, eqs. (4) were chosen because they are a par- ticularly straightforward example of this. (This point will *1 For the purposes of the proof of eq. (3) and for those of be extensively discussed in ref. [9]). Cantor set basin bound- the present letter, eq. (1) should be regarded as shorthand aries due to horseshoes axe not curves (i.e., they cannot be for a = lime~olnf/ln e. In addition, we note that, if the represented as x = x(u), 0 = O(u) with x(u) and O(u) con- basins are unbounded regions, then we restrict attention tinuous functions of u), and have been known for a long to a bounded region of phase space (e.g., the rectangular time (see, for example ref. [6]). Fractal basin boundaries region of fig. 1) for the purposes of calculating f in (1) which are curves (unlike that for fig. 2) have been dis- and N(6) in (2). cussed by Grebogi et al. [7]. 416 Volume 99A, number 9 PHYSICS LETTERS 26 December 1983 We now wish to explore the consequences for pre- for e = 0.125, 59% of the initial conditions are uncer- diction of this infinitely fine scaled structure. To do tain; for e = 0.002, 26%; for e = 3 × 10 -5, 12%. Thus this, consider an initial condition (00, x0). We ask, even apparently small uncertainty in the initial con- what is the effect of a small change e in the x-coordi- ditions yieMs substantial fractions of the phase space nate? Thus we iterate the initial conditions (00, x0), which are uncertain as to which final state is eventually (00, x 0 + e), and (00, x 0 - e) until they approach one attained. The implication is that extraordinarily high of the attractors. If either or both of the perturbed accuracy of initial conditions may sometimes be nec- initial conditions yield orbits which do not approach essary for the reliable prediction of the eventual final the same attractor as the unperturbed initial condi- state. tion, we say that (00, x0) is uncertain. Now say that Furthermore, the example just discussed is by no we randomly choose a large number of initial condi- means extreme or unusual. As evidence for this we tions in the rectangle shown in fig. 2, and let fdenote cite two examples: the fraction of these which we find to be uncertain. (1) The Lorenz system.