Ordinary Differential Equations/Maximum domain of solution

From Wikibooks, open books for an open world
Jump to: navigation, search

Even if a differential equation satisfies the Picard–Lindelöf theorem, it may still not have a solution for all of \mathbb{R} that is uniquely determined by a single initial condition. Here we study the maximum interval to which a solution may be extended.

Uniqueness of solution over intervals[edit]

Theorem Local uniqueness implies global uniqueness over intervals

Hypothesis

  • y is a solution to an IVP
  • y is locally unique (by the Picard–Lindelöf theorem for example)
  • the domain of y is an interval (which contains x_0, otherwise the initial condition makes no sense)

y(x) is the only solution to that IVP of the which has that domain.

Maximal solution[edit]

Definition[edit]

Definition Maximal solution and maximal domain of solution

The maximal solution y_{max} of an IVP which is locally unique is the solution such that

  • has an interval as domain (which contains x_0, otherwise the initial condition makes no sense)
  • is not a restriction of any other solution whose domain is a larger interval

Dom(y_{max}) is called the maximal domain of solution and is denoted ]x_-,x_+[


Unicity[edit]

Theorem Every IVP has an unique maximal solution

Theorem If there is a solution y_\mathbb{R} with domain \mathbb{R} then it is the maximal solution.

We must verify that:

  • Dom(y_\mathbb{R}) is an interval. Obviously true.
  • y_{\mathbb{R}} is not a restriction of any other solution. Obvious because the domain of y_{\mathbb{R}} is already all of \mathbb{R}

Behavior at boundary[edit]

Theorem Behavior at the boundary of a maximal solution

Hypothesis

The domain of a maximal solution is smaller than infinity (x_\pm \neq \pm \infty)

Atx_\pm one or two of the following happen:

  • explosion in finite time: \lim_{x \to x_\pm} ||y(x)|| = \infty
  • y falls out of the domain of F: \lim_{x \to x_\pm} \in \overline{Dom(F)}

Sufficient condition for Dom(y)=\mathbb{R}[edit]

Theorem Growth at most linear implies Dom(y)=\mathbb{R}

Hypothesis

F grows at most linearly with y

 ||F(x,y)|| \leq a||y||^2 + b

The domain of the maximal solution is all of \mathbb{R}

Examples[edit]

y'=y,\, y(0)=1[edit]

Example y'=y,\, y(0)=1 has an infinity of solutions

F(x,y)=y , which is C^1 in both x and y, and therefore satisfies the Picard–Lindelöf theorem on all of its domain, so that solutions are locally unique.

All of the following are solutions to this IVP:

 \begin{align}y_1 : \, & ]-2,2[ \, \to \, ]e^{-2},e^{2}[ \\ & x \mapsto e^x \end{align}

 \begin{align}y_2 : \, & ]-1,3[ \, \to \, ]e^{-1},e^{3}[ \\ & x \mapsto e^x \end{align}

 \begin{align}y_\mathbb{R} : \, & \mathbb{R} , \to \, \mathbb{R}_+ \\& x \mapsto e^x \end{align}

 \begin{align}y_3 : & Dom(y_3) = ]-1,1[ \, \cup \, ]2,4[ \\ & x \mapsto \begin{cases}e^x & x \in ]-1,1[ \\y(3)e^{(x-3)} & ]2,4[\end{cases}\end{align} Different domains mean completely different solutions

A function is a set of ordered pairs {(x,f(x))}, with x \in Dom(f). If the domains are different, the sets are completely different, and therefore the solutions are completely different too.

y_1,y_2 and y_\mathbb{R} are the only solutions on their respective domains

This is a conclusion of the theorem on the uniqueness on intervals, because their domains are all intervals.

If the domain is not an interval there is generally no uniqueness

For the fixed domain of y_3

 ]-1,1[ \, \cup \, ]2,4[

which is clearly not an interval, then for any value y(3) we choose we have a different solution.

Therefore, even if the domain is fixed but not an interval and there is local uniqueness, there may not be global uniqueness. This is the major reason why we restrict ourselves to a maximal domain that is an interval, even if there may be larger domains with non-unique solutions.

To determine the solution uniquely the initial value at 0 is not enough. We would have to set a value on ]2,4[ such as y(3). This happens because the initial condition at 0 is separated from the interval ]2,4[. We need to set the a value inside of ]2,4[, such as y(3) in order to fully determine the solution.

y_\mathbb{R} is the maximal solution y_{max}

Since the domain is all of \mathbb{R}, it must be the maximal solution by the proposition above.

We could immediately see that Dom(y_{max})=\mathbb{R} without solving it because the growth of F is linear, as stated on the theorem of sublinear growth.

We note also that as stated on the definition, any solution that has a domain that is an interval such as y_1 and y_2 are restrictions of y_\mathbb{R}.

Solutions which are not defined on intervals may not be a restriction of the the maximal solution y_{max}

y_3 is not a a restriction of y_\mathbb{R} unless y(3) \neq e^3. However, y_3 is not unique, and in unacceptable in a physical situation since there is no causality between the two intervals of its domain. So it is not very serious if it is not a restriction of y_\mathbb{R}, which is generally the 'best' one.


y' = y^2[edit]

Example y' = y^2

F(x,y)=y^2, which is C^1 \, \forall (x,y) \in \mathbb{R}^2 and therefore Lipschitz continuous, satisfying the Picard–Lindelöf theorem.

  • the maximal domain may depend on the initial condition
  • the maximal domain may be different than \mathbb{R}
  • at the borders of the maximal domain the function may go to infinity

The maximal solutions are

 \begin{cases} \mathbb{R} & y_0 = 0 \\]x_0-\frac{1}{y_0},+\infty[ & y_0 > 0 \\]-\infty,x_0-\frac{1}{y_0}[ & y_0 < 0 \end{cases}

which clearly depend on the initial condition, and are not all of \mathbb{R} unless for the trivial solution.

At x_0-\frac{1}{y_0} the maximal solution tends to infinity.

Extending the domain other side of the singularity leads to non uniqueness

Fix y(1)=1. One might want to extend the domain of y_{max} to \mathbb{R}^*. But then there is no more uniqueness since for any a>0 the following is a solution

 \begin{align}y_a : \, & Dom(y_1) = \mathbb{R}^* \\ & x \mapsto \begin{cases}\frac{1}{x} & x>0 \\\frac{1}{x-a} & x<0\end{cases} \end{align}

This is no surprise since \mathbb{R}^* is not an interval.

y_a above are not maximal solutions

Actually any solution that is an interval is a restriction of those solutions. But this is not a maximal solution since its domain is not an interval.

We wouldn't like those to be maximal solutions since they would not be unique.


y' = \frac{-xy}{ln(y)}[edit]

Example y' = \frac{-xy}{ln(y)},\, y(0) = e^2

F(x,y)=\frac{-xy}{ln(y)}, which is C^1 and therefore Lipschitz continuous \forall y > 0, satisfying the Picard–Lindelöf theorem. If y=0, f is undefined, so that

 \Omega=\mathbb{R} \times \mathbb{R}^*

The maximum solution is

 \begin{align}y : & \, Dom(y) = \, ]-2,2[ \\ & x \mapsto \exp(\sqrt{4-x^2}) \end{align} At the border the solution may fall out of the domain of F

The solution is only defined for x \in ]-2,2[, because if x=\pm 2 then y=0, which is not in the domain of F. The maximal solution ends at those points, which is one of the two possibilities for finite domains. The solution does not explode in this case.