The Nonlinear Pendulum

Size: px
Start display at page:

Download "The Nonlinear Pendulum"

Transcription

1 The Nonlinear Pendulum Evan Sheridan Feburary 18th 013 Abstract Both the non-linear linear pendulum are investigated compared using the pendulum.c program that utilizes the trapezoid method for solving differential equations. The damped pendulum was then analysed using the improved Runge-Kutta method after time it was found the non-linear pendulum approximated the Linear pendulum. Finally, the damped driven oscillator was analysed signs of chaotic behaviour were exhibited for A > 1.50m. 1

2 1 Aims To use the pendulum.c program to determine the motion of the simple pendulum To alter the pendulum.c program so that it uses the Runge-Kutta Method analyse the motion of the damped driven pendulum Use Gnuplot as an aid to interpret data. Backround Theory In general for a simple pendulum we have the equation of motion: d θ dt = g sinθ) 1.1) l If we make a small angle approximation, i.e sinθ) θ then we get : d θ dt = g l θ 1.) We call 1.) the equation of motion for the linear pendulum 1.1) the equation for the nonlinear pendulum. Now we wish to solve for the solution of these second order differential equations. In general we use a computational method to do this it can be done by 3 of the the following methods: The Simple Euler Method, The Trapezoid Method The Runge-Kutta Method. The Simple Euler Method In this method we take a second order differential equation reduce the problem to solving to first order differential equations. We do this by introducing the a new variable. In the specific case of the non-linear pendulum we have the following: ωt) = dθ dt fθ, ω, t) = dω dt Now if we wish to solve these equations numerically we need to begin with our initial values θ 0, ω 0 use a Taylor Expansion in order to derive θ 1,...,n, ω 1,...,n. The Taylor Expansion is given by: θt + t) = θt) + dθt) dt t + d θt) t d + O t 3 ) t! ωt + t) = ωt) + dω dt t + d ω d t which when writing in iterative notation becomes: t! + O t 3 ) θ n+1 = θ n + ω n t + O t ) 1

3 ω n+1 = ω n θ n t + O t ) Thus, now using these iterative functions we having the initial values we can solve for the motion trajectory. Obviously, in order to get any idea of what the motion looks like we would have to do quite a few iterations this is why computational methods are used for it would be far too tedious to this by h. The Simple Euler Method lays the foundation the next two methods are improvements on it. The Trapezoid Method This method makes a better approximation for the iterative functions. Instead we find that the area under a trapezoid for a θt) curve is gives the following result: θ n+1 θ n + t dθt) dθt + t + dt dt By explicitly writing out the Taylor Expansions for both θ n, ω n we get the following: k1a + k b θ n+1 = θ n + k1b+ k b ω n+1 = ω n + with k 1a = ω n t; k a = ω n + fθ n, ω n, t) t) t; k b = tfθ n+1, ω n+1 + k 1b, t n+1 ) All in all, the trapezoid method is a better approximation for solving the second order differential equation than the Simple Euler Method. The Runge-Kutta Method The equation of the Non-Linear Damped Driven Pendulum is given by: fθ, ω, t) = β sinθ) kω + A cosωt) In order to solve for the motion trajectory in this instance we use the 4th Order Runge-Kutta Method. Using the Runge-Kutta method we get the follwing: k1a + k a + k 3a + k 4a θt + h) = θt) + 6 k1b + k b + k 3b + k 4b ωt + h) = ωt) + 6 with k a = h ω + k ) 1b k 3a = h ω + k ) b k 1a = hω k 1b = hfθ, ω, t) k b = hf k 3b = hf θ + k 1a, ω + k 1b, t + h ) θ + k a, ω + k b, t + h ) k 4a = hω + k 3b ) k 4b = hfθ + k 3a, ω + k 3b, t + h)

4 The Experiment Exercises 1 1i) In this instance we analysed equation 1.) for a set of different starting values. We got the following plots: In each case the pendulum is simply oscillating about the it s equilibrium position as expected. 1ii) In this instance we analyses equation 1.1) against 1.) plotted θt) got the following results: As we can see for θ = 0 at the beginning the Non-Linear Linear case are quite similar as we expect because of the small angle approximation. However, as θ gets larger the trajectories differ as expected. Similar behaviour for θ = 0.. The dramatic difference occurs for θ = 3.14 which is evident from the plot. Essentially, this says that we must be very careful when applying the small angle approximation the approximation must really be quite small. 3

5 Out of the curiosity the Phase Portraits of the Non-Linear Undamped Pendulum were plotted: The red phase path corresponding to θ = 3.14) is identified as the separatrix. It is the boundary between oscillations of the nonlinear pendulum rotations. Exercise i) The pendulum.c program was edited to apply the Runge-Kutta method for β = 1, θ = 3.14, ω = 0 it was demonstrated with the below illustration that we get equivalent results as the trapezoid method: 4

6 iii) For values of β = 1, k = 0.5 for the damped pendulum we compared the linear the nonlinear cases. We got the following plot: As we can see, after a sufficient amount of time the behaviour of the damped linear pendulum approximates very well to the non-linear case which is as expected. iv) The Runge-Kutta program was changed to analyse the damped-driven oscillator. The first 500 iterations were neglected as it took time to settle into the attractor when we plotted the phase portraits. This was achieved by adding a simple if statement into the for loop. The Phase Portraits are plotted below for certain values of A, θ, ω 5

7 Conclusions The seperatrix was successfully found in the case of the Non-Linear undamped pendulum. This is a very hy tool for it tells us at what value of θ we should expect to transform from oscillitory motion to rotary motion. As expected this should be in around a value of π which is where the separatrix was found to be. After a sufficient amount of time, roughly 0 seconds in our case, the linearly damped pendulum had approximately the same behaviour of the nonlinearly damped pendulum. We expect this as a result of the damping which results in θ decreasing when θ becomes small enough i.e after 0 seconds) we can apply the small angle approximation sinθ) θ resulting in equivalent behaviour. Analysing the behaviour of the damped-driven oscillator using phase portraits we found that at A = 1.50 the system exhibited chaotic behaviour. Increasing the number of iterations the program was run through it was found that the plot became quite a bit denser which indicate that the path in phase space was never repeated which is an indication of chaotic behaviour. Period doubling was witnessed at A = 1.07 which also would indicate that we should see chaotic behaviour for larger values of A. Finally, it was noted that that the initial values of θ, ω settled into the attractor as long as A was not in the chaotic regime. 6

The Nonlinear Pendulum

The Nonlinear Pendulum The Nonlinear Pendulum - Pádraig Ó Conbhuí - 08531749 TP Monday 1. Abstract This experiment was performed to examine the effects that linearizing equations has on the accuracy of results and to find ways

More information

Mechanical Resonance and Chaos

Mechanical Resonance and Chaos Mechanical Resonance and Chaos You will use the apparatus in Figure 1 to investigate regimes of increasing complexity. Figure 1. The rotary pendulum (from DeSerio, www.phys.ufl.edu/courses/phy483l/group_iv/chaos/chaos.pdf).

More information

Computational Physics (6810): Session 8

Computational Physics (6810): Session 8 Computational Physics (6810): Session 8 Dick Furnstahl Nuclear Theory Group OSU Physics Department February 24, 2014 Differential equation solving Session 7 Preview Session 8 Stuff Solving differential

More information

Linear and Nonlinear Oscillators (Lecture 2)

Linear and Nonlinear Oscillators (Lecture 2) Linear and Nonlinear Oscillators (Lecture 2) January 25, 2016 7/441 Lecture outline A simple model of a linear oscillator lies in the foundation of many physical phenomena in accelerator dynamics. A typical

More information

Oscillatory Motion. Simple pendulum: linear Hooke s Law restoring force for small angular deviations. small angle approximation. Oscillatory solution

Oscillatory Motion. Simple pendulum: linear Hooke s Law restoring force for small angular deviations. small angle approximation. Oscillatory solution Oscillatory Motion Simple pendulum: linear Hooke s Law restoring force for small angular deviations d 2 θ dt 2 = g l θ small angle approximation θ l Oscillatory solution θ(t) =θ 0 sin(ωt + φ) F with characteristic

More information

Oscillatory Motion. Simple pendulum: linear Hooke s Law restoring force for small angular deviations. Oscillatory solution

Oscillatory Motion. Simple pendulum: linear Hooke s Law restoring force for small angular deviations. Oscillatory solution Oscillatory Motion Simple pendulum: linear Hooke s Law restoring force for small angular deviations d 2 θ dt 2 = g l θ θ l Oscillatory solution θ(t) =θ 0 sin(ωt + φ) F with characteristic angular frequency

More information

A plane autonomous system is a pair of simultaneous first-order differential equations,

A plane autonomous system is a pair of simultaneous first-order differential equations, Chapter 11 Phase-Plane Techniques 11.1 Plane Autonomous Systems A plane autonomous system is a pair of simultaneous first-order differential equations, ẋ = f(x, y), ẏ = g(x, y). This system has an equilibrium

More information

Nonlinear Oscillators: Free Response

Nonlinear Oscillators: Free Response 20 Nonlinear Oscillators: Free Response Tools Used in Lab 20 Pendulums To the Instructor: This lab is just an introduction to the nonlinear phase portraits, but the connection between phase portraits and

More information

Chaotic motion. Phys 420/580 Lecture 10

Chaotic motion. Phys 420/580 Lecture 10 Chaotic motion Phys 420/580 Lecture 10 Finite-difference equations Finite difference equation approximates a differential equation as an iterative map (x n+1,v n+1 )=M[(x n,v n )] Evolution from time t

More information

Transitioning to Chaos in a Simple Mechanical Oscillator

Transitioning to Chaos in a Simple Mechanical Oscillator Transitioning to Chaos in a Simple Mechanical Oscillator Hwan Bae Physics Department, The College of Wooster, Wooster, Ohio 69, USA (Dated: May 9, 8) We vary the magnetic damping, driver frequency, and

More information

Midterm EXAM PHYS 401 (Spring 2012), 03/20/12

Midterm EXAM PHYS 401 (Spring 2012), 03/20/12 Midterm EXAM PHYS 401 (Spring 2012), 03/20/12 Name: Signature: Duration: 75 minutes Show all your work for full/partial credit! In taking this exam you confirm to adhere to the Aggie Honor Code: An Aggie

More information

CHAOS -SOME BASIC CONCEPTS

CHAOS -SOME BASIC CONCEPTS CHAOS -SOME BASIC CONCEPTS Anders Ekberg INTRODUCTION This report is my exam of the "Chaos-part" of the course STOCHASTIC VIBRATIONS. I m by no means any expert in the area and may well have misunderstood

More information

CDS 101 Precourse Phase Plane Analysis and Stability

CDS 101 Precourse Phase Plane Analysis and Stability CDS 101 Precourse Phase Plane Analysis and Stability Melvin Leok Control and Dynamical Systems California Institute of Technology Pasadena, CA, 26 September, 2002. mleok@cds.caltech.edu http://www.cds.caltech.edu/

More information

Oscillations. Simple Harmonic Motion of a Mass on a Spring The equation of motion for a mass m is attached to a spring of constant k is

Oscillations. Simple Harmonic Motion of a Mass on a Spring The equation of motion for a mass m is attached to a spring of constant k is Dr. Alain Brizard College Physics I (PY 10) Oscillations Textbook Reference: Chapter 14 sections 1-8. Simple Harmonic Motion of a Mass on a Spring The equation of motion for a mass m is attached to a spring

More information

Numerical Methods for the Solution of Differential Equations

Numerical Methods for the Solution of Differential Equations Numerical Methods for the Solution of Differential Equations Markus Grasmair Vienna, winter term 2011 2012 Analytical Solutions of Ordinary Differential Equations 1. Find the general solution of the differential

More information

Lecture V : Oscillatory motion and spectral analysis

Lecture V : Oscillatory motion and spectral analysis Lecture V : Oscillatory motion and spectral analysis I. IDEAL PENDULUM AND STABILITY ANALYSIS Let us remind ourselves of the equation of motion for the pendulum. Remembering that the external torque applied

More information

Theoretical physics. Deterministic chaos in classical physics. Martin Scholtz

Theoretical physics. Deterministic chaos in classical physics. Martin Scholtz Theoretical physics Deterministic chaos in classical physics Martin Scholtz scholtzzz@gmail.com Fundamental physical theories and role of classical mechanics. Intuitive characteristics of chaos. Newton

More information

Lecture IV: Time Discretization

Lecture IV: Time Discretization Lecture IV: Time Discretization Motivation Kinematics: continuous motion in continuous time Computer simulation: Discrete time steps t Discrete Space (mesh particles) Updating Position Force induces acceleration.

More information

PHYSICS 210 SOLUTION OF THE NONLINEAR PENDULUM EQUATION USING FDAS

PHYSICS 210 SOLUTION OF THE NONLINEAR PENDULUM EQUATION USING FDAS PHYSICS 210 SOLUTION OF THE NONLINEAR PENDULUM EQUATION USING FDAS 1. PHYSICAL & MATHEMATICAL FORMULATION O θ L r T m W 1 1.1 Derivation of the equation of motion O Consider idealized pendulum: Mass of

More information

Checking the Radioactive Decay Euler Algorithm

Checking the Radioactive Decay Euler Algorithm Lecture 2: Checking Numerical Results Review of the first example: radioactive decay The radioactive decay equation dn/dt = N τ has a well known solution in terms of the initial number of nuclei present

More information

Chapter 2 PARAMETRIC OSCILLATOR

Chapter 2 PARAMETRIC OSCILLATOR CHAPTER- Chapter PARAMETRIC OSCILLATOR.1 Introduction A simple pendulum consists of a mass m suspended from a string of length L which is fixed at a pivot P. When simple pendulum is displaced to an initial

More information

Simple Harmonic Motion

Simple Harmonic Motion Simple Harmonic Motion (FIZ 101E - Summer 2018) July 29, 2018 Contents 1 Introduction 2 2 The Spring-Mass System 2 3 The Energy in SHM 5 4 The Simple Pendulum 6 5 The Physical Pendulum 8 6 The Damped Oscillations

More information

4. Complex Oscillations

4. Complex Oscillations 4. Complex Oscillations The most common use of complex numbers in physics is for analyzing oscillations and waves. We will illustrate this with a simple but crucially important model, the damped harmonic

More information

Lecture V: The game-engine loop & Time Integration

Lecture V: The game-engine loop & Time Integration Lecture V: The game-engine loop & Time Integration The Basic Game-Engine Loop Previous state: " #, %(#) ( #, )(#) Forces -(#) Integrate velocities and positions Resolve Interpenetrations Per-body change

More information

MATH 415, WEEKS 7 & 8: Conservative and Hamiltonian Systems, Non-linear Pendulum

MATH 415, WEEKS 7 & 8: Conservative and Hamiltonian Systems, Non-linear Pendulum MATH 415, WEEKS 7 & 8: Conservative and Hamiltonian Systems, Non-linear Pendulum Reconsider the following example from last week: dx dt = x y dy dt = x2 y. We were able to determine many qualitative features

More information

Physics 106b: Lecture 7 25 January, 2018

Physics 106b: Lecture 7 25 January, 2018 Physics 106b: Lecture 7 25 January, 2018 Hamiltonian Chaos: Introduction Integrable Systems We start with systems that do not exhibit chaos, but instead have simple periodic motion (like the SHO) with

More information

Chapter 14. PowerPoint Lectures for University Physics, Thirteenth Edition Hugh D. Young and Roger A. Freedman. Lectures by Wayne Anderson

Chapter 14. PowerPoint Lectures for University Physics, Thirteenth Edition Hugh D. Young and Roger A. Freedman. Lectures by Wayne Anderson Chapter 14 Periodic Motion PowerPoint Lectures for University Physics, Thirteenth Edition Hugh D. Young and Roger A. Freedman Lectures by Wayne Anderson Exam 3 results Class Average - 57 (Approximate grade

More information

STABILITY. Phase portraits and local stability

STABILITY. Phase portraits and local stability MAS271 Methods for differential equations Dr. R. Jain STABILITY Phase portraits and local stability We are interested in system of ordinary differential equations of the form ẋ = f(x, y), ẏ = g(x, y),

More information

P321(b), Assignement 1

P321(b), Assignement 1 P31(b), Assignement 1 1 Exercise 3.1 (Fetter and Walecka) a) The problem is that of a point mass rotating along a circle of radius a, rotating with a constant angular velocity Ω. Generally, 3 coordinates

More information

2.3 Damping, phases and all that

2.3 Damping, phases and all that 2.3. DAMPING, PHASES AND ALL THAT 107 2.3 Damping, phases and all that If we imagine taking our idealized mass on a spring and dunking it in water or, more dramatically, in molasses), then there will be

More information

Physics 101 Discussion Week 12 Explanation (2011)

Physics 101 Discussion Week 12 Explanation (2011) Physics 101 Discussion Week 12 Eplanation (2011) D12-1 Horizontal oscillation Q0. This is obviously about a harmonic oscillator. Can you write down Newton s second law in the (horizontal) direction? Let

More information

2nd Year Computational Physics Week 4 (standard): Ordinary differential equations

2nd Year Computational Physics Week 4 (standard): Ordinary differential equations 2nd Year Computational Physics Week 4 (standard): Ordinary differential equations Last compiled September 28, 2017 1 Contents 1 Introduction 4 2 Prelab Questions 4 3 Simple Differential Equations 5 3.1

More information

The dynamics of the Forced Damped Pendulum. John Hubbard Cornell University and Université de Provence

The dynamics of the Forced Damped Pendulum. John Hubbard Cornell University and Université de Provence The dynamics of the Forced Damped Pendulum John Hubbard Cornell University and Université de Provence Three ways to view the pendulum Three ways to view the pendulum 1. As a physical object Three ways

More information

Nonlinear Dynamic Systems Homework 1

Nonlinear Dynamic Systems Homework 1 Nonlinear Dynamic Systems Homework 1 1. A particle of mass m is constrained to travel along the path shown in Figure 1, which is described by the following function yx = 5x + 1x 4, 1 where x is defined

More information

MAS212 Assignment #2: The damped driven pendulum

MAS212 Assignment #2: The damped driven pendulum MAS Assignment #: The damped driven pendulum Dr. Sam Dolan (s.dolan@sheffield.ac.uk) Introduction: In this assignment you will use Python to investigate a non-linear differential equation which models

More information

Parametric Resonance and Elastic Pendulums

Parametric Resonance and Elastic Pendulums Parametric Resonance and Elastic Pendulums Ravitej Uppu Abstract In this I try to extend the theoretical conception of Elastic Pendulum that can be explained by the Driven Pendulums that I presented during

More information

Nonlinear Oscillations and Chaos

Nonlinear Oscillations and Chaos CHAPTER 4 Nonlinear Oscillations and Chaos 4-. l l = l + d s d d l l = l + d m θ m (a) (b) (c) The unetended length of each spring is, as shown in (a). In order to attach the mass m, each spring must be

More information

Wave Phenomena Physics 15c

Wave Phenomena Physics 15c Wave Phenomena Physics 15c Lecture Harmonic Oscillators (H&L Sections 1.4 1.6, Chapter 3) Administravia! Problem Set #1! Due on Thursday next week! Lab schedule has been set! See Course Web " Laboratory

More information

Euler s Method, cont d

Euler s Method, cont d Jim Lambers MAT 461/561 Spring Semester 009-10 Lecture 3 Notes These notes correspond to Sections 5. and 5.4 in the text. Euler s Method, cont d We conclude our discussion of Euler s method with an example

More information

ACM/CMS 107 Linear Analysis & Applications Fall 2016 Assignment 4: Linear ODEs and Control Theory Due: 5th December 2016

ACM/CMS 107 Linear Analysis & Applications Fall 2016 Assignment 4: Linear ODEs and Control Theory Due: 5th December 2016 ACM/CMS 17 Linear Analysis & Applications Fall 216 Assignment 4: Linear ODEs and Control Theory Due: 5th December 216 Introduction Systems of ordinary differential equations (ODEs) can be used to describe

More information

Computer Problems for Methods of Solving Ordinary Differential Equations

Computer Problems for Methods of Solving Ordinary Differential Equations Computer Problems for Methods of Solving Ordinary Differential Equations 1. Have a computer make a phase portrait for the system dx/dt = x + y, dy/dt = 2y. Clearly indicate critical points and separatrices.

More information

Class 4: More Pendulum results

Class 4: More Pendulum results Class 4: More Pendulum results The pendulum is a mechanical problem with a long and interesting history, from Galileo s first ansatz that the period was independent of the amplitude based on watching priests

More information

APPM 2460 CHAOTIC DYNAMICS

APPM 2460 CHAOTIC DYNAMICS APPM 2460 CHAOTIC DYNAMICS 1. Introduction Today we ll continue our exploration of dynamical systems, focusing in particular upon systems who exhibit a type of strange behavior known as chaos. We will

More information

Differential Equations: Homework 8

Differential Equations: Homework 8 Differential Equations: Homework 8 Alvin Lin January 08 - May 08 Section.6 Exercise Find a general solution to the differential equation using the method of variation of parameters. y + y = tan(t) r +

More information

28. Pendulum phase portrait Draw the phase portrait for the pendulum (supported by an inextensible rod)

28. Pendulum phase portrait Draw the phase portrait for the pendulum (supported by an inextensible rod) 28. Pendulum phase portrait Draw the phase portrait for the pendulum (supported by an inextensible rod) θ + ω 2 sin θ = 0. Indicate the stable equilibrium points as well as the unstable equilibrium points.

More information

The object of this experiment is to study systems undergoing simple harmonic motion.

The object of this experiment is to study systems undergoing simple harmonic motion. Chapter 9 Simple Harmonic Motion 9.1 Purpose The object of this experiment is to study systems undergoing simple harmonic motion. 9.2 Introduction This experiment will develop your ability to perform calculations

More information

Chasing Chaos With a Magnetic Pendulum

Chasing Chaos With a Magnetic Pendulum Chasing Chaos With a Magnetic Pendulum PHY 300 - Junior Phyics Laboratory Hassan Bukhari Roll no: 2012-10-0152 Department of Physcis LUMS SSE Saturday, October, 23, 2010 1 Introduction Chaos expresses

More information

Coordinate Curves for Trajectories

Coordinate Curves for Trajectories 43 The material on linearizations and Jacobian matrices developed in the last chapter certainly expanded our ability to deal with nonlinear systems of differential equations Unfortunately, those tools

More information

Periodic Skeletons of Nonlinear Dynamical Systems in the Problems of Global Bifurcation Analysis

Periodic Skeletons of Nonlinear Dynamical Systems in the Problems of Global Bifurcation Analysis Periodic Skeletons of Nonlinear Dynamical Systems in the Problems of Global Bifurcation Analysis M Zakrzhevsky, I Schukin, A Klokov and E Shilvan PERIODIC SKELETONS OF NONLINEAR DYNAMICAL SYSTEMS IN THE

More information

Transverse Linearization for Controlled Mechanical Systems with Several Passive Degrees of Freedom (Application to Orbital Stabilization)

Transverse Linearization for Controlled Mechanical Systems with Several Passive Degrees of Freedom (Application to Orbital Stabilization) Transverse Linearization for Controlled Mechanical Systems with Several Passive Degrees of Freedom (Application to Orbital Stabilization) Anton Shiriaev 1,2, Leonid Freidovich 1, Sergey Gusev 3 1 Department

More information

Assignments VIII and IX, PHYS 301 (Classical Mechanics) Spring 2014 Due 3/21/14 at start of class

Assignments VIII and IX, PHYS 301 (Classical Mechanics) Spring 2014 Due 3/21/14 at start of class Assignments VIII and IX, PHYS 301 (Classical Mechanics) Spring 2014 Due 3/21/14 at start of class Homeworks VIII and IX both center on Lagrangian mechanics and involve many of the same skills. Therefore,

More information

Chaos and Liapunov exponents

Chaos and Liapunov exponents PHYS347 INTRODUCTION TO NONLINEAR PHYSICS - 2/22 Chaos and Liapunov exponents Definition of chaos In the lectures we followed Strogatz and defined chaos as aperiodic long-term behaviour in a deterministic

More information

Basic Theory of Differential Equations

Basic Theory of Differential Equations page 104 104 CHAPTER 1 First-Order Differential Equations 16. The following initial-value problem arises in the analysis of a cable suspended between two fixed points y = 1 a 1 + (y ) 2, y(0) = a, y (0)

More information

AIMS Exercise Set # 1

AIMS Exercise Set # 1 AIMS Exercise Set #. Determine the form of the single precision floating point arithmetic used in the computers at AIMS. What is the largest number that can be accurately represented? What is the smallest

More information

Chaos. Lendert Gelens. KU Leuven - Vrije Universiteit Brussel Nonlinear dynamics course - VUB

Chaos. Lendert Gelens. KU Leuven - Vrije Universiteit Brussel   Nonlinear dynamics course - VUB Chaos Lendert Gelens KU Leuven - Vrije Universiteit Brussel www.gelenslab.org Nonlinear dynamics course - VUB Examples of chaotic systems: the double pendulum? θ 1 θ θ 2 Examples of chaotic systems: the

More information

221B Lecture Notes on Resonances in Classical Mechanics

221B Lecture Notes on Resonances in Classical Mechanics 1B Lecture Notes on Resonances in Classical Mechanics 1 Harmonic Oscillators Harmonic oscillators appear in many different contexts in classical mechanics. Examples include: spring, pendulum (with a small

More information

Nonlinear Autonomous Systems of Differential

Nonlinear Autonomous Systems of Differential Chapter 4 Nonlinear Autonomous Systems of Differential Equations 4.0 The Phase Plane: Linear Systems 4.0.1 Introduction Consider a system of the form x = A(x), (4.0.1) where A is independent of t. Such

More information

I. THE PENDULUM. y 1 (t) = θ(t) (2) y 2 (t) = dθ(t) dt

I. THE PENDULUM. y 1 (t) = θ(t) (2) y 2 (t) = dθ(t) dt 1 I. THE PENDULUM In this section, we will explore using Mathematica as a tool, several aspects of the pendulum that you might have done as anexperiment inthe lab. We begin by looking at the forces that

More information

Nonlinear System Analysis

Nonlinear System Analysis Nonlinear System Analysis Lyapunov Based Approach Lecture 4 Module 1 Dr. Laxmidhar Behera Department of Electrical Engineering, Indian Institute of Technology, Kanpur. January 4, 2003 Intelligent Control

More information

1 Pushing your Friend on a Swing

1 Pushing your Friend on a Swing Massachusetts Institute of Technology MITES 017 Physics III Lecture 05: Driven Oscillations In these notes, we derive the properties of both an undamped and damped harmonic oscillator under the influence

More information

Nonlinear Dynamics: Synchronisation

Nonlinear Dynamics: Synchronisation Nonlinear Dynamics: Synchronisation Bristol Centre for Complexity Sciences Ian Ross BRIDGE, School of Geographical Sciences, University of Bristol October 19, 2007 1 / 16 I: Introduction 2 / 16 I: Fireflies

More information

I ml. g l. sin. l 0,,2,...,n A. ( t vs. ), and ( vs. ). d dt. sin l

I ml. g l. sin. l 0,,2,...,n A. ( t vs. ), and ( vs. ). d dt. sin l Advanced Control Theory Homework #3 Student No.597 Name : Jinseong Kim. The simple pendulum dynamics is given by g sin L A. Derive above dynamic equation from the free body diagram. B. Find the equilibrium

More information

Numerical solution of ODEs

Numerical solution of ODEs Numerical solution of ODEs Arne Morten Kvarving Department of Mathematical Sciences Norwegian University of Science and Technology November 5 2007 Problem and solution strategy We want to find an approximation

More information

Introduction. Prediction MATH February 2017

Introduction. Prediction MATH February 2017 21 February 2017 Predicting the future is very difficult, especially if it s about the future. Niels Bohr Can we say what is going to happen: in the next minute? tomorrow? next year? Predicting the future

More information

1.11 Some Higher-Order Differential Equations

1.11 Some Higher-Order Differential Equations page 99. Some Higher-Order Differential Equations 99. Some Higher-Order Differential Equations So far we have developed analytical techniques only for solving special types of firstorder differential equations.

More information

More Details Fixed point of mapping is point that maps into itself, i.e., x n+1 = x n.

More Details Fixed point of mapping is point that maps into itself, i.e., x n+1 = x n. More Details Fixed point of mapping is point that maps into itself, i.e., x n+1 = x n. If there are points which, after many iterations of map then fixed point called an attractor. fixed point, If λ

More information

Chaotic motion. Phys 750 Lecture 9

Chaotic motion. Phys 750 Lecture 9 Chaotic motion Phys 750 Lecture 9 Finite-difference equations Finite difference equation approximates a differential equation as an iterative map (x n+1,v n+1 )=M[(x n,v n )] Evolution from time t =0to

More information

TP 3:Runge-Kutta Methods-Solar System-The Method of Least Squares

TP 3:Runge-Kutta Methods-Solar System-The Method of Least Squares TP :Runge-Kutta Methods-Solar System-The Method of Least Squares December 8, 2009 1 Runge-Kutta Method The problem is still trying to solve the first order differential equation dy = f(y, x). (1) dx In

More information

Mathematical Model of Forced Van Der Pol s Equation

Mathematical Model of Forced Van Der Pol s Equation Mathematical Model of Forced Van Der Pol s Equation TO Tsz Lok Wallace LEE Tsz Hong Homer December 9, Abstract This work is going to analyze the Forced Van Der Pol s Equation which is used to analyze the

More information

B2.III Revision notes: quantum physics

B2.III Revision notes: quantum physics B.III Revision notes: quantum physics Dr D.M.Lucas, TT 0 These notes give a summary of most of the Quantum part of this course, to complement Prof. Ewart s notes on Atomic Structure, and Prof. Hooker s

More information

A645/A445: Exercise #1. The Kepler Problem

A645/A445: Exercise #1. The Kepler Problem A645/A445: Exercise #1 The Kepler Problem Due: 2017 Sep 26 1 Numerical solution to the one-degree-of-freedom implicit equation The one-dimensional implicit solution is: t = t o + x x o 2(E U(x)). (1) The

More information

The Harmonic Oscillator

The Harmonic Oscillator The Harmonic Oscillator Math 4: Ordinary Differential Equations Chris Meyer May 3, 008 Introduction The harmonic oscillator is a common model used in physics because of the wide range of problems it can

More information

Wave Phenomena Physics 15c

Wave Phenomena Physics 15c Wave Phenomena Physics 15c Lecture Harmonic Oscillators (H&L Sections 1.4 1.6, Chapter 3) Administravia! Problem Set #1! Due on Thursday next week! Lab schedule has been set! See Course Web " Laboratory

More information

LECTURE 8: DYNAMICAL SYSTEMS 7

LECTURE 8: DYNAMICAL SYSTEMS 7 15-382 COLLECTIVE INTELLIGENCE S18 LECTURE 8: DYNAMICAL SYSTEMS 7 INSTRUCTOR: GIANNI A. DI CARO GEOMETRIES IN THE PHASE SPACE Damped pendulum One cp in the region between two separatrix Separatrix Basin

More information

LINEAR RESPONSE THEORY

LINEAR RESPONSE THEORY MIT Department of Chemistry 5.74, Spring 5: Introductory Quantum Mechanics II Instructor: Professor Andrei Tokmakoff p. 8 LINEAR RESPONSE THEORY We have statistically described the time-dependent behavior

More information

Chapter 16: Oscillations

Chapter 16: Oscillations Chapter 16: Oscillations Brent Royuk Phys-111 Concordia University Periodic Motion Periodic Motion is any motion that repeats itself. The Period (T) is the time it takes for one complete cycle of motion.

More information

Nonlinear Dynamics Semi-classical Model of Quantum Spin

Nonlinear Dynamics Semi-classical Model of Quantum Spin 1 Nonlinear Dynamics Semi-classical Model of Quantum Spin Joshua J. Heiner David R. Thayer Department of Physics and Astronomy Department of Physics and Astronomy University of Wyoming University of Wyoming

More information

Introduction to Dynamical Systems Basic Concepts of Dynamics

Introduction to Dynamical Systems Basic Concepts of Dynamics Introduction to Dynamical Systems Basic Concepts of Dynamics A dynamical system: Has a notion of state, which contains all the information upon which the dynamical system acts. A simple set of deterministic

More information

Oscillations Simple Harmonic Motion

Oscillations Simple Harmonic Motion Oscillations Simple Harmonic Motion Lana Sheridan De Anza College Dec 1, 2017 Overview oscillations simple harmonic motion (SHM) spring systems energy in SHM pendula damped oscillations Oscillations and

More information

Chap. 20: Initial-Value Problems

Chap. 20: Initial-Value Problems Chap. 20: Initial-Value Problems Ordinary Differential Equations Goal: to solve differential equations of the form: dy dt f t, y The methods in this chapter are all one-step methods and have the general

More information

Chapter 2: Complex numbers

Chapter 2: Complex numbers Chapter 2: Complex numbers Complex numbers are commonplace in physics and engineering. In particular, complex numbers enable us to simplify equations and/or more easily find solutions to equations. We

More information

Maps and differential equations

Maps and differential equations Maps and differential equations Marc R. Roussel November 8, 2005 Maps are algebraic rules for computing the next state of dynamical systems in discrete time. Differential equations and maps have a number

More information

Appendix 3B MATLAB Functions for Modeling and Time-domain analysis

Appendix 3B MATLAB Functions for Modeling and Time-domain analysis Appendix 3B MATLAB Functions for Modeling and Time-domain analysis MATLAB control system Toolbox contain the following functions for the time-domain response step impulse initial lsim gensig damp ltiview

More information

Figure 1: Schematic of ship in still water showing the action of bouyancy and weight to right the ship.

Figure 1: Schematic of ship in still water showing the action of bouyancy and weight to right the ship. MULTI-DIMENSIONAL SYSTEM: In this computer simulation we will explore a nonlinear multi-dimensional system. As before these systems are governed by equations of the form x 1 = f 1 x 2 = f 2.. x n = f n

More information

Mathematics for chemical engineers. Numerical solution of ordinary differential equations

Mathematics for chemical engineers. Numerical solution of ordinary differential equations Mathematics for chemical engineers Drahoslava Janovská Numerical solution of ordinary differential equations Initial value problem Winter Semester 2015-2016 Outline 1 Introduction 2 One step methods Euler

More information

Exam in Systems Engineering/Process Control

Exam in Systems Engineering/Process Control Department of AUTOMATIC CONTROL Exam in Systems Engineering/Process Control 27-6-2 Points and grading All answers must include a clear motivation. Answers may be given in English or Swedish. The total

More information

Second quantization: where quantization and particles come from?

Second quantization: where quantization and particles come from? 110 Phys460.nb 7 Second quantization: where quantization and particles come from? 7.1. Lagrangian mechanics and canonical quantization Q: How do we quantize a general system? 7.1.1.Lagrangian Lagrangian

More information

Numerical Methods for Initial Value Problems; Harmonic Oscillators

Numerical Methods for Initial Value Problems; Harmonic Oscillators Lab 1 Numerical Methods for Initial Value Problems; Harmonic Oscillators Lab Objective: Implement several basic numerical methods for initial value problems (IVPs), and use them to study harmonic oscillators.

More information

PH36010: Numerical Methods - Evaluating the Lorenz Attractor using Runge-Kutta methods Abstract

PH36010: Numerical Methods - Evaluating the Lorenz Attractor using Runge-Kutta methods Abstract PH36010: Numerical Methods - Evaluating the Lorenz Attractor using Runge-Kutta methods Mr. Benjamen P. Reed (110108461) IMPACS, Aberystwyth University January 31, 2014 Abstract A set of three coupled ordinary

More information

In the presence of viscous damping, a more generalized form of the Lagrange s equation of motion can be written as

In the presence of viscous damping, a more generalized form of the Lagrange s equation of motion can be written as 2 MODELING Once the control target is identified, which includes the state variable to be controlled (ex. speed, position, temperature, flow rate, etc), and once the system drives are identified (ex. force,

More information

Lecture XXVI. Morris Swartz Dept. of Physics and Astronomy Johns Hopkins University November 5, 2003

Lecture XXVI. Morris Swartz Dept. of Physics and Astronomy Johns Hopkins University November 5, 2003 Lecture XXVI Morris Swartz Dept. of Physics and Astronomy Johns Hopins University morris@jhu.edu November 5, 2003 Lecture XXVI: Oscillations Oscillations are periodic motions. There are many examples of

More information

Ratchet Baryogenesis during Reheating

Ratchet Baryogenesis during Reheating Ratchet Baryogenesis during Reheating Neil D. Barrie Kavli IPMU (WPI) February 17, 2018 K. Bamba, N. B., A. Sugamoto, T. Takeuchi and K. Yamashita, arxiv:1610.03268, and upcoming work. N. D. Barrie (Kavli

More information

Question: Total. Points:

Question: Total. Points: MATH 308 May 23, 2011 Final Exam Name: ID: Question: 1 2 3 4 5 6 7 8 9 Total Points: 0 20 20 20 20 20 20 20 20 160 Score: There are 9 problems on 9 pages in this exam (not counting the cover sheet). Make

More information

Solution Set Five. 2 Problem #2: The Pendulum of Doom Equation of Motion Simple Pendulum Limit Visualization...

Solution Set Five. 2 Problem #2: The Pendulum of Doom Equation of Motion Simple Pendulum Limit Visualization... : Solution Set Five Northwestern University, Classical Mechanics Classical Mechanics, Third Ed.- Goldstein November 4, 2015 Contents 1 Problem #1: Coupled Mass Oscillator System. 2 1.1 Normal Modes.......................................

More information

Chapter 11 ORDINARY DIFFERENTIAL EQUATIONS

Chapter 11 ORDINARY DIFFERENTIAL EQUATIONS Chapter 11 ORDINARY DIFFERENTIAL EQUATIONS The general form of a first order differential equations is = f(x, y) with initial condition y(a) = y a We seek the solution y = y(x) for x > a This is shown

More information

CALIFORNIA INSTITUTE OF TECHNOLOGY Control and Dynamical Systems

CALIFORNIA INSTITUTE OF TECHNOLOGY Control and Dynamical Systems CDS 101 1. Åström and Murray, Exercise 1.3 2. Åström and Murray, Exercise 1.4 3. Åström and Murray, Exercise 2.6, parts (a) and (b) CDS 110a 1. Åström and Murray, Exercise 1.4 2. Åström and Murray, Exercise

More information

Consistency and Convergence

Consistency and Convergence Jim Lambers MAT 77 Fall Semester 010-11 Lecture 0 Notes These notes correspond to Sections 1.3, 1.4 and 1.5 in the text. Consistency and Convergence We have learned that the numerical solution obtained

More information

MAS212 Assignment #2: The damped driven pendulum

MAS212 Assignment #2: The damped driven pendulum MAS Assignment #: The damped driven pendulum Sam Dolan (January 8 Introduction In this assignment we study the motion of a rigid pendulum of length l and mass m, shown in Fig., using both analytical and

More information

1 (2n)! (-1)n (θ) 2n

1 (2n)! (-1)n (θ) 2n Complex Numbers and Algebra The real numbers are complete for the operations addition, subtraction, multiplication, and division, or more suggestively, for the operations of addition and multiplication

More information

Do not turn over until you are told to do so by the Invigilator.

Do not turn over until you are told to do so by the Invigilator. UNIVERSITY OF EAST ANGLIA School of Mathematics Main Series UG Examination 216 17 INTRODUCTION TO NUMERICAL ANALYSIS MTHE612B Time allowed: 3 Hours Attempt QUESTIONS 1 and 2, and THREE other questions.

More information