For the pendulum equation \(\ddot{\theta}+(g / l) \sin \theta=0\) find the equations of the trajectories in the phase plane and sketch the phase portrait of the system. Show that on the separatrices $$ \theta(t)=(2 n+1) \pi \pm 4 \arctan [\exp (\omega t+\alpha)] $$ where \(n\) is an integer, \(\alpha=\) constant and \(\omega=\sqrt{g / l}\). [Hint: Use the substitution \(\left.u=\tan \frac{1}{4}\\{\theta-(2 n+1) \pi\\} .\right]\)

Short Answer

Expert verified
Based on the given pendulum equation, we have transformed it into a first-order system of equations using the suggested substitution and obtained the trajectories in the phase plane. After analyzing the behavior of the system, we have sketched the phase portrait by identifying the fixed points and understanding the direction of the trajectories. Finally, we have shown the given equation for the separatrices using the results obtained and analyzing the behavior of the system at the separatrices' trajectories.

Step by step solution

01

Apply the suggested substitution

Using the Hint given in the problem, let's apply the suggested substitution: $$ u = \tan\left(\frac{1}{4}\left(\theta - (2n + 1)\pi\right)\right) $$ Next, we will need the derivatives of \(\theta\) to transform the given equation into a first-order equation: $$ \dot{\theta} = \frac{4\dot{u}}{1 + u^2} $$ and $$ \ddot{\theta} = \frac{4}{(1 + u^2)^2}(1 - u^2)\ddot{u}-\frac{8u}{(1 + u^2)^2}\dot{u}^2 $$ Substitute the derivatives of \(\theta\) into the given equation: $$ \frac{4}{(1 + u^2)^2}(1 - u^2)\ddot{u}-\frac{8u}{(1 + u^2)^2}\dot{u}^2+(g / l) \sin\left(4\arctan{u}+(2 n+1)\pi\right)=0 $$
02

Find the equation of trajectories in the phase plane

Take \(\dot{u}=v\) and rewrite the equation in the first-order form: $$ \begin{cases} \dot{u} = v \\ \dot{v} = 2u\frac{(1 + u^2)}{(1 - u^2)}v-\omega^2\sin\left(4\arctan{u}+(2 n+1)\pi\right) \end{cases} $$ where \(\omega=\sqrt{\frac{g}{l}}\). Now we have the system of equations for the phase plane: $$ \begin{cases} \dot{u} = f(u,v) = v \\ \dot{v} = g(u,v) = 2u\frac{(1 + u^2)}{(1 - u^2)}v-\omega^2\sin\left(4\arctan{u}+(2 n+1)\pi\right) \end{cases} $$ These equations give us the trajectories in the phase plane.
03

Sketch the phase portrait of the system

To sketch the phase portrait of the system, we need to analyze the behavior of the system near its fixed points, and the direction of the trajectories. Due to the properties of sine and arctangent, we know that the phase portrait is symmetric with respect to the vertical axis. Also, the pendulum equation is periodical. Thus, we only need to analyze one area of the phase plane. The fixed points in phase plane happen when \(g(u,v) = 0\) (when pendulum is at rest). They are: 1. \((u,v) = (0,0)\) (unstable focus) 2. \((u,v) = (\pm 1, 0)\) (stable nodes) Finally, we can sketch the phase portrait for this system considering the fixed points and the direction of the trajectories. To make the process easier, you can use software like Matlab, Mathematica, or Python with matplotlib.
04

Show the equation for the separatrices

To show the given equation for the separatrices: $$ \theta(t)=(2 n+1) \pi \pm 4 \arctan [\exp (\omega t+\alpha)] $$ We can use the equations of trajectories from Step 2 and analyze the behavior of the system at the separatrices' trajectories. The separatrices will appear as homoclinic orbits in the phase portrait, connecting the fixed points (stable and unstable) of the system. Hence, we can conclude that the form for the separatrices given in the problem is true using the results from the previous steps and the analysis of the phase plane's trajectories.

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

Show that Rayleigh's equation in the form (13.28) has a single critical point at \((0,0)\) and that this is always unstable. Making use of substitution \(\dot{x}=v / \sqrt{3}\) show that \(v\) satisfies the Van der Pol equation (13.29).

For the Arms-Race model system (13.26) with all parameters positive show that there is an asymptotically stable coexistence or a runaway escalation according as \(c_{1} c_{2}>a_{1} a_{2}\) or \(c_{1} c_{2}

The simple SIR model equations for the transmission of a disease are $$ \begin{aligned} \dot{S} &=-a S I \\ \dot{I} &=a S I-b I \\ \dot{R} &=b I \end{aligned} $$ where \(S(t), I(t), R(t)\) are respectively susceptibles, infectives, removed/recovered and \(a, b\) are positive constants. (a) Show that the overall population \(N=S+I+R\) remains constant, so that we may consider \((S, I)\) in a projected phase plane. Hence show that a trajectory with initial values \(\left(S_{0}, I_{0}\right)\) has equation \(I(S)=\) \(I_{0}+S_{0}-S+(b / a) \ln \left(S / S_{0}\right)\) (b) Using the function \(I(S)\) show that an epidemic can occur only if the number of susceptibles \(S_{0}\) in the population exceeds the threshold level \(b / a\) and that the disease stops spreading through lack of infectives rather than through lack of susceptibles. (c) For the trajectory which corresponds to \(S_{0}=(b / a)+\delta, I_{0}=\epsilon\) with \(\delta, \epsilon\) small and positive, show that, to a good approximation, there are \((b / a)-\delta\) susceptibles who escape infection [the KermackMcKendrick theorem of epidemiology \((1926 / 27)]\).

For the Rikitake dynamo system (13.35): (a) Show that there are two real critical points at \(\left(\pm k, \pm 1 / k, \mu k^{2}\right)\) in the \(\left(X_{1}, X_{2}, Y\right)\) phase space, where \(k\) is given by \(A=\mu\left(k^{2}-1 / k^{2}\right)\). (b) Show that the stability of these critical points is determined by eigenvalues \(\lambda\) satisfying the cubic $$ (\lambda+2 \mu)\left[\lambda^{2}+\left(k^{2}+\frac{1}{k^{2}}\right)\right]=0 $$ so that the points are not asymptotically stable in this approximation. (For the full system they are actually unstable.) (c) Show that the divergence of the phase-space flow velocity is negative, so that the flow causes volume to contract. (d) Given that \(\bar{Y}=\sqrt{2}(Y-A / 2)\), show that $$ \frac{1}{2} \frac{\mathrm{d}}{\mathrm{d} t}\left(X_{1}^{2}+X_{2}^{2}+\bar{Y}^{2}\right)=-\mu\left(X_{1}^{2}+X_{2}^{2}\right)+\sqrt{2} \bar{Y} $$ and use this result to determine in which region of the space the trajectories all have a positive inward component towards \(X_{1}=\) \(X_{2}=\bar{Y}=0\) on surfaces \(X_{1}^{2}+X_{2}^{2}+\bar{Y}^{2}=\) constant

Consider the gradient system (13.17) in the case \(U(x, y)=x^{2}(x-1)^{2}+y^{2}\). Find the critical points and their character. Sketch the phase portrait for the system.

See all solutions

Recommended explanations on Physics 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