数理方法-解偏微分方程（施工中）

习题选讲

验证 u(x,t) = F(x + at) + G(x - at) 是振动方程 u_{tt} - a^2 u_{xx} = 0 的通解，其中 F,G 是任意的二阶可微函数。

解：由链式法则

\begin{aligned} u_t &= F'(x + at) \cdot a + G'(x - at) \cdot (-a) \\ u_{tt} &= F''(x + at) \cdot a^2 + G''(x - at) \cdot (-a)^2 \\ u_x &= F'(x + at) + G'(x - at) \\ u_{xx} &= F''(x + at) + G''(x - at) \end{aligned}

容易得到 u_{tt} - a^2 u_{xx} = 0，表明 u(x,t) = F(x + at) + G(x - at) 是振动方程的解。

下面证明是通解：

令 (x,t) \to (\xi, \eta) 满足 \xi = x + at, \eta = x - at，由链式法则：

\begin{aligned} u_x &= u_\xi + u_\eta \\ u_{xx} &= u_{\xi\xi} + 2u_{\xi\eta} + u_{\eta\eta} \\ u_t &= u_\xi \cdot a + u_\eta \cdot (-a) \\ u_{tt} &= (u_{\xi\xi} - 2u_{\xi\eta} + u_{\eta\eta}) \cdot a^2 \end{aligned}

代入方程 u_{tt} - a^2 u_{xx} = 0，得到 -4a^2 u_{\xi\eta} = 0 或 u_{\xi\eta} = 0。

上述推导表明，u_{tt} - a^2 u_{xx} = 0 \Leftrightarrow u_{\xi\eta} = 0，于是开始积分：

\begin{aligned} u_\xi &= C_1(\xi) \\ u &= \int C_1(\xi) d\xi + C_2(\eta) \end{aligned}

改写原方程为 u(\xi, \eta) = F(\xi) + G(\eta)，是 u_{\xi\eta} = 0 的通解。

那么 u(x,t) = F(x + at) + G(x - at) 也是振动方程 u_{tt} - a^2 u_{xx} = 0 的通解，其中 F(x + at), G(x - at) 分别是左行波和右行波，a 是行波速度。

二阶线性偏微分方程的化简与分类

一般形式

A u_{xx} + 2B u_{xy} + C u_{yy} + D u_x + E u_y + F u = 0 或 f

其中 A,B,C,\ldots 是 (x,y) 的已知函数或者常数，且 A,B,C 不全为 0，u = u(x,y) 是未知函数。

通常情况下，称二阶导数部分为"主要部分"，其余部分简记为 H。于是方程改写为 A u_{xx} + 2B u_{xy} + C u_{yy} = H。

判别式 \Delta

定义 \Delta = (2B)^2 - 4AC = 4(B^2 - AC)，可用"特征线方法"简化方程。

利用理论，方程具有以下形式的特征方程：A (dy)^2 - 2B dx dy + C (dx)^2 = 0。

方程的解（曲线）称为特征线。利用"特征线"可以适当对方程进行简化，并进行分类：

\Delta > 0，双曲型
\Delta < 0，椭圆型
\Delta = 0，抛物型

分类与化简

常系数情形

a_{11} u_{xx} + 2a_{12} u_{xy} + a_{22} u_{yy} = H，其中 a_{11}, a_{12}, a_{22} \in \mathbb{R}，H 是低阶项部分。

令 A = \begin{pmatrix} a_{11} & a_{12} \\ a_{12} & a_{22} \end{pmatrix}，是一个对称矩阵，对应一个二次型 \left( \frac{\partial}{\partial x}, \frac{\partial}{\partial y} \right) A \begin{pmatrix} \frac{\partial}{\partial x} \\ \frac{\partial}{\partial y} \end{pmatrix} u = H，存在正交变换 (x,y) \to (\xi,\eta)，满足

\begin{pmatrix} \frac{\partial}{\partial x} \\ \frac{\partial}{\partial y} \end{pmatrix} = Q \begin{pmatrix} \frac{\partial}{\partial \xi} \\ \frac{\partial}{\partial \eta} \end{pmatrix} \Leftrightarrow \begin{pmatrix} \xi \\ \eta \end{pmatrix} = Q \begin{pmatrix} x \\ y \end{pmatrix}

于是在这个变换下，u(\xi,\eta) 满足方程 \lambda_1 u_{\xi\xi} + \lambda_2 u_{\eta\eta} = \tilde{H}，其中 \lambda_1, \lambda_2 是矩阵 A 的特征值。有以下三种情况：

\lambda_1 \lambda_2 > 0，这时 \Delta < 0，令 \alpha = \frac{1}{\sqrt{\lambda_1}} \xi, \beta = \frac{1}{\sqrt{\lambda_2}} \eta，则方程可以简化为 u_{\alpha\alpha} + u_{\beta\beta} = \tilde{\tilde{H}}，其主部是拉普拉斯方程，对应椭圆型。
\lambda_1 \lambda_2 = 0，这时 \Delta = 0，令 \alpha = \frac{1}{\sqrt{\lambda_1}} \xi, \beta = \eta，则方程可以简化为 u_{\alpha\alpha} = \tilde{\tilde{H}} \triangleq b_1 u_\alpha + b_2 u_\beta + c u，其主部是热传导方程。进一步可以简化为 u_\beta - a^2 u_{\alpha\alpha} = 0 或 f，左边称为导热算子，或热传导的标准型。
\lambda_1 \lambda_2 < 0，这时 \Delta > 0，不妨设 \alpha = \frac{1}{\sqrt{\lambda_1}} \xi, \beta = \frac{1}{\sqrt{-\lambda_2}} \eta，则方程简化为 u_{\alpha\alpha} - u_{\beta\beta} = \tilde{\tilde{H}}，其主部是波动方程。

变系数情形

a_{11} u_{xx} + 2a_{12} u_{xy} + a_{22} u_{yy} = H，其中 a_{11}, a_{12}, a_{22} 是关于 x,y 的函数。

这时，矩阵 A 是函数矩阵，不再适用"二次型正交变换"，因为经过链式法则，得到的是同样复杂、形式相近的式子。比如令 \xi = \xi(x,y), \eta = \eta(x,y)，方程可化为 \bar{a}_{11} u_{\xi\xi} + 2\bar{a}_{12} u_{\xi\eta} + \bar{a}_{22} u_{\eta\eta} = \bar{H}，其中

\begin{aligned} \bar{a}_{11} &= a_{11} \xi_x^2 + 2a_{12} \xi_x \xi_y + a_{22} \xi_y^2 \\ \bar{a}_{12} &= a_{11} \xi_x \eta_x + a_{12} (\xi_x \eta_y + \xi_y \eta_x) + a_{22} \xi_y \eta_y \\ \bar{a}_{22} &= a_{11} \eta_x^2 + 2a_{12} \eta_x \eta_y + a_{22} \eta_y^2 \end{aligned}

如果取一阶偏微分方程 a_{11} \varphi_x^2 + 2a_{12} \varphi_x \varphi_y + a_{22} \varphi_y^2 = 0 的一个特解作为新自变量 \xi = \xi(x,y)，则 \bar{a}_{11} = 0，类似的，如果取另一个特解作为自变量 \eta，则 \bar{a}_{22} = 0，从而方程得以化简。

注意到 a_{11} \varphi_x^2 + 2a_{12} \varphi_x \varphi_y + a_{22} \varphi_y^2 = 0 的求解可以化为常微分方程求解：将方程变形为 a_{11} \left( -\frac{\varphi_x}{\varphi_y} \right)^2 - 2a_{12} \left( -\frac{\varphi_x}{\varphi_y} \right) + a_{22} = 0，由隐函数导数理论，可以得到 a_{11} \left( \frac{dy}{dx} \right)^2 - 2a_{12} \left( \frac{dy}{dx} \right) + a_{22} = 0，或写为 a_{11} dy^2 - 2a_{12} dx dy + a_{22} dx^2 = 0，称为二阶线性偏微分方程的特征方程，特征方程的解称为特征线。