§1

偏微分方程(Partial Differential Equation, PDE)是含有未知多元函数及其偏导数的方程。它刻画"场"随空间与时间的演化,是连续介质物理的基石。

一、什么是 PDE:阶、线性、拟线性

设 $u(x,y)$ 或 $u(x,t)$ 为未知函数,含偏导数如 $u_x,u_{xx},u_{xy}$ 的方程即 PDE。

  • 阶:方程中出现的最高阶偏导数的阶。例如 $u_{xx}+u_{yy}=0$ 是二阶。
  • 线性 PDE:未知函数 $u$ 及其各偏导数均为一次,且系数只依赖于自变量。如

$$a(x,y)u_{xx}+b(x,y)u_{xy}+c(x,y)u_{yy}+d(x,y)u_x+e(x,y)u_y+f(x,y)u=g(x,y).$$

其中 $u$ 不乘 $u$,也不出现 $u_x^2$、$u u_x$ 等非线性项。

  • 拟线性 PDE:最高阶项系数只依赖于自变量与低阶项(但低阶项可非线性)。如一阶拟线性

$$a(x,y,u)u_x+b(x,y,u)u_y=c(x,y,u).$$

  • 半线性 / 完全非线性:依最高阶项是否线性进一步细分(如 $u_{xx}u_{yy}-u_{xy}^2=0$ 完全非线性)。

> 本模块聚焦线性二阶 PDE,因为它们有最完整的理论(叠加原理、分离变量、傅里叶方法)。

常见方程一览

方程型物理含义
$u_t=\alpha u_{xx}$抛物热扩散、浓度扩散
$u_{tt}=c^2 u_{xx}$双曲波传播、振动
$u_{xx}+u_{yy}=0$椭圆稳态温度、静电势
$u_t+u u_x=0$拟线性一阶交通流、浅水波

二、两个自变量的二阶线性 PDE 与分类

最一般形式:

$$A\,u_{xx}+2B\,u_{xy}+C\,u_{yy}+D\,u_x+E\,u_y+F\,u=G,$$

其中 $A,B,C$ 为 $x,y$ 的函数(为简单先视为常数)。其类型由判别式

$$\Delta=B^2-A\,C$$

决定:

判别式类型典型方程物理特征
$\Delta>0$双曲型$u_{tt}=c^2 u_{xx}$(波动)有特征线、信息传播
$\Delta=0$抛物型$u_t=\alpha u_{xx}$(热)单向演化、扩散、平滑
$\Delta<0$椭圆型$u_{xx}+u_{yy}=0$(拉普拉斯)无时间维、稳态、调和

分类的几何直觉:二次多项式 $A\,X^2+2B\,XY+C\,Y^2$ 的定型性,对应方程"沿某些方向有信息传播(双曲)""只沿一个方向演化(抛物)""各向同性无演化(椭圆)"。

例. 判断 $u_{xx}+2u_{xy}+u_{yy}=0$ 的类型。

$A=1,B=1,C=1\Rightarrow \Delta=1-1=0$,为抛物型(退化情形,对应热类)。

判别式的来源(特征线视角)

分类的本质是看方程沿哪些方向"传播信息"。对二阶主部

$$A\,u_{xx}+2B\,u_{xy}+C\,u_{yy}=0,$$

作自变量的可逆变换 $(x,y)\mapsto(\xi,\eta)$,使主部只余 $u_{\xi\xi}$、$u_{\xi\eta}$ 或 $u_{\eta\eta}$。方向 $dy/dx$ 满足的特征方程

$$A\left(\frac{dy}{dx}\right)^2-2B\left(\frac{dy}{dx}\right)+C=0.$$

其判别式正是 $B^2-AC$:

  • $\Delta>0$:有两个相异实根方向 → 两条特征线 → 双曲型;
  • $\Delta=0$:一个重根方向(退化) → 抛物型;
  • $\Delta<0$:无实特征方向 → 椭圆型。

这说明"分类"不是人为约定,而是方程内在的几何结构。

三、定解条件与适定性

PDE 单独不能定解,还需定解条件。对一个含时间 $t$ 的问题,常见组合:

  • 初值问题(Cauchy 问题):给定 $t=0$ 时的状态 $u(x,0)=\varphi(x)$、$u_t(x,0)=\psi(x)$。波动方程需要初位移+初速度。
  • 边值问题:在区域边界 $\partial\Omega$ 上给条件:
  • Dirichlet 条件:$u|_{\partial\Omega}=g$(给定边界函数值,如边界温度);
  • Neumann 条件:$\frac{\partial u}{\partial n}\big|_{\partial\Omega}=h$(给定边界法向导数,如绝热/绝缘);
  • Robin 条件:混合 $a u+b\,u_n=c$。
  • 混合问题:既有初值又有边界(如有限长杆的热传导:$u(x,0)=\varphi(x)$ 且两端 $u(0,t),u(L,t)$ 给定)。

三类边界条件对比

条件数学形式物理意义典型场景
Dirichlet$u=g$边界量已知边界恒温、固定位移
Neumann$u_n=h$边界通量已知绝热、绝缘、自由端
Robin$a u+b\,u_n=c$量与通量耦合对流换热、弹性支撑

一个"好"的定解问题应具备适定性(Hadamard):解存在、唯一、且连续依赖于定解数据。椭圆型多配边值,抛物/双曲型多配初值或混合条件——这正是分类的物理意义。

四、叠加原理

对线性齐次 PDE,若 $u_1,u_2$ 是解,则任意线性组合 $c_1u_1+c_2u_2$ 也是解;若右端非齐次 $Lu=G$,且 $G=G_1+G_2$,可分别求特解再相加。叠加原理是分离变量法与傅里叶级数能"拼"出一般解的根本原因。

叠加原理的两个层次

1. 齐次方程:解的线性组合仍是解。这是傅里叶级数"逐模态叠加"合法性的基础。

2. 非齐次方程:若 $Lu_1=G_1,\ Lu_2=G_2$,则 $u_1+u_2$ 满足 $L(u_1+u_2)=G_1+G_2$。因此非齐次问题可拆成"齐次解 + 特解"。

五、例题与解答

例题1(分类). 分类 $u_{xx}+u_{yy}=u_t$(一维热+源的空间项)。

把它写成对 $(x,t)$ 的二阶方程:$u_{xx}+0\cdot u_{xt}-u_t$ 中时间只出现一阶,故对 $(x,t)$ 是"二阶对 $x$、一阶对 $t$"。标准二阶两变量形式取变量 $(x,t)$:$A=1,B=0,C=0\Rightarrow \Delta=0$,抛物型——正是热方程类。

例题2(定解条件匹配). 一根长 $L$ 的均匀杆,一端恒温 $0$、一端恒温 $T_0$,初始温度分布 $\varphi(x)$。写出定解问题。

$$\begin{cases}u_t=\alpha u_{xx},&00,\\ u(0,t)=0,\quad u(L,t)=T_0,&t>0\ \text{(Dirichlet 边界)},\\ u(x,0)=\varphi(x),&0

这是典型的抛物型混合问题。

例题3(分类与特征). 判断 $u_{tt}-4u_{xx}=0$ 的类型并写出特征线。

视变量为 $(t,x)$,则 $A=-4,\ B=0,\ C=1$(对应 $C u_{tt}+2B u_{tx}+A u_{xx}$),或用 $(x,t)$ 记 $u_{xx}-4u_{tt}=0$ 得 $A=1,B=0,C=-4$,$\Delta=0-(-4)=4>0$,双曲型。特征方程 $dx/dt=\pm 2$,即 $x-2t=\text{const}$ 与 $x+2t=\text{const}$ 两族特征线,波速 $c=2$。

例题4(Neumann 条件的物理). 一根绝热长杆,两端绝热,初值 $\varphi(x)$。应选何种边界条件?

绝热意味着边界无热流,即法向导数为 $0$:$u_x(0,t)=0,\ u_x(L,t)=0$(Neumann 条件)。由于无热量进出,总热量守恒,温度最终趋于初始平均温度 $\frac1L\int_0^L\varphi(x)\,dx$。

练习

1. 判断下列方程类型(双曲/抛物/椭圆):(a) $u_{tt}=9u_{xx}$;(b) $u_t=u_{xx}+u_{yy}$;(c) $u_{xx}-2u_{xy}+u_{yy}=0$。

2. 一根长 $L$ 的杆,左端恒温 $0$、右端绝热,初值 $\varphi(x)$。写出定解问题(注明每个条件类型)。

3. 设 $u_1$ 满足 $u_t=\alpha u_{xx}$,$u_2$ 满足 $u_t=\alpha u_{xx}+f(x,t)$。利用叠加原理,说明如何由它们构造 $u_t=\alpha u_{xx}+2f(x,t)$ 的解。

4. 对 $u_{xx}+u_{yy}=0$,为什么只需边界条件、不需要初值?

参考答案与提示

1. (a) 双曲($\Delta=9>0$);(b) 抛物(时间一阶、空间二阶,主部判别式为 $0$);(c) 抛物($\Delta=(-1)^2-1\cdot1=0$)。

2. $\begin{cases}u_t=\alpha u_{xx},&00,\\ u(x,0)=\varphi(x),&0

3. 由线性性,$u=u_1+u_2$ 满足 $u_t=\alpha u_{xx}+f$;再叠加一次得 $u_1+2u_2$ 满足右端 $2f$。

4. 椭圆型无时间维,描述"处处已达平衡"的稳态场,只由边界值唯一确定(极值原理/最大值原理保证唯一性),故无需、也不存在初值。

本章小结

  • PDE 按最高阶偏导的线性结构分类;两变量二阶线性由 $\Delta=B^2-AC$ 判为双曲/抛物/椭圆。
  • 三类对应热(扩散)、波(传播)、拉普拉斯(稳态)三大物理过程。
  • 解 PDE 须配定解条件;线性齐次方程满足叠加原理,是后续分离变量法的基础。
  • 分类由特征线结构决定:双曲有两条实特征线、抛物一条、椭圆无实特征线。

互动演示

拖动下方控件观察动态过程。

二阶 PDE 分类:判别式 Δ=B²−A·C。拖动 A/B/C 观察类型变化。