1651 字
8 分钟
泛函分析 16|变分方法与微分方程的弱解

数学总导航 · 本课程目录 · ← 上一章 · 下一章 →

两端固定的晾衣绳,挂上重量后弯成什么形状?#

将绳子拉紧,两端固定,均匀挂上轻小负载。用 u(t)u(t) 表示向下位移,一个归一化的小挠度模型是 −u′′=1-u''=1、u(0)=u(1)=0u(0)=u(1)=0。

解是 u(t)=t(1−t)/2u(t)=t(1-t)/2,中点下垂最多,为 1/81/8。这个数依赖我们选择的单位与归一化张力;它只是用来展示计算机制,不是任何实绳的直接工程预测。

如果重量分布粗糙甚至分段变化,直接要求每一点都存在二阶导数会很别扭。弱形式把要求转成整体平衡:让每一种小的试探位移都满足相应的积分关系。接下来会逐步说明这个转写过程,以及为什么它保证唯一解。

从绳子的方程走到积分平衡,中间不能省#

先假设 u,vu,v 足够光滑,且试探位移 vv 在两端为零。把 −u′′=f-u''=f 乘 vv 再积分:−∫u′′v=∫fv-\int u''v=\int fv。分部积分后,左边为 −[u′v]01+∫u′v′-[u'v]_0^1+\int u'v'。

因为 vv 端点为零,边界项消失,留下 ∫u′v′=∫fv\int u'v'=\int fv。现在两边只要求一阶弱导数与函数平方可积,就可以把这条关系延伸为弱解定义。

这不是说所有经典方程都可无条件倒推回来;从弱解得到更强的经典正则性还需额外条件。这里先保证弱意义下问题完整且稳定。

把能量增量完整展开#

能量 E(v)=12∫∣v′∣2−∫fvE(v)=\tfrac12\int|v'|^2-\int fv。在候选 uu 周围改变 εh\varepsilon h,有

E(u+εh)−E(u)=ε(∫u′h′−∫fh)+ε22∫∣h′∣2.E(u+\varepsilon h)-E(u) =\varepsilon\left(\int u'h'-\int fh\right) +\frac{\varepsilon^2}{2}\int|h'|^2.

若 uu 为最小点,一阶系数对所有 hh 必须为零,便得到弱形式。反过来,一旦一阶系数全零,余下二次项非负,说明 uu 真是全局最小点。固定边界又排除非零常数方向,使最小点唯一。

用一个形状模板就能算出答案#

先只允许 uN(t)=c t(1−t)u_N(t)=c\,t(1-t)。代入能量,∫01(1−2t)2dt=1/3\int_0^1(1-2t)^2dt=1/3,∫01t(1−t)dt=1/6\int_0^1t(1-t)dt=1/6;当 f=1f=1,

E(c)=c26−c6=16(c−1/2)2−124.E(c)=\frac{c^2}{6}-\frac c6 =\frac16(c-1/2)^2-\frac1{24}.

因此 c=1/2c=1/2 最优。这个形状模板恰好包含精确解,所以一次就找到答案。一般载荷下,一个模板不够,就增加更多基函数,得到 Galerkin 线性系统。

强制性到底防住了什么#

强制性要求能量二次部分至少按 α∥u∥2\alpha\|u\|^2 增长,防止解逃到很远但代价不增加。连续性则控制相互作用不会突然放大。

Lax–Milgram 把这两种控制结合起来,给存在、唯一和稳定性。Céa 引理继续告诉你:数值解与真解的差,不会比所选有限维空间能做到的最佳近似差太多。理论并不替你选择网格,却能明确指出改善近似应往哪里努力。

带你做:载荷加倍,解为什么也加倍#

对弱问题 a(u,v)=F(v)a(u,v)=F(v),如果载荷变为 2F2F,线性性使 2u2u 满足新方程;唯一性保证它就是新解。因此在线性小挠度模型中,重量加倍,下垂也加倍。

若载荷发生小扰动 ΔF\Delta F,两解之差满足 a(Δu,v)=ΔF(v)a(\Delta u,v)=\Delta F(v)。取 v=Δuv=\Delta u,由强制性得

α∥Δu∥2≤∥ΔF∥∥Δu∥,∥Δu∥≤∥ΔF∥α.\alpha\|\Delta u\|^2\le\|\Delta F\|\|\Delta u\|, \quad\|\Delta u\|\le\frac{\|\Delta F\|}{\alpha}.

这份估计直接回答“称重误差会带来多少形状误差”,而不是只保证一个抽象解存在。

真实绳子在大变形时可能不再线性,材料也可能改变张力。模型中的加倍规律依赖线性前提;数学推导清楚,才能知道实际应用时要核验哪一种近似是否仍合理。

回头整理:定义、工具与练习#

不直接追问曲线的二阶导数#

考虑一根固定两端的弦:

−u′′=f在 (0,1),u(0)=u(1)=0.-u''=f\quad\text{在 }(0,1),\qquad u(0)=u(1)=0.

如果 ff 不连续,要求处处存在经典二阶导数可能太严格。把方程乘测试函数后积分,转移一次导数,得到

∫01u′v′ dt=∫01fv dt对所有 v∈H01(0,1).\int_0^1u'v'\,dt=\int_0^1fv\,dt \quad\text{对所有 }v\in H_0^1(0,1).

现在只需一阶弱导数平方可积。这叫弱形式,边界条件通过函数空间编码。它不是随意放宽,而是保留积分意义下的全部平衡关系。

Lax–Milgram 定理#

以下使用实 Hilbert 空间版本。设 a:H×H→Ra:H\times H\to\mathbb R 为双线性形式,满足

∣a(u,v)∣≤M∥u∥∥v∥,a(v,v)≥α∥v∥2,α>0.|a(u,v)|\le M\|u\|\|v\|,\qquad a(v,v)\ge\alpha\|v\|^2,\quad\alpha>0.

后者称强制性。对每个有界线性泛函 FF,存在唯一 u∈Hu\in H 使 a(u,v)=F(v)a(u,v)=F(v) 对全部 vv 成立,并有 ∥u∥≤∥F∥/α\|u\|\le\|F\|/\alpha。

证明路线:Riesz 表示把 aa 写为 ⟨Au,v⟩\langle Au,v\rangle。强制性给 ∥Au∥≥α∥u∥\|Au\|\ge\alpha\|u\|,故 AA 单射且像闭。若 yy 垂直于像,取 v=yv=y 得 a(y,y)=0a(y,y)=0,强制性迫使 y=0y=0,故像还稠密,从而满射。范数界由把 v=uv=u 代回得到。

复数形式应改用半双线性并要求 Re⁡a(v,v)≥α∥v∥2\operatorname{Re}a(v,v)\ge\alpha\|v\|^2,本章不混用两种约定。

核验 Poisson 问题的每个条件#

取 H=H01(0,1)H=H_0^1(0,1),用 ∥u∥H=∥u′∥2\|u\|_H=\|u'\|_2,由 Poincaré 不等式知它与原 H1H^1 范数等价,故完备。

令 a(u,v)=∫u′v′a(u,v)=\int u'v',则有界常数和强制常数都可取 1。若 f∈L2f\in L^2,F(v)=∫fvF(v)=\int fv 满足

∣F(v)∣≤∥f∥2∥v∥2≤∥f∥2∥v′∥2.|F(v)|\le\|f\|_2\|v\|_2\le\|f\|_2\|v'\|_2.

因此弱解存在唯一,而且 ∥u′∥2≤∥f∥2\|u'\|_2\le\|f\|_2。这同时给了输入扰动的稳定性估计。

当 f=1f=1 时,显式解为 u(t)=t(1−t)/2u(t)=t(1-t)/2。把它代入弱形式可验证;对一般 L2L^2 输入,即使没有漂亮公式,定理仍给结论。

能量最小化#

当 aa 对称时,同一个 uu 最小化

E(v)=12a(v,v)−F(v).E(v)=\frac12a(v,v)-F(v).

展开可得 E(u+h)−E(u)=a(h,h)/2≥α∥h∥2/2E(u+h)-E(u)=a(h,h)/2\ge\alpha\|h\|^2/2。所以驻点不仅局部最优,而且唯一全局最优。若 aa 不对称,Lax–Milgram 仍可成立,但不能用这个能量泛函直接产生原方程。

更一般的直接法:证明能量有下界并强制;取极小化序列;用自反性提取弱收敛子列;保证可行集弱闭、能量弱下半连续;最后用严格凸性得唯一。每一步有独立条件。

Galerkin:回到有限个未知数#

选有限维子空间 HN=span⁡{ϕ1,…,ϕN}H_N=\operatorname{span}\{\phi_1,\ldots,\phi_N\},设 uN=∑cjϕju_N=\sum c_j\phi_j,要求 a(uN,ϕi)=F(ϕi)a(u_N,\phi_i)=F(\phi_i)。得到刚度矩阵 Aij=a(ϕj,ϕi)A_{ij}=a(\phi_j,\phi_i)。

Galerkin 正交性为 a(u−uN,vN)=0a(u-u_N,v_N)=0。Céa 引理给出

∥u−uN∥H≤Mαinf⁡vN∈HN∥u−vN∥H.\|u-u_N\|_H\le\frac M\alpha\inf_{v_N\in H_N}\|u-v_N\|_H.

证明把误差的强制性界与正交性结合,再用双线性形式的有界性即可。它告诉我们:离散解的误差与所选空间的最佳近似误差同阶,泛函分析真正连接了存在性与数值算法。

练习与解答#

  1. 对 f=1f=1、H1=span⁡{t(1−t)}H_1=\operatorname{span}\{t(1-t)\},求 Galerkin 系数。
  2. 若删去零边界条件,只用导数范数会怎样?

解答:1. ∫(1−2t)2dt=1/3\int(1-2t)^2dt=1/3,∫t(1−t)dt=1/6\int t(1-t)dt=1/6,所以系数为 1/21/2,恰好恢复精确解。2. 非零常数的导数范数为零,不再是范数,强制性与唯一性都需要重新处理。


数学总导航 · 本课程目录 · ← 上一章 · 下一章 →

泛函分析 16|变分方法与微分方程的弱解
https://f1yingwhite.github.io/posts/数学/泛函分析/16-变分方法与微分方程的弱解/
作者
F1yingWhite
发布于
2026-11-03
许可协议
CC BY-NC-SA 4.0