（低阶）偏微分方程的分类与特征线

                     

贡献者： int256; addis

　　 对于低阶偏微分方程而言，我们可以找到一族曲线，偏微分方程的初值信息沿着这条曲线传播（体现为在这条曲线上值不变或衰减），这条曲线就被称为特征线，而由于解沿这曲线传递，我们就可以通过解常微分方程的方法求出这条曲线（将在下面的例子中看到这一点）。最直白的说，就是偏微分方程这个系统随时间的演化，沿着特征线将变为 “常微分方程”。

　　 我们可以假想初值是一物理系统的初始状态，从初值情况开始，解将进行传播并产生波，描述这个 “波” 传递的过程就是特征线法，故又称行波法。

　　 同时，特征线法为二阶 PDE 提供了一个分类方法。

　　 本文默认 PDE 的范围是 $-\mathrm{\infty }<x,y,\cdots <+\mathrm{\infty },t>0$，用 $u_t^{\prime }$ 表示 $\frac{\partial u}{\partial t}$，用 $u_x^{\prime }$ 表示 $\frac{\partial u}{\partial x}$。默认二阶导连续，即 $u_{xy}^{″}=u_{yx}^{″}=\frac{\partial }{\partial y}\left(\frac{\partial u}{\partial x}\right)=\frac{\partial }{\partial x}\left(\frac{\partial u}{\partial y}\right)$。

　　 特征线有如下性质，

1. 对于 $n$ 维（这里包含时间维）的 PDE，其特征线总是 $n-1$ 维的。
2. 解的 “间断” 性质通过、且仅能通过特征线传播。（也就是说，不同时刻，同一特征线上的解相同，但有时数值会沿特征线衰减。）
3. 特征线上 PDE 的解是一样的（或衰减），特征线的斜率是解 “传播” 的速度。

　　 其中性质 $2$ 对应着如果不存在特征线（特征线不是实的，类似一元二次方程没有实数根称为根不存在），那么 PDE 的解连续（这就是 elliptic PDE）。

　　 同时求出特征线还有助于解决诸如影响区域、依赖区域一类的问题（这是由解信息沿特征线传播决定的）。

1. 特征线

　　 考虑一个一阶拟线性方程的柯西问题，一般形式为

　　 这有一个直观的几何解释：在 $(x,t,u)$ 的三维空间中，曲面上任意一点 $P_0(x_0,t_0,u_0)$ 的法向量为 $\overrightarrow{n}=(u_x,u_t,-1){|}_{P_0}$。在曲面 $u=u(x,t)$ 上，曲线

在点 $P_0$ 处的切向量为 显然，$\overrightarrow{n}$ 与 $\overrightarrow{\tau }$ 在 $P_0$ 互相垂直。

　　 当考虑矢量 $\overrightarrow{\mu }=(a(x,t,u),b(x,t,u),c(x,t,u)){|}_{P_0}$，则微分方程恰好表示法向量 $\overrightarrow{n}$ 与 $\overrightarrow{\mu }$ 在 $P_0$ 点处互相垂直，从而在曲面 $u=u(x,t)$ 上以 $\overrightarrow{\mu }$ 为切矢量且与曲线 $\mathrm{\Gamma }$ 相交，交点是 $(\alpha ,0,\phi (\alpha ))$ 的曲线方程为：

　　 其中 $\alpha$ 是曲线的参数，这组微分方程的解就被称为特征线。这条曲线在 $(x,t)$ 平面的投影就是一般意义的特征线，给出原微分方程的解曲面。

2. 一维一阶 PDE

　　 特征线法为一维一阶 PDE 提供了一个解法：先求出特征线，后代入求出初值在特征线上的传播情况即可。

　　 下面利用两个一维一阶的 PDE 举例来说明如何求解 PDE 的特征线。

　　 另一个方法则是考虑特征线下的坐标系 $(s,\tau )$，其中 $\tau$ 仍是某个常数，由选择的 “哪条” 特征线而定，而 $s\ge 0$ 是代表在某条特征线上的位置。

　　 同时，对于二元二阶线性的 PDE，其特征线是 2 维的，并且我们将会看到其体现出圆锥曲线的性质，同时各类圆锥曲线的特征线对应的 PDE 有不同的性质。下面我们来讨论特征线与 PDE 的分类。

3. 二阶线性 PDE 的分类

　　 在一般的数理方程或介绍 PDE 的书籍中，会把二阶线性 PDE 将之主要分为以下三类，分别是

1. 椭圆类（elliptic PDE），例如泊松方程 ${\mathbf{\nabla }}^{2}u=f$；
2. 抛物类（parabolic PDE），例如热传导方程 $k{\mathbf{\nabla }}^{2}T=\frac{\partial T}{\partial t}$；
3. 双曲类（hyperbolic PDE），例如波动方程 ${\mathbf{\nabla }}^{2}w=\frac{1}{c^2}\frac{{\partial }^{2}w}{\partial t^2}$。

　　 这三类的解信息沿特征线传播，有性质与对应的物理问题：

1. 椭圆类：解的信息以无限速度传播，且解不应当存在间断。适用于平衡问题（静态或准静态，例如无源力场）。
2. 抛物类：解的信息以无限速度传播，适用于扩散问题（例如热传导）。
3. 双曲类：解的信息以有限速度传播，且能保留间断信息，适用于波动问题。

　　 另外还有一类 PDE 被称为 ultrahyperbolic PDE（超双曲类）。下面将讲解是如何分类的。

二元二阶线性 PDE

　　 首先保持二阶导连续的假设，即 $\frac{\partial }{\partial x}\frac{\partial u}{\partial y}=\frac{\partial }{\partial y}\frac{\partial u}{\partial x}$。则偏微分方程具有一般形式： $$A{u}_{xx}^{″}+B{u}_{xy}^{″}+C{u}_{yy}^{″}+D{u}_x^{\prime }+E{u}_y^{\prime }+F=0,$$ 与一阶线性 PDE 对比，可以得到特征线应当满足方程： $$\{\begin{array}{rl}\mathrm{d}(u_x^{\prime })& =u_{xx}^{″}\mathrm{d}x+u_{xy}^{″}\mathrm{d}y,\\ \mathrm{d}(u_y^{\prime })& =u_{xy}^{″}\mathrm{d}x+u_{yy}^{″}\mathrm{d}y.\end{array}$$ 联立这三式，令有三个矩阵：

则这三式可转化为 $\mathbf{Mp}=\mathbf{N}$。 显然，$u$ 的各个二阶偏导项不是唯一确定的，所以 $det\mathbf{M}=0$，也就是需要： 也就是： 发现是一个关于 $\mathrm{d}y/\mathrm{d}x$ 的二次方程，判别式 $\mathrm{\Delta }=B^2-4AC$。对比一般利用二次表达式的圆锥曲线定义，可以将特征线分类如下。

1. $\mathrm{\Delta }<0$，存在虚特征线（即无特征线），称这种为椭圆类 PDE；
2. $\mathrm{\Delta }=0$，存在一族特征线，称这种为抛物类 PDE；
3. $\mathrm{\Delta }>0$，存在两族不同的特征线，称这种为双曲类 PDE。

　　 这就是分类的依据。

　　 值得注意的是，PDE 可能在不同区域有不同分类。例如 $y{u}_{xx}^{″}-u_{yy}^{″}=0$，在 $y>0$ 时是双曲类；在 $y=0$ 时是抛物类，在 $y<0$ 时是椭圆类。

多元二阶线性 PDE

　　 考虑线性微分算符 $\hat{L}$：

忽略低阶项（一阶的 ${\mathbf{B}}_i$ 与 $F$），考虑系数矩阵 ${\mathbf{A}}_{i,j}$ 的特征值，我们将 $\hat{L}u=0$ 分类如下。

1. 当特征值均大于 $0$ 或均小于 $0$，即特征值全部同号且都非零，归类为椭圆类；
2. 当特征值有一个为 $0$，其余均大于 $0$ 或均小于 $0$，即特征值除有一 $0$ 外均同号，归类为抛物类；
3. 当特征值一个为正数，其他为负数；或一个为负数，其他为正数。即特征值均非 $0$，并且有仅有一非 $0$ 特征值与其他非 $0$ 特征值符号相反，归类为双曲类；
4. 正特征值和负特征值的个数都均大于一，且特征值均非 $0$，归类为超双曲类。

　　 一般来说，对于一个有物理意义的 PDE，特征值有一个为 $0$ 的归类为抛物类，其余全同号为椭圆类，否则大概率为双曲类。

---

[1] ^ Stanley J. Farlow, Partial Differential Equations for Scientists and Engineers. 1993.