This preview shows pages 1–3. Sign up to view the full content.
This preview has intentionally blurred sections. Sign up to view the full version.
View Full Document
Unformatted text preview: Solution to Homework #2, 36754 7 February 2006 Exercise 5.3 (The Logistic Map as a Measure Preserving Transformation) The logistic map with a = 4 is a measurepreserving transfor mation, and the measure it preserves has the density 1 / p x (1 x ) (on the unit interval. 1. Verify that this density is invariant under the action of the logistic map. 2. Simulate the logistic map with uniformly distributed X . What happens to the density of X t as t ? a There are a number of ways to do this; here is one which will feed into later material. It does not use the fact that the invariant density is the (1 / 2 , 1 / 2) distribution, or that it has a closedform integral, Z x dy p y (1 y ) = 2 sin 1 x though both of these are fine facts to have handy. Lets write the mapping as F ( x ) = 4 x (1 x ). Solving a simple quadratic equation gives us the fact that F 1 ( x ) is the set 1 2 ( 1 1 x ) , 1 2 ( 1 + 1 x ) . Notice, for later use, that the two solutions add up to 1. Notice also that F 1 ([0 ,x ]) = 0 , 1 2 ( 1 1 x ) 1 2 ( 1 + 1 x ) , 1 . Now we consider P ( X n +1 x ), the cumulative distribution function of X n +1 . P ( X n +1 x ) = P ( X n +1 [0 ,x ]) = P ( X n F 1 ([0 ,x ]) ) = P X n , 1 2 ( 1 1 x ) 1 2 ( 1 + 1 x ) , 1 = Z 1 2 ( 1 1 x ) n ( y ) dy + Z 1 1 2 ( 1+ 1 x ) n ( y ) dy 1 where n is the density of X n . So we have an integral equation for the evolution of the density, Z x n +1 ( y ) dy = Z 1 2 ( 1 1 x ) n ( y ) dy + Z 1 1 2 ( 1+ 1 x ) n ( y ) dy This sort of integral equation is complicated to solve directly. Instead, take the derivative of both sides with respect to x ; we can do this through the fundamental theorem of calculus. On the left hand side, this will just give n +1 ( x ), the density we want. n +1 ( x ) = d dx Z 1 2 ( 1 1 x ) n ( y ) dy + d dx Z 1 1 2 ( 1+ 1 x ) n ( y ) dy = n 1 2 ( 1 1 x ) d dx 1 2 ( 1 1 x ) n 1 2 ( 1 + 1 x ) d dx 1 2 ( 1 + 1 x ) = 1 4 1 x n 1 2 ( 1 1 x ) + n 1 2 ( 1 + 1 x ) Notice that this defines a linear operator taking densities to densities. (YouNotice that this defines a linear operator taking densities to densities....
View Full
Document
 Spring '06
 Schalizi

Click to edit the document details