Ordinary differential equation

From Wikibase
Jump to navigation Jump to search

Languages: English · français · Esperanto

An ordinary differential equation (ODE) is a differential equation whose unknown function depends on a single independent variable. An equation whose unknown depends on several independent variables is a partial differential equation instead.

General form

An ordinary differential equation of order $n$ has the general form

$$F\bigl(x,y,y',y'',\ldots,y^{(n)}\bigr)=0$$

Example:

A mass $m$ attached to a spring of stiffness $k$, displaced from its rest position and released, obeys the second-order equation

$$\frac{d^{2}x}{dt^{2}}+\omega_{0}^{2}x=0$$

where $\omega_{0}=\sqrt{k/m}$ is the angular frequency of the oscillation.[1]

Reduction to a first-order system

For an ordinary differential equation, considering successive derivatives of the unknown function as new unknown functions may reduce the order of the original equation without losing any information:

$$F\bigl(x,y,y',y'',\ldots,y^{(n)}\bigr)=0 \Longrightarrow \begin{cases} y_2=y_1' \\ y_3=y_2' \\ \vdots \\ y_n=y_{n-1}' \\ y_n'=f(x,y_1,y_2,\ldots,y_n) \end{cases}$$

Example:

A pendulum obeys the equation of motion

$$\frac{d^{2}\theta}{dt^{2}}+\frac{g}{L}\sin\theta=0$$

which involves the second derivative.

To reduce the equation to the first degree, introduce the angular velocity $\omega=\frac{d\theta}{dt}$ as a second unknown:

$$\begin{cases} \omega=\frac{d\theta}{dt} \\ \frac{d\omega}{dt}=-\frac{g}{L}\sin\theta \end{cases}$$

From equation 2 it can be inferred that equilibria is reached at $\theta=0$ or $\theta=\pi$

Linear equations and superposition

An equation is linear when its terms in the unknown are of the first degree; collecting them on the left gives

$$L[y]\equiv y^{(n)}+a_{n-1}(x)\,y^{(n-1)}+\cdots+a_{1}(x)\,y'+a_{0}(x)\,y=g(x)$$

where the coefficients $a_{0}$, …, $a_{n-1}$ depend on $x$ only. The operator $L$ is linear, $L[c_{1}y_{1}+c_{2}y_{2}]=c_{1}L[y_{1}]+c_{2}L[y_{2}]$, so solutions combine:

$$ \begin{cases} L[y_{1}]=L[y_{2}]=0 \ \Longrightarrow\ L[c_{1}y_{1}+c_{2}y_{2}]=c_{1}\cdot 0+c_{2}\cdot 0=0 & \text{(superposition)} \\ L[y_{h}]=0,\quad L[y_{p}]=g \ \Longrightarrow\ L[y_{h}+y_{p}]=0+g=g & \text{(homogeneous + particular)} \end{cases} $$

Hence the general solution of a linear equation is $y=y_{h}+y_{p}$: the homogeneous solution $y_{h}$ carries the $n$ arbitrary constants, and $y_{p}$ is any one solution of the full equation. Neither rule survives for a nonlinear equation, such as $y'=y^{2}$:

$$y_{1}'=y_{1}^{2},\quad y_{2}'=y_{2}^{2} \ \Longrightarrow\ (y_{1}+y_{2})'=y_{1}^{2}+y_{2}^{2}\neq (y_{1}+y_{2})^{2}$$

so the sum of two solutions is not a solution.[2]

Solutions and initial conditions

A solution is a function that satisfies the equation identically, as $\cos\omega_{0}t$ satisfies the harmonic oscillator equation:

$$x=\cos\omega_{0}t \ \Longrightarrow\ \ddot{x}=-\omega_{0}^{2}\cos\omega_{0}t=-\omega_{0}^{2}x \ \Longrightarrow\ \ddot{x}+\omega_{0}^{2}x=0$$

Together with $\sin\omega_{0}t$ it gives the general solution of that equation, with one arbitrary constant for each order:

$$x(t)=C_{1}\cos\omega_{0}t+C_{2}\sin\omega_{0}t$$

Fixing the constants requires as many conditions as the order. For a mass pulled to the position $x_{0}$ and released from rest at $t=0$,

$$x(0)=C_{1}=x_{0} \ \Longrightarrow\ \dot{x}(t)=-C_{1}\omega_{0}\sin\omega_{0}t+C_{2}\omega_{0}\cos\omega_{0}t \ \Longrightarrow\ \dot{x}(0)=C_{2}\omega_{0}=0$$

so $C_{2}=0$, and the particular solution is the oscillation

$$x(t)=x_{0}\cos\omega_{0}t$$

Imposed at a single point, the $n$ conditions

$$y(x_{0})=y_{0},\qquad y'(x_{0})=y_{1},\qquad \ldots,\qquad y^{(n-1)}(x_{0})=y_{n-1}$$

form an initial value problem; at the first order the single condition is


Conditions imposed at two different points are boundary conditions, and the problem is a boundary value problem, which may have no solution or infinitely many.[1]

See also

References

Further reading