Warning: foreach() argument must be of type array|object, bool given in /var/www/html/web/app/themes/studypress-core-theme/template-parts/header/mobile-offcanvas.php on line 20

For certain \(r\) intervals, or windows, the Lorenz equations exhibit a period- doubling property similar to that of the logistic difference equation discussed in Section \(2.9 .\) Careful calculations may reveal this phenomenon. (a) One period-doubling window contains the value \(r=100 .\) Let \(r=100\) and plot the trajectory starting at \((5,5,5)\) or some other initial point of your choice. Does the solution appear to be periodic? What is the period? (b) Repeat the calculation in part (a) for slightly smaller values of \(r .\) When \(r \cong 99.98\), you may be able to observe that the period of the solution doubles. Try to observe this result by performing calculations with nearby values of \(r\). (c) As \(r\) decreases further, the period of the solution doubles repeatedly. The next period doubling occurs at about \(r=99.629 .\) Try to observe this by plotting trajectories for nearby values of \(r .\)

Short Answer

Expert verified
Answer: As the value of r decreases in specific r intervals, the period of the solution of the Lorenz equations doubles, demonstrating the period-doubling property.

Step by step solution

Achieve better grades quicker with Premium

  • Unlimited AI interaction
  • Study offline
  • Say goodbye to ads
  • Export flashcards

Over 22 million students worldwide already upgrade their learning with Vaia!

01

Understanding the Lorenz equations

The Lorenz equations are given by: dx/dt = σ(y-x) dy/dt = rx - y - xz dz/dt = xy - bz Where x, y, z are the dependent variables, t is the independent variable (time), and σ, r, b are constants. We will use these equations to study the period-doubling phenomenon with different values of r.
02

Plotting the trajectory for r=100

Given r=100 and the initial point (5,5,5), we need to solve the Lorenz equations numerically and plot the trajectory of the solution. This can be done using numerical methods like Euler's method, Runge-Kutta method, or any other suitable numerical solver. After solving and plotting the trajectory, we should observe whether the solution appears to be periodic or not.
03

Determine the period for r=100

To find the period of the solution for r=100, we should analyze the plotted trajectory and look for repeating patterns. More precisely, we should look for the time it takes for the function to repeat itself, which is the period.
04

Plotting the trajectory for a smaller value of r

Next, let's consider r≈99.98 and plot the trajectory of the solution starting from the same initial point (5,5,5). We can follow the same process as in Step 2, using a numerical solver for the Lorenz equations.
05

Observe the period doubling

Now we should analyze the plotted trajectory for r≈99.98 and determine the period. At this point, we should observe that the period has doubled as compared to the solution with r=100.
06

Plotting the trajectory for r≈99.629

Finally, we will work with r≈99.629. We should solve the Lorenz equations and plot the trajectory of the solution, just like we did in Steps 2 and 4.
07

Observe the period doubling again

Now, let's analyze the plotted trajectory for r≈99.629 and find the period of the solution. We should be able to see that the period has doubled again compared to the solution with r≈99.98. In conclusion, we have demonstrated the period-doubling property of the Lorenz equations for certain r intervals by observing the changes in the period of the solution as r decreases.

Key Concepts

These are the key concepts you need to understand to accurately answer the question.

Period-Doubling
Period-doubling is a fascinating phenomenon often encountered in dynamical systems, which can be described as the process where the period of a system's oscillation doubles as a parameter, such as 'r' in the Lorenz equations, is changed. This intriguing behavior is a hallmark of chaos, indicative of a system transitioning towards chaotic behavior as the bifurcation parameter passes through a critical value.

Imagine a predictable oscillation completing one cycle in a given time. As the parameter changes, new cycles emerge which take twice as long to complete, hence the term 'period-doubling'. As we continue to adjust the parameter, further period-doublings can occur, leading to an exponential increase in the number of periods. Relating to the exercise, by adjusting the parameter 'r', students are able to witness the emergence and evolution of period-doubling in the Lorenz system, a clear illustration of complex non-linear dynamics at work.
Numerical Methods
Numerical methods are algorithms used to solve mathematical problems that are difficult or impossible to tackle analytically. In the context of differential equations, such as the Lorenz equations, numerical solutions become indispensable due to the complexity of the system.

When dealing with the Lorenz equations, the common numerical methods include Euler's Method, Runge-Kutta Method, and others. These are iterative techniques which generate approximations to the solution over time. They are vital in exploring dynamic systems like the Lorenz attractor, as they enable students to compute the trajectories and observe the evolution of the system over time without requiring a closed-form solution. The reliability and ease with which numerical methods can be implemented, often via computer software, make them an essential tool in both academic studies and practical applications.
Trajectory Plotting
Trajectory plotting is a visual representation of the path that a dynamic system follows through its phase space over time. For the Lorenz equations, the phase space is three-dimensional, spanned by the variables x, y, and z.

It's a tool that can provide profound insights into the behavior of complex systems. By plotting the trajectory of solutions to the Lorenz equations at different values of 'r', students can visually grasp the behavior changes from periodic to chaotic. With initial conditions set, such as \( (5,5,5) \) in the exercise, and after running simulations using numerical methods, the generated data points can be plotted to reveal the intricate structure of the Lorenz attractor. Students can then compare these structures to identify period-doubling transitions and other dynamical phenomena.
Differential Equations
Differential equations are mathematical equations that relate a function to its derivatives, expressing how the rate of change of a quantity depends on other quantities.

For example, the Lorenz equations are a set of three coupled, non-linear differential equations that model certain types of atmospheric convection. They are a fundamental piece of the chaos theory puzzle and a great example of how simple rules can lead to complex and unpredictable behavior. Understanding these equations involves studying the rates at which variables change concerning one another, which is a key component in a multitude of scientific disciplines. When students approach these differential equations with an eye toward both their mathematical intricacy and their representation of real-world phenomena, they gain a deeper appreciation for the complex interplay between parameters and system behavior.

One App. One Place for Learning.

All the tools & learning materials you need for study success - in one app.

Get started for free

Most popular questions from this chapter

(a) A special case of the Lienard equation of Problem 8 is $$ \frac{d^{2} u}{d t^{2}}+\frac{d u}{d t}+g(u)=0 $$ where \(g\) satisfies the conditions of Problem 6 . Letting \(x=u, y=d u / d t,\) show that the origin is a critical point of the resulting system. This equation can be interpreted as describing the motion of a spring-mass system with damping proportional to the velocity and a nonlinear restoring force. Using the Liapunov function of Problem \(6,\) show that the origin is a stable critical point, but note that even with damping we cannot conclude asymptotic stability using this Liapunov function. (b) Asymptotic stability of the critical point \((0,0)\) can be shown by constructing a better Liapunov function as was done in part (d) of Problem 7 . However, the analysis for a general function \(g\) is somewhat sophisticated and we only mention that appropriate form for \(V\) is $$ V(x, y)=\frac{1}{2} y^{2}+A y g(x)+\int_{0}^{x} g(s) d s $$ where \(A\) is a positive constant to be chosen so that \(V\) is positive definite and \(\hat{V}\) is negative definite. For the pendulum problem \([g(x)=\sin x]\) use \(V\) as given by the preceding equation with \(A=\frac{1}{2}\) to show that the origin is asymptotically stable. Hint: Use \(\sin x=x-\alpha x^{3} / 3 !\) and \(\cos x=1-\beta x^{2} / 2 !\) where \(\alpha\) and \(\beta\) depend on \(x,\) but \(0<\alpha<1\) and \(0<\beta<1\) for \(-\pi / 2

show that the given system has no periodic solutions other than constant solutions. $$ d x / d t=x+y+x^{3}-y^{2}, \quad d y / d t=-x+2 y+x^{2} y+y^{3} / 3 $$

The equation of motion of an undamped pendulum is \(d^{2} \theta / d t^{2}+\omega^{2} \sin \theta=0,\) where \(\omega^{2}=g / L .\) Let \(x=\theta, y=d \theta / d t\) to obtain the system of equations $$ d x / d t=y, \quad d y / d t=-\omega^{2} \sin x $$ (a) Show that the critical points are \((\pm n \pi, 0), n=0,1,2, \ldots,\) and that the system is almost lincar in the neighborhood of cach critical point. (b) Show that the critical point \((0,0)\) is a (stable) center of the corresponding linear system. Using Theorem 9.3.2 what can be said about the nonlinear system? The situation is similar at the critical points \((\pm 2 n \pi, 0), n=1,2,3, \ldots\) What is the physical interpretation of these critical points? (c) Show that the critical point \((\pi, 0)\) is an (unstable) saddle point of the corresponding linear system. What conclusion can you draw about the nonlinear system? The situation is similar at the critical points \([\pm(2 n-1) \pi, 0], n=1,2,3, \ldots\) What is the physical interpretation of these critical points? (d) Choose a value for \(\omega^{2}\) and plot a few trajectories of the nonlinear system in the neighborhood of the origin. Can you now draw any further conclusion about the nature of the critical point at \((0,0)\) for the nonlinear system? (e) Using the value of \(\omega^{2}\) from part (d) draw a phase portrait for the pendulum. Compare your plot with Figure 9.3 .5 for the damped pendulum.

(a) Find an equation of the form \(H(x, y)=c\) satisfied by the trajectories. (b) Plot several level curves of the function \(H\). These are trajectories of the given system. Indicate the direction of motion on each trajectory. $$ d x / d t=2 y, \quad d y / d t=8 x $$

(a) Determine all critical points of the given system of equations. (b) Find the corresponding linear system near each critical point. (c) Find the eigenalues of each linear system. What conclusions can you then draw about the nonlinear system? (d) Draw a phase portrait of the nonlinear system to confirm your conclusions, or to extend them in those cases where the linear system does not provide definite information about the nonlinear system. $$ d x / d t=-(x-y)(1-x-y), \quad d y / d t=x(2+y) $$

See all solutions

Recommended explanations on Math Textbooks

View all explanations

What do you think about this solution?

We value your feedback to improve our textbook solutions.

Study anywhere. Anytime. Across all devices.

Sign-up for free