« View all new features in Mathematica 9  previous  |  next 
New in Mathematica 9Parametric Differential Equations

Sensitivity of the Lorenz Equations 

Visualize the sensitivity of the Lorenz equations with respect to a parameter.

In[1]:=
Click for copyable input
X
In[2]:=
Click for copyable input
X

A plot of the solution shows a part of the classic Lorenz attractor. The color changes from red to blue as time increases.

In[3]:=
Click for copyable input
X
Out[3]=

The derivative with respect to shows how sensitive the solution is to changes in . The colors correspond to the same time values as in the plot above.

In[4]:=
Click for copyable input
X
Out[4]=

Two solutions with nearby deviate quickly. The large sensitivity is an indication that this is likely.

In[5]:=
Click for copyable input
X
Out[5]=

Another way to see the sensitivities is to consider a coordinate basis local to a point on the trajectory, so that one coordinate points along the trajectory and two are orthogonal. The orthogonal directions are chosen so that one component is in the direction of the nearest equilibrium point, since around these points the trajectory stays near a plane.

In[6]:=
Click for copyable input
X
Out[6]=
In[7]:=
Click for copyable input
X
In[8]:=
Click for copyable input
X
Out[8]=

The largest part of the sensitivity is in the direction along the trajectory.

Using the decomposition of the sensitivity above, a plot can be made that shows the scaled magnitude of the sensitivity in the directions orthogonal to the trajectory.

In[9]:=
Click for copyable input
X
Out[9]=