Try the fastest way to create flashcards
Question

Either by hand or using a computer software package, sketch the phase plane diagrams for the given system. Discuss the stability of each limit cycle and critical point. Here x=rcosθx=r \cos \theta and y=rsinθy=r \sin \theta.
dxdt=y+xsin(1/r)\frac{d x}{d t}=y+x \sin (1 / r),

dydt=x+ysin(1/r)\frac{d y}{d t}=-x+y \sin (1 / r)

Solution

Verified
Answered 6 months ago
Answered 6 months ago
Step 1
1 of 4

Multiplying the first equation by xx, the second one by yy and adding, we get

xdxdt+ydydt=x(y+xsin(1/r))+y(x+ysin(1/r))=xy+x2sin(1/r)xy+y2sin(1/r)Simplify=x2sin(1/r)+y2sin(1/r)Take out common factor=(x2+y2)sin(1/r)Use the fact that x2+y2=r2=r2sin(1/r)\begin{align*} x\frac{dx}{dt}+y\frac{dy}{dt}&=x\big(y+x\sin(1/r)\big)+y\big(-x+y\sin(1/r)\big)\\ &=xy+x^2\sin(1/r)-xy+y^2\sin(1/r)\\ &\qquad\text{Simplify}\\ &=x^2\sin(1/r)+y^2\sin(1/r)\\ &\qquad\text{Take out common factor}\\ &=(x^2+y^2)\sin(1/r)\\ &\qquad\text{Use the fact that $x^2+y^2=r^2$}\\ &=r^2\sin(1/r) \end{align*}

Using the identity

xdxdt+ydydt=rdrdt,x\frac{dx}{dt}+y\frac{dy}{dt}=r\frac{dr}{dt},

we have

rdrdt=r2sin(1/r).r\frac{dr}{dt}=r^2\sin(1/r).

Simplify

drdt=rsin(1/r).\frac{dr}{dt}=r\sin(1/r).

We have that dr/dt=0dr/dt=0 when sin(1/r)=0    1/r=nπ    r=1/(nπ), n=1,2,\sin(1/r)=0\implies 1/r=n\pi\implies r=1/(n\pi),~n=1,2,\ldots. Thus, the origin (r=0r=0) is not an isolated critical point. Note that dr/dt>0dr/dt>0 (sin(1/r)>0\sin(1/r)>0) for 2nπ<1/r<(2n+1)π2n\pi<1/r<(2n+1)\pi or 1/[(2n+1)π]<r<1/(2nπ)1\big/[(2n+1)\pi]<r<1\big/(2n\pi), and dr/dt<0dr/dt<0 for 1/(2nπ)<r<1/[(2n1)π]1\big/(2n\pi)<r<1\big/[(2n-1)\pi].

Consequently, the trajectories spiral into the limit cycles r=1/(2nπ)r=1/(2n\pi) (hence, stable) and away from the limit cycles r=1/[(2n+1)π]r=1\big/[(2n+1)\pi] (hence, unstable). {\bf The origin is a stable critical point}. To realize this, notice that for every ϵ>0\epsilon>0 there exists a positive integer nϵn_\epsilon such that 1/(2nϵπ)<ϵ1/(2n_\epsilon\pi)<\epsilon, so we can simply choose δϵ=1/(2nϵπ)\delta_\epsilon=1/(2n_\epsilon\pi). Then, the fact that r=1/(2nϵπ)r=1/(2n_\epsilon\pi) is a stable limit cycle guarantees that any trajectory that begins within δϵ\delta_\epsilon of the origin remains within ϵ\epsilon of the origin.

Create a free account to view solutions

Create a free account to view solutions

Recommended textbook solutions

Fundamentals of Differential Equations 9th Edition by Arthur David Snider, Edward B. Saff, R. Kent Nagle

Fundamentals of Differential Equations

9th EditionISBN: 9780321977069 (5 more)Arthur David Snider, Edward B. Saff, R. Kent Nagle
2,119 solutions
Fundamentals of Differential Equations and Boundary Value Problems 7th Edition by Arthur David Snider, Edward B. Saff, R. Kent Nagle

Fundamentals of Differential Equations and Boundary Value Problems

7th EditionISBN: 9780321977106Arthur David Snider, Edward B. Saff, R. Kent Nagle
2,485 solutions
A First Course in Differential Equations with Modeling Applications 10th Edition by Dennis G. Zill

A First Course in Differential Equations with Modeling Applications

10th EditionISBN: 9781111827052 (2 more)Dennis G. Zill
2,369 solutions
Differential Equations with Boundary-Value Problems 9th Edition by Dennis G. Zill

Differential Equations with Boundary-Value Problems

9th EditionISBN: 9781305965799 (3 more)Dennis G. Zill
3,184 solutions

More related questions

1/4

1/7