好的,我们这次来讲解 “偏微分方程”(Partial Differential Equations,简称 PDEs)。
我会从最基础的概念开始,逐步深入到它的分类、解法思路和一些重要例子。
1. 什么是偏微分方程?
偏微分方程 是包含未知函数的偏导数的方程。
未知函数通常是多元函数,例如 \(u(x, y)\) 或 \(u(t, x)\)。
这与常微分方程(ODE)不同——ODE 的未知函数是一元函数,只涉及常导数。
例子(一维热方程):
\[\frac{\partial u}{\partial t} = k \frac{\partial^2 u}{\partial x^2} \]
这里 \(u(t, x)\) 是未知函数(比如表示温度),\(t\) 是时间,\(x\) 是空间坐标,\(k\) 是热扩散系数。方程表示温度随时间的变化与空间二阶导数(曲率)成正比。
2. 偏微分方程的阶数
方程的阶 是方程中出现的最高阶偏导的阶数。
例如:
\[\frac{\partial^2 u}{\partial x^2} + \frac{\partial^2 u}{\partial y^2} = 0 \]
是二阶 PDE(拉普拉斯方程)。
3. 线性与非线性
- 线性 PDE:未知函数及其所有偏导以线性形式出现(系数可以是自变量的函数)。
例如:
\[ a(x,y) \frac{\partial u}{\partial x} + b(x,y) \frac{\partial u}{\partial y} + c(x,y) u = d(x,y) \]
齐次线性 PDE:\(d(x,y) = 0\)。
- 半线性 PDE:最高阶偏导是线性的,但系数可依赖于未知函数的低阶偏导。
例如:
\[ \frac{\partial^2 u}{\partial x^2} + \frac{\partial^2 u}{\partial y^2} + u \frac{\partial u}{\partial x} = 0 \]
- 拟线性 PDE:最高阶偏导的系数依赖于未知函数及其低阶偏导。
例如:
\[ u_x u_{xx} + u_{yy} = 0 \]
这里最高阶项 \(u_{xx}\) 的系数 \(u_x\) 依赖于未知函数的一阶导。
- 完全非线性 PDE:最高阶偏导非线性出现。
例如:
\[ (u_{xx})^2 + u_{yy} = 0 \]
4. 三类经典线性二阶 PDE
在两个自变量情形下,一般二阶线性 PDE 可写为:
\[A u_{xx} + 2B u_{xy} + C u_{yy} + D u_x + E u_y + F u = G \]
其中 \(A, B, C, \dots\) 可以是 \(x, y\) 的函数。
在固定点处,根据判别式 \(\Delta = B^2 - AC\) 分类:
-
椭圆型(Elliptic):\(\Delta < 0\)
典型例子:拉普拉斯方程 \(u_{xx} + u_{yy} = 0\)
描述稳态现象(无时间演化),如平衡温度分布、势场。
边值问题(Boundary Value Problem)。 -
抛物型(Parabolic):\(\Delta = 0\)
典型例子:热方程 \(u_t = k u_{xx}\)
描述扩散、时间演化且具有“平滑化”效应。
初边值问题(Initial-Boundary Value Problem)。 -
双曲型(Hyperbolic):\(\Delta > 0\)
典型例子:波动方程 \(u_{tt} = c^2 u_{xx}\)
描述振动、波传播,信息以有限速度传播。
初值问题(Cauchy 问题)或初边值问题。
这种分类与二次曲线类似,并且影响到方程的解的性质(如信息传播域、对边界数据的依赖性)。
5. 定解条件
PDE 的通解通常包含任意函数(类似于 ODE 的任意常数),所以要得到特定解,需要附加条件:
- 边界条件(Boundary Conditions):在空间区域的边界上给定函数值(狄利克雷条件)或法向导数(诺伊曼条件)等。
- 初始条件(Initial Conditions):在初始时刻 \(t=0\) 给出未知函数(及时间导数)的值。
例如热方程:给出初始温度分布 \(u(0, x) = f(x)\) 和边界条件 \(u(t, 0) = 0\)。
6. 解法思路简介
没有通用公式解所有 PDE,但有一些经典方法:
-
分离变量法
假设解可写成 \(u(t, x) = T(t) X(x)\),代入 PDE 得到分别关于 \(t\) 和 \(x\) 的常微分方程,求解后再叠加(傅里叶级数)。 -
特征线法(适用于一阶 PDE 和双曲型 PDE)
沿某些曲线(特征曲线)将 PDE 化为 ODE。 -
积分变换法
用傅里叶变换或拉普拉斯变换消去某些自变量(特别是当区域是无界时),将 PDE 化为 ODE。 -
格林函数法 / 基本解
求点源产生的解,然后利用线性叠加原理(卷积)得到一般解。 -
变分法 / 有限元法
将 PDE 转化为等价的变分问题(极小化某个能量泛函),然后数值求解。 -
数值方法:有限差分法、有限元法、有限体积法等。
7. 重要例子
-
拉普拉斯方程 \(\nabla^2 u = 0\)
解的性质:平均值性质、极值原理、光滑性。 -
热方程 \(u_t = k \nabla^2 u\)
基本解是高斯核,体现无穷传播速度与平滑效应。 -
波动方程 \(u_{tt} = c^2 \nabla^2 u\)
达朗贝尔公式(一维情形)显示波以速度 \(c\) 传播,保持初始间断。 -
薛定谔方程(线性)
\(i\hbar \psi_t = -\frac{\hbar^2}{2m} \psi_{xx} + V \psi\)
量子力学的基本方程。 -
非线性例子:
- 伯格斯方程 \(u_t + u u_x = \nu u_{xx}\)(模拟激波)
- KdV 方程 \(u_t + u u_x + \delta u_{xxx} = 0\)(描述孤立子)
- 纳维-斯托克斯方程(流体力学)
8. 适定性概念
一个 PDE 问题如果是适定的(Hadamard),需满足:
- 解存在(Existence)
- 解唯一(Uniqueness)
- 解连续依赖于初始/边界数据(Stability)
许多 PDE 研究围绕证明特定问题的适定性展开。
9. 进一步方向
- 索伯列夫空间:提供 PDE 解的弱解(广义解)的合适函数空间。
- 正则性理论:弱解是否光滑?
- 几何分析:用 PDE 研究几何问题(如里奇流)。
- 自由边界问题:相变、渗流等。
希望这个循序渐进的概述能帮你对“偏微分方程”有一个清晰且系统的初步认识。如果你想深入某个具体类别或方法,我可以继续展开。