微分方程的前世今生

q1338 2020-02-08

展开全文

微分方程（英语：Differential equation，DE）是一种数学方程，用来描述某一类函数与其导数之间的关系。微分方程的解是一个符合方程的函数。而在初等数学的代数方程里，其解是常数值。

微分方程的应用十分广泛，可以解决许多与导数有关的问题。物理中许多涉及变力的运动学、动力学问题，如空气的阻力为速度函数的落体运动等问题，很多可以用微分方程求解。此外，微分方程在化学、工程学、经济学和人口统计等领域都有应用。

数学领域对微分方程的研究着重在几个不同的面向，但大多数都是关心微分方程的解。只有少数简单的微分方程可以求得解析解。不过即使没有找到其解析解，仍然可以确认其解的部分性质。在无法求得解析解时，可以利用数值分析的方式，利用电脑来找到其数值解。动力系统理论强调对于微分方程系统的量化分析，而许多数值方法可以计算微分方程的数值解，且有一定的准确度。

分类

微分方程可分为以下几类，而随着微分方程种类的不同，其相关研究的方式也会随之不同。

常微分方程及偏微分方程

常微分方程（ODE）是指一微分方程的未知数是单一自变量的函数。最简单的常微分方程，未知数是一个实数或是复数的函数，但未知数也可能是一个向量函数或是矩阵函数，后者可对应一个由常微分方程组成的系统。微分方程的表达通式是：

{\displaystyle f\left(x,{\frac {d^{n}y}{dx^{n}}},{\frac {d^{(n-1)}y}{dx^{(n-1)}}},\cdots ,{\frac {dy}{dx}},y\right)=0}f\left(x,{\frac {d^{n}y}{dx^{n}}},{\frac {d^{{(n-1)}}y}{dx^{{(n-1)}}}},\cdots ,{\frac {dy}{dx}},y\right)=0

常微分方程常依其阶数分类，阶数是指自变量导数的最高阶数[1]:p.3，最常见的二种为一阶微分方程及二阶微分方程。例如以下的贝塞尔方程：

{\displaystyle x^{2}{\frac {d^{2}y}{dx^{2}}}+x{\frac {dy}{dx}}+(x^{2}-\alpha ^{2})y=0}x^{2}{\frac {d^{2}y}{dx^{2}}}+x{\frac {dy}{dx}}+(x^{2}-\alpha ^{2})y=0

（其中y为因变量）为二阶微分方程，其解为贝塞尔函数。

偏微分方程（PDE）是指一微分方程的未知数是多个自变量的函数[2]，且方程中有未知数对自变量的偏微分。偏微分方程的阶数定义类似常微分方程，但更细分为椭圆型、双曲线型及抛物线型的偏微分方程，尤其在二阶偏微分方程中上述的分类更是重要。有些偏微分方程在整个自变量的值域中无法归类在上述任何一种型式中，这种偏微分方程则称为混合型。像以下的方程就是偏微分方程：

{\displaystyle {\frac {\partial u}{\partial t}}+t{\frac {\partial u}{\partial x}}=0.}{\frac {\partial u}{\partial t}}+t{\frac {\partial u}{\partial x}}=0.

线性及非线性编辑

常微分方程及偏微分方程都可以分为线性及非线性二类。

若微分方程中没有出现应变数及其微分项的乘积，此微分方程为线性微分方程，否则即为非线性微分方程。

齐次线性微分方程是线性微分方程中更细的分类，微分方程的解乘上一系数或是与另一个解相加后的结果仍为微分方程的解。

若线性微分方程的系数均为常数，则为常系数线性微分方程。常系数线性微分方程可以利用拉氏转换转换为代数方程[1]:p.315-316，因此简化求解的过程。

针对非线性的微分方程，只有相当少数的方法可以求得微分方程的解析解，而且这些方法需要微分方程有特别的对称性。长时间时非线性微分方程可能会出现非常复杂的特性，也可能会有混沌现象。有关非线性微分方程的一些基本问题，例如解的存在性、唯一性及初始值非线性微分方程的适定性问题，以及边界值非线性微分方程都是相当难的问题，甚至针对特定非线性微分方程的上述基本问题都被视为是数学理论的一大突破。例如2000年提出的7个千禧年大奖难题中，其中一个是纳维-斯托克斯存在性与光滑性，都是探讨纳维－斯托克斯方程其解的数学性质[3]，截至2018年8月此问题仍尚未被证明。

线性微分方程常常用来近似非线性微分方程，不过只在特定的条件下才能近似。例如单摆的运动方程为非线性的微分方程，但在小角度时可以近似为线性的微分方程。

举例编辑

以下是常微分方程的一些例子，其中{\displaystyle u}u为未知的函数，自变量为{\displaystyle x}x，{\displaystyle c}c及{\displaystyle \omega }\omega 均为常数。

非齐次一阶常系数线性微分方程：

{\displaystyle {\frac {du}{dx}}=cu+x^{2}.}{\frac {du}{dx}}=cu+x^{2}.

齐次二阶线性微分方程：

{\displaystyle {\frac {d^{2}u}{dx^{2}}}-x{\frac {du}{dx}}+u=0.}{\frac {d^{2}u}{dx^{2}}}-x{\frac {du}{dx}}+u=0.

描述谐振子的齐次二阶常系数线性微分方程：

{\displaystyle {\frac {d^{2}u}{dx^{2}}}+\omega ^{2}u=0.}{\frac {d^{2}u}{dx^{2}}}+\omega ^{2}u=0.

非齐次一阶非线性微分方程：

{\displaystyle {\frac {du}{dx}}=u^{2}+1.}{\frac {du}{dx}}=u^{2}+1.

描述长度为{\displaystyle L}L的单摆的二阶非线性微分方程：

{\displaystyle L{\frac {d^{2}u}{dx^{2}}}+g\sin u=0.}L{\frac {d^{2}u}{dx^{2}}}+g\sin u=0.

以下是偏微分方程的一些例子，其中{\displaystyle u}u为未知的函数，自变量为{\displaystyle x}x及{\displaystyle t}t或者是{\displaystyle x}x及{\displaystyle y}y。

齐次一阶线性偏微分方程：

{\displaystyle {\frac {\partial u}{\partial t}}+t{\frac {\partial u}{\partial x}}=0.}{\frac {\partial u}{\partial t}}+t{\frac {\partial u}{\partial x}}=0.

拉普拉斯方程，是椭圆型的齐次二阶常系数线性偏微分方程：

{\displaystyle {\frac {\partial ^{2}u}{\partial x^{2}}}+{\frac {\partial ^{2}u}{\partial y^{2}}}=0.}{\frac {\partial ^{2}u}{\partial x^{2}}}+{\frac {\partial ^{2}u}{\partial y^{2}}}=0.

KdV方程，是三阶的非线性偏微分方程：

{\displaystyle {\frac {\partial u}{\partial t}}=6u{\frac {\partial u}{\partial x}}-{\frac {\partial ^{3}u}{\partial x^{3}}}.}{\frac {\partial u}{\partial t}}=6u{\frac {\partial u}{\partial x}}-{\frac {\partial ^{3}u}{\partial x^{3}}}.

性质

普遍性的数学描述

许多物理或是化学的基本定律都可以写成微分方程的形式。在生物学及经济学中，微分方程用来作为复杂系统的数学模型。微分方程的数学理论最早是和方程对应的科学领域一起出现，而微分方程的解就可以用在该领域中。不过有时二个截然不同的科学领域会形成相同的微分方程，此时微分方程对应的数学理论可以看到不同现象后面一致的原则。

例如考虑光和声音在空气中的传播，以及池塘水面上的波动，这些都可以用同一个二阶的偏微分方程来描述，此方程即为波动方程，因此可以将光和声音视为一种波，和水面上的水波有些类似之处。约瑟夫·傅里叶所发展的热传导理论，其统御方程是另一个二阶偏微分方程－热传导方程，扩散作用看似和热传导不同，但也适用同一个统御方程，而经济学中的布莱克-休斯方程也和热传导方程有关。

微分方程的解编辑

微分方程的解通常是一个函数表达式{\displaystyle y=f(x)\,}y=f(x)\,（含一个或多个待定常数，由初始条件确定）。例如：

{\displaystyle {\frac {dy}{dx}}=\sin x}{\frac {dy}{dx}}=\sin x，

的解是

{\displaystyle y=-\cos x+C}y=-\cos x+C，

其中{\displaystyle C}C是待定常数；

例如，如果知道

{\displaystyle y=f(\pi )=2}y=f(\pi )=2，

则可推出

{\displaystyle C=1}C=1，

而可知 {\displaystyle y=-\cos x+1}y=-\cos x+1，

简易微分方程的求解方法编辑

一阶线性常微分方程编辑

对于一阶线性常微分方程，常用的方法是常数变易法：

对于方程：{\displaystyle y'+p(x)y+q(x)=0}y'+p(x)y+q(x)=0

可知其通解：{\displaystyle y=C(x)e^{-\int p(x)\,dx}}y=C(x)e^{{-\int p(x)\,dx}}

然后将这个通解代回到原式中，即可求出{\displaystyle C(x)}C(x)的值

二阶常系数齐次常微分方程编辑

对于二阶常系数齐次常微分方程，常用方法是求出其特征方程的解

对于方程：{\displaystyle y''+py'+qy=0}y''+py'+qy=0

其特征方程：{\displaystyle r^{2}+pr+q=0}r^{2}+pr+q=0

根据其特征方程，判断根的分布情况，然后得到方程的通解： {\displaystyle y=c_{1}y_{1}+c_{2}y_{2}}y=c_{1}y_{1}+c_{2}y_{2}

一般的通解形式为

（在{\displaystyle r_{1}=r_{2}}{\displaystyle r_{1}=r_{2}}的情况下）:{\displaystyle y=(C_{1}+C_{2}x)e^{rx}}y=(C_{1}+C_{2}x)e^{{rx}}

（在{\displaystyle r_{1}\neq r_{2}}{\displaystyle r_{1}\neq r_{2}}的情况下）:{\displaystyle y=C_{1}e^{r_{1}x}+C_{2}e^{r_{2}x}}y=C_{1}e^{{r_{1}x}}+C_{2}e^{{r_{2}x}}

（在共轭复数根的情况下）：{\displaystyle y=e^{\alpha x}(C_{1}\cos(\beta x)+C_{2}\sin(\beta x))}y=e^{{\alpha x}}(C_{1}\cos(\beta x)+C_{2}\sin(\beta x))

约束条件编辑

微分方程的约束条件是指其解需符合的条件，依常微分方程及偏微分方程的不同，有不同的约束条件。

常微分方程常见的约束条件是函数在特定点的值，若是高阶的微分方程，会加上其各阶导数的值，有这类约束条件的常微分方程称为初值问题。

若是二阶的常微分方程，也可能会指定函数在二个特定点的值，此时的问题即为边界值问题。若边界条件指定二点数值，称为狄利克雷边界条件（第一类边值条件），此外也有指定二个特定点上导数的边界条件，称为诺伊曼边界条件（第二类边值条件）等。

偏微分方程常见的问题以边界值问题为主，不过边界条件则是指定一特定超曲面的值或导数需符定特定条件。

解的存在性及唯一性编辑

存在性是指给定一微分方程及约束条件，判断其解是否存在。唯一性是指在上述条件下，是否只存在一个解。

针对常微分方程的初值问题，皮亚诺存在性定理可判别解的存在性，柯西-利普希茨定理则可以判别解的存在性及唯一性。

针对偏微分方程，柯西-克瓦列夫斯基定理（英语：Cauchy–Kowalevski theorem）可以判别解的存在性及唯一性。皮亚诺存在性定理可以判断常微分方程初值问题的解是否存在。

历史编辑

微分方程的起源约在十七世纪末，为了解决物理及天文学问题而产生，大约和微积分的发展同时。惠更斯在1693年的《教师学报》中提到常微分方程，雅各布·白努利在1691年建立悬链线的微分方程，并求得其函数。微分方程在十八世纪中期成为一个独立的学科[4]，而微分方程也带动许多当时的科学发展，例如海王星的发现就和微分方程的分析有关[5]。

偏微分方程是由傅里叶开始的，他在1822年发表《热的解析理论》，提出热传导方程的偏微分方程，并且利用分离变量法求得级数解，并且开始有关傅里叶级数的研究。另外在十九世纪有关拉普拉斯方程的研究也是偏微分方程的重要发展。拉普拉斯和泊松都有许多的贡献，后来乔治·格林提出了相关格林函数及格林公式等概念，并带动斯托克斯、麦克斯韦及后来电磁学相关的研究。而流体力学的纳维-斯托克斯方程及弹性介质的柯西方程也是在十九世纪提出的偏微分方程。[5]。后来许多的理论都是以偏微分方程的形式出现，量子力学的基础方程薛定谔方程也是偏微分方程，广义相对论中的爱因斯坦重力场方程也有类似偏微分的协变导数。