DOCSLIB.ORG

Chapter 3, Dynamic Behavior

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
Chapter 3, Dynamic Behavior Chapter 3 Dynamic Behavior In this chapter we give a broad discussion of the behavior of dynamical sys- tems. We focus on systems modeled by di®erential equations, but consider the general nonlinear case. This allows us to discuss equilibrium points, stability, limit cycles and other key concepts of dynamical systems. We also introduce some measures of performance for input/output systems that provide additional tools for characterizing dynamic behavior. 3.1 Solving Di®erential Equations In the last chapter, we saw that one of the methods of modeling dynamical systems is through the use of ordinary di®erential equations (ODEs). A state space, input/output system has the form x_ = f(x; u) (3.1) y = g(x); where x Rn is the state, u Rp is the input, and y Rq is the output. The 2 2 2 smooth maps f : Rn Rp Rn and g : Rn Rq represent the dynamics £ ! ! and measurements for the system. We will focus in this text on single input, single output (SISO) systems, for which p = q = 1. We begin by investigating systems in which the input has been set to a function of the state, u = ®(x). This is one of the simplest types of feedback, in which the system regulates its own behavior. The di®erential equations in this case become x_ = f(x; ®(x)) = F (x): (3.2) As before, we write x = (x ; : : : ; x ) Rn for the state vector. 1 n 2 61 62 CHAPTER 3. DYNAMIC BEHAVIOR In order to understand the dynamic behavior of this system, we need to analyze the features of the solutions of this equation. While in some simple situations we can write down the solutions in analytical form, more often we must rely on computational approaches. We begin by describing the class of solutions for this problem. Initial Value Problems We say that x(t) is a solution of the di®erential equation (3.2) on the time interval t R to t R if 0 2 f 2 x_(t) = F (x(t)) for all t t < t . 0 · f A given di®erential equation may have many solutions. We will most often be interested in the initial value problem, where x(t) is prescribed at a given time t R and we wish to ¯nd a solution valid for all future time, t > t . 0 2 0 We say that x(t) is a solution of the di®erential equation (3.2) with initial value x Rn at t R if 0 2 0 2 x(t ) = x and x_(t) = F (x(t)) for all t t < t . 0 0 0 · f For most di®erential equations we will encounter, there is a unique solution that is de¯ned for t t t . The solution may de¯ned for all time t > t , 0 · · f 0 in which case we take t = . Because we will primarily be interested in f 1 solutions of the initial value problem for ODEs, we will often refer to this simply as the solution of an ODE. We will usually assume that t0 is equal to 0. In the case when F is indepen- Ä dent of time (as in equation (3.2)), we can do so without loss of generality by choosing a new independent (time) variable, ¿ = t t . ¡ 0 Example 3.1 (Damped oscillator). Consider a damped, linear oscillator, in- troduced in Example 2.3. The equations of motion for the system are mqÄ + bq_ + kq = 0; where q is the displacement of the oscillator from its rest position. We assume that b2 < 4km, corresponding to a lightly damped system (the reason for this particular choice will become clear later). We can rewrite this in state space form by setting x1 = q and x2 = q_, giving x_ 1 = x2 k b x_ = x x : 2 ¡m 1 ¡ m 2 3.1. SOLVING DIFFERENTIAL EQUATIONS 63 1 0.5 2 , x 0 1 x −0.5 −1 0 5 10 15 20 time (sec) Figure 3.1: Response of the damped oscillator to the initial condition x0 = (1; 0). In vector form, the right hand side can be written as x F (x) = 2 : k x b x ·¡ m 1 ¡ m 2¸ The solution to the initial value problem can be written in a number of di®erent ways and will be explored in more detail in Chapter 4. Here we simply assert that the solution can be written as bt x20 x (t) = e¡ 2m x cos !t + (x + ) sin !t 1 10 10 ! bt ³ b ´ b x20 x (t) = e¡ 2m x + x (! ) cos !t + x ! + (x + ) sin !t 2 20 10 ¡ 2 10 2 10 ! µ ¶ ¡ ¢ ¡ ¢ where x = (x ; x ) is the initial condition and ! = p4km b2=2m. This 0 10 20 ¡ solution can be veri¯ed by substituting it into the di®erential equation. We see that the solution is explicitly dependent on the initial condition and it can be shown that this solution is unique. A plot of the initial condition response is shown in Figure 3.1. We note that this form of the solution only holds for b2 4km < 0, corresponding to an \underdamped" oscillator. ¡ 64 CHAPTER 3. DYNAMIC BEHAVIOR Numerical Solutions One of the bene¯ts of the computer revolution that we can bene¯t from is that it is very easy to obtain a numerical solution of the di®erential equation when the initial condition is given. A nice consequence of this is as soon as we have a model in the form of (3.2), it is straightforward to generate the behavior of x for di®erent initial conditions, as we saw briey in the previous chapter. Modern computing environments allow simulation of di®erential equa- tions as a basic operation. In particular, MATLAB provides several tools for representing, simulating, and analyzing ordinary di®erential equations of the form in equation (3.2). To de¯ne an ODE in MATLAB, we de¯ne a function representing the right hand side of equation (3.2): function dxdt = sysname(t, x) dxdt = [ F1(x); F2(x); ... Fn(x); ]; Each function Fi(x) takes a (column) vector x and returns the ith element of the di®erential equation. The ¯rst argument to the function sysname, t, represents the current time and allows for the possibility of time-varying dif- ferential equations, in which the right hand side of the ODE in equation (3.2) depends explicitly on time. ODEs de¯ne in this fashion can be simulated by using the MATLAB ode45 command: ode45('file', [0,T], [x10, x20, ..., xn0]) The ¯rst argument is the name of the ¯le containing the ODE declaration, the second argument gives the time interval over which the simulation should be performed and the ¯nal argument gives the vector of initial conditions. The default action of the ode45 command is to plot the time response of each of the states of the system. Example 3.2 (Balance system). Consider the balance system given in Exam- ple 2.1 and reproduced Figure 3.2a. Suppose that a coworker has designed 3.1. SOLVING DIFFERENTIAL EQUATIONS 65 m 1.2 p theta 1 0.8 0.6 p 0.4 0.2 position (m), angle (deg) 0 F M −0.2 −0.4 0 2 4 6 8 10 12 14 16 18 20 time (sec) (a) (b) Figure 3.2: Balance system: (a) simpli¯ed diagram and (b) initial condition response. a control law that will hold the position of the system steady in the upright position at p = 0. The form of the control law is F = Kx; where x = (p; ; p_; _) R4 is the state of the system, F is the input, and 2 K = (k1; k2; k3; k4) is the vector of \gains" for the control law. The equations of motion for the system, in state space form, are p p_ d _ 2 3 = 2 1 3 dt p_ M + m ml cos ¡ bx_ + ml sin _2 + Kx 6_7 6 2 ¡ 7 6 7 6 J + ml ml cos mgl sin 7 4 5 4· ¸ · ¡ ¸5 p y = : · ¸ We use the following parameters for the system (corresponding roughly to a human being balanced on a stabilizing cart): M = 10 kg m = 80 kg b = 0.1 Ns/m J = 100 kg m2/s2 l = 1 m g = 9.8 m/s2 K = 1 120 4 20 ¡ ¡ This system can now be sim£ulated using MATLAB¤ or a similar numerical tool. The results are shown in Figure 3.2b, with initial condition x0 = (1; 0; 0; 0). We see from the plot that after an initial transient, the angle and position of the system return to zero (and remain there). 66 CHAPTER 3. DYNAMIC BEHAVIOR 100 100 80 80 60 60 x x 40 40 20 20 PSfrag replacements 0 PSfrag replacements 0 0 0.2 0.4 0.6 0.8 1 0 2 4 6 8 10 t t (a) (b) Figure 3.3: Solutions to the di®erential equation (3.3) (a) and (3.4) (b). Ä Existence and Uniqueness Without imposing some conditions on the function F the di®erential equa- tion (3.2) may not have a solution for all t, and there is no guarantee that the solution is unique. We illustrate this with two examples. Example 3.3 (Finite escape time). Let x R and consider the di®erential 2 equation dx = x2 (3.3) dt with initial condition x(0) = 1. By di®erentiation we can verify that the function 1 x(t) = (3.4) 1 t ¡ satis¯es the di®erential equation and it also satis¯es the initial condition.
Load more

Recommended publications
  • Chapter 8 Nonlinear Systems
    Chapter 8 Nonlinear systems 8.1 Linearization, critical points, and equilibria Note: 1 lecture, §6.1–§6.2 in [EP], §9.2–§9.3 in [BD] Except for a few brief detours in chapter 1, we considered mostly linear equations. Linear equations suffice in many applications, but in reality most phenomena require nonlinear equations. Nonlinear equations, however, are notoriously more difficult to understand than linear ones, and many strange new phenomena appear when we allow our equations to be nonlinear. Not to worry, we did not waste all this time studying linear equations. Nonlinear equations can often be approximated by linear ones if we only need a solution “locally,” for example, only for a short period of time, or only for certain parameters. Understanding linear equations can also give us qualitative understanding about a more general nonlinear problem. The idea is similar to what you did in calculus in trying to approximate a function by a line with the right slope. In § 2.4 we looked at the pendulum of length L. The goal was to solve for the angle θ(t) as a function of the time t. The equation for the setup is the nonlinear equation L g θ�� + sinθ=0. θ L Instead of solving this equation, we solved the rather easier linear equation g θ�� + θ=0. L While the solution to the linear equation is not exactly what we were looking for, it is rather close to the original, as long as the angleθ is small and the time period involved is short. You might ask: Why don’t we just solve the nonlinear problem? Well, it might be very difficult, impractical, or impossible to solve analytically, depending on the equation in question.
    [Show full text]
  • Introduction to Mathematical Modeling Difference Equations, Differential Equations, & Linear Algebra (The First Course of a Two-Semester Sequence)
    Introduction to Mathematical Modeling Difference Equations, Differential Equations, & Linear Algebra (The First Course of a Two-Semester Sequence) Dr. Eric R. Sullivan [email protected] Department of Mathematics Carroll College, Helena, MT Content Last Updated: January 8, 2018 1 ©Eric Sullivan. Some Rights Reserved. This work is licensed under a Creative Commons Attribution-NonCommercial-ShareAlike 4.0 International License. You may copy, distribute, display, remix, rework, and per- form this copyrighted work, but only if you give credit to Eric Sullivan, and all deriva- tive works based upon it must be published under the Creative Commons Attribution- NonCommercial-Share Alike 4.0 United States License. Please attribute this work to Eric Sullivan, Mathematics Faculty at Carroll College, [email protected]. To view a copy of this license, visit https://creativecommons.org/licenses/by-nc-sa/4.0/ or send a letter to Creative Commons, 171 Second Street, Suite 300, San Francisco, Cali- fornia, 94105, USA. Contents 0 To the Student and the Instructor5 0.1 An Inquiry Based Approach...........................5 0.2 Online Texts and Other Resources........................7 0.3 To the Instructor..................................7 1 Fundamental Notions from Calculus9 1.1 Sections from Active Calculus..........................9 1.2 Modeling Explorations with Calculus...................... 10 2 An Introduction to Linear Algebra 16 2.1 Why Linear Algebra?............................... 16 2.2 Matrix Operations and Gaussian Elimination................. 20 2.3 Gaussian Elimination: A First Look At Solving Systems........... 23 2.4 Systems of Linear Equations........................... 28 2.5 Linear Combinations............................... 35 2.6 Inverses and Determinants............................ 38 2.6.1 Inverses................................... 40 2.6.2 Determinants..............................
    [Show full text]
  • Mathematical Modeling of Control Systems
    OGATA-CH02-013-062hr 7/14/09 1:51 PM Page 13 2 Mathematical Modeling of Control Systems 2–1 INTRODUCTION In studying control systems the reader must be able to model dynamic systems in math- ematical terms and analyze their dynamic characteristics.A mathematical model of a dy- namic system is defined as a set of equations that represents the dynamics of the system accurately, or at least fairly well. Note that a mathematical model is not unique to a given system.A system may be represented in many different ways and, therefore, may have many mathematical models, depending on one’s perspective. The dynamics of many systems, whether they are mechanical, electrical, thermal, economic, biological, and so on, may be described in terms of differential equations. Such differential equations may be obtained by using physical laws governing a partic- ular system—for example, Newton’s laws for mechanical systems and Kirchhoff’s laws for electrical systems. We must always keep in mind that deriving reasonable mathe- matical models is the most important part of the entire analysis of control systems. Throughout this book we assume that the principle of causality applies to the systems considered.This means that the current output of the system (the output at time t=0) depends on the past input (the input for t<0) but does not depend on the future input (the input for t>0). Mathematical Models. Mathematical models may assume many different forms. Depending on the particular system and the particular circumstances, one mathemati- cal model may be better suited than other models.
    [Show full text]
  • Mathematical Modeling and Ordinary Differential Equations
    MATHEMATICAL MODELING AND ORDINARY DIFFERENTIAL EQUATIONS I-Liang Chern Department of Mathematics National Taiwan University 2007, 2015 January 6, 2016 2 Contents 1 First-Order Single Differential Equations 1 1.1 What is mathematical modeling? . .1 1.2 Relaxation and Equilibria . .2 1.3 Modeling population dynamics of single species . .6 1.4 Techniques to solve single first-order equations . 17 1.4.1 Linear first-order equation . 17 1.4.2 Separation of variables . 20 1.4.3 Other special classes that are solvable . 21 1.5 Vector Fields and Family of Curves . 22 1.5.1 Vector Fields . 22 1.5.2 Family of curves and Orthogonal trajectories . 27 1.5.3 *Envelop. 29 1.5.4 *An example from thermodynamics – existence of entropy . 31 1.6 Existence and Uniqueness . 33 1.7 *Numerical Methods: First Order Difference Equations . 37 1.7.1 Euler method . 37 1.7.2 First-order difference equation . 38 1.8 Historical Note . 39 2 Second Order Linear Equations 41 2.1 Models for linear oscillators . 41 2.1.1 The spring-mass system . 41 2.1.2 Electric circuit system . 42 2.2 Methods to solve homogeneous equations . 43 2.2.1 Homogeneous equations (complex case) . 43 2.2.2 Homogeneous equation (real case) . 45 2.3 Methods to solve Inhomogeneous equations . 47 2.3.1 Method of underdetermined coefficients . 48 2.3.2 Method of Variation of Constants . 49 2.4 Linear oscillators . 52 2.4.1 Harmonic oscillators . 52 2.4.2 Damping . 53 3 4 CONTENTS 2.4.3 Forcing and Resonance .
    [Show full text]
  • In Order to Derive His Equations, Lorenz (1963)
    THE LORENZ EQUATIONS MARIANNA KATERINA PELEKANOS MSc MATHEMATICS SCHOOL OF MATHEMATICAL SCIENCES QUEEN MARY UNIVERSITY September 10, 2009 1 ABSTRACT The project aims to analyse the work of Edward Lorenz. The way dynamical systems with complex behaviour, such as the weather behave after a long period of time. The project will mention the Lorenz equations and what they represent. What happens to the system if we vary one parameter having the other parameters fixed will also be discussed here. In addition it will be discussed about what happens to the dynamical system if we slightly change the initial conditions. Finally it will talk about what is meant by sensitive dependence on initial conditions and the butterfly effect and how these to concepts are related to chaos theory. 2 TABLE OF CONTENTS Acknowledgements 4 1. Introduction 5 2. Derivation of the Lorenz Equations 7 2.1 The Bernard Rayleigh Experiment 7 2.2 Derivation of the Lorenz Equations 9 3. The Parameter r 18 3.1 Solutions of the Lorenz Equations 18 3.2 Different Values of r 19 4. Nontrivial Fixed Points 21 4.1 The Region r>1 21 4.2 Time Dependent non Stationary Solutions 23 4.3 Period Doubling Scenarios 24 5. The Butterfly Effect 27 6. Conclusion 29 Bibliography 30 3 ACKNOWLEDGEMENTS This dissertation was an interesting experience for me, but would not have been possible without the support of my supervisor Doctor Wolfram Just which I thank very much for all the help and support he have given me. My special thanks to my parents, who made this Masters possible and also to my friends for their support.
    [Show full text]
  • Variations of the Lorenz-96 Model
    Variations of the Lorenz-96 Model Master Thesis in Mathematics: Mathematics and Complex Dynamical Systems Name: Anouk F.G. Pelzer, s2666545 First supervisor: Dr. A.E. Sterk Second supervisor: Prof. Dr. H. Waalkens Date: 29 June 2020 Abstract. In this thesis we will investigate the differences between the monoscale Lorenz-96 model and its modifications. We obtain these modifications by changing the structure of the nonlinear terms in the original Lorenz-96 model. We will treat the gen- eral modified model and three specific systems and compare these with the Lorenz-96 model. To determine the dynamics of the models, we will study the eigenvalues of the Jacobian matrix at the trivial equilibrium, Lyapunov coefficients to distinguish between sub- and supercritical Hopf bifurcations, Lyapunov exponents to determine when there is chaos for example etc. Some of the results are that in the Lorenz-96 model there are no escaping orbits, while this is only true under some conditions for the modified systems. Furthermore there are differences in bifurcations. The external forcing F ap- pearing in the equations of all these models can be used as a bifurcation parameter. The trivial equilibrium is always stable when F is zero for every model. For F < 0 some modifications has a stable trivial equilibrium, while this is not the case for the original model. For this model there only occur first Hopf bifurcations when F is positive, but for its modifications there can be other (first) bifurcations too, like a Pitchfork bifurcation. There are even bifurcations occuring for these modified systems, while these don't ap- pear in the original Lorenz-96 model case, such as degenerate bifurcations meaning that suddenly there appear simultaneously more than one equilibrium and two eigenvalues of the Jacobian matrix are zero at the same bifurcation parameter F .
    [Show full text]
  • Analysis of Discovery of Chaos: Social and Cognitive Aspects
    DOCUMENT RESUME ED 421 340 SE 061 533 AUTHOR Kim, J. B. TITLE Analysis of Discovery of Chaos: Social and Cognitive Aspects. PUB DATE 1998-04-00 NOTE 25p.; Paper presented at the Annual Meeting of the American Educational Research Association (San Diego, CA, April 13-17, 1998). PUB TYPE Reports Evaluative (142)-- Speeches/Meeting Papers (150) EDRS PRICE MF01/PC01 Plus Postage. DESCRIPTORS Cognitive Processes; *Computer Simulation; Computers; *Discovery Processes; Higher Education; Knowledge Representation; Physics; *Science and Society; *Scientific Concepts; *Scientific Principles; Scientists; Technology IDENTIFIERS *Chaos; *Chaos Theory ABSTRACT The purpose of this study was to examine Edward Lorenz's psychological processes and other environmental aspects in the discovery of chaos at that time. The general concept of chaos is discussed based on relations with previous scientific theories such as Newtonian physics and quantum mechanics. The constraints of discovery in terms of available technology at the time of discovery are explained. The following general arguments are made in this paper:(1) the idea of chaos is discovered utilizing the computer and it is hard to investigate without computer technology;(2) Lorenz's accumulated beliefs (knowledge) create a perceptual problem space which is a path toward the discovery of chaos;(3) knowledge can be a tool to guide discovery;(4) the approval of the discovery is a generalized and abstract process; and (5) visual representations give a general picture of the effects of chaos simulation. (Contains 36 references.) (DDR) ******************************************************************************** * Reproductions supplied by EDRS are the best that can be made * * from the original document. * ******************************************************************************** Discovery of Chaos 1 Analysis of Discovery of Chaos: Social and CognitiveAspects Kim, J.-B.
    [Show full text]
  • Chapter 3 One-Dimensional Systems
    Chapter 3 One-Dimensional Systems In this chapter we describe geometrical methods of analysis of one-dimensional dynam- ical systems, i.e., systems having only one variable. An example of such a system is the space-clamped membrane having Ohmic leak current IL ˙ C V = ¡gL(V ¡ EL) : (3.1) Here the membrane voltage V is a time-dependent variable, and the capacitance C, leak conductance gL and leak reverse potential EL are constant parameters described in the previous chapter. We use this and other one-dimensional neural models to introduce and illustrate the most important concepts of dynamical system theory: equilibrium, stability, attractor, phase portrait, and bifurcation. 3.1 Electrophysiological Examples The Hodgkin-Huxley description of dynamics of membrane potential and voltage-gated conductances can be reduced to a one-dimensional system when all transmembrane conductances have fast kinetics. For the sake of illustration, let us consider a space- clamped membrane having leak current and a fast voltage-gated current Ifast having only one gating variable p, Leak I I z }| L { z }|fast { ˙ C V = ¡ gL(V ¡ EL) ¡ g p (V ¡ E) (3.2) p˙ = (p1(V ) ¡ p)=¿(V ) (3.3) with dimensionless parameters C = 1, gL = 1, and g = 1. Suppose that the gating kinetic (3.3) is much faster than the voltage kinetic (3.2), which means that the voltage- sensitive time constant ¿(V ) is very small, i.e. ¿(V ) ¿ 1, in the entire biophysical voltage range. Then, the gating process may be treated as being instantaneous, and the asymptotic value p = p1(V ) may be used
    [Show full text]
  • Investigation of Lorenz Equation System with Variable Step Size Strategy
    Proceedings of International Conference on Mathematics and Mathematics Education (ICMME 2019) Turk. J. Math. Comput. Sci. 11(Special Issue)(2019) 16–20 c MatDer https://dergipark.org.tr/tjmcs http://tjmcs.matder.org.tr MATDER Investigation of Lorenz Equation System with Variable Step Size Strategy Saniye I˙nce Polat1 ,Gulnur¨ C¸ elik Kızılkan2;∗ , 1S¸ehit Mustafa C¸uhadar Secondary School, 42250, Sel¸cuklu,Konya. 2Department of Mathematics-Computer, Faculty of Science, Necmettin Erbakan University, 42090, Meram, Konya. Received: 24-09-2019 • Accepted: 25-11-2019 Abstract. In this study, variable step size strategy has been considered to analyze the numerical solution of the Lorenz system with chaotic structure. Phase portraits have been obtained for this chaotic system. The effectiveness of the variable step size strategy for the solution of this chaotic system has been discussed. 2010 AMS Classification: 37N30, 65P20. Keywords: Lorenz systems, numerical solution, variable step size. 1. Introduction As a meteorologist, Edward Norton Lorenz studied to analyze the weather forecasts in 1963 and he developed a mathematical model consisting of three non-linear differential equations. In his paper “Deterministic Nonperiodic Flow”, Lorenz explained that these equations are sensitive to the initial conditions. Thus, the idea of “deterministic chaos” was discovered. In later years, chaos was thought to be very useful for many technological disciplines and was studied extensively in science and engineering ( [1,6, 13, 14]). Therefore, numerical and approximate solutions of the Lorenz system are analyzed in many studies ( [2,7,8, 11, 15]). The mathematical model of the Lorenz system has given as following: dx dy dz = σ(y − x); = rx − y − z; = xy − bz (1.1) dt dt dt where x, y and z are state variables; σ is the prandtl number, b and r are positive constant parameters.
    [Show full text]
  • Edward Norton Lorenz Discovererofchaos
    GENERAL ARTICLE Edward Norton Lorenz DiscovererofChaos V Krishnamurthy Edward Lorenz discovered nonperiodic behavior in deter- ministic nonlinear systems and laid the foundation of chaos theory. He showed that chaos exhibits sensitive depen- dence on initial conditions implying that long-range weather prediction is difficult because of errors in the observations used as initial conditions. Lorenz described the intricate structure of chaotic attractors and quantified predictability. His important contributions in meteorology include energy V Krishnamurthy is at cycle, slow manifold and general circulation. the Center for Ocean- Land-Atmosphere Discovery of Chaos Studies, George Mason University, USA, and Weather forecasts routinely issued nowadays by major prediction has worked at MIT, University of Maryland centers in the world are prepared by integrating global-scale and Abdus Salam ICTP. numerical models on supercomputers. Weather prediction models He was a doctoral are based on dynamical equations governing the atmosphere, student of Edward ocean, land and other components. The first dynamical weather Lorenz. His interests include chaos, monsoon forecast was reported by the project led by Jule Charney and John variability and climate von Neumann at the Institute for Advanced Study in Princeton in change . 1950 using the pioneering Electronic Numerical Integrator and Computer (ENIAC). However, during the 1950s, the weather forecasts were usually made by statistical models that were primarily linear methods relying on past observed data. In 1955, the Department of Meteorology at the Massachusetts Institute of Technology (MIT) appointed Edward Lorenz as a new faculty member to lead the on-going statistical forecasting project. Lorenz examined numerous statistical schemes and convinced himself that the statistical predictions were similar to subjective predictions and that even one-day forecasts were mediocre.
    [Show full text]
  • Intro to Chaos Theory and Its Application to Not-So-Sciency-Things
    Intro to Chaos Theory and its Application to Not-So-Sciency-Things Wes Brown In the 1960‟s, an assistant professor at MIT popularized one of the most catchy phrases in science: the Butterfly Effect. First mentioned in A Sound of Thunder, a short story by Ray Bradbury, the seemingly inconsequential death of a butterfly in the past makes a realistic impact in the future. As a literary device used to entertain, the concept suggests that our present conditions can be dramatically altered by the most insignificant change in the past. Unlike Bradbury though, Edward Norton Lorenz used the literary device in a paper he wrote in 1972: Predictability: Does the Flap of a Butterfly’s Wings in Brazil Set Off a Tornado in Texas? As a mathematician and meteorologist, Lorenz was skeptical about the linear systems being developed to predict the weather. As such, he began creating models of his own, including a 12-variable computer model used for weather prediction. While doing so, Lorenz noticed that small variations in his initial variables produced profound changes in his results. Further pursuing the idea, the meteorologist published a paper describing his findings in the Journal of the Atmospheric Sciences. His 1963 paper, Deterministic Nonperiodic Flow, is credited for laying the foundation for Chaos Theory. Often associated with physics and mathematics, Chaos Theory has now become a field of study that examines the properties of systems that are considered “chaotic.” In short, these systems are characterized as being anything that is highly sensitive to initial conditions. These initial conditions have rich interactions with each other, operate far from equilibrium, and self- organize into patterns that can develop into emergent behavior.
    [Show full text]
  • Nonlinear Ordinary Differential Equations
    Nonlinear Ordinary Differential Equations by Peter J. Olver University of Minnesota 1. Introduction. These notes are concerned with initial value problems for systems of ordinary dif- ferential equations. Here our emphasis will be on nonlinear phenomena and properties, particularly those with physical relevance. Finding a solution to a differential equation may not be so important if that solution never appears in the physical model represented by the system, or is only realized in exceptional circumstances. Thus, equilibrium solu- tions, which correspond to configurations in which the physical system does not move, only occur in everyday situations if they are stable. An unstable equilibrium will not ap- pear in practice, since slight perturbations in the system or its physical surroundings will immediately dislodge the system far away from equilibrium. Of course, very few nonlinear systems can be solved explicitly, and so one must typ- ically rely on a numerical scheme to accurately approximate the solution. Basic methods for initial value problems, beginning with the simple Euler scheme, and working up to the extremely popular Runge–Kutta fourth order method, will be the subject of the final section of the chapter. However, numerical schemes do not always give accurate results, and we briefly discuss the class of stiff differential equations, which present a more serious challenge to numerical analysts. Without some basic theoretical understanding of the nature of solutions, equilibrium points, and stability properties, one would not be able to understand when numerical so- lutions (even those provided by standard well-used packages) are to be trusted. Moreover, when testing a numerical scheme, it helps to have already assembled a repertoire of nonlin- ear problems in which one already knows one or more explicit analytic solutions.
    [Show full text]