Characteristic equation (calculus)

Last updated November 03, 2023

In mathematics, the characteristic equation (or auxiliary equation^[1]) is an algebraic equation of degree $n$ upon which depends the solution of a given $n$ th- order differential equation ^[2] or difference equation.^[3]^[4] The characteristic equation can only be formed when the differential or difference equation is linear and homogeneous, and has constant coefficients.^[1] Such a differential equation, with $y$ as the dependent variable, superscript $(n)$ denoting n^th-derivative, and $a n, a n - 1, ..., a 1, a 0$ as constants,

The characteristic roots (roots of the characteristic equation) also provide qualitative information about the behavior of the variable whose evolution is described by the dynamic equation. For a differential equation parameterized on time, the variable's evolution is stable if and only if the real part of each root is negative. For difference equations, there is stability if and only if the modulus of each root is less than 1. For both types of equation, persistent fluctuations occur if there is at least one pair of complex roots.

The method of integrating linear ordinary differential equations with constant coefficients was discovered by Leonhard Euler, who found that the solutions depended on an algebraic 'characteristic' equation.^[2] The qualities of the Euler's characteristic equation were later considered in greater detail by French mathematicians Augustin-Louis Cauchy and Gaspard Monge.^[2]^[6]

Derivation

Starting with a linear homogeneous differential equation with constant coefficients $a n, a n - 1, ..., a 1, a 0$ ,

a_{n}y^{(n)}+a_{n-1}y^{(n-1)}+\cdots +a_{1}y^{\prime }+a_{0}y=0,

it can be seen that if $y (x) = e rx$ , each term would be a constant multiple of $e rx$ . This results from the fact that the derivative of the exponential function $e rx$ is a multiple of itself. Therefore, $y' = re rx$ , $y ″ = r 2 e rx$ , and $y (n) = r n e rx$ are all multiples. This suggests that certain values of $r$ will allow multiples of $e rx$ to sum to zero, thus solving the homogeneous differential equation.^[5] In order to solve for $r$ , one can substitute $y = e rx$ and its derivatives into the differential equation to get

a_{n}r^{n}e^{rx}+a_{n-1}r^{n-1}e^{rx}+\cdots +a_{1}re^{rx}+a_{0}e^{rx}=0

Since $e rx$ can never equal zero, it can be divided out, giving the characteristic equation

a_{n}r^{n}+a_{n-1}r^{n-1}+\cdots +a_{1}r+a_{0}=0.

By solving for the roots, $r$ , in this characteristic equation, one can find the general solution to the differential equation.^[1]^[6] For example, if $r$ has roots equal to 3, 11, and 40, then the general solution will be $y(x)=c_{1}e^{3x}+c_{2}e^{11x}+c_{3}e^{40x}$ , where $c_{1}$ , $c_{2}$ , and $c_{3}$ are arbitrary constants which need to be determined by the boundary and/or initial conditions.

Formation of the general solution

Solving the characteristic equation for its roots, $r 1, ..., r n$ , allows one to find the general solution of the differential equation. The roots may be real or complex, as well as distinct or repeated. If a characteristic equation has parts with distinct real roots, $h$ repeated roots, or $k$ complex roots corresponding to general solutions of $y D (x)$ , $y R 1 (x), ..., y R h (x)$ , and $y C 1 (x), ..., y C k (x)$ , respectively, then the general solution to the differential equation is

y(x)=y_{\mathrm {D} }(x)+y_{\mathrm {R} _{1}}(x)+\cdots +y_{\mathrm {R} _{h}}(x)+y_{\mathrm {C} _{1}}(x)+\cdots +y_{\mathrm {C} _{k}}(x)

Example

The linear homogeneous differential equation with constant coefficients

y^{(5)}+y^{(4)}-4y^{(3)}-16y''-20y'-12y=0

has the characteristic equation

r^{5}+r^{4}-4r^{3}-16r^{2}-20r-12=0

By factoring the characteristic equation into

(r-3)(r^{2}+2r+2)^{2}=0

one can see that the solutions for $r$ are the distinct single root $r 1 = 3$ and the double complex roots $r 2,3,4,5 = 1 \pm i$ . This corresponds to the real-valued general solution

y(x)=c_{1}e^{3x}+e^{x}(c_{2}\cos x+c_{3}\sin x)+xe^{x}(c_{4}\cos x+c_{5}\sin x)

with constants $c 1, ..., c 5$ .

Distinct real roots

The superposition principle for linear homogeneous differential equations says that if $u 1, ..., u n$ are $n$ linearly independent solutions to a particular differential equation, then $c 1 u 1 + \dots + c n u n$ is also a solution for all values $c 1, ..., c n$ .^[1]^[7] Therefore, if the characteristic equation has distinct real roots $r 1, ..., r n$ , then a general solution will be of the form

y_{\mathrm {D} }(x)=c_{1}e^{r_{1}x}+c_{2}e^{r_{2}x}+\cdots +c_{n}e^{r_{n}x}

Repeated real roots

If the characteristic equation has a root $r 1$ that is repeated $k$ times, then it is clear that $y p (x) = c 1 e r 1 x$ is at least one solution.^[1] However, this solution lacks linearly independent solutions from the other $k - 1$ roots. Since $r 1$ has multiplicity $k$ , the differential equation can be factored into^[1]

\left({\frac {d}{dx}}-r_{1}\right)^{k}y=0.

The fact that $y p (x) = c 1 e r 1 x$ is one solution allows one to presume that the general solution may be of the form $y (x) = u (x) e r 1 x$ , where $u (x)$ is a function to be determined. Substituting $ue r 1 x$ gives

\left({\frac {d}{dx}}-r_{1}\right)\!ue^{r_{1}x}={\frac {d}{dx}}\left(ue^{r_{1}x}\right)-r_{1}ue^{r_{1}x}={\frac {d}{dx}}(u)e^{r_{1}x}+r_{1}ue^{r_{1}x}-r_{1}ue^{r_{1}x}={\frac {d}{dx}}(u)e^{r_{1}x}

when $k = 1$ . By applying this fact $k$ times, it follows that

\left({\frac {d}{dx}}-r_{1}\right)^{k}ue^{r_{1}x}={\frac {d^{k}}{dx^{k}}}(u)e^{r_{1}x}=0.

By dividing out $e r 1 x$ , it can be seen that

{\frac {d^{k}}{dx^{k}}}(u)=u^{(k)}=0.

Therefore, the general case for $u (x)$ is a polynomial of degree $k - 1$ , so that $u (x) = c 1 + c 2 x + c 3 x 2 + \dots + c k x k -1$ .^[6] Since $y (x) = ue r 1 x$ , the part of the general solution corresponding to $r 1$ is

y_{\mathrm {R} }(x)=e^{r_{1}x}\!\left(c_{1}+c_{2}x+\cdots +c_{k}x^{k-1}\right).

Complex roots

If a second-order differential equation has a characteristic equation with complex conjugate roots of the form $r 1 = a + bi$ and $r 2 = a - bi$ , then the general solution is accordingly $y (x) = c 1 e (a + bi) x + c 2 e (a - bi) x$ . By Euler's formula, which states that $e iθ = cos θ + i sin θ$ , this solution can be rewritten as follows:

{\begin{aligned}y(x)&=c_{1}e^{(a+bi)x}+c_{2}e^{(a-bi)x}\\&=c_{1}e^{ax}(\cos bx+i\sin bx)+c_{2}e^{ax}(\cos bx-i\sin bx)\\&=\left(c_{1}+c_{2}\right)e^{ax}\cos bx+i(c_{1}-c_{2})e^{ax}\sin bx\end{aligned}}

where $c 1$ and $c 2$ are constants that can be non-real and which depend on the initial conditions.^[6] (Indeed, since $y (x)$ is real, $c 1 - c 2$ must be imaginary or zero and $c 1 + c 2$ must be real, in order for both terms after the last equals sign to be real.)

For example, if $c1 = c2 = .mw-parser-output .sfrac{white-space:nowrap}.mw-parser-output .sfrac.tion,.mw-parser-output .sfrac .tion{display:inline-block;vertical-align:-0.5em;font-size:85%;text-align:center}.mw-parser-output .sfrac .num,.mw-parser-output .sfrac .den{display:block;line-height:1em;margin:0 0.1em}.mw-parser-output .sfrac .den{border-top:1px solid}.mw-parser-output .sr-only{border:0;clip:rect(0,0,0,0);height:1px;margin:-1px;overflow:hidden;padding:0;position:absolute;width:1px}1/2$ , then the particular solution $y 1 (x) = e ax cos bx$ is formed. Similarly, if $c 1 = 1 / 2 i$ and $c 2 = - 1 / 2 i$ , then the independent solution formed is $y 2 (x) = e ax sin bx$ . Thus by the superposition principle for linear homogeneous differential equations, a second-order differential equation having complex roots $r = a \pm bi$ will result in the following general solution:

y_{\mathrm {C} }(x)=e^{ax}(C_{1}\cos bx+C_{2}\sin bx)

This analysis also applies to the parts of the solutions of a higher-order differential equation whose characteristic equation involves non-real complex conjugate roots.

Related Research Articles

In mathematics, a complex number is an element of a number system that extends the real numbers with a specific element denoted $i$ , called the imaginary unit and satisfying the equation $; every complex number can be expressed in the form, where a and b are real numbers. Because no real number satisfies the above equation, i was called an imaginary number by René Descartes. For the complex number, a is called the real part, and b is called the imaginary part . The set of complex numbers is denoted by either of the symbols or C . Despite the historical nomenclature "imaginary", complex numbers are regarded in the mathematical sciences as just as "real" as the real numbers and are fundamental in many aspects of the scientific description of the natural world.$

In algebra, a quadratic equation is any equation that can be rearranged in standard form as

In mathematics, the trigonometric functions are real functions which relate an angle of a right-angled triangle to ratios of two side lengths. They are widely used in all sciences that are related to geometry, such as navigation, solid mechanics, celestial mechanics, geodesy, and many others. They are among the simplest periodic functions, and as such are also widely used for studying periodic phenomena through Fourier analysis.

In mathematics and physics, Laplace's equation is a second-order partial differential equation named after Pierre-Simon Laplace, who first studied its properties. This is often written as

<span class="mw-page-title-main">Navier–Stokes equations</span> Equations describing the motion of viscous fluid substances

The Navier–Stokes equations are partial differential equations which describe the motion of viscous fluid substances, named after French engineer and physicist Claude-Louis Navier and Irish physicist and mathematician George Gabriel Stokes. They were developed over several decades of progressively building the theories, from 1822 (Navier) to 1842–1850 (Stokes).

In mathematics, the Lambert $W$ function, also called the omega function or product logarithm, is a multivalued function, namely the branches of the converse relation of the function $f (w) = we w$ , where $w$ is any complex number and $e w$ is the exponential function.

In mathematics, the Laplace operator or Laplacian is a differential operator given by the divergence of the gradient of a scalar function on Euclidean space. It is usually denoted by the symbols $, (where is the nabla operator), or . In a Cartesian coordinate system, the Laplacian is given by the sum of second partial derivatives of the function with respect to each independent variable. In other coordinate systems, such as cylindrical and spherical coordinates, the Laplacian also has a useful form. Informally, the Laplacian Δ f (p) of a function f at a point p measures by how much the average value of f over small spheres or balls centered at p deviates from f (p) .$

<span class="mw-page-title-main">Cubic equation</span> Polynomial equation of degree 3

In algebra, a cubic equation in one variable is an equation of the form

<span class="mw-page-title-main">Quadratic function</span> Polynomial function of degree two

In mathematics, a quadratic polynomial is a polynomial of degree two in one or more variables. A quadratic function is the polynomial function defined by a quadratic polynomial. Before the 20th century, the distinction was unclear between a polynomial and its associated polynomial function; so "quadratic polynomial" and "quadratic function" were almost synonymous. This is still the case in many elementary courses, where both terms are often abbreviated as "quadratic".

In mathematics, a Green's function is the impulse response of an inhomogeneous linear differential operator defined on a domain with specified initial conditions or boundary conditions.

In mathematics, a linear differential equation is a differential equation that is defined by a linear polynomial in the unknown function and its derivatives, that is an equation of the form

In physics, the Hamilton–Jacobi equation, named after William Rowan Hamilton and Carl Gustav Jacob Jacobi, is an alternative formulation of classical mechanics, equivalent to other formulations such as Newton's laws of motion, Lagrangian mechanics and Hamiltonian mechanics.

In mathematics, an Euler–Cauchy equation, or Cauchy–Euler equation, or simply Euler's equation is a linear homogeneous ordinary differential equation with variable coefficients. It is sometimes referred to as an equidimensional equation. Because of its particularly simple equidimensional structure, the differential equation can be solved explicitly.

In mathematics, the Euler–Tricomi equation is a linear partial differential equation useful in the study of transonic flow. It is named after mathematicians Leonhard Euler and Francesco Giacomo Tricomi.

In mathematics, delay differential equations (DDEs) are a type of differential equation in which the derivative of the unknown function at a certain time is given in terms of the values of the function at previous times. DDEs are also called time-delay systems, systems with aftereffect or dead-time, hereditary systems, equations with deviating argument, or differential-difference equations. They belong to the class of systems with the functional state, i.e. partial differential equations (PDEs) which are infinite dimensional, as opposed to ordinary differential equations (ODEs) having a finite dimensional state vector. Four points may give a possible explanation of the popularity of DDEs:

Aftereffect is an applied problem: it is well known that, together with the increasing expectations of dynamic performances, engineers need their models to behave more like the real process. Many processes include aftereffect phenomena in their inner dynamics. In addition, actuators, sensors, and communication networks that are now involved in feedback control loops introduce such delays. Finally, besides actual delays, time lags are frequently used to simplify very high order models. Then, the interest for DDEs keeps on growing in all scientific areas and, especially, in control engineering.
Delay systems are still resistant to many classical controllers: one could think that the simplest approach would consist in replacing them by some finite-dimensional approximations. Unfortunately, ignoring effects which are adequately represented by DDEs is not a general alternative: in the best situation, it leads to the same degree of complexity in the control design. In worst cases, it is potentially disastrous in terms of stability and oscillations.
Voluntary introduction of delays can benefit the control system.
In spite of their complexity, DDEs often appear as simple infinite-dimensional models in the very complex area of partial differential equations (PDEs).

The gradient theorem, also known as the fundamental theorem of calculus for line integrals, says that a line integral through a gradient field can be evaluated by evaluating the original scalar field at the endpoints of the curve. The theorem is a generalization of the second fundamental theorem of calculus to any curve in a plane or space rather than just the real line.

In mathematics, the spectral theory of ordinary differential equations is the part of spectral theory concerned with the determination of the spectrum and eigenfunction expansion associated with a linear ordinary differential equation. In his dissertation, Hermann Weyl generalized the classical Sturm–Liouville theory on a finite closed interval to second order differential operators with singularities at the endpoints of the interval, possibly semi-infinite or infinite. Unlike the classical case, the spectrum may no longer consist of just a countable set of eigenvalues, but may also contain a continuous part. In this case the eigenfunction expansion involves an integral over the continuous part with respect to a spectral measure, given by the Titchmarsh–Kodaira formula. The theory was put in its final simplified form for singular differential equations of even degree by Kodaira and others, using von Neumann's spectral theorem. It has had important applications in quantum mechanics, operator theory and harmonic analysis on semisimple Lie groups.

<span class="mw-page-title-main">Vibrations of a circular membrane</span>

A two-dimensional elastic membrane under tension can support transverse vibrations. The properties of an idealized drumhead can be modeled by the vibrations of a circular membrane of uniform thickness, attached to a rigid frame. Due to the phenomenon of resonance, at certain vibration frequencies, its resonant frequencies, the membrane can store vibrational energy, the surface moving in a characteristic pattern of standing waves. This is called a normal mode. A membrane has an infinite number of these normal modes, starting with a lowest frequency one called the fundamental mode.

In mathematics, a linear recurrence with constant coefficients sets equal to 0 a polynomial that is linear in the various iterates of a variable—that is, in the values of the elements of a sequence. The polynomial's linearity means that each of its terms has degree 0 or 1. A linear recurrence denotes the evolution of some variable over time, with the current time period or discrete moment in time denoted as $t$ , one period earlier denoted as $t - 1$ , one period later as $t + 1$ , etc.

In mathematics, the exponential response formula (ERF), also known as exponential response and complex replacement, is a method used to find a particular solution of a non-homogeneous linear ordinary differential equation of any order. The exponential response formula is applicable to non-homogeneous linear ordinary differential equations with constant coefficients if the function is polynomial, sinusoidal, exponential or the combination of the three. The general solution of a non-homogeneous linear ordinary differential equation is a superposition of the general solution of the associated homogeneous ODE and a particular solution to the non-homogeneous ODE. Alternative methods for solving ordinary differential equations of higher order are method of undetermined coefficients and method of variation of parameters.

References

1 2 3 4 5 6 7 Edwards, C. Henry; Penney, David E. (2008). "Chapter 3". Differential Equations: Computing and Modeling. David Calvis. Upper Saddle River, New Jersey: Pearson Education. pp. 156–170. ISBN 978-0-13-600438-7.
1 2 3 Smith, David Eugene. "History of Modern Mathematics: Differential Equations". University of South Florida.
↑ Baumol, William J. (1970). Economic Dynamics (3rd ed.). p. 172.
↑ Chiang, Alpha (1984). Fundamental Methods of Mathematical Economics (3rd ed.). pp. 578, 600. ISBN 9780070107809.
1 2 Chu, Herman; Shah, Gaurav; Macall, Tom. "Linear Homogeneous Ordinary Differential Equations with Constant Coefficients". eFunda. Retrieved 1 March 2011.
1 2 3 4 5 Cohen, Abraham (1906). An Elementary Treatise on Differential Equations. D. C. Heath and Company.
↑ Dawkins, Paul. "Differential Equation Terminology". Paul's Online Math Notes. Retrieved 2 March 2011.

This page is based on this Wikipedia article
Text is available under the CC BY-SA 4.0 license; additional terms may apply.
Images, videos and audio are available under their respective licenses.

[edwards-1] 1 2 3 4 5 6 7 Edwards, C. Henry; Penney, David E. (2008). "Chapter 3". Differential Equations: Computing and Modeling. David Calvis. Upper Saddle River, New Jersey: Pearson Education. pp. 156–170. ISBN 978-0-13-600438-7.

[smith-2] 1 2 3 Smith, David Eugene. "History of Modern Mathematics: Differential Equations". University of South Florida.

[3] Baumol, William J. (1970). Economic Dynamics (3rd ed.). p. 172.

[4] Chiang, Alpha (1984). Fundamental Methods of Mathematical Economics (3rd ed.). pp. 578, 600. ISBN 9780070107809.

[eFunda-5] 1 2 Chu, Herman; Shah, Gaurav; Macall, Tom. "Linear Homogeneous Ordinary Differential Equations with Constant Coefficients". eFunda. Retrieved 1 March 2011.

[cohen-6] 1 2 3 4 5 Cohen, Abraham (1906). An Elementary Treatise on Differential Equations. D. C. Heath and Company.

[dawkins-7] Dawkins, Paul. "Differential Equation Terminology". Paul's Online Math Notes. Retrieved 2 March 2011.

[1]

[2]

[3]

[4]

[5]

[6]

[7]