Movatterモバイル変換


[0]ホーム

URL:


Jump to content
WikipediaThe Free Encyclopedia
Search

Ordinary differential equation

From Wikipedia, the free encyclopedia
Differential equation containing derivatives with respect to only one variable
parabolic projectile motion showing velocity vector
Thetrajectory of aprojectile launched from acannon follows a curve determined by an ordinary differential equation that is derived from Newton's second law.
Differential equations
Scope
Classification
Solution
People

In mathematics, anordinary differential equation (ODE) is adifferential equation (DE) dependent on only a single independentvariable. As with any other DE, its unknown(s) consists of one (or more)function(s) and involves thederivatives of those functions.[1] The term "ordinary" is used in contrast withpartial differential equations (PDEs) which may be with respect tomore than one independent variable,[2] and, less commonly, in contrast withstochastic differential equations (SDEs) where the progression is random.[3]

Differential equations

[edit]

Alinear differential equation is a differential equation that is defined by alinear polynomial in the unknown function and its derivatives, that is anequation of the form

a0(x)y+a1(x)y+a2(x)y++an(x)y(n)+b(x)=0,{\displaystyle a_{0}(x)y+a_{1}(x)y'+a_{2}(x)y''+\cdots +a_{n}(x)y^{(n)}+b(x)=0,}

wherea0(x),,an(x){\displaystyle a_{0}(x),\ldots ,a_{n}(x)} andb(x){\displaystyle b(x)} are arbitrarydifferentiable functions that do not need to be linear, andy,,y(n){\displaystyle y',\ldots ,y^{(n)}} are the successive derivatives of the unknown functiony{\displaystyle y} of the variablex{\displaystyle x}.[4]

Among ordinary differential equations, linear differential equations play a prominent role for several reasons. Mostelementary andspecial functions that are encountered inphysics andapplied mathematics are solutions of linear differential equations (seeHolonomic function). When physical phenomena are modeled with non-linear equations, they are generally approximated by linear differential equations for an easier solution. The few non-linear ODEs that can be solved explicitly are generally solved by transforming the equation into an equivalent linear ODE (see, for exampleRiccati equation).[5]

Some ODEs can be solved explicitly in terms of known functions andintegrals. When that is not possible, the equation for computing theTaylor series of the solutions may be useful. For applied problems,numerical methods for ordinary differential equations can supply an approximation of the solution[6].

Background

[edit]

Ordinary differential equations (ODEs) arise in many contexts ofmathematics andsocial andnatural sciences. Mathematical descriptions of change use differentials and derivatives. Various differentials, derivatives, and functions become related via equations, such that a differential equation is a result that describes dynamically changing phenomena, evolution, and variation. Often, quantities are defined as the rate of change of other quantities (for example, derivatives of displacement with respect to time), orgradients of quantities, which is how they enter differential equations.[7]

Specific mathematical fields includegeometry andanalytical mechanics. Scientific fields include much ofphysics andastronomy (celestial mechanics),meteorology (weather modeling),chemistry (reaction rates),[8]biology (infectious diseases, genetic variation),ecology andpopulation modeling (population competition),economics (stock trends, interest rates and the market equilibrium price changes).

Many mathematicians have studied differential equations and contributed to the field, includingNewton,Leibniz, theBernoulli family,Riccati,Clairaut,d'Alembert, andEuler.

A simple example isNewton's second law of motion—the relationship between the displacementx{\displaystyle x} and the timet{\displaystyle t} of an object under the forceF{\displaystyle F}, is given by the differential equation

md2x(t)dt2=F(x(t)){\displaystyle m{\frac {\mathrm {d} ^{2}x(t)}{\mathrm {d} t^{2}}}=F(x(t))\,}

which constrains themotion of a particle of constant massm{\displaystyle m}. In general,F{\displaystyle F} is a function of the positionx(t){\displaystyle x(t)} of the particle at timet{\displaystyle t}. The unknown functionx(t){\displaystyle x(t)} appears on both sides of the differential equation, and is indicated in the notationF(x(t)){\displaystyle F(x(t))}.[9][10][11][12]

Definitions

[edit]

In what follows,y{\displaystyle y} is adependent variable representing an unknown functiony=f(x){\displaystyle y=f(x)} of theindependent variablex{\displaystyle x}. Thenotation for differentiation varies depending upon the author and upon which notation is most useful for the task at hand. In this context, theLeibniz's notationdydx,d2ydx2,,dnydxn{\displaystyle {\frac {dy}{dx}},{\frac {d^{2}y}{dx^{2}}},\ldots ,{\frac {d^{n}y}{dx^{n}}}} is more useful for differentiation andintegration, whereasLagrange's notationy,y,,y(n){\displaystyle y',y'',\ldots ,y^{(n)}} is more useful for representinghigher-order derivatives compactly, andNewton's notation(y˙,y¨,y...){\displaystyle ({\dot {y}},{\ddot {y}},{\overset {...}{y}})} is often used in physics for representing derivatives of low order with respect to time.

General definition

[edit]
See also:Order of differential equation

GivenF{\displaystyle F}, a function ofx{\displaystyle x},y{\displaystyle y}, and derivatives ofy{\displaystyle y}. Then an equation of the form

F(x,y,y,,y(n1))=y(n){\displaystyle F\left(x,y,y',\ldots ,y^{(n-1)}\right)=y^{(n)}}

is called anexplicit ordinary differential equation of ordern{\displaystyle n}.[13][14]

More generally, animplicit ordinary differential equation of ordern{\displaystyle n} takes the form:[15]

F(x,y,y,y, , y(n))=0{\displaystyle F\left(x,y,y',y'',\ \ldots ,\ y^{(n)}\right)=0}

There are further classifications:

Autonomous
A differential equation isautonomous if it does not depend on the variablex.
Linear
A differential equation islinear ifF{\displaystyle F} can be written as alinear combination of the derivatives ofy{\displaystyle y}; that is, it can be rewritten as
y(n)=i=0n1ai(x)y(i)+r(x){\displaystyle y^{(n)}=\sum _{i=0}^{n-1}a_{i}(x)y^{(i)}+r(x)}
whereai(x){\displaystyle a_{i}(x)} andr(x){\displaystyle r(x)} are continuous functions ofx{\displaystyle x}.[13][16][17]The functionr(x){\displaystyle r(x)} is called thesource term, leading to further classification.[16][18]
Homogeneous
A linear differential equation ishomogeneous ifr(x)=0{\displaystyle r(x)=0}. In this case, there is always the "trivial solution"y=0{\displaystyle y=0}.
Nonhomogeneous (or inhomogeneous)
A linear differential equation isnonhomogeneous ifr(x)0{\displaystyle r(x)\neq 0}.
Non-linear
A differential equation that is not linear.

System of ODEs

[edit]
Main article:System of differential equations

A number of coupled differential equations form a system of equations. Ify{\displaystyle \mathbf {y} } is a vector whose elements are functions;y(x)=[y1(x),y2(x),,ym(x)]{\displaystyle \mathbf {y} (x)=[y_{1}(x),y_{2}(x),\ldots ,y_{m}(x)]}, andF{\displaystyle \mathbf {F} } is avector-valued function ofy{\displaystyle \mathbf {y} } and its derivatives, then

y(n)=F(x,y,y,y,,y(n1)){\displaystyle \mathbf {y} ^{(n)}=\mathbf {F} \left(x,\mathbf {y} ,\mathbf {y} ',\mathbf {y} '',\ldots ,\mathbf {y} ^{(n-1)}\right)}

is anexplicit system of ordinary differential equations ofordern{\displaystyle n} anddimensionm{\displaystyle m}. Incolumn vector form:

(y1(n)y2(n)ym(n))=(f1(x,y,y,y,,y(n1))f2(x,y,y,y,,y(n1))fm(x,y,y,y,,y(n1))){\displaystyle {\begin{pmatrix}y_{1}^{(n)}\\y_{2}^{(n)}\\\vdots \\y_{m}^{(n)}\end{pmatrix}}={\begin{pmatrix}f_{1}\left(x,\mathbf {y} ,\mathbf {y} ',\mathbf {y} '',\ldots ,\mathbf {y} ^{(n-1)}\right)\\f_{2}\left(x,\mathbf {y} ,\mathbf {y} ',\mathbf {y} '',\ldots ,\mathbf {y} ^{(n-1)}\right)\\\vdots \\f_{m}\left(x,\mathbf {y} ,\mathbf {y} ',\mathbf {y} '',\ldots ,\mathbf {y} ^{(n-1)}\right)\end{pmatrix}}}

These are not necessarily linear. Theimplicit analogue is:

F(x,y,y,y,,y(n))=0{\displaystyle \mathbf {F} \left(x,\mathbf {y} ,\mathbf {y} ',\mathbf {y} '',\ldots ,\mathbf {y} ^{(n)}\right)={\boldsymbol {0}}}

where0=(0,0,,0){\displaystyle {\boldsymbol {0}}=(0,0,\ldots ,0)} is thezero vector. In matrix form

(f1(x,y,y,y,,y(n))f2(x,y,y,y,,y(n))fm(x,y,y,y,,y(n)))=(000){\displaystyle {\begin{pmatrix}f_{1}(x,\mathbf {y} ,\mathbf {y} ',\mathbf {y} '',\ldots ,\mathbf {y} ^{(n)})\\f_{2}(x,\mathbf {y} ,\mathbf {y} ',\mathbf {y} '',\ldots ,\mathbf {y} ^{(n)})\\\vdots \\f_{m}(x,\mathbf {y} ,\mathbf {y} ',\mathbf {y} '',\ldots ,\mathbf {y} ^{(n)})\end{pmatrix}}={\begin{pmatrix}0\\0\\\vdots \\0\end{pmatrix}}}

For a system of the formF(x,y,y)=0{\displaystyle \mathbf {F} \left(x,\mathbf {y} ,\mathbf {y} '\right)={\boldsymbol {0}}}, some sources also require that theJacobian matrixF(x,u,v)v{\displaystyle {\frac {\partial \mathbf {F} (x,\mathbf {u} ,\mathbf {v} )}{\partial \mathbf {v} }}} benon-singular in order to call this an implicit ODE [system]; an implicit ODE system satisfying this Jacobian non-singularity condition can be transformed into an explicit ODE system. In the same sources, implicit ODE systems with a singular Jacobian are termeddifferential algebraic equations (DAEs). This distinction is not merely one of terminology; DAEs have fundamentally different characteristics and are generally more involved to solve than (nonsingular) ODE systems.[19][20][21] Presumably for additional derivatives, theHessian matrix and so forth are also assumed non-singular according to this scheme,[citation needed] although note thatany ODE of order greater than one can be (and usually is) rewritten as system of ODEs of first order,[22] which makes the Jacobian singularity criterion sufficient for this taxonomy to be comprehensive at all orders.

The behavior of a system of ODEs can be visualized through the use of aphase portrait.

Solutions

[edit]

Given a differential equation

F(x,y,y,,y(n))=0{\displaystyle F\left(x,y,y',\ldots ,y^{(n)}\right)=0}

a functionu:IRR{\displaystyle u:I\subset \mathbb {R} \to \mathbb {R} }, whereI{\displaystyle I} is an interval, is called asolution orintegral curve forF{\displaystyle F}, ifu{\displaystyle u} isn{\displaystyle n}-times differentiable onI{\displaystyle I}, and

F(x,u,u, , u(n))=0xI.{\displaystyle F(x,u,u',\ \ldots ,\ u^{(n)})=0\quad x\in I.}

Given two solutionsu:JRR{\displaystyle u:J\subset \mathbb {R} \to \mathbb {R} } andv:IRR{\displaystyle v:I\subset \mathbb {R} \to \mathbb {R} },u{\displaystyle u} is called anextension ofv{\displaystyle v} ifIJ{\displaystyle I\subset J} and

u(x)=v(x)xI.{\displaystyle u(x)=v(x)\quad x\in I.\,}

A solution that has no extension is called amaximal solution. A solution defined on all ofR{\displaystyle \mathbb {R} } is called aglobal solution.

Ageneral solution of ann{\displaystyle n}th-order equation is a solution containingn{\displaystyle n} arbitrary independentconstants of integration. Aparticular solution is derived from the general solution by setting the constants to particular values, often chosen to fulfill set 'initial conditions orboundary conditions'.[23] Asingular solution is a solution that cannot be obtained by assigning definite values to the arbitrary constants in the general solution.[24]

In the context of linear ODE, the terminologyparticular solution can also refer to any solution of the ODE (not necessarily satisfying the initial conditions), which is then added to thehomogeneous solution (a general solution of the homogeneous ODE), which then forms a general solution of the original ODE. This is the terminology used in theguessing method section in this article, and is frequently used when discussing themethod of undetermined coefficients andvariation of parameters.

Solutions of finite duration

[edit]

For non-linear autonomous ODEs it is possible under some conditions to develop solutions of finite duration,[25] meaning here that from its own dynamics, the system will reach the value zero at an ending time and stays there in zero forever after. These finite-duration solutions can't be analytical functions on the whole real line, and because they will be non-Lipschitz functions at their ending time, they are not included in the uniqueness theorem of solutions of Lipschitz differential equations.

As example, the equation:

y=sgn(y)|y|,y(0)=1{\displaystyle y'=-{\text{sgn}}(y){\sqrt {|y|}},\,\,y(0)=1}

Admits the finite duration solution:

y(x)=14(1x2+|1x2|)2{\displaystyle y(x)={\frac {1}{4}}\left(1-{\frac {x}{2}}+\left|1-{\frac {x}{2}}\right|\right)^{2}}

Theories

[edit]

Singular solutions

[edit]

The theory ofsingular solutions of ordinary andpartial differential equations was a subject of research from the time of Leibniz, but only since the middle of the nineteenth century has it received special attention. A valuable but little-known work on the subject is that of Houtain (1854).Darboux (from 1873) was a leader in the theory, and in the geometric interpretation of these solutions he opened a field worked by various writers, notablyCasorati andCayley. To the latter is due (1872) the theory of singular solutions of differential equations of the first order as accepted circa 1900.

Reduction to quadratures

[edit]

The primitive attempt in dealing with differential equations had in view a reduction toquadratures, that is, expressing the solutions in terms of known function and their integrals. This is possible for linear equations with constant coefficients, it appeared in the 19th century that this is generally impossible in other cases. Hence, analysts began the study (for their own) of functions that are solutions of differential equations, thus opening a new and fertile field.Cauchy was the first to appreciate the importance of this view. Thereafter, the real question was no longer whether a solution is possible by quadratures, but whether a given differential equation suffices for the definition of a function, and, if so, what are the characteristic properties of such functions.

Fuchsian theory

[edit]
Main article:Frobenius method

Two memoirs byFuchs[26] inspired a novel approach, subsequently elaborated by Thomé andFrobenius. Collet was a prominent contributor beginning in 1869. His method for integrating a non-linear system was communicated to Bertrand in 1868.Clebsch (1873) attacked the theory along lines parallel to those in his theory ofAbelian integrals. As the latter can be classified according to the properties of the fundamental curve that remains unchanged under a rational transformation, Clebsch proposed to classify the transcendent functions defined by differential equations according to the invariant properties of the corresponding surfacesf=0{\displaystyle f=0} under rational one-to-one transformations.

Lie's theory

[edit]

From 1870,Sophus Lie's work put the theory of differential equations on a better foundation. He showed that the integration theories of the older mathematicians can, usingLie groups, be referred to a common source, and that ordinary differential equations that admit the sameinfinitesimal transformations present comparable integration difficulties. He also emphasized the subject oftransformations of contact.

Lie's group theory of differential equations has been certified, namely: (1) that it unifies the many ad hoc methods known for solving differential equations, and (2) that it provides powerful new ways to find solutions. The theory has applications to both ordinary and partial differential equations.[27]

A general solution approach uses the symmetry property of differential equations, the continuousinfinitesimal transformations of solutions to solutions (Lie theory). Continuousgroup theory,Lie algebras, anddifferential geometry are used to understand the structure of linear and non-linear (partial) differential equations for generating integrable equations, to find itsLax pairs, recursion operators,Bäcklund transform, and finally finding exact analytic solutions to DE.

Symmetry methods have been applied to differential equations that arise in mathematics, physics, engineering, and other disciplines.

Sturm–Liouville theory

[edit]
Main article:Sturm–Liouville theory

Sturm–Liouville theory is a theory of a special type of second-order linear ordinary differential equation. Their solutions are based oneigenvalues and correspondingeigenfunctions of linear operators defined via second-orderhomogeneous linear equations. The problems are identified as Sturm–Liouville problems (SLP) and are named afterJ. C. F. Sturm andJ. Liouville, who studied them in the mid-1800s. SLPs have an infinite number of eigenvalues, and the corresponding eigenfunctions form a complete, orthogonal set, which makes orthogonal expansions possible. This is a key idea in applied mathematics, physics, and engineering.[28] SLPs are also useful in the analysis of certain partial differential equations.

Existence and uniqueness of solutions

[edit]

There are several theorems that establish existence and uniqueness of solutions toinitial value problems involving ODEs both locally and globally. The two main theorems are

TheoremAssumptionConclusion
Peano existence theoremF{\displaystyle F}continuouslocal existence only
Picard–Lindelöf theoremF{\displaystyle F}Lipschitz continuouslocal existence and uniqueness

In their basic form both of these theorems only guarantee local results, though the latter can be extended to give a global result, for example, if the conditions ofGrönwall's inequality are met.

Also, uniqueness theorems like the Lipschitz one above do not apply toDAE systems, which may have multiple solutions stemming from their (non-linear) algebraic part alone.[29]

Local existence and uniqueness theorem simplified

[edit]

The theorem can be stated simply as follows.[30] For the equation and initial value problem:y=F(x,y),y0=y(x0){\displaystyle y'=F(x,y)\,,\quad y_{0}=y(x_{0})}ifF{\displaystyle F} andF/y{\displaystyle \partial F/\partial y} are continuous in a closed rectangleR=[x0a,x0+a]×[y0b,y0+b]{\displaystyle R=[x_{0}-a,x_{0}+a]\times [y_{0}-b,y_{0}+b]}in thexy{\displaystyle x-y} plane, wherea{\displaystyle a} andb{\displaystyle b} arereal (symbolically:a,bR{\displaystyle a,b\in \mathbb {R} }) and×{\displaystyle \times } denotes theCartesian product, square brackets denoteclosed intervals, then there is an intervalI=[x0h,x0+h][x0a,x0+a]{\displaystyle I=[x_{0}-h,x_{0}+h]\subset [x_{0}-a,x_{0}+a]}for somehR{\displaystyle h\in \mathbb {R} } wherethe solution to the above equation and initial value problem can be found. That is, there is a solution and it is unique. Since there is no restriction onF{\displaystyle F} to be linear, this applies to non-linear equations that take the formF(x,y){\displaystyle F(x,y)}, and it can also be applied to systems of equations.

Global uniqueness and maximum domain of solution

[edit]

When the hypotheses of the Picard–Lindelöf theorem are satisfied, then local existence and uniqueness can be extended to a global result. More precisely:[31]

For each initial condition(x0,y0){\displaystyle (x_{0},y_{0})} there exists a unique maximum (possibly infinite) open interval

Imax=(x,x+),x±R{±},x0Imax{\displaystyle I_{\max }=(x_{-},x_{+}),x_{\pm }\in \mathbb {R} \cup \{\pm \infty \},x_{0}\in I_{\max }}

such that any solution that satisfies this initial condition is arestriction of the solution that satisfies this initial condition with domainImax{\displaystyle I_{\max }}.

In the case thatx±±{\displaystyle x_{\pm }\neq \pm \infty }, there are exactly two possibilities

whereΩ{\displaystyle \Omega } is the open set in whichF{\displaystyle F} is defined, andΩ¯{\displaystyle \partial {\bar {\Omega }}} is its boundary.

Note that the maximum domain of the solution

Example.
y=y2{\displaystyle y'=y^{2}}

This means thatF(x,y)=y2{\displaystyle F(x,y)=y^{2}}, which isC1{\displaystyle C^{1}} and therefore locally Lipschitz continuous, satisfying the Picard–Lindelöf theorem.

Even in such a simple setting, the maximum domain of solution cannot be allR{\displaystyle \mathbb {R} } since the solution is

y(x)=y0(x0x)y0+1{\displaystyle y(x)={\frac {y_{0}}{(x_{0}-x)y_{0}+1}}}

which has maximum domain:

{Ry0=0(,x0+1y0)y0>0(x0+1y0,+)y0<0{\displaystyle {\begin{cases}\mathbb {R} &y_{0}=0\\[4pt]\left(-\infty ,x_{0}+{\frac {1}{y_{0}}}\right)&y_{0}>0\\[4pt]\left(x_{0}+{\frac {1}{y_{0}}},+\infty \right)&y_{0}<0\end{cases}}}

This shows clearly that the maximum interval may depend on the initial conditions. The domain ofy{\displaystyle y} could be taken as beingR(x0+1/y0),{\displaystyle \mathbb {R} \setminus (x_{0}+1/y_{0}),} but this would lead to a domain that is not an interval, so that the side opposite to the initial condition would be disconnected from the initial condition, and therefore not uniquely determined by it.

The maximum domain is notR{\displaystyle \mathbb {R} } because

limxx±y(x),{\displaystyle \lim _{x\to x_{\pm }}\|y(x)\|\to \infty ,}

which is one of the two possible cases according to the above theorem.

Reduction of order

[edit]

Differential equations are usually easier to solve if theorder of the equation can be reduced.

Reduction to a first-order system

[edit]

Any explicit differential equation of ordern{\displaystyle n},

F(x,y,y,y, , y(n1))=y(n){\displaystyle F\left(x,y,y',y'',\ \ldots ,\ y^{(n-1)}\right)=y^{(n)}}

can be written as a system ofn{\displaystyle n} first-order differential equations by defining a new family of unknown functions

yi=y(i1).{\displaystyle y_{i}=y^{(i-1)}.\!}

fori=1,2,,n{\displaystyle i=1,2,\ldots ,n}. Then{\displaystyle n}-dimensional system of first-order coupled differential equations is then

y1=y2y2=y3yn1=ynyn=F(x,y1,,yn).{\displaystyle {\begin{array}{rcl}y_{1}'&=&y_{2}\\y_{2}'&=&y_{3}\\&\vdots &\\y_{n-1}'&=&y_{n}\\y_{n}'&=&F(x,y_{1},\ldots ,y_{n}).\end{array}}}

more compactly in vector notation:

y=F(x,y){\displaystyle \mathbf {y} '=\mathbf {F} (x,\mathbf {y} )}

where

y=(y1,,yn),F(x,y1,,yn)=(y2,,yn,F(x,y1,,yn)).{\displaystyle \mathbf {y} =(y_{1},\ldots ,y_{n}),\quad \mathbf {F} (x,y_{1},\ldots ,y_{n})=(y_{2},\ldots ,y_{n},F(x,y_{1},\ldots ,y_{n})).}

Summary of exact solutions

[edit]

Some differential equations have solutions that can be written in an exact and closed form. Several important classes are given here.

In the table below,P(x){\displaystyle P(x)},Q(x){\displaystyle Q(x)},P(y){\displaystyle P(y)},Q(y){\displaystyle Q(y)}, andM(x,y){\displaystyle M(x,y)},N(x,y){\displaystyle N(x,y)} are anyintegrable functions ofx{\displaystyle x},y{\displaystyle y};b{\displaystyle b} andc{\displaystyle c} are real given constants;C1,C2,{\displaystyle C_{1},C_{2},\ldots } are arbitrary constants (complex in general). The differential equations are in their equivalent and alternative forms that lead to the solution through integration.

In the integral solutions,λ{\displaystyle \lambda } andε{\displaystyle \varepsilon } are dummy variables of integration (the continuum analogues of indices insummation), and the notationxF(λ)dλ{\displaystyle \int ^{x}F(\lambda )\,d\lambda } just means to integrateF(λ){\displaystyle F(\lambda )} with respect toλ{\displaystyle \lambda }, thenafter the integration substituteλ=x{\displaystyle \lambda =x}, without adding constants (explicitly stated).

Separable equations

[edit]
Differential equationSolution methodGeneral solution
First-order, separable inx{\displaystyle x} andy{\displaystyle y} (general case, see below for special cases)[32]

P1(x)Q1(y)+P2(x)Q2(y)dydx=0P1(x)Q1(y)dx+P2(x)Q2(y)dy=0{\displaystyle {\begin{aligned}P_{1}(x)Q_{1}(y)+P_{2}(x)Q_{2}(y)\,{\frac {dy}{dx}}&=0\\P_{1}(x)Q_{1}(y)\,dx+P_{2}(x)Q_{2}(y)\,dy&=0\end{aligned}}}

Separation of variables (divide byP2Q1{\displaystyle P_{2}Q_{1}}).xP1(λ)P2(λ)dλ+yQ2(λ)Q1(λ)dλ=C{\displaystyle \int ^{x}{\frac {P_{1}(\lambda )}{P_{2}(\lambda )}}\,d\lambda +\int ^{y}{\frac {Q_{2}(\lambda )}{Q_{1}(\lambda )}}\,d\lambda =C}
First-order, separable inx{\displaystyle x}[30]

dydx=F(x)dy=F(x)dx{\displaystyle {\begin{aligned}{\frac {dy}{dx}}&=F(x)\\dy&=F(x)\,dx\end{aligned}}}

Direct integration.y=xF(λ)dλ+C{\displaystyle y=\int ^{x}F(\lambda )\,d\lambda +C}
First-order, autonomous, separable iny{\displaystyle y}[30]

dydx=F(y)dy=F(y)dx{\displaystyle {\begin{aligned}{\frac {dy}{dx}}&=F(y)\\dy&=F(y)\,dx\end{aligned}}}

Separation of variables (divide byF{\displaystyle F}).x=ydλF(λ)+C{\displaystyle x=\int ^{y}{\frac {d\lambda }{F(\lambda )}}+C}
First-order, separable inx{\displaystyle x} andy{\displaystyle y}[30]

P(y)dydx+Q(x)=0P(y)dy+Q(x)dx=0{\displaystyle {\begin{aligned}P(y){\frac {dy}{dx}}+Q(x)&=0\\P(y)\,dy+Q(x)\,dx&=0\end{aligned}}}

Integrate throughout.yP(λ)dλ+xQ(λ)dλ=C{\displaystyle \int ^{y}P(\lambda )\,d\lambda +\int ^{x}Q(\lambda )\,d\lambda =C}

General first-order equations

[edit]
Differential equationSolution methodGeneral solution
First-order, homogeneous[30]

dydx=F(yx){\displaystyle {\frac {dy}{dx}}=F\left({\frac {y}{x}}\right)}

Sety = ux, then solve by separation of variables inu andx.ln(Cx)=y/xdλF(λ)λ{\displaystyle \ln(Cx)=\int ^{y/x}{\frac {d\lambda }{F(\lambda )-\lambda }}}
First-order, separable[32]

yM(xy)+xN(xy)dydx=0yM(xy)dx+xN(xy)dy=0{\displaystyle {\begin{aligned}yM(xy)+xN(xy)\,{\frac {dy}{dx}}&=0\\yM(xy)\,dx+xN(xy)\,dy&=0\end{aligned}}}

Separation of variables (divide byxy{\displaystyle xy}).

ln(Cx)=xyN(λ)dλλ[N(λ)M(λ)]{\displaystyle \ln(Cx)=\int ^{xy}{\frac {N(\lambda )\,d\lambda }{\lambda [N(\lambda )-M(\lambda )]}}}

IfN=M{\displaystyle N=M}, the solution isxy=C{\displaystyle xy=C}.

Exact differential, first-order[30]

M(x,y)dydx+N(x,y)=0M(x,y)dy+N(x,y)dx=0{\displaystyle {\begin{aligned}M(x,y){\frac {dy}{dx}}+N(x,y)&=0\\M(x,y)\,dy+N(x,y)\,dx&=0\end{aligned}}}

whereMy=Nx{\displaystyle {\frac {\partial M}{\partial y}}={\frac {\partial N}{\partial x}}}

Integrate throughout.F(x,y)=xM(λ,y)dλ+yY(λ)dλ=yN(x,λ)dλ+xX(λ)dλ=C{\displaystyle {\begin{aligned}F(x,y)&=\int ^{x}M(\lambda ,y)\,d\lambda +\int ^{y}Y(\lambda )\,d\lambda \\&=\int ^{y}N(x,\lambda )\,d\lambda +\int ^{x}X(\lambda )\,d\lambda =C\end{aligned}}}

whereY(y)=N(x,y)yxM(λ,y)dλ{\displaystyle Y(y)=N(x,y)-{\frac {\partial }{\partial y}}\int ^{x}M(\lambda ,y)\,d\lambda } andX(x)=M(x,y)xyN(x,λ)dλ{\displaystyle X(x)=M(x,y)-{\frac {\partial }{\partial x}}\int ^{y}N(x,\lambda )\,d\lambda }

Inexact differential, first-order[30]

M(x,y)dydx+N(x,y)=0M(x,y)dy+N(x,y)dx=0{\displaystyle {\begin{aligned}M(x,y){\frac {dy}{dx}}+N(x,y)&=0\\M(x,y)\,dy+N(x,y)\,dx&=0\end{aligned}}}

whereMyNx{\displaystyle {\frac {\partial M}{\partial y}}\neq {\frac {\partial N}{\partial x}}}

Integration factorμ(x,y){\displaystyle \mu (x,y)} satisfying

(μM)y=(μN)x{\displaystyle {\frac {\partial (\mu M)}{\partial y}}={\frac {\partial (\mu N)}{\partial x}}}

Ifμ(x,y){\displaystyle \mu (x,y)} can be found in a suitable way, then

F(x,y)=xμ(λ,y)M(λ,y)dλ+yY(λ)dλ=yμ(x,λ)N(x,λ)dλ+xX(λ)dλ=C{\displaystyle {\begin{aligned}F(x,y)=&\int ^{x}\mu (\lambda ,y)M(\lambda ,y)\,d\lambda +\int ^{y}Y(\lambda )\,d\lambda \\=&\int ^{y}\mu (x,\lambda )N(x,\lambda )\,d\lambda +\int ^{x}X(\lambda )\,d\lambda =C\end{aligned}}}

whereY(y)=N(x,y)yxμ(λ,y)M(λ,y)dλ{\displaystyle Y(y)=N(x,y)-{\frac {\partial }{\partial y}}\int ^{x}\mu (\lambda ,y)M(\lambda ,y)\,d\lambda }andX(x)=M(x,y)xyμ(x,λ)N(x,λ)dλ{\displaystyle X(x)=M(x,y)-{\frac {\partial }{\partial x}}\int ^{y}\mu (x,\lambda )N(x,\lambda )\,d\lambda }

General second-order equations

[edit]
Differential equationSolution methodGeneral solution
Second-order, autonomous[33]

d2ydx2=F(y){\displaystyle {\frac {d^{2}y}{dx^{2}}}=F(y)}

Multiply both sides of equation by2dydx{\displaystyle 2{\frac {dy}{dx}}}, substitute2dydxd2ydx2=ddx(dydx)2=2dydxF(y){\displaystyle 2{\frac {dy}{dx}}{\frac {d^{2}y}{dx^{2}}}={\frac {d}{dx}}\left({\frac {dy}{dx}}\right)^{2}=2{\frac {dy}{dx}}F(y)}, then integrate twice.x=±ydλ2λF(ε)dε+C1+C2{\displaystyle x=\pm \int ^{y}{\frac {d\lambda }{\sqrt {2\int ^{\lambda }F(\varepsilon )\,d\varepsilon +C_{1}}}}+C_{2}}

Linear to thenth order equations

[edit]
Differential equationSolution methodGeneral solution
First-order, linear, inhomogeneous, function coefficients[30]

dydx+P(x)y=Q(x){\displaystyle {\frac {dy}{dx}}+P(x)y=Q(x)}

Integrating factor:exP(λ)dλ.{\displaystyle e^{\int ^{x}P(\lambda )\,d\lambda }.}y=exP(λ)dλ[xeλP(ε)dεQ(λ)dλ+C]{\displaystyle y=e^{-\int ^{x}P(\lambda )\,d\lambda }\left[\int ^{x}e^{\int ^{\lambda }P(\varepsilon )\,d\varepsilon }Q(\lambda )\,d\lambda +C\right]}
Second-order, linear, inhomogeneous, function coefficients

d2ydx2+2p(x)dydx+(p(x)2+p(x))y=q(x){\displaystyle {\frac {d^{2}y}{dx^{2}}}+2p(x){\frac {dy}{dx}}+\left(p(x)^{2}+p'(x)\right)y=q(x)}

Integrating factor:exP(λ)dλ{\displaystyle e^{\int ^{x}P(\lambda )\,d\lambda }}y=exP(λ)dλ[x(ξeλP(ε)dεQ(λ)dλ)dξ+C1x+C2]{\displaystyle y=e^{-\int ^{x}P(\lambda )\,d\lambda }\left[\int ^{x}\left(\int ^{\xi }e^{\int ^{\lambda }P(\varepsilon )\,d\varepsilon }Q(\lambda )\,d\lambda \right)d\xi +C_{1}x+C_{2}\right]}
Second-order, linear, inhomogeneous, constant coefficients[34]

d2ydx2+bdydx+cy=r(x){\displaystyle {\frac {d^{2}y}{dx^{2}}}+b{\frac {dy}{dx}}+cy=r(x)}

Complementary functionyc{\displaystyle y_{c}}: assumeyc(x)=eαx{\displaystyle y_{c}(x)=e^{\alpha x}}, substitute and solve polynomial inα{\displaystyle \alpha }, to find thelinearly independent functionseαjx{\displaystyle e^{\alpha _{j}x}}.

Particular integralyp{\displaystyle y_{p}}: in general themethod of variation of parameters, though for very simpler(x){\displaystyle r(x)} inspection may work.[30]

y=yc+yp{\displaystyle y=y_{c}+y_{p}}

Ifb2>4c{\displaystyle b^{2}>4c}, then

yc=C1ex2(b+b24c)+C2ex2(bb24c){\displaystyle y_{c}=C_{1}e^{-{\frac {x}{2}}\,\left(b+{\sqrt {b^{2}-4c}}\right)}+C_{2}e^{-{\frac {x}{2}}\,\left(b-{\sqrt {b^{2}-4c}}\right)}}

Ifb2=4c{\displaystyle b^{2}=4c}, then

yc=(C1x+C2)ebx2{\displaystyle y_{c}=(C_{1}x+C_{2})e^{-{\frac {bx}{2}}}}

Ifb2<4c{\displaystyle b^{2}<4c}, then

yc=ebx2[C1sin(x4cb22)+C2cos(x4cb22)]{\displaystyle y_{c}=e^{-{\frac {bx}{2}}}\left[C_{1}\sin \left(x\,{\frac {\sqrt {4c-b^{2}}}{2}}\right)+C_{2}\cos \left(x\,{\frac {\sqrt {4c-b^{2}}}{2}}\right)\right]}

n{\displaystyle n}th-order, linear, inhomogeneous, constant coefficients[34]

j=0nbjdjydxj=r(x){\displaystyle \sum _{j=0}^{n}b_{j}{\frac {d^{j}y}{dx^{j}}}=r(x)}

Complementary functionyc{\displaystyle y_{c}}: assumeyc(x)=eαx{\displaystyle y_{c}(x)=e^{\alpha x}}, substitute and solve polynomial inα{\displaystyle \alpha }, to find thelinearly independent functionseαjx{\displaystyle e^{\alpha _{j}x}}.

Particular integralyp{\displaystyle y_{p}}: in general themethod of variation of parameters, though for very simpler(x){\displaystyle r(x)} inspection may work.[30]

y=yc+yp{\displaystyle y=y_{c}+y_{p}}

Sinceαj{\displaystyle \alpha _{j}} are the solutions of thepolynomial ofdegreen{\displaystyle n}:j=1n(ααj)=0{\textstyle \prod _{j=1}^{n}(\alpha -\alpha _{j})=0}, then:forαj{\displaystyle \alpha _{j}} all different,yc=j=1nCjeαjx{\displaystyle y_{c}=\sum _{j=1}^{n}C_{j}e^{\alpha _{j}x}}for each rootαj{\displaystyle \alpha _{j}} repeatedkj{\displaystyle k_{j}} times,yc=j=1n(=1kjCj,x1)eαjx{\displaystyle y_{c}=\sum _{j=1}^{n}\left(\sum _{\ell =1}^{k_{j}}C_{j,\ell }x^{\ell -1}\right)e^{\alpha _{j}x}}for someαj{\displaystyle \alpha _{j}} complex, then settingαj=χj+iγj{\displaystyle \alpha _{j}=\chi _{j}+i\gamma _{j}}, and usingEuler's formula, allows some terms in the previous results to be written in the formCjeαjx=Cjeχjxcos(γjx+φj){\displaystyle C_{j}e^{\alpha _{j}x}=C_{j}e^{\chi _{j}x}\cos(\gamma _{j}x+\varphi _{j})}whereφj{\displaystyle \varphi _{j}} is an arbitrary constant (phase shift).

Guessing solutions

[edit]
See also:Method of undetermined coefficients
This sectiondoes notcite anysources. Please helpimprove this section byadding citations to reliable sources. Unsourced material may be challenged andremoved.(January 2020) (Learn how and when to remove this message)

When all other methods for solving an ODE fail, or in the cases where we have some intuition about what the solution to a DE might look like, it is sometimes possible to solve a DE simply by guessing the solution and validating it is correct. To use this method, we simply guess a solution to the differential equation, and then plug the solution into the differential equation to verify whether it satisfies the equation. If it does then we have a particular solution to the DE, otherwise we start over again and try another guess. For instance we could guess that the solution to a DE has the form:y=Aeiαt{\displaystyle y=Ae^{i\alpha t}} since this is a very common solution that physically behaves in a sinusoidal way.

In the case of a first order ODE that is non-homogeneous we need to first find a solution to the homogeneous portion of the DE, otherwise known as the associated homogeneous equation, and then find a solution to the entire non-homogeneous equation by guessing. Finally, we add both of these solutions together to obtain the general solution to the ODE, that is:

general solution=general solution of the associated homogeneous equation+particular solution{\displaystyle {\text{general solution}}={\text{general solution of the associated homogeneous equation}}+{\text{particular solution}}}

Software for ODE solving

[edit]
  • Maxima, an open-sourcecomputer algebra system.
  • COPASI, a free (Artistic License 2.0) software package for the integration and analysis of ODEs.
  • MATLAB, a technical computing application (MATrix LABoratory)
  • GNU Octave, a high-level language, primarily intended for numerical computations.
  • Scilab, an open source application for numerical computation.
  • Maple, a proprietary application for symbolic calculations.
  • Mathematica, a proprietary application primarily intended for symbolic calculations.
  • SymPy, a Python package that can solve ODEs symbolically
  • Julia (programming language), a high-level language primarily intended for numerical computations.
  • SageMath, an open-source application that uses a Python-like syntax with a wide range of capabilities spanning several branches of mathematics.
  • SciPy, a Python package that includes an ODE integration module.
  • Chebfun, an open-source package, written inMATLAB, for computing with functions to 15-digit accuracy.
  • GNU R, an open source computational environment primarily intended for statistics, which includes packages for ODE solving.

See also

[edit]

Notes

[edit]
  1. ^Dennis G. Zill (15 March 2012).A First Course in Differential Equations with Modeling Applications. Cengage Learning.ISBN 978-1-285-40110-2.Archived from the original on 17 January 2020. Retrieved11 July 2019.
  2. ^"What is the origin of the term "ordinary differential equations"?".hsm.stackexchange.com.Stack Exchange. Retrieved2016-07-28.
  3. ^Karras, Tero; Aittala, Miika; Aila, Timo; Laine, Samuli (2022). "Elucidating the Design Space of Diffusion-Based Generative Models".arXiv:2206.00364 [cs.CV].
  4. ^Butcher, J. C. (2000-12-15)."Numerical methods for ordinary differential equations in the 20th century".Journal of Computational and Applied Mathematics. Numerical Analysis 2000. Vol. VI: Ordinary Differential Equations and Integral Equations.125 (1):1–29.Bibcode:2000JCoAM.125....1B.doi:10.1016/S0377-0427(00)00455-6.ISSN 0377-0427.
  5. ^Greenberg, Michael D. (2012).Ordinary differential equations. Hoboken, N.J: Wiley.ISBN 978-1-118-23002-2.
  6. ^Acton, Forman S. (1990).Numerical methods that work. Spectrum. Washington, D.C: Mathematical Association of America.ISBN 978-1-4704-5727-3.
  7. ^Denis, Byakatonda (2020-12-10). "An Overview of Numerical and Analytical Methods for solving Ordinary Differential Equations".arXiv:2012.07558 [math.HO].
  8. ^Mathematics for Chemists, D.M. Hirst,Macmillan Press, 1976, (No ISBN) SBN: 333-18172-7
  9. ^Kreyszig (1972, p. 64)
  10. ^Simmons (1972, pp. 1, 2)
  11. ^Halliday & Resnick (1977, p. 78)
  12. ^Tipler (1991, pp. 78–83)
  13. ^abHarper (1976, p. 127)
  14. ^Kreyszig (1972, p. 2)
  15. ^Simmons (1972, p. 3)
  16. ^abKreyszig (1972, p. 24)
  17. ^Simmons (1972, p. 47)
  18. ^Harper (1976, p. 128)
  19. ^Kreyszig (1972, p. 12)
  20. ^Ascher & Petzold (1998, p. 12)
  21. ^Achim Ilchmann; Timo Reis (2014).Surveys in Differential-Algebraic Equations II. Springer. pp. 104–105.ISBN 978-3-319-11050-9.
  22. ^Ascher & Petzold (1998, p. 5)
  23. ^Kreyszig (1972, p. 78)
  24. ^Kreyszig (1972, p. 4)
  25. ^Vardia T. Haimo (1985). "Finite Time Differential Equations".1985 24th IEEE Conference on Decision and Control. pp. 1729–1733.doi:10.1109/CDC.1985.268832.S2CID 45426376.
  26. ^Crelle, 1866, 1868
  27. ^Dresner (1999, p. 9)
  28. ^Logan, J. (2013). Applied mathematics (4th ed.).
  29. ^Ascher & Petzold (1998, p. 13)
  30. ^abcdefghijElementary Differential Equations and Boundary Value Problems (4th Edition), W.E. Boyce, R.C. Diprima, Wiley International, John Wiley & Sons, 1986,ISBN 0-471-83824-1
  31. ^Boscain; Chitour 2011, p. 21
  32. ^abMathematical Handbook of Formulas and Tables (3rd edition), S. Lipschutz, M. R. Spiegel, J. Liu, Schaum's Outline Series, 2009, ISC_2N 978-0-07-154855-7
  33. ^Further Elementary Analysis, R. Porter, G.Bell & Sons (London), 1978,ISBN 0-7135-1594-5
  34. ^abMathematical methods for physics and engineering, K.F. Riley, M.P. Hobson, S.J. Bence, Cambridge University Press, 2010, ISC_2N 978-0-521-86153-3

References

[edit]

Bibliography

[edit]

External links

[edit]
Wikibooks has a book on the topic of:Calculus/Ordinary differential equations
Wikimedia Commons has media related toOrdinary differential equations.
Classification
Operations
Attributes of variables
Relation to processes
Solutions
Existence/uniqueness
Solution topics
Solution methods
Examples
Mathematicians
Major topics inmathematical analysis
Authority control databases: NationalEdit this at Wikidata
Retrieved from "https://en.wikipedia.org/w/index.php?title=Ordinary_differential_equation&oldid=1317664412"
Categories:
Hidden categories:

[8]ページ先頭

©2009-2025 Movatter.jp