从弱导数到Green公式：Sobolev空间、嵌入定理与有限元

泛函分析学到 Sobolev 空间这一章进度条基本就走到整门课最能落地的地方了。前面从赋范空间、Hahn-Banach、一致有界原理一路推过去抽象归抽象但总有点悬在半空的感觉到了 Sobolev 空间和 Green 公式这里所有抽象工具突然全落到了偏微分方程的弱解、变分问题和有限元离散上。这篇笔记主要想顺三件事Sobolev 空间 $W^{k,p}$ 到底是被什么需求逼出来的嵌入定理和迹定理凭什么成立以及 Green 公式从经典形式推广到弱形式时边界条件究竟发生了什么变化。如果你正在啃 Evans、Brezis 或者国内泛函教材的对应章节被一堆 $1/p - k/n$ 的指数关系绕得头晕那下面这些推导过程和踩坑记录应该能帮你把线捋顺。1. 从导数不够用说起$W^{k,p}$ 是逼出来的任何一本教材讲到 Sobolev 空间开头几乎都是直接甩定义。但说实话如果不知道它在解决什么问题定义看完就忘。所以我想先回到痛点我们本来有一堆好用的函数空间为什么还要再造一个1.1 经典可微空间在 $L^p$ 语境下的短板先看 $C^1([0,1])$ 配上范数 $|f| \sup|f| \sup|f|$这其实是个 Banach 空间看上去挺好。但只要把它放到用积分度量的 $L^p$ 视角里问题立刻冒出来$C^1$ 在 $L^p$ 范数下不完备。举个最扎心的例子取一串光滑函数去逼近折线 $f(x)|x-\tfrac12|$。这串函数本身在 $L^2$ 意义下是 Cauchy 列但它们几乎处处收敛到的那个折线函数在 $x\tfrac12$ 处导数根本不存在导数序列的极限只是一个跳跃量。这意味着如果我们老老实实待在 $C^1$ 里做极小化很可能极小化序列的极限跑出去了。这正是变分法最怕的局面——你构造了一个能量递减的序列却不能保证它有极限落在同一个空间里。经典解的存在性证明老是失败根子就在这里。1.2 弱导数把求导的负担转移给测试函数解决办法非常巧妙既然逐点可导这件事太苛刻那就换个判据。分部积分公式 $\int_\Omega u,\varphi,dx -\int_\Omega u,\varphi,dx$ 对光滑 $u$ 和紧支光滑测试函数 $\varphi\in C_c^\infty(\Omega)$ 成立。我们干脆把这个公式当成定义来用。具体说设 $u\in L^1_{\text{loc}}(\Omega)$如果存在 $v\in L^1_{\text{loc}}(\Omega)$ 使得$$\int_\Omega u,D^\alpha\varphi,dx (-1)^{|\alpha|}\int_\Omega v,\varphi,dx,\quad \forall \varphi\in C_c^\infty(\Omega)$$对所有多重指标 $|\alpha|\le k$ 成立就称 $vD^\alpha u$ 是 $u$ 的 $\alpha$ 阶弱导数。注意这里的精妙之处测试函数无穷可微且紧支所以右边的积分永远有意义求导的脏活被丢给了这些测试函数。弱导数如果存在在同一 $L^1_{\text{loc}}$ 等价类里是唯一的这一点用变分基本引理就能证。提示弱导数不是几乎处处导数。有些函数几乎处处不可导却有弱导数反过来也有处处可导但弱导数不存在的情况。千万别把两者混为一谈。1.3 完备性才是 Sobolev 空间的立身之本有了弱导数Sobolev 空间的定义就顺理成章$$W^{k,p}(\Omega){u\in L^p(\Omega): D^\alpha u\in L^p(\Omega),\ \forall |\alpha|\le k}$$配上范数 $|u|{W^{k,p}}\big(\sum{|\alpha|\le k}|D^\alpha u|_p^p\big)^{1/p}$。$p2$ 时记 $H^kW^{k,2}$是 Hilbert 空间可以做内积。关键定理是$W^{k,p}$ 在这个范数下完备因此是 Banach 空间。为什么这最重要因为变分法、Lax-Milgram、Riesz 表示这一整套工具全都建立在完备性上。没有完备性你构造的近似解序列可能压根不收敛到空间内的元素一切存在性证明就都是空谈。我一开始只把它当技术性引理后来回头再看才意识到这才是 Sobolev 空间真正的价值内核。2. 范数、半范数与几个容易记混的记号定义看完之后实操里最容易栽跟头的地方不是定理本身而是记号。$W^{k,p}$、$W_0^{k,p}$、$H^s$ 这几个来回出现初学时几乎必混淆一次。2.1 半范数在估计里反而更好用$W^{k,p}$ 里最常用的不是完整范数而是只取最高阶导数的半范数$$|u|{W^{k,p}}\Big(\sum{|\alpha|k}|D^\alpha u|_p^p\Big)^{1/p}$$半范数的问题在于它不是范数——常数函数或次数低于 $k$ 的多项式的 $k$ 阶半范数为零却不等于零函数。但正因为这样它特别适合用来做模掉低阶项的论证。在证明 Poincaré 不等式、紧嵌入这些结论时真正起作用的就是这个最高阶项。我踩过的坑是一开始总想着半范数不严谨还是写完整范数保险结果证明越写越臃肿最后发现用半范数加一个低阶范数做等价估计三行就写完了。2.2 $W_0^{k,p}$ 到底零在哪里$W_0^{k,p}(\Omega)$ 的定义是 $C_c^\infty(\Omega)$ 在 $|\cdot|_{W^{k,p}}$ 下的闭包。注意闭包这两个字它意味着这里的边值为零不是逐点点态意义上的零而是可以被光滑紧支函数逼近出零边值。这是初学者最反直觉的一点$W_0^{1,p}$ 里的函数未必在边界上取零但它确实是零边值条件的正确函数类。记法上的建议$H_0^1$ 用得最多它对应齐次 Dirichlet 边值问题的解空间。看到 $H_0^1$脑子里应该自动浮现边界上被钉死了这个图像。2.3 $H^s$ 与 $W^{k,p}$ 的关系别搞反教材里 $H^s$ 有时指 $W^{s,2}$$s$ 可以是实数用 Fourier 变换定义有时又专门指整数阶的 $W^{k,2}$。这个记号在不同书里不完全一致是考试和读文献时的高频混淆点。我的做法是读论文时先看它开头的记号表确认 $H^s$ 是整数阶还是分数阶别硬背。如果 $s$ 为实数且用 Fourier 定义那么 $H^s$ 是分数阶 Sobolev 空间如果只说 $H^k$$k$ 整数通常就是 $W^{k,2}$。这俩在 $s$ 为非整数时不相等差别就在于光滑性的粒度。3. 嵌入定理Sobolev 空间里最值钱的部分嵌入定理是我认为 Sobolev 理论最漂亮也最实用的结论。它干的事很简单告诉你有多少阶、多少可积性的导数就能换来多少光滑性或多少可积性。有了它很多先验估计才能升级成真正有用的结论。3.1 维数决定命运$p$ 与 $n$ 的三种关系设 $\Omega\subset\mathbb{R}^n$ 是有界区域边界适当光滑。对 $W^{1,p}(\Omega)$结果完全取决于 $p$ 和维数 $n$ 的大小关系情形嵌入结果直觉$p n$$W^{1,p}\hookrightarrow L^q$$\frac1q\frac1p-\frac1n$可积性被抬高到 $qp$$p n$$W^{1,p}\hookrightarrow L^q$任意有限 $q$临界情形几乎任意可积$p n$$W^{1,p}\hookrightarrow C^{0,\gamma}$$\gamma1-\frac np$直接换来 Hölder 连续性这张表要背下来。它的物理直觉是导数可积性越强、维数越低函数就越规则。三维空间里 $H^1W^{1,2}$ 落在 $pn$ 那一栏所以 $H^1\hookrightarrow L^6$$n3$ 时 $q6$。这个 $L^6$ 在非线性 PDE 里就是著名的临界指数。3.2 紧嵌入与 Rellich-Kondrachov光有连续嵌入还不够变分法真正需要的是紧嵌入。Rellich-Kondrachov 定理说的是对有界 Lipschitz 区域$W^{1,p}\hookrightarrow\hookrightarrow L^q$$1\le q\frac{np}{n-p}$$pn$ 时是紧的。紧性带来一个巨大好处——有界序列必有强收敛子列这就让极小化序列能真正抽出收敛的子列。注意紧嵌入对 $q$ 有严格上界到临界指数 $q\frac{np}{n-p}$ 时紧性就丢了只剩连续嵌入。这正是临界指数问题特别难的根本原因也是很多非线性分析里紧性缺失讨论的源头。3.3 那个反复出现的 $\frac np$ 到底什么意思$\frac np$ 这个数在 Sobolev 理论里出现的频率高得吓人。它其实是一个尺度分析的结果当你在 $\mathbb{R}^n$ 里做伸缩 $u_\lambda(x)u(\lambda x)$$L^p$ 范数和导数范数按不同的幂次缩放要让两者在同一个尺度下平衡指数就得这么组合。理解这一点之后再看到嵌入定理里那些指数关系就不会觉得是凭空变出来的而是尺度不变性的必然产物。我的经验是碰到记不住的指数关系就在草稿纸上做一次伸缩变换验算两分钟就能推出来比死记强。4. 迹定理为什么边界值能被抓住在经典框架里函数在边界取什么值是天经地义的直接代入就行。但到了 $L^p$ 空间函数只是等价类边界是零测集谈边界值根本没有意义。这是从经典 PDE 跨到弱解理论时最大的观念障碍迹定理就是来填这个坑的。4.1 边界值在 $L^p$ 里根本控制不住$L^p(\Omega)$ 里的函数可以任意改动一个零测集上的取值而不改变等价类而 $\partial\Omega$ 恰恰是零测集。所以$u$ 在边界等于某个函数这句话在 $L^p$ 里是没有定义的。但物理问题里边界条件又必须谈怎么办答案是只要 $u$ 有一阶可积导数即 $u\in W^{1,p}$光滑性就足够让我们在边界上读出一个良定义的值。4.2 迹算子的连续性与它的核对适当的区域Lipschitz 边界即可存在一个连续线性算子$$\text{tr}: W^{1,p}(\Omega)\to W^{1-1/p,p}(\partial\Omega)$$它把内点函数映到边界函数并且在光滑函数上与通常的边界限制一致。注意像空间是分数阶 Sobolev 空间 $W^{1-1/p,p}$指数又掉了一层 $1/p$这个掉一层其实就是光滑性在跨越边界时损失掉的部分跟前面尺度分析是同一套逻辑。更漂亮的是$\text{tr}$ 的核恰好是 $W_0^{1,p}(\Omega)$。这句话把两件事完美地串了起来——被光滑紧支函数逼近和边界迹为零是等价的。我第一次看到这个结论时觉得特别清爽因为之前对 $W_0^{1,p}$ 的理解一直停在闭包这种操作定义上有了迹定理它才真正有了几何意义。4.3 有限元里边界条件凭什么能直接代进去做实... 做数值计算的同学可能更有感受有限元组装刚度矩阵时Dirichlet 边界条件可以直接删行删列或者强制置值。为什么合法本质上就是迹定理在背后保证——因为 $H_0^1$ 里的函数迹为零所以齐次边值问题的解空间就是 $H_0^1$。而自然边界条件Neumann 型则会自动出现在弱形式里不需要强加。这个本质边界条件 vs 自然边界条件的区分理解透彻之后整个有限元的边界处理逻辑就通了。我当年学有限元时是硬记套路直到补了迹定理才明白为什么两类条件处理方式完全不同。5. Green 公式的推广从经典到弱形式终于到 Green 公式了。它在微弱意义下其实就是一个高维分部积分但它把边界条件、弱解定义、对称性这几件事全串在了一起。5.1 经典 Green 公式回顾对光滑函数 $u,v$ 和光滑边界 $\partial\Omega$经典的 Green 第二公式是$$\int_\Omega (u,\Delta v - v,\Delta u),dx \int_{\partial\Omega}\left(u\frac{\partial v}{\partial n}-v\frac{\partial u}{\partial n}\right)dS$$还有第一公式分部积分形式$$\int_\Omega \nabla u\cdot\nabla v,dx -\int_\Omega u,\Delta v,dx \int_{\partial\Omega} u\frac{\partial v}{\partial n},dS$$这两个公式是弱形式推导的全部弹药。经典的推导要求 $u,v\in C^2$、边界 $C^1$条件苛刻。5.2 分部积分在 Sobolev 框架下的严格版本推广的关键在于只要 $u,v\in H^1$公式右边那些量就都有意义了。左边 $\nabla u\cdot\nabla v$ 是 $L^2$ 里的可积函数没问题。右边的面积分怎么办这时迹定理出场——$u$ 在边界上的迹属于 $H^{1/2}$法向导数 $\partial u/\partial n$ 则需要 $u\in H^2$ 才有 $H^{1/2}$ 级别的迹。这解释了为什么不同边值问题的弱形式对函数光滑性的要求不一样二阶问题的弱形式通常只需要 $H^1$ 就够正是因为最高阶导数被分摊掉了一次。更严格地说Sobolev 框架下的分部积分可以写成对 $u\in H^1$、$\mathbf{F}\in H(\text{div})$ 有$$\int_\Omega u,\nabla\cdot\mathbf{F},dx -\int_\Omega \nabla u\cdot\mathbf{F},dx \langle u,\mathbf{F}\cdot n,, 1\rangle_{\partial\Omega}$$其中边界项是一个对偶配对。这种写法在处理混合有限元、间断 Galerkin 时特别有用因为它不需要预先假设边界足够光滑。5.3 弱解与两类边界条件的本质区别以 Poisson 问题 $-\Delta u f$在 $\Omega$ 内、$u0$在 $\partial\Omega$ 上为例。两边同时乘以测试函数 $v\in H_0^1$分部积分$$-\int_\Omega \Delta u, v,dx \int_\Omega \nabla u\cdot\nabla v,dx - \int_{\partial\Omega}\frac{\partial u}{\partial n}v,dS$$因为 $v\in H_0^1$它在边界上的迹为零最后那个面积分自动消失。于是弱形式就是求 $u\in H_0^1$ 使得$$\int_\Omega \nabla u\cdot\nabla v,dx \int_\Omega f v,dx,\quad \forall v\in H_0^1$$这个方程用 Lax-Milgram 直接就有唯一解因为双线性型 $a(u,v)\int\nabla u\cdot\nabla v$ 的强制性由 Poincaré 不等式保证。两类边界条件的区别在这里体现得淋漓尽致Dirichlet 条件是本质条件进了试探空间和测试空间Neumann 条件是自然条件它不走空间约束而是自动出现在边界积分项里。如果把 Neumann 条件 $\partial u/\partial ng$ 加上面积分不消失反而成为方程的一部分。我踩的最大一个坑就是一开始把 Neumann 条件也想塞进空间结果空间定义不下