Estimating Distance Between Solutions to ODE

  • Thread starter Thread starter AKG
  • Start date Start date
  • Tags Tags
    Ode
AI Thread Summary
The discussion centers on proving that for the ODE x'(t) = cos(t) - 2x(t) - x(t)^{2005}, the distance between two solutions x_a and x_b with initial values a and b satisfies the inequality |x_a(t) - x_b(t)| < |a - b| exp(-2t) for t > 0. It is established that the solutions do not cross, allowing the removal of absolute values under the assumption that a > b. The Mean Value Theorem is applied to derive a relationship involving the flow associated with the ODE, leading to the need to show that the derivative of the flow is less than exp(-2t). The discussion highlights challenges in demonstrating the positivity of certain integrals and the need for a clearer approach to the problem, particularly in estimating the behavior of the solutions over time. The thread concludes with a request for guidance on correcting the approach to prove the desired inequality.
AKG
Science Advisor
Homework Helper
Messages
2,559
Reaction score
4
Consider the ODE:

x&#039;(t) = f(t, x(t)) = \cos (t) - 2x(t) - x(t)^{2005}

For any two solutions x_a,\ x_b with initial values x_a(0) = a,\ x_b(0) = b, prove that for any t > 0, the following inequality holds:

|x_a(t) - x_b(t)| &lt; |a - b|\exp (-2t)

Since the solutions can't cross each other, we know that if a > b, then x_a(t) &gt; x_b(t) for all t, so we can stipulate that a > b and remove the absolute value bars. We have:

\frac{x_a(t) - x_b(t)}{a - b} &lt; \exp (-2t)

Now, consider the time-t mapping \phi _t defined by:

\phi _t (x_0) = x_{x_0}(t) = \phi (t, x_0)

where \phi is the flow associated with the ODE. The above inequality can be re-written:

\frac{\phi _t (a) - \phi _t (b)}{a - b} &lt; \exp (-2t)

By Mean Value Theorem:

\frac{\phi _t (a) - \phi _t (b)}{a - b} \geq \exp (-2t) \Rightarrow (\exists c \in (a, b))(\phi _t &#039;(c) \geq \exp (-2t))

(\forall c \in (a,b))(\phi _t &#039;(c) &lt; \exp(-2t)) \Rightarrow \frac{\phi _t (a) - \phi _t (b)}{a - b} &lt; \exp (-2t)

Since this must hold for arbitrary a, b, we it suffices to show that:

(\forall c)(\phi _t&#039;(c) &lt; \exp(-2t))*

\phi _t&#039;(c) = \frac{\partial \phi}{\partial \x_0}(t, c)

My book gives that:

\frac{\partial \phi}{\partial x_0}(t, c) = \exp \left (\int _0 ^t \frac{\partial f}{\partial x_0}(s, \phi (s, c))\, ds\right )

I want to show that this is less than exp(-2t), which means:

\int _0 ^t \frac{\partial f}{\partial x_0}(s, \phi (s, c))\, ds &lt; -2t

\int _0 ^t \frac{\partial}{\partial x_0}\left (\cos (s) - 2\phi (s, c) - \phi (s, c)^{2005}\right )\, ds &lt; -2t

\int _0 ^t -2\frac{\partial \phi}{\partial x_0} (s, c) - 2005\phi (s, c)^{2004}\frac{\partial \phi}{\partial x_0}(s, c)\, ds &lt; \int _0 ^t -2\, ds

\int _0 ^t \frac{\partial \phi}{\partial x_0}(s, c)\left (1 + \frac{2005}{2}\phi (s, c)^{2004}\right ) - 1\, ds &gt; 0

It would be enough to show that the integrand is positive:

\frac{\partial \phi}{\partial x_0}(s, c)\left (1 + \frac{2005}{2}\phi (s, c)^{2004}\right ) - 1 &gt; 0*

\frac{1}{\frac{\partial \phi}{\partial x_0}(s, c)} &lt; 1 + \frac{2005}{2}\phi (s, c)^{2004}

And it would suffice to show that:

\frac{1}{\frac{\partial \phi}{\partial x_0}(s, c)} &lt; 1*

\frac{\partial \phi}{\partial x_0}(s, c) &gt; 1

\exp \left (\int _0 ^s \frac{\partial f}{\partial x_0}(r, \phi (r, c))\, dr\right ) &gt; \exp (0)

\int _0 ^s \frac{\partial f}{\partial x_0}(r, \phi (r, c))\, dr &gt; 0

\int _0 ^s -2\frac{\partial \phi}{\partial x_0} (r, c) - 2005\phi (r, c)^{2004}\frac{\partial \phi}{\partial x_0}(r, c)\, dr &gt; 0

\int _0 ^s \frac{\partial \phi}{\partial x_0} (r, c)\left (1 + \frac{2005}{2}\phi (r, c)^{2004}\right )\, dr &lt; 0

But 1 is positive, anything to the exponent 2004 is positive, and the parital derivative of \phi is positive, otherwise two solutions would cross. So I cannot show this last inequality, but the steps indicated by * were not necessary step, just sufficient, so those steps are where I went wrong. But how do I correct it? I can't seem to figure the problem out. My guess is that the whole approach above is wrong, so how do I do this problem, especially in general? Thanks.
 
Last edited:
Physics news on Phys.org
AKG said:
By Mean Value Theorem:

\frac{\phi _t (a) - \phi _t (b)}{a - b} \geq \exp (-2t) \Rightarrow (\exists c \in (a, b))(\phi _t &#039;(c) \geq \exp (-2t))

(([/color] ~(\exists c \in (a, b))(\phi _t &#039;(c) \geq \exp (-2t)) \Rightarrow \frac{\phi _t (a) - \phi _t (b)}{a - b} &lt; \exp (-2t) ))[/color]

(\forall c \in (a,b))(\phi _t &#039;(c) &lt; \exp(-2t)) \Rightarrow \frac{\phi _t (a) - \phi _t (b)}{a - b} &lt; \exp (-2t)

The third line is the contrapositive of the first.

However, the middle line (in big red quotation marks) is unnecessary and incorrect anyways.
 
The problem with the second line was that the "~" symbol didn't show up. Anyways, any ideas as to how to do the prolbem?
 
I'm thinking about it. :smile:
 
Last edited by a moderator:
When I asked someone else this, they told me that, by mean value theorem:

\exists !u\ \mbox{such that} \ |x_1(a, t) - x_2(b, t)| = \frac{dx}{da} (u(t), t) |a - b|

Although I didn't find out what x_1,\ x_2,\ x were. They said that, using this, estimate u to get the answer. This doesn't make any sense to me, but maybe this gives you an idea as to how to do this?
 
TL;DR Summary: I came across this question from a Sri Lankan A-level textbook. Question - An ice cube with a length of 10 cm is immersed in water at 0 °C. An observer observes the ice cube from the water, and it seems to be 7.75 cm long. If the refractive index of water is 4/3, find the height of the ice cube immersed in the water. I could not understand how the apparent height of the ice cube in the water depends on the height of the ice cube immersed in the water. Does anyone have an...
Thread 'Variable mass system : water sprayed into a moving container'
Starting with the mass considerations #m(t)# is mass of water #M_{c}# mass of container and #M(t)# mass of total system $$M(t) = M_{C} + m(t)$$ $$\Rightarrow \frac{dM(t)}{dt} = \frac{dm(t)}{dt}$$ $$P_i = Mv + u \, dm$$ $$P_f = (M + dm)(v + dv)$$ $$\Delta P = M \, dv + (v - u) \, dm$$ $$F = \frac{dP}{dt} = M \frac{dv}{dt} + (v - u) \frac{dm}{dt}$$ $$F = u \frac{dm}{dt} = \rho A u^2$$ from conservation of momentum , the cannon recoils with the same force which it applies. $$\quad \frac{dm}{dt}...
Back
Top