Dynamical Systems: how to find equation for Poincare map?

Click For Summary
The discussion revolves around deriving the Poincare map for a dynamical system described in polar coordinates, specifically focusing on the equations \(\frac{dr}{dt} = r(a - r)\) and \(\frac{d\theta}{dt} = -1\). The solution to these equations indicates that all trajectories converge to a limit cycle as \(t\) approaches infinity for \(a > 0\). The user seeks clarification on how to express the Poincare map \(P(r)\) in terms of the previously derived solution, particularly how to relate \(r_n\) to \(P(r_{n-1})\). They note that the Poincare map captures intersections of trajectories with a defined plane, and they emphasize the need to evaluate \(r(2\pi)\) based on the previous state \(r_{n-1}\). The discussion highlights the challenge of connecting the derived expression to the Poincare map definition.
Master1022
Messages
590
Reaction score
116
Homework Statement
Using the solution given in the previous part, show that for ## a > 0 ## one can define the Poincare map: ## P: \mathbb{R}^{+} \righarrow \mathbb{R}^{+} ## by: (expression below)
Relevant Equations
Differential equations
Hi,

I was attempting a question on the dynamical systems topic of Poincare maps, and was struggling to understand a certain part of it.

Knowledge from prior parts of the questions:
There was a system which we converted to polar coordinates to get: (## a ## is an arbitrary real constant)
\frac{dr}{dt} = r (a - r)
\frac{d\theta}{dt} = -1

and the first equation was solved with initial condition: ## r(t = 0) = r_0 ## and ## \theta(t = 0 ) = 0 ## to give:
r(t) = \frac{a r_0}{r_0 + (a - r_0) e^{-at}}
\theta(t) = -t
We can see that all solutions tend to a limit cycle (## r = a ##) as ## t \rightarrow \infty## for ## a > 0 ##.

Question I am stuck on: Using the solution given in the previous part, show that for ## a > 0 ## one can define the Poincare map: ## P: \mathbb{R}^{+} \rightarrow \mathbb{R}^{+} ## by ## r_n = P(r_{n - 1}) ## with:
P(r) = \frac{a r}{r + (a - r) e^{-2 a \pi}}

Attempt:
This should be quite a simple thing to do, but I don't quite understand each step in the logic to get there. So I know the Poincare map can be thought of as a plane where we mark the intersections of the trajectory with the plane. This means that a fixed point on a Poincare map corresponds to a limit cycle for the full trajectory.

This plane will always be at the same angle ## \pm 2 n \pi ## (the x-axis). Are we interested in all the times ## t ## when ## = - 2 n \pi ##? These times ## t ## can be found by doing: ## -t = -2 n \pi \rightarrow t = 2 n \pi ##.

Substituting this into the expression for ## r(t) ## yields:
r(t) = \frac{a r_0}{r_0 + (a - r_0) e^{-2a n \pi}} \rightarrow r_n = \frac{a r_0}{r_0 + (a - r_0) e^{-2a n \pi}}

However, I am not sure how to show that the expression is equal to ## P(r)##. What is the reasoning behind that?

Thanks in advance
 
Physics news on Phys.org
Don't solve for r_n = r(2n\pi) subject to r(0) = r_0. Instead solve for r(2\pi) subject to r(0) = r_{n-1}.
 
First, I tried to show that ##f_n## converges uniformly on ##[0,2\pi]##, which is true since ##f_n \rightarrow 0## for ##n \rightarrow \infty## and ##\sigma_n=\mathrm{sup}\left| \frac{\sin\left(\frac{n^2}{n+\frac 15}x\right)}{n^{x^2-3x+3}} \right| \leq \frac{1}{|n^{x^2-3x+3}|} \leq \frac{1}{n^{\frac 34}}\rightarrow 0##. I can't use neither Leibnitz's test nor Abel's test. For Dirichlet's test I would need to show, that ##\sin\left(\frac{n^2}{n+\frac 15}x \right)## has partialy bounded sums...