Delayed differential equations – conditions for a steady-state solution

calculusdelay differential equationsordinary differential equations

Is there a way of finding the condition required for a delayed differential equation to have a steady state solution? By 'steady state' I mean that $\frac{dx}{dt}=0$, or alternatively, $x(t)$ approaching some constant as $t \rightarrow\infty$.

If we have an ordinary differential equation:

$\frac{d}{dt}x(t) = ax(t) + c$

then we can solve it using the integrating factor method to get the following solution:

$x(t) = e^{at}x(0) + -\frac{c}{a}[1-e^{at}]$

assuming we have the initial condition $x(0)$. Call this equation (*). We can now see that a steady state solution exists only if $a<0$, so all of the exponential terms become $0$ as $t \rightarrow \infty$.

My question is this: is there an equivalent condition for a steady state solution to exist for a delayed differential equation of the form:

$\frac{d}{dt}x(t) = a x(t) + bx(t-\tau) + c$ ?

My hand-waving solution is that the condition is $(a +b)<0 $ for steady state to exist, since when $t \rightarrow\infty$, there isn't a difference between $t$ and $t-\tau$. This also recovers the correct answer in the case where $\tau \rightarrow 0$. I was hoping someone might be able to confirm or dis-confirm this in a slightly more rigourous way. Thanks!

Best Answer

To have a steady state response we need a stable system. Now calling $u(t) = x(t)-\frac{c}{a+b}$ (assuming $a +b \ne 0$ etc.) we have

$$ \dot{u}(t)=au(t)+bu(t-\tau) $$

and after Laplace transforming

$$ U(s) = \frac{u(0)}{s-a-be^{-\tau s}} $$

so $u(t)$ is stable as far as the poles for $U(s)$ are left complex plane included.

As $s-a-be^{-\tau s} = x-a-b e^{-\tau x}\cos(\tau y) + i\left(y+be^{-\tau x}\sin(\tau y)\right)$ we can construct some graphics to show the poles dependence on $a,b,\tau$

Follows a plot showing the case $a = 1, b = -2, \tau = 1$

enter image description here

and the case with $a = -2, b = 1, \tau = 1$

enter image description here

The poles are at the intersection between the blue and red curves. In the first case we can verify that the response is unstable and in the second case the response is stable. In delayed DE appears infinite poles and we should assure that no one invades the positive complex plane.

It seems that for this differential equation it is sufficient that $a + b <0$, to be stable although this would have to be demonstrated.

Related Question