在下文中, y {\displaystyle y} 是因变量,代表独立变量 x {\displaystyle x} 的未知函数 y = f ( x ) {\displaystyle y=f(x)} 。导数的记号因作者不同而异,也取决于具体任务中最便于使用的记号。其中,莱布尼茨记号 d y d x , d 2 y d x 2 , … , d n y d x n {\displaystyle {\frac {dy}{dx}},{\frac {d^{2}y}{dx^{2}}},\ldots ,{\frac {d^{n}y}{dx^{n}}}} 在求导 和积分 时更为便利;而拉格朗日记号 y ′ , y ″ , … , y ( n ) {\displaystyle y',y'',\ldots ,y^{(n)}} 在紧凑表示高阶导数时更有优势;牛顿记号 ( y ˙ , y ¨ , y . . . ) {\displaystyle ({\dot {y}},{\ddot {y}},{\overset {...}{y}})} 则常用于物理学中表示相对于时间的低阶导数。
一般定义 给定一个函数 F {\displaystyle F} ,它依赖于 x {\displaystyle x} 、y {\displaystyle y} 以及 y {\displaystyle y} 的导数。那么形如
F ( x , y , y ′ , … , y ( n − 1 ) ) = y ( n ) {\displaystyle F\left(x,y,y',\ldots ,y^{(n-1)}\right)=y^{(n)}} 的方程称为 n {\displaystyle n} 阶显式常微分方程 。
更一般地,n {\displaystyle n} 阶隐式常微分方程的形式为:
F ( x , y , y ′ , y ″ , … , y ( n ) ) = 0 {\displaystyle F\left(x,y,y',y'',\ \ldots ,\ y^{(n)}\right)=0} 进一步分类如下:
自治的 如果微分方程不显式依赖于变量 x ,则称其为自治(驻定)的。 线性的 如果 F {\displaystyle F} 可以写成 y {\displaystyle y} 的各阶导数的线性组合,即方程可以改写为 y ( n ) = ∑ i = 0 n − 1 a i ( x ) y ( i ) + r ( x ) {\displaystyle y^{(n)}=\sum _{i=0}^{n-1}a_{i}(x)y^{(i)}+r(x)} 其中 a i ( x ) {\displaystyle a_{i}(x)} 和 r ( x ) {\displaystyle r(x)} 是 x {\displaystyle x} 的连续函数 ,则称该微分方程为线性的。 函数 r ( x ) {\displaystyle r(x)} 称为源项(source term),从而引出进一步的分类。 齐次的 如果 r ( x ) = 0 {\displaystyle r(x)=0} ,则线性微分方程是齐次的。此时总存在“平凡解” y = 0 {\displaystyle y=0} 。 非齐次的 如果 r ( x ) ≠ 0 {\displaystyle r(x)\neq 0} ,则线性微分方程是非齐次的。 非线性的 不是线性的微分方程。
常微分方程组 多个相互耦合的微分方程构成一个方程组。如果 y {\displaystyle \mathbf {y} } 是一个向量,其分量为函数
y ( x ) = [ y 1 ( x ) , y 2 ( x ) , … , y m ( x ) ] {\displaystyle \mathbf {y} (x)=[y_{1}(x),y_{2}(x),\ldots ,y_{m}(x)]} 且 F {\displaystyle \mathbf {F} } 是 y {\displaystyle \mathbf {y} } 及其导数的向量值函数,则
y ( n ) = F ( x , y , y ′ , y ″ , … , y ( n − 1 ) ) {\displaystyle \mathbf {y} ^{(n)}=\mathbf {F} \left(x,\mathbf {y} ,\mathbf {y} ',\mathbf {y} '',\ldots ,\mathbf {y} ^{(n-1)}\right)} 称为 n {\displaystyle n} 阶 m {\displaystyle m} 维显式常微分方程组 。用列向量形式表示为:
( y 1 ( n ) y 2 ( n ) ⋮ y m ( n ) ) = ( f 1 ( x , y , y ′ , y ″ , … , y ( n − 1 ) ) f 2 ( x , y , y ′ , y ″ , … , y ( n − 1 ) ) ⋮ f m ( x , y , y ′ , y ″ , … , y ( n − 1 ) ) ) {\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}}} 这些方程组不一定是线性的。其隐式对应形式为:
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}}} 其中 0 = ( 0 , 0 , … , 0 ) {\displaystyle {\boldsymbol {0}}=(0,0,\ldots ,0)} 是零向量。用矩阵 形式表示为:
( f 1 ( x , y , y ′ , y ″ , … , y ( n ) ) f 2 ( x , y , y ′ , y ″ , … , y ( n ) ) ⋮ f m ( x , y , y ′ , y ″ , … , y ( n ) ) ) = ( 0 0 ⋮ 0 ) {\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}}} 对于形如 F ( x , y , y ′ ) = 0 {\displaystyle \mathbf {F} \left(x,\mathbf {y} ,\mathbf {y} '\right)={\boldsymbol {0}}} 的方程组,某些文献还要求雅可比矩阵 ∂ F ( x , u , v ) ∂ v {\displaystyle {\frac {\partial \mathbf {F} (x,\mathbf {u} ,\mathbf {v} )}{\partial \mathbf {v} }}} 非奇异 ,才将其称为隐式常微分方程组;满足这一雅可比非奇异条件的隐式方程组可以转化为显式方程组。在同一文献中,雅可比矩阵奇异的隐式方程组被称为微分代数方程 (DAE)。这一区分不仅仅是术语上的差异;微分代数方程具有本质上不同的特性,通常比(非奇异的)常微分方程组更难求解。
据推测,对于更高阶导数的情形,按照这一方案也假设黑塞矩陣等相应矩阵非奇异[來源請求] ,不过需要注意,任何高于一阶的常微分方程都可以(且通常会)改写为一阶常微分方程组,这使得雅可比奇异性判别准则足以全面覆盖所有阶数的分类。
常微分方程组的性态可以通过相图来进行可视化。
常微分方程的解 给定微分方程
F ( x , y , y ′ , … , y ( n ) ) = 0 {\displaystyle F\left(x,y,y',\ldots ,y^{(n)}\right)=0} 一个函数 u : I ⊂ R → R {\displaystyle u:I\subset \mathbb {R} \to \mathbb {R} } ,其中 I {\displaystyle I} 是区间,如果 u {\displaystyle u} 在 I {\displaystyle I} 上 n {\displaystyle n} 次可微 ,且
F ( x , u , u ′ , … , u ( n ) ) = 0 x ∈ I . {\displaystyle F(x,u,u',\ \ldots ,\ u^{(n)})=0\quad x\in I.} 则称 u {\displaystyle u} 为该方程的解 或积分曲线 。
给定两个解 u : J ⊂ R → R {\displaystyle u:J\subset \mathbb {R} \to \mathbb {R} } 和 v : I ⊂ R → R {\displaystyle v:I\subset \mathbb {R} \to \mathbb {R} } ,如果 I ⊂ J {\displaystyle I\subset J} 且
u ( x ) = v ( x ) x ∈ I . {\displaystyle u(x)=v(x)\quad x\in I.\,} 则称 u {\displaystyle u} 是 v {\displaystyle v} 的延拓 。
没有延拓的解称为极大解 (maximal solution)。在整个 R {\displaystyle \mathbb {R} } 上定义的解称为全局解 (global solution)。
n {\displaystyle n} 阶方程的通解 是指含有 n {\displaystyle n} 个任意独立積分常數的解。特解 是指通过给通解中的常数赋以特定值而得到的,通常选择这些值以满足给定的初值条件或边值条件。奇异解 是指不能通过给通解中的任意常数赋以特定值而得到的解。
解的存在唯一性 不止一个定理可在局部和整体意义下证实常微分方程初值问题解的存在唯一性,其中两个主要定理为皮亚诺存在性定理与柯西-利普希茨定理。这两个定理的基本形式均只保证局部结果,但后者可推广至整体结果,如满足格朗沃尔不等式的条件时。
此外,上述利普希茨型唯一性定理不适用于微分代数方程组 ,因其(非线性)代数部分本身可能导致多个解。
局部存在唯一性 该定理可简述如下:
对于方程及初值问题
y ′ = F ( x , y ) , y 0 = y ( x 0 ) {\displaystyle y'=F(x,y)\,,\quad y_{0}=y(x_{0})} 若 F {\displaystyle F} 及 ∂ F / ∂ y {\displaystyle \partial F/\partial y} 在 x − y {\displaystyle x-y} 平面上的闭矩形
R = [ x 0 − a , x 0 + a ] × [ y 0 − b , y 0 + b ] {\displaystyle R=[x_{0}-a,x_{0}+a]\times [y_{0}-b,y_{0}+b]} 内连续(a , b ∈ R {\displaystyle a,b\in \mathbb {R} } 且符号 × {\displaystyle \times } 表示笛卡尔积,方括号表示闭区间),则存在某个h ∈ R {\displaystyle h\in \mathbb {R} } ,使得区间
I = [ x 0 − h , x 0 + h ] ⊂ [ x 0 − a , x 0 + a ] {\displaystyle I=[x_{0}-h,x_{0}+h]\subset [x_{0}-a,x_{0}+a]} 上的上述方程及初值问题存在解,且该解唯一。由于未要求 F {\displaystyle F} 为线性函数,故此结论适用于形如 F ( x , y ) {\displaystyle F(x,y)} 的非线性方程,亦可应用于方程组。
解的全局唯一性和最大定义域 满足柯西-利普希茨定理的假设条件时,局部存在唯一性可以扩展至全局。更确切地说:
对于每个初值条件 ( x 0 , y 0 ) {\displaystyle (x_{0},y_{0})} ,存在唯一的最大(可能无限)开区间
I max = ( x − , x + ) , x ± ∈ R ∪ { ± ∞ } , x 0 ∈ I max {\displaystyle I_{\max }=(x_{-},x_{+}),x_{\pm }\in \mathbb {R} \cup \{\pm \infty \},x_{0}\in I_{\max }} 使得任何满足该初值条件的解,都是定义在 I max {\displaystyle I_{\max }} 上且满足该初值条件的解的限制。
在 x ± ≠ ± ∞ {\displaystyle x_{\pm }\neq \pm \infty } 的情形下,只有以下两种可能:
有限时间爆炸:lim sup x → x ± ‖ y ( x ) ‖ → ∞ {\displaystyle \limsup _{x\to x_{\pm }}\|y(x)\|\to \infty } 离开定义域:lim x → x ± y ( x ) ∈ ∂ Ω ¯ {\displaystyle \lim _{x\to x_{\pm }}y(x)\in \partial {\bar {\Omega }}} 其中 Ω {\displaystyle \Omega } 是 F {\displaystyle F} 定义的开集,∂ Ω ¯ {\displaystyle \partial {\bar {\Omega }}} 是其闭包的边界。
需要注意的是:
解的最大定义域总是开区间(以保证唯一性); 解的最大定义域可能严格小于 R {\displaystyle \mathbb {R} } ; 解的最大定义域可能依赖于具体的初值条件 ( x 0 , y 0 ) {\displaystyle (x_{0},y_{0})} 。
一些精确封闭形式的解 一些微分方程有精确封闭形式的解,这里给出几个重要的类型。
在下表中,P ( x ) , Q ( x ) ; P ( y ) , Q ( y ) {\displaystyle P(x),Q(x);P(y),Q(y)} 和M ( x , y ) , N ( x , y ) {\displaystyle M(x,y),N(x,y)} 是任意关于x , y {\displaystyle x,y} 的可积函数,b , c {\displaystyle b,c} 是给定的实常数,C , C 1 , C 2 … {\displaystyle C,C_{1},C_{2}\ldots } 是任意常数(一般为复数)。这些微分方程的等价或替代形式通过积分可以得到解。
在积分解中,λ {\displaystyle \lambda } 和ϵ {\displaystyle \epsilon } 是积分变量(求和下标的连续形式),记号∫ x F ( λ ) d λ {\displaystyle \int ^{x}F(\lambda )\mathrm {d} \lambda } 只表示F ( λ ) {\displaystyle F(\lambda )} 对λ {\displaystyle \lambda } 积分,在积分以后λ = x {\displaystyle \lambda {}=x} 替换,无需加常数(明确说明)。
微分方程 解法 通解 可分离微分方程 一阶,变量x {\displaystyle x} 和y {\displaystyle y} 均可分离(一般情况,下面有特殊情况)P 1 ( x ) Q 1 ( y ) + P 2 ( x ) Q 2 ( y ) d y d x = 0 {\displaystyle P_{1}(x)Q_{1}(y)+P_{2}(x)Q_{2}(y)\,{\frac {\mathrm {d} y}{\mathrm {d} x}}=0\,\!}
P 1 ( x ) Q 1 ( y ) d x + P 2 ( x ) Q 2 ( y ) d y = 0 {\displaystyle P_{1}(x)Q_{1}(y)\,\mathrm {d} x+P_{2}(x)Q_{2}(y)\,\mathrm {d} y=0\,\!}
分离变量(除以P 2 Q 1 {\displaystyle P_{2}Q_{1}} )。 ∫ x P 1 ( λ ) P 2 ( λ ) d λ + ∫ y Q 2 ( λ ) Q 1 ( λ ) d λ = C {\displaystyle \int ^{x}{\frac {P_{1}(\lambda )}{P_{2}(\lambda )}}\,\mathrm {d} \lambda +\int ^{y}{\frac {Q_{2}(\lambda )}{Q_{1}(\lambda )}}\,\mathrm {d} \lambda =C\,\!} 一阶,变量x {\displaystyle x} 可分离d y d x = F ( x ) {\displaystyle {\frac {\mathrm {d} y}{\mathrm {d} x}}=F(x)\,\!}
d y = F ( x ) d x {\displaystyle \mathrm {d} y=F(x)\,\mathrm {d} x\,\!}
直接积分。 y = ∫ x F ( λ ) d λ + C {\displaystyle y=\int ^{x}F(\lambda )\,\mathrm {d} \lambda +C\,\!} 一阶自治,变量y {\displaystyle y} 可分离d y d x = F ( y ) {\displaystyle {\frac {\mathrm {d} y}{\mathrm {d} x}}=F(y)\,\!}
d y = F ( y ) d x {\displaystyle \mathrm {d} y=F(y)\,\mathrm {d} x\,\!}
分离变量(除以F {\displaystyle F} )。 x = ∫ y d λ F ( λ ) + C {\displaystyle x=\int ^{y}{\frac {\mathrm {d} \lambda }{F(\lambda )}}+C\,\!} 一阶,变量x {\displaystyle x} 和y {\displaystyle y} 均可分离P ( y ) d y d x + Q ( x ) = 0 {\displaystyle P(y){\frac {\mathrm {d} y}{\mathrm {d} x}}+Q(x)=0\,\!}
P ( y ) d y + Q ( x ) d x = 0 {\displaystyle P(y)\,\mathrm {d} y+Q(x)\,\mathrm {d} x=0\,\!}
整个积分。 ∫ y P ( λ ) d λ + ∫ x Q ( λ ) d λ = C {\displaystyle \int ^{y}P(\lambda )\,{\mathrm {d} \lambda }+\int ^{x}Q(\lambda )\,\mathrm {d} \lambda =C\,\!} 一般一阶微分方程 一阶,齐次d y d x = F ( y x ) {\displaystyle {\frac {\mathrm {d} y}{\mathrm {d} x}}=F\left({\frac {y}{x}}\right)\,\!}
令y = u x {\displaystyle y=ux} ,然后通过分离变量u {\displaystyle u} 和x {\displaystyle x} 求解。 ln ( C x ) = ∫ y x d λ F ( λ ) − λ {\displaystyle \ln(Cx)=\int ^{\frac {y}{x}}{\frac {\mathrm {d} \lambda }{F(\lambda )-\lambda }}\,\!} 一阶,可分离变量y M ( x y ) + x N ( x y ) d y d x = 0 {\displaystyle yM(xy)+xN(xy)\,{\frac {\mathrm {d} y}{\mathrm {d} x}}=0\,\!}
y M ( x y ) d x + x N ( x y ) d y = 0 {\displaystyle yM(xy)\,\mathrm {d} x+xN(xy)\,\mathrm {d} y=0\,\!}
分离变量(除以x y {\displaystyle xy} )。 ln ( C x ) = ∫ x y N ( λ ) d λ λ [ N ( λ ) − M ( λ ) ] {\displaystyle \ln(Cx)=\int ^{xy}{\frac {N(\lambda )\,\mathrm {d} \lambda }{\lambda [N(\lambda )-M(\lambda )]}}\,\!}
如果N = M {\displaystyle N=M} ,解为x y = C {\displaystyle xy=C} 。
正合微分,一阶M ( x , y ) d y d x + N ( x , y ) = 0 {\displaystyle M(x,y){\frac {\mathrm {d} y}{\mathrm {d} x}}+N(x,y)=0\,\!}
M ( x , y ) d y + N ( x , y ) d x = 0 {\displaystyle M(x,y)\,\mathrm {d} y+N(x,y)\,\mathrm {d} x=0\,\!}
其中∂ M ∂ x = ∂ N ∂ y {\displaystyle {\frac {\partial M}{\partial x}}={\frac {\partial N}{\partial y}}\,\!}
全部積分 F ( x , y ) = ∫ y M ( x , λ ) d λ + ∫ x N ( λ , y ) d λ + Y ( y ) + X ( x ) = C {\displaystyle {\begin{aligned}F(x,y)&=\int ^{y}M(x,\lambda )\,\mathrm {d} \lambda +\int ^{x}N(\lambda ,y)\,\mathrm {d} \lambda \\&+Y(y)+X(x)=C\end{aligned}}\,\!} 其中Y ( y ) {\displaystyle Y(y)} 和X ( x ) {\displaystyle X(x)} 是积分出来的函数而不是常数,将它们列在这里以使最终函数F ( x , y ) {\displaystyle F(x,y)} 满足初始条件。
非正合微分 ,一阶M ( x , y ) d y d x + N ( x , y ) = 0 {\displaystyle M(x,y){\frac {\mathrm {d} y}{\mathrm {d} x}}+N(x,y)=0\,\!}
M ( x , y ) d y + N ( x , y ) d x = 0 {\displaystyle M(x,y)\,\mathrm {d} y+N(x,y)\,\mathrm {d} x=0\,\!}
其中∂ M ∂ x ≠ ∂ N ∂ y {\displaystyle {\frac {\partial M}{\partial x}}\neq {\frac {\partial N}{\partial y}}\,\!}
积分因子μ ( x , y ) {\displaystyle \mu (x,y)} 满足 ∂ ( μ M ) ∂ x = ∂ ( μ N ) ∂ y {\displaystyle {\frac {\partial (\mu M)}{\partial x}}={\frac {\partial (\mu N)}{\partial y}}\,\!}
如果可以得到μ ( x , y ) {\displaystyle \mu (x,y)} : F ( x , y ) = ∫ y μ ( x , λ ) M ( x , λ ) d λ + ∫ x μ ( λ , y ) N ( λ , y ) d λ + Y ( y ) + X ( x ) = C {\displaystyle {\begin{aligned}F(x,y)&=\int ^{y}\mu (x,\lambda )M(x,\lambda )\,\mathrm {d} \lambda +\int ^{x}\mu (\lambda ,y)N(\lambda ,y)\,\mathrm {d} \lambda \\&+Y(y)+X(x)=C\\\end{aligned}}\,\!}
一般二阶微分方程 二阶,自治d 2 y d x 2 = F ( y ) {\displaystyle {\frac {\mathrm {d} ^{2}y}{\mathrm {d} x^{2}}}=F(y)\,\!}
原方程乘以2 d y d x {\displaystyle {\frac {2\mathrm {d} y}{\mathrm {d} x}}} ,代换2 d y d x d 2 y d x 2 = d d x ( d y d x ) 2 {\displaystyle 2{\frac {\mathrm {d} y}{\mathrm {d} x}}{\frac {\mathrm {d} ^{2}y}{dx^{2}}}={\frac {\mathrm {d} }{\mathrm {d} x}}\left({\frac {\mathrm {d} y}{\mathrm {d} x}}\right)^{2}\,\!} ,然后两次积分。 x = ± ∫ y d λ 2 ∫ λ F ( ϵ ) d ϵ + C 1 + C 2 {\displaystyle x=\pm \int ^{y}{\frac {\mathrm {d} \lambda }{\sqrt {2\int ^{\lambda }F(\epsilon )\,\mathrm {d} \epsilon +C_{1}}}}+C_{2}\,\!} 线性微分方程(最高到n {\displaystyle n} 阶) 一阶线性,非齐次的函数系数d y d x + P ( x ) y = Q ( x ) {\displaystyle {\frac {\mathrm {d} y}{\mathrm {d} x}}+P(x)y=Q(x)\,\!}
积分因子:e ∫ x P ( λ ) d λ {\displaystyle e^{\int ^{x}P(\lambda )\,d\lambda }} 。 y = e − ∫ x P ( λ ) d λ [ ∫ x e ∫ λ P ( ϵ ) d ϵ Q ( λ ) d λ + C ] {\displaystyle y=e^{-\int ^{x}P(\lambda )\,\mathrm {d} \lambda }\left[\int ^{x}e^{\int ^{\lambda }P(\epsilon )\,\mathrm {d} \epsilon }Q(\lambda )\,{\mathrm {d} \lambda }+C\right]} 二阶线性,非齐次的常系数d 2 y d x 2 + b d y d x + c y = r ( x ) {\displaystyle {\frac {\mathrm {d} ^{2}y}{\mathrm {d} x^{2}}}+b{\frac {\mathrm {d} y}{\mathrm {d} x}}+cy=r(x)\,\!}
余函数y c {\displaystyle y_{c}} :设y c = e α x {\displaystyle y_{c}=\mathrm {e} ^{\alpha x}} ,代换并解出α {\displaystyle \alpha } 中的多项式,求出线性无关函数e α j x {\displaystyle e^{\alpha _{j}x}} 。 特解y p {\displaystyle y_{p}} :一般运用常数变易法 ,虽然对于非常容易的r ( x ) {\displaystyle r(x)} 可以直观判断。
y = y c + y p {\displaystyle y=y_{c}+y_{p}} 如果b 2 > 4 c {\displaystyle b^{2}>4c} ,则:
y c = C 1 e ( − b + b 2 − 4 c ) x 2 + C 2 e − ( b + b 2 − 4 c ) x 2 {\displaystyle y_{c}=C_{1}e^{\left(-b+{\sqrt {b^{2}-4c}}\right){\frac {x}{2}}}+C_{2}e^{-\left(b+{\sqrt {b^{2}-4c}}\right){\frac {x}{2}}}\,\!}
如果b 2 = 4 c {\displaystyle b^{2}=4c} ,则:
y c = ( C 1 x + C 2 ) e − b x 2 {\displaystyle y_{c}=(C_{1}x+C_{2})e^{-{\frac {bx}{2}}}\,\!}
如果b 2 < 4 c {\displaystyle b^{2}<4c} ,则:
y c = e − b x 2 [ C 1 sin ( | b 2 − 4 c | x 2 ) + C 2 cos ( | b 2 − 4 c | x 2 ) ] {\displaystyle y_{c}=e^{-{\frac {bx}{2}}}\left[C_{1}\sin {\left({\sqrt {\left|b^{2}-4c\right|}}{\frac {x}{2}}\right)}+C_{2}\cos {\left({\sqrt {\left|b^{2}-4c\right|}}{\frac {x}{2}}\right)}\right]\,\!}
n {\displaystyle n} 阶线性,非齐次常系数∑ j = 0 n b j d j y d x j = r ( x ) {\displaystyle \sum _{j=0}^{n}b_{j}{\frac {\mathrm {d} ^{j}y}{\mathrm {d} x^{j}}}=r(x)\,\!}
余函数y c {\displaystyle y_{c}} :设y c = e α x {\displaystyle y_{c}=\mathrm {e} ^{\alpha x}} ,代换并解出α {\displaystyle \alpha } 中的多项式,求出线性无关函数e α j x {\displaystyle e^{\alpha _{j}x}} 。 特解y p {\displaystyle y_{p}} :一般运用常数变易法 ,虽然对于非常容易的r ( x ) {\displaystyle r(x)} 可以直观判断。
y = y c + y p {\displaystyle y=y_{c}+y_{p}} 由于α j {\displaystyle \alpha _{j}} 为n {\displaystyle n} 阶多项式 的解: ∏ j = 1 n ( α − α j ) = 0 {\displaystyle \prod _{j=1}^{n}\left(\alpha -\alpha _{j}\right)=0\,\!} ,于是:
对于各不相同的α j {\displaystyle \alpha _{j}} ,
y c = ∑ j = 1 n C j e α j x {\displaystyle y_{c}=\sum _{j=1}^{n}C_{j}e^{\alpha _{j}x}\,\!}
每个根α j {\displaystyle \alpha _{j}} 重复k j {\displaystyle k_{j}} 次,
y c = ∑ j = 1 n ( ∑ ℓ = 1 k j C ℓ x ℓ − 1 ) e α j x {\displaystyle y_{c}=\sum _{j=1}^{n}\left(\sum _{\ell =1}^{k_{j}}C_{\ell }x^{\ell -1}\right)e^{\alpha _{j}x}\,\!}
对于一些复数值的αj ,令α = χj + i γj ,使用欧拉公式,前面结果中的一些项就可以写成
C j e α j x = C j e χ j x cos ( γ j x + ϕ j ) {\displaystyle C_{j}e^{\alpha _{j}x}=C_{j}e^{\chi _{j}x}\cos(\gamma _{j}x+\phi _{j})\,\!} 的形式,其中ϕj 为任意常量(相移)。
有限持续时间(Finite Duration)的解 对于非线性自治常微分方程,在某些条件下可以构造有限持续时间的解(Solutions of Finite Duration),这里指的是系统在自身动力学的驱动下,在某个有限的终止时刻到达零值,并在之后永远停留在零。
这些有限持续时间的解不可能是整个实直线上的解析函数,而且由于它们在终止时刻处是非利普希茨函数,因此不被利普希茨型微分方程解的唯一性定理所涵盖。
举个例子,方程
y ′ = − sgn ( y ) | y | , y ( 0 ) = 1 {\displaystyle y'=-{\text{sgn}}(y){\sqrt {|y|}},\,\,y(0)=1} 具有有限持续时间的解
y ( x ) = 1 4 ( 1 − x 2 + | 1 − x 2 | ) 2 {\displaystyle y(x)={\frac {1}{4}}\left(1-{\frac {x}{2}}+\left|1-{\frac {x}{2}}\right|\right)^{2}} 该解在 x > 2 {\displaystyle x>2} 时取值为零,并在 x = 2 {\displaystyle x=2} 处平滑连接(右导数为零),但右端函数在 y = 0 {\displaystyle y=0} 处不满足利普希茨条件,从而允许这种“有限时间熄灭”并停留的行为。