Determine the stability of a fixed point of oscillations

Click For Summary
The discussion focuses on determining the stability of fixed points in a system of coupled differential equations representing chemical reactions. The user successfully observed oscillatory behavior through numerical solutions using Euler's Method but seeks to analyze the stability of the fixed points analytically. The key approach involves using Taylor expansion to evaluate the derivative f'(x*) and its implications for stability based on whether |f'(x*)| is greater or less than 1. Suggestions include examining the steady state and analyzing the signs of the derivatives near this point, as well as sketching isoclines for better visualization. The conversation highlights the challenge of transitioning from numerical to analytical methods in this context.
peadar2211
Messages
1
Reaction score
0

Homework Statement


I have a system of coupled differential equations representing chemical reactions and given certain initial conditions for the equations I can observe oscillation behaviour when I solved the equations numerically using Euler's Method. However, then it asks to investigate the stability of the fixed point of each equation analytically and compare them to the numerically determined values. The exact wording is:
The steady state solution for the rate equations may be termed a fixed point(x*) since xn+1=f(xn) =xn or x*=f(x*) where f(x) represents the left side of the expression for the Euler algorithm.
To determine the stability of a fixed point we let xn=x*+en and xn+1=x*+en+1 where |ei|<<1. Using a Taylor expansion this can be written as x*+enf'(x*). Substituting for xn+1 we get that the derivative f'(x*)=en+1/en. If |f'(x*)|>1 the trajectory will diverge from x* since en+1>en. The opposite is true for |f'(x*)|<1 and the fixed point is stable. Investigate the stability of the fixed point analytically and compare with your numerical results.

Homework Equations


Euler's Method: xn+1=xn+Δt*dxn/dt
The differential equations I solved, although I don't think they're relevant to this part are
X'=A-(B+1)X+X2Y and Y'=BX-X2Y

The Attempt at a Solution


I know that I need to investigate f'(x*) so I've tried differentiating Euler's method with respect to xn, but then I have 1+(something) where the something is either really complicated because it would be d/dx(dx/dt)Δt, or it would be 0. Either way, it doesn't seem like something I could use to either graph or solve analytically to see when it has a value above 1 (As it is constant in one of those cases).
 
Physics news on Phys.org
As far as I can make out you have a problem with just two variables, X and Y, and are quoting a text elaborating the formulation for n variables, xn.

IMO you would do better to just considered two for now.

You have been able to solve the differential equations, exceptional situation, lucky you! However even without solving them you can quite easily obtain what the steady state is. Or are, because in a non-linear case there may be more than one. In this case?

When you know a stationary state the next thing Is to see how the differential equation looks like just near it - in first place what are signs of X' and Y' near the s.p.? Helpful to sketch the isoclines (X' = 0 and Y' = 0) and put in little arrows showing which way sample points in each sector are moving.
 
Last edited:
When student posts a problem in his first and only post and then never comes back even to the site in over a week we can count the problem as solved? - what is the practice here?
 
Question: A clock's minute hand has length 4 and its hour hand has length 3. What is the distance between the tips at the moment when it is increasing most rapidly?(Putnam Exam Question) Answer: Making assumption that both the hands moves at constant angular velocities, the answer is ## \sqrt{7} .## But don't you think this assumption is somewhat doubtful and wrong?

Similar threads

  • · Replies 2 ·
Replies
2
Views
2K
  • · Replies 12 ·
Replies
12
Views
2K
  • · Replies 2 ·
Replies
2
Views
1K
Replies
10
Views
3K
  • · Replies 18 ·
Replies
18
Views
3K
Replies
3
Views
2K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 7 ·
Replies
7
Views
2K
  • · Replies 2 ·
Replies
2
Views
2K
Replies
2
Views
1K