Differential equation, limit cycles, stability, phase plane. By 11/12

In summary: For part (c), you need to determine the region of initial positions (x(0), y(0)) such that the trajectory eventually ends up at (0,0). This can be done by solving the system of equations for x and y as functions of t and then finding the values of (x(0), y(0)) such that x(t) = 0 and y(t) = 0 for some t > 0. The resulting region would be a subset of the (x,y) plane.
  • #1
debian
7
0

Homework Statement



The trajectory of an arrow in space obeys the following system of equations:

[itex]\dot{x}[/itex] = y+[itex](x^2+y^2-3)^2 (x^3-x+xy^2)[/itex]

[itex]\dot{y}[/itex] = y+[itex](x^2+y^2-3)^2 (y^3-y+x^2y)[/itex]

1. Questions

a) Derive an ODE for the radial coordiante r(t) = [itex]\sqrt[]{x^2(t)+y^2(t)}[/itex]

b) Show that the system admits two limit cycles and classify their stability

c) The target is standing at the origin (0,0) of the phase plane. Let Po=(x(0),y(0))[itex]\neq[/itex] (0,0) denote the initial position of the arrow at the time t = 0. Determine and plot the region in the (x,y) plane containing all the points Po such that the arrow will hit the target at some instant of time t > 0


The Attempt at a Solution



I derived ODE in radial coordiante and got:

[itex]\dot{r}[/itex]=[itex](r^2-3)^2(r^2-1)r[/itex]

So we have 5 points to check the stability: [itex]0, 1,-1,\sqrt{3},-\sqrt{3} [/itex]

Than I calculate the derivative of [itex]\dot{r}[/itex] and evaluate in those 5 points: I get five zeros so they are undecided.

Then I daw the plot of [itex]\dot{r}[/itex] and if the [itex]\dot{r}[/itex]>0 it does mean that arrow on the phase line is to the right and if [itex]\dot{r}[/itex]<0 the arrow is to the left.

So The phase line looks like this: <- (-\sqrt{3}) <- (-1) -> (0) <- (1) -> (\sqrt{3}) ->

Is that correct? What should I do next?
 
Physics news on Phys.org
  • #2
debian said:

Homework Statement



The trajectory of an arrow in space obeys the following system of equations:

[itex]\dot{x}[/itex] = y+[itex](x^2+y^2-3)^2 (x^3-x+xy^2)[/itex]

[itex]\dot{y}[/itex] = y+[itex](x^2+y^2-3)^2 (y^3-y+x^2y)[/itex]

Are you sure these are correct? Because they tell me that
[tex]r\dot r = x\dot x + y\dot y[/tex]
will not be independent of [itex]\theta[/itex], because the linear terms give me
[tex]xy + y^2 = r^2(\cos\theta\sin\theta + \sin^2\theta)[/tex].

Assuming that you are correct that
[tex]\dot r = r(r^2 - 3)^2(r^2 - 1) = f(r)[/tex]
then since by convention [itex]r \geq 0[/itex], there is a fixed point at the origin and two orbits, at [itex]r = \sqrt 3[/itex] and [itex]r = 1[/itex].

To determine the stability, you look at the sign of [itex]f'(r)[/itex]. If it's negative then the point is stable, and if it's positive then the point is unstable. If it's zero then it's inconclusive and you have to look at higher derivatives.

Alternatively, in one dimension r is increasing if [itex]\dot r > 0[/itex] and decreasing if [itex]\dot r < 0[/itex]. By inspection, [itex]\dot r < 0[/itex] if [itex]0 < r < 1[/itex] and [itex]\dot r > 0[/itex] when [itex]1 < r < \sqrt 3[/itex] and [itex]r > \sqrt 3[/itex]. So the origin is stable, the orbit at 1 is unstable, and the orbit at [itex]\sqrt 3[/itex] is stable from below and unstable from above.

For part (c) you need to think about what trajectories eventually end up at the origin.
 
  • #3
There is a typo: shoul be -x instead of one of y

[itex]\dot{x}[/itex] = y+[itex](x^2+y^2-3)^2 (x^3-x+xy^2)[/itex]

[itex]\dot{y}[/itex] = -x+[itex](x^2+y^2-3)^2 (y^3-y+x^2y)[/itex]
 
  • #4
debian said:
There is a typo: shoul be -x instead of one of y

[itex]\dot{x}[/itex] = y+[itex](x^2+y^2-3)^2 (x^3-x+xy^2)[/itex]

[itex]\dot{y}[/itex] = -x+[itex](x^2+y^2-3)^2 (y^3-y+x^2y)[/itex]

Then your result for [itex]\dot r[/itex] is correct.
 

1. What is a differential equation?

A differential equation is a mathematical equation that describes the relationship between a function and its derivatives. It is often used to model physical systems and predict their behavior over time.

2. What are limit cycles in differential equations?

Limit cycles are closed curves in the phase plane of a differential equation that represent stable oscillatory behavior. This means that the system will continue to cycle between certain values, rather than approaching a fixed point.

3. How do you determine the stability of a system in a phase plane?

The stability of a system in a phase plane can be determined by analyzing the direction and magnitude of the vectors in the phase plane. If the vectors point towards a fixed point, the system is stable. If they point away from the fixed point, the system is unstable.

4. What is the significance of phase planes in differential equations?

Phase planes provide a visual representation of the behavior of a system over time. They allow us to see how the variables of a system change and interact with each other, and can help us determine the stability and behavior of the system.

5. How are differential equations used in real-world applications?

Differential equations are used in a wide range of fields, including physics, engineering, economics, and biology. They are used to model and understand complex systems and phenomena, and can be used to make predictions and solve problems in these fields.

Similar threads

  • Calculus and Beyond Homework Help
Replies
1
Views
284
Replies
7
Views
521
  • Calculus and Beyond Homework Help
Replies
7
Views
284
  • Calculus and Beyond Homework Help
Replies
21
Views
840
  • Calculus and Beyond Homework Help
Replies
8
Views
802
  • Calculus and Beyond Homework Help
Replies
0
Views
166
  • Calculus and Beyond Homework Help
Replies
2
Views
733
  • Calculus and Beyond Homework Help
Replies
2
Views
513
  • Calculus and Beyond Homework Help
Replies
3
Views
571
  • Calculus and Beyond Homework Help
Replies
6
Views
760
Back
Top