To show that the limit set for the Poincare map \({{\bf{x}}_{\bf{n}}}{\bf{ = x(2\pi n),}}{{\bf{v}}_{\bf{n}}}{\bf{ = x'(2\pi n)}}\) where x(t) is a solution to equation (6), is an ellipse and that this ellipse is the same for any initial values\({\bf{(}}{{\bf{x}}_{\bf{o}}}{\bf{,}}{{\bf{\nu }}_{\bf{o}}}{\bf{)}}\), do the following:

(a)Argue that since the initial values affect only the transient solution to (6), the limit set for the Poincare map is independent of the initial values.

(b)Now show that for nlarge\({{\bf{x}}_{\bf{n}}} \approx {\bf{asin(2}}\sqrt {\bf{2}} {\bf{\pi n + }}\psi {\bf{)}},{{\bf{v}}_{\bf{n}}} \approx {\bf{c + }}\sqrt {\bf{2}} {\bf{acos(2}}\sqrt {\bf{2}} {\bf{\pi n + }}\psi {\bf{)}}\),

Where\({\bf{a = (1 + 2(0}}{\bf{.22}}{{\bf{)}}^{\bf{2}}}{{\bf{)}}^{\frac{{{\bf{ - 1}}}}{{\bf{2}}}}}{\bf{,c = (0}}{\bf{.22}}{{\bf{)}}^{{\bf{ - 1}}}}\), and\(\psi {\bf{ = arctan}}\left\{ {{\bf{ - [(0}}{\bf{.22)}}\sqrt {\bf{2}} {{\bf{]}}^{{\bf{ - 1}}}}} \right\}\).

(c)Use the result of part (b) to conclude that the ellipse\({{\bf{x}}^{\bf{2}}}{\bf{ + }}\frac{{{{{\bf{(v - c)}}}^{\bf{2}}}}}{{\bf{2}}}{\bf{ = }}{{\bf{a}}^{\bf{2}}}\)contains the limit set of the Poincare map.

Short Answer

Expert verified
  1. As \(t \to \infty \)the results tend to zero. So the points \(\left( {x\left( t \right),v\left( t \right)} \right)\) are independent of the initial conditions
  2. The values are

\(\begin{array}{l}{x_n} = a\sin \left( {\sqrt 2 \left( {2\pi n} \right) + \varphi } \right)\\{v_n} = \frac{1}{{0.22}} + \frac{{\sqrt 2 }}{{\sqrt {1 + 2{{(0.22)}^2}} }}\cos \left( {\sqrt 2 \left( {2\pi n} \right) + \varphi } \right)\end{array}\)

  1. \({x^2}_n + \frac{{{{\left( {{v_n} - c} \right)}^2}}}{2} = {a^2}\)having center\(\left( {0,c} \right)\).

Step by step solution

01

(a)Step 1: Show that the map is independent of initial conditions

The solution to equation (6) is given by \(x(t) = {x_h}(t) + {x_p}(t)\).where

\(\begin{array}{c}{x_h} = A{e^{ - 0.11t}}\sin \left( {\sqrt {9879} t + \phi } \right)\\{x_p} = \frac{1}{{0.22}}\sin t + \frac{1}{{\sqrt {1 + 2(0.22)} }}\sin \left( {\sqrt 2 t + \phi } \right)\\\tan \varphi = - \frac{1}{{0.22\sqrt 2 }}\end{array}\)

Is a steady state term.

Now, differentiate \(x\left( t \right)\)then:

\(\begin{array}{c}v\left( t \right) = x{'_h} + x{'_p}\\v\left( t \right) = x{'_h} + \frac{1}{{0.22}}\cos t + \frac{{\sqrt 2 }}{{\sqrt {1 + 2{{\left( {0.22} \right)}^2}} }}\cos \left( {\sqrt 2 t + \phi } \right)\end{array}\)

The steady solution doesn't depend upon the initial conditions. These values affect only constants A and \(\phi \) .

Hence,\(t \to \infty \), the results tend to be zero. So the points \(\left( {x\left( t \right),v\left( t \right)} \right)\) are independent of the initial conditions.

02

(b) Step 2: Find the value of \({\bf{(}}{{\bf{x}}_{\bf{n}}}{\bf{,}}{{\bf{v}}_{\bf{n}}}{\bf{)}}\)

Put \(t = 2\pi n\) then

\(\begin{array}{l}{x_n} = x\left( {2\pi n} \right)\\{x_n} = {x_h}\left( {2\pi n} \right) + \frac{1}{{\sqrt {1 + 2{{\left( {0.22} \right)}^2}} }}\sin \left( {\sqrt 2 \left( {2\pi n} \right) + \varphi } \right)\\{v_n} = v\left( {2\pi n} \right)\\{v_n} = x{'_h}\left( {2\pi n} \right) + \frac{1}{{0.22}} + \frac{{\sqrt 2 }}{{\sqrt {1 + 2{{\left( {0.22} \right)}^2}} }}\cos \left( {\sqrt 2 \left( {2\pi n} \right) + \varphi } \right)\end{array}\)

As \(n \to \infty ,{x_h}(2\pi n) \to 0,x{'_h}(2\pi n) \to 0\)

\(\begin{array}{l}{x_n} = \frac{1}{{\sqrt {1 + 2{{\left( {0.22} \right)}^2}} }}\sin \left( {\sqrt 2 \left( {2\pi n} \right) + \varphi } \right)\\{x_n} = a\sin \left( {\sqrt 2 \left( {2\pi n} \right) + \varphi } \right)\\{v_n} = \frac{1}{{0.22}} + \frac{{\sqrt 2 }}{{\sqrt {1 + 2{{\left( {0.22} \right)}^2}} }}\cos \left( {\sqrt 2 \left( {2\pi n} \right) + \varphi } \right)\\{v_n} = c + \sqrt 2 a\cos \left( {\sqrt 2 \left( {2\pi n} \right) + \varphi } \right)\end{array}\)

So, the values are

\(\begin{array}{l}{x_n} = a\sin \left( {\sqrt 2 \left( {2\pi n} \right) + \varphi } \right)\\{v_n} = \frac{1}{{0.22}} + \frac{{\sqrt 2 }}{{\sqrt {1 + 2{{\left( {0.22} \right)}^2}} }}\cos \left( {\sqrt 2 \left( {2\pi n} \right) + \varphi } \right)\end{array}\)

03

(c) Step 3: Verify the given result.

From the part (b) for large n.

\(\begin{array}{c}{x^2}_n = {a^2}{\sin ^2}\left( {\sqrt 2 \left( {2\pi n} \right) + \varphi } \right)\\{\left( {{v_n} - c} \right)^2} = 2{a^2}{\cos ^2}\left( {\sqrt 2 \left( {2\pi n} \right) + \varphi } \right)\\{x^2}_n + \frac{{{{\left( {{v_n} - c} \right)}^2}}}{2} = {a^2}\left[ {{{\sin }^2}\left( {\sqrt 2 \left( {2\pi n} \right) + \varphi } \right) + {{\cos }^2}\left( {\sqrt 2 \left( {2\pi n} \right) + \varphi } \right)} \right]\\{x^2}_n + \frac{{{{\left( {{v_n} - c} \right)}^2}}}{2} = {a^2}\end{array}\)

Therefore, the results are\({x^2}_n + \frac{{{{\left( {{v_n} - c} \right)}^2}}}{2} = {a^2}\) having a center\(\left( {0,c} \right)\).

Unlock Step-by-Step Solutions & Ace Your Exams!

  • Full Textbook Solutions

    Get detailed explanations and key concepts

  • Unlimited Al creation

    Al flashcards, explanations, exams and more...

  • Ads-free access

    To over 500 millions flashcards

  • Money-back guarantee

    We refund you if you fail your exam.

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

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

In Problems 3 – 18, use the elimination method to find a general solution for the given linear system, where differentiation is with respect to t.

x'+y'-x=5,x'+y'+y=1

In Problems 29 and 30, determine the range of values (if any) of the parameter that will ensure all solutions x(t), and y(t) of the given system remain bounded as t+.

dxdt=λx-y,dydt=3x+y

In Problems 19–24, convert the given second-order equation into a first-order system by setting v=y’. Then find all the critical points in the yv-plane. Finally, sketch (by hand or software) the direction fields, and describe the stability of the critical points (i.e., compare with Figure 5.12).

y''(t)+y(t)-y(t)3=0

In Problems 19 – 21, solve the given initial value problem.

d2xdt2=y;x0=3,x'0=1,d2ydt2=x;y0=1,y'0=-1

Sticky Friction. An alternative for the damping friction model F = -by′ discussed in Section 4.1 is the “sticky friction” model. For a mass sliding on a surface as depicted in Figure 5.18, the contact friction is more complicated than simply -by′. We observe, for example, that even if the mass is displaced slightly off the equilibrium location y = 0, it may nonetheless remain stationary due to the fact that the spring force -ky is insufficient to break the static friction’s grip. If the maximum force that the friction can exert is denoted by m, then a feasible model is given by

\({{\bf{F}}_{{\bf{friction}}}}{\bf{ = }}\left\{ \begin{array}{l}{\bf{ky,if}}\left| {{\bf{ky}}} \right|{\bf{ < }}\mu {\bf{andy' = 0}}\\\mu {\bf{sign(y),if}}\left| {{\bf{ky}}} \right| \ge {\bf{0andy' = 0}}\\ - \mu {\bf{sign(y'),ify'}} \ne 0.\end{array} \right.\)

(The function sign (s) is +1 when s 7 0, -1 when s 6 0, and 0 when s = 0.) The motion is governed by the equation (16) \({\bf{m}}\frac{{{{\bf{d}}^{\bf{2}}}{\bf{y}}}}{{{\bf{d}}{{\bf{t}}^{\bf{2}}}}}{\bf{ = - ky + }}{{\bf{F}}_{{\bf{friction}}}}\)Thus, if the mass is at rest, friction balances the spring force if \(\left| {\bf{y}} \right|{\bf{ < }}\frac{\mu }{{\bf{k}}}\)but simply opposes it with intensity\(\mu \)if\(\left| {\bf{y}} \right| \ge \frac{\mu }{{\bf{k}}}\). If the mass is moving, friction opposes the velocity with the same intensity\(\mu \).

  1. Taking m =\(\mu \) = k = 1, convert (16) into the firstorder system y′ = v (17)\({\bf{v' = }}\left\{ \begin{array}{l}{\bf{0,if}}\left| {\bf{y}} \right|{\bf{ < 1andv = 0}}{\bf{.}}\\{\bf{ - y + sign(y),if}}\left| {\bf{y}} \right| \ge {\bf{1andv = 0}}\\{\bf{ - y - sign(v),ifv}} \ne 0\end{array} \right.\) ,
  2. Form the phase plane equation for (17) when v ≠ 0 and solve it to derive the solutions\({{\bf{v}}^{\bf{2}}}{\bf{ + (y \pm 1}}{{\bf{)}}^{\bf{2}}}{\bf{ = c}}\).where the plus sign prevails for v>0 and the minus sign for v<0.
  3. Identify the trajectories in the phase plane as two families of concentric semicircles. What is the centre of the semicircles in the upper half-plane? The lower half-plane?
  4. What are the critical points for (17)?
  5. Sketch the trajectory in the phase plane of the mass released from rest at y = 7.5. At what value for y does the mass come to rest?
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