三大经典偏微分方程的分离变量法
2023年2月7日
本文作为 UM-SIN2004 课程内容的归纳总结,只讨论二阶偏微分方程边值问题的分离变量法。当然,除了分离变量法还有极坐标系、格林函数法等等。
我们先来看一下这三类方程模型和应用场景:
方程类型 | 标准形式 | 典型应用场景 |
---|
波动方程 (双曲型) | ∂t2∂2u=c2∇2u | 声波、光波、弦振动 |
热传导方程 (抛物型) | ∂t∂u=α∇2u | 热扩散、粒子布朗运动 |
拉普拉斯方程 (椭圆型) | ∇2u=0 | 静电势、稳态温度场 |
1. 分离变量法通用框架
1.1 核心思想
分离变量法(Separation of Variables)通过假设解为单变量函数的乘积,将偏微分方程转化为常微分方程组:
u(x,t)=X(x)Y(y)⋯T(t)
适用条件:
- 线性齐次PDE
- 边界条件可分离
- 定义域规则(矩形、圆形等)
1.2 通用求解步骤
- 变量分离:设解为乘积形式,代入原方程
- 常微分方程:整理得到关于各变量的ODE
- 本征值问题:通过边界条件求解空间部分
- 时间演化:根据PDE类型确定时间部分解
- 叠加解:线性组合所有本征模
- 系数确定:利用初始/边界条件展开
2. 波动方程求解(一维弦振动)
2.1 定解问题
⎩⎨⎧∂t2∂2u=c2∂x2∂2u,u(0,t)=u(L,t)=0u(x,0)=ϕ(x),ut(x,0)=ψ(x)0<x<L,t>0
2.2 求解过程
-
分离变量:设 u(x,t)=X(x)T(t),得:
c2TT′′=XX′′=−λ
-
空间部分:
X′′+λX=0⇒Xn(x)=sin(Lnπx),λn=(Lnπ)2
-
时间部分:
Tn′′+c2λnTn=0⇒Tn(t)=Ancos(Lnπct)+Bnsin(Lnπct)
-
通解:
u(x,t)=n=1∑∞[Ancos(Lnπct)+Bnsin(Lnπct)]sin(Lnπx)
-
系数确定:
- An 来自 ϕ(x) 的傅里叶正弦展开
- Bn 来自 ψ(x) 的傅里叶正弦展开
3. 热传导方程求解(一维杆)
3.1 定解问题
⎩⎨⎧∂t∂u=α∂x2∂2u,u(0,t)=u(L,t)=0u(x,0)=f(x)0<x<L,t>0
3.2 求解过程
-
分离变量:设 u(x,t)=X(x)T(t),得:
αTT′=XX′′=−λ
-
空间部分(同波动方程):
Xn(x)=sin(Lnπx)
-
时间部分:
Tn(t)=Cne−α(Lnπ)2t
-
通解:
u(x,t)=n=1∑∞Cne−α(Lnπ)2tsin(Lnπx)
-
系数确定:
- Cn 来自 f(x) 的傅里叶正弦展开
4. 拉普拉斯方程求解(矩形域)
4.1 定解问题
⎩⎨⎧∂x2∂2u+∂y2∂2u=0,u(0,y)=u(a,y)=0u(x,0)=0,u(x,b)=g(x)0<x<a,0<y<b
4.2 求解过程
-
分离变量:设 u(x,y)=X(x)Y(y),得:
XX′′=−YY′′=−λ
-
x方向:
Xn(x)=sin(anπx),λn=(anπ)2
-
y方向:
Yn(y)=Dnsinh(anπy)
-
通解:
u(x,y)=n=1∑∞Dnsin(anπx)sinh(anπy)
-
系数确定:
- 通过 g(x) 的傅里叶正弦展开求 Dn
参考
