How can I find linear functions that will tend to equilibrium in a given system?

  • Context: Graduate 
  • Thread starter Thread starter Mynock
  • Start date Start date
  • Tags Tags
    Linear Systems
Click For Summary
SUMMARY

The discussion focuses on finding linear functions f(r) and h(r) for the system defined by the equations dp/dt = A*f(r) and dr/dt = -B*h(r), where A and B are positive constants. The solution involves using linear functions of the form f(r) = α·r and h(r) = β·r, leading to the stability condition β > 0. The derived solutions indicate that the system will reach equilibrium at p = γ and r = 0, provided that the condition for stability is met.

PREREQUISITES
  • Understanding of linear differential equations
  • Familiarity with stability analysis in dynamical systems
  • Knowledge of exponential decay functions
  • Basic calculus, particularly differentiation and integration
NEXT STEPS
  • Study linear differential equations and their solutions
  • Research stability criteria in dynamical systems
  • Explore the implications of equilibrium points in mathematical models
  • Learn about the role of constants in differential equations
USEFUL FOR

Mathematicians, physicists, and engineers interested in dynamical systems, particularly those analyzing stability and equilibrium in linear models.

Mynock
Messages
1
Reaction score
0
Hey, I'm trying to find two linear functions, f(r) and h(r) for the following system:

dp/dt = A*f(r)

dr/dt = -B*h(r)

where A and B are constants greater than zero. I'm trying to find linear functions that will tend to equillibrium, and also where

limit df(r)/dt = 0
r->0

I have been trying various linear functions and have been unable to come up with a solution. Is there a solution? and if so, what would be one? It's probalby something simple that I'm overlooking. Any help would be appreciated. Thanks.
 
Physics news on Phys.org
The linear (homogeneous) functions in one variable, 'r' in this case, have the general form:

[tex]f(r) \ = \ \alpha\cdot r,\quad g(r) \ = \ \beta\cdot r, \qquad \alpha,\beta\in\mathbb{R}[/tex]

So this system can be trivially solved, since the variables are separated. On one hand, you have:

[tex]\dot{r}\ = \ -B\beta\, r[/tex]

so

[tex]r(t) \ = \ r_0e^{-B\beta t}[/tex]

The condition for stability is obviously [tex]\beta > 0[/tex].

On the other hand:

[tex]\dot{p}\quad = \quad A\alpha\, r \quad = \quad A\alpha r_0 e^{-B\beta t}[/tex]

Then,

[tex]p(t) \quad = \quad \gamma \ - \ \frac{A\alpha r_0}{B\beta} e^{-B\beta t}[/tex]

where [tex]\gamma[/tex] is an integration constant. This means that, as long as [tex]\beta > 0[/tex], your system will be stable, with equilibrium point [tex]p = \gamma, r = 0[/tex].
 

Similar threads

  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 2 ·
Replies
2
Views
1K
  • · Replies 2 ·
Replies
2
Views
3K
  • · Replies 4 ·
Replies
4
Views
3K
  • · Replies 0 ·
Replies
0
Views
4K
  • · Replies 3 ·
Replies
3
Views
3K
  • · Replies 4 ·
Replies
4
Views
2K
  • · Replies 1 ·
Replies
1
Views
3K
  • · Replies 1 ·
Replies
1
Views
3K
  • · Replies 8 ·
Replies
8
Views
3K