本页目录

基础衔接 18 · 热核正则化与随机分布:在哪个空间里真的收敛?

先修:分布与弱导数、Sobolev 空间、Wick 平方与截断极限。Wick 平方讲只让一个随机对象与常数测试函数配对;本讲把整个 Fourier 场放入 \(L^2(\Omega;H^s)\),给出精确正则性阈值,再证明热核平滑在同一概率空间上回到原随机分布。

每个光滑测试函数都看见极限,为什么还要指定 Sobolev 空间?

1. 先把圆周、基和范数的归一化钉牢

把圆周写成 \(\mathbb T=[0,2\pi]\),使用归一化测度

\[ \frac{dx}{2\pi}. \]

实三角基

\[ e_k^c(x)=\sqrt2\cos(kx), \qquad e_k^s(x)=\sqrt2\sin(kx), \qquad k\ge1 \]

在 \(L^2(\mathbb T,dx/(2\pi))\) 中正交归一。本讲去掉常数模;这让 \(\alpha=0\) 的例子成为均值为零的空间白噪声,不是包含独立零模的完整白噪声。

令 \(\xi_k^c,\xi_k^s\) 全部相互独立且服从 \(N(0,1)\),固定 \(\alpha\ge0\),定义

\[ X_N(x)=\sum_{k=1}^N k^{-\alpha} \bigl(\xi_k^c e_k^c(x)+\xi_k^s e_k^s(x)\bigr). \]

所有截断共用同一组随机系数。从 \(N\) 到 \(M\) 只打开新频率,不重新抽取旧频率。

若 \(f\) 的实 Fourier 坐标为 \(\widehat f_k^c,\widehat f_k^s\),本讲采用

\[ \|f\|_{H^s}^2 =\sum_{k\ge1}(1+k^2)^s \left(|\widehat f_k^c|^2+|\widehat f_k^s|^2\right) \]

作为去零模部分的 Sobolev 范数。\(s<0\) 不是“负的范数”,而是高频受到较小权重的较粗空间;其完备化包含普通点值未必存在的分布。

2. 不只看一个配对:直接计算整个空间范数

正交性和 \(E[(\xi_k^c)^2]=E[(\xi_k^s)^2]=1\) 给出

\[ \boxed{ E\|X_N\|_{H^s}^2 =2\sum_{k=1}^N(1+k^2)^s k^{-2\alpha}. } \]

前面的因子 2 来自每个正频率的余弦与正弦两份实坐标。若改用复指数基,必须同时处理正负频率与实值共轭约束;不能只换符号而保留错误计数。

对 \(M>N\),同一耦合使前 \(N\) 个频率逐项消掉,因此

\[ \boxed{ E\|X_M-X_N\|_{H^s}^2 =2\sum_{k=N+1}^M(1+k^2)^s k^{-2\alpha}. } \]

当 \(k\) 很大时,被加项与 \(2k^{2s-2\alpha}\) 同阶。\(p\) 级数判据于是给出

\[ \sum_{k\ge1}(1+k^2)^s k^{-2\alpha}<\infty \quad\Longleftrightarrow\quad 2s-2\alpha<-1. \]

也就是精确阈值

\[ \boxed{s<\alpha-\frac12.} \]

在这个范围内,\((X_N)\) 是 \(L^2(\Omega;H^s)\) Cauchy 序列。因为 \(L^2(\Omega;H^s)\) 完备,存在同一概率空间上的 \(H^s\) 值随机变量 \(X\),使

\[ X_N\longrightarrow X \qquad\text{于 }L^2(\Omega;H^s). \]

这句话同时指定了随机变量拓扑和空间拓扑,比“每个频率看起来稳定”强得多。

3. 临界值不能靠有限曲线蒙混

在临界值 \(s=\alpha-1/2\),取倍增截断 \(M=2N\):

\[ E\|X_{2N}-X_N\|_{H^{\alpha-1/2}}^2 =2\sum_{k=N+1}^{2N} (1+k^2)^{\alpha-1/2}k^{-2\alpha}. \]

单项可以写成

\[ \frac2k\left(1+\frac1{k^2}\right)^{\alpha-1/2}, \]

所以

\[ \boxed{ E\|X_{2N}-X_N\|_{H^{\alpha-1/2}}^2 \longrightarrow2\log2\ne0. } \]

序列在临界空间不是 Cauchy。若 \(s>\alpha-1/2\),倍增壳层更不会趋零。图上某几个有限 \(N\) 的曲线变平,不能覆盖这个无限频率反例。

这里的二阶矩发散已经否定 \(L^2(\Omega;H^s)\) 收敛,但单凭“期望无穷”还不能断言每个样本的范数无穷。要得到几乎必然的临界失败,再设

\[w_k=(1+k^2)^{\alpha-1/2}k^{-2\alpha}\sim k^{-1},\qquad Z_k=(\xi_k^c)^2+(\xi_k^s)^2.\]

\(Z_k\) 独立,\(E Z_k=2\)、\(\operatorname{Var}Z_k=4\)。因为 \(\sum_k4w_k^2<\infty\),独立中心随机变量级数 \(\sum_k w_k(Z_k-2)\) 几乎必然收敛(可由其部分和是有界于 \(L^2\) 的鞅及鞅收敛定理得到)。而 \(\sum_k2w_k=\infty\),故 \(\sum_k w_kZ_k=\infty\) 几乎必然。这才证明 \(X\) 几乎必然不属于临界空间;更高阶空间嵌入临界空间,因而也被排除。低于临界时范数平方期望有限,则范数几乎必然有限。

在收敛范围内还能给出可交卷的尾界。令

\[ p=2\alpha-2s>1, \qquad C_s=2^{\max(s,0)}. \]

对 \(k\ge1\),有 \((1+k^2)^s\le C_sk^{2s}\),故

\[ \boxed{ E\|X-X_N\|_{H^s}^2 \le\frac{2C_s}{p-1}N^{1-p}. } \]

这是均方 \(H^s\) 误差界;开平方才是 \(L^2(\Omega;H^s)\) 范数界。

4. 一个测试函数的收敛,究竟少了什么?

对光滑测试函数 \(\varphi\),记

\[ \widehat\varphi_k^c=\langle\varphi,e_k^c\rangle, \qquad \widehat\varphi_k^s=\langle\varphi,e_k^s\rangle, \]

其中配对使用 \(dx/(2\pi)\)。则

\[ \langle X_N,\varphi\rangle =\sum_{k=1}^Nk^{-\alpha} \left(\xi_k^c\widehat\varphi_k^c +\xi_k^s\widehat\varphi_k^s\right), \]

并且

\[ \boxed{ \operatorname{Var}\langle X_N,\varphi\rangle =\sum_{k=1}^Nk^{-2\alpha} \left(|\widehat\varphi_k^c|^2 +|\widehat\varphi_k^s|^2\right). } \]

光滑函数的 Fourier 系数比任意固定负幂衰减得更快,所以这份标量方差级数收敛。对每个固定 \(\varphi\),配对因此在 \(L^2(\Omega)\) 中有极限。

但只检查一个 \(\varphi\) 完全可能漏掉高频。例如 \(\varphi(x)=\cos x\) 时,\(N\ge1\) 后

\[ \langle X_N,\cos x\rangle =\frac1{\sqrt2}\xi_1^c, \]

与所有 \(k\ge2\) 的坐标无关。这个标量从第一步就“完美收敛”,却无法判断整个场是否在 \(H^{-1/2}\)、\(L^2\) 或任何其他指定空间收敛。

要声称随机分布在某个拓扑中收敛,需要一组控制测试函数族或空间范数的一致估计,并处理随机律的紧性或强耦合。本模型用上一节的 \(L^2(\Omega;H^s)\) 尾和直接完成了这件事;单个配对只是必要的观察窗口,不是充分证书。

5. 为什么白噪声不是“每一点抽一个普通随机数”?

取 \(\alpha=0\)。对任意固定 \(x\),截断场的点方差为

\[ E[X_N(x)^2] =2\sum_{k=1}^N \left(\cos^2kx+\sin^2kx\right) =2N. \]

所以 \(X_N(x)\) 不会在 \(L^2(\Omega)\) 中收敛。事实上它服从 \(N(0,2N)\),对每个有限 \(M\),\(P(|X_N(x)|\le M)\to0\),连有限随机点值的依概率极限也不存在。更本质地,上一节给出的极限只属于

\[ H^s\qquad(s<-1/2), \]

而点值评价并不是这些负正则性空间上的连续泛函。白噪声的基本对象是 \(\varphi\mapsto\langle X,\varphi\rangle\),不是一张由相互独立点值组成的普通函数表。

由于本讲去掉零模,\(\alpha=0\) 时对一般测试函数有

\[ \operatorname{Var}\langle X,\varphi\rangle =\|\varphi-\langle\varphi,1\rangle\|_{L^2}^2. \]

若把独立常数模补回,右边才变成完整的 \(\|\varphi\|_{L^2}^2\)。这个差别不影响高频正则性阈值,但影响协方差的精确公式。

6. 热核怎样平滑,又在哪个拓扑回去?

圆周 Laplacian 满足

\[ \Delta e_k^c=-k^2e_k^c, \qquad \Delta e_k^s=-k^2e_k^s. \]

因此对 \(\epsilon>0\),定义

\[ X_\epsilon=e^{\epsilon\Delta}X =\sum_{k\ge1}e^{-\epsilon k^2}k^{-\alpha} \bigl(\xi_k^ce_k^c+\xi_k^se_k^s\bigr). \]

热核不是“把所有频率减去同一个常数”,而是给第 \(k\) 个频率乘 \(e^{-\epsilon k^2}\)。高频按 \(k^2\) 指数压低。对任意 \(r\in\mathbb R\),

\[ E\|X_\epsilon\|_{H^r}^2 =2\sum_{k\ge1}(1+k^2)^r k^{-2\alpha}e^{-2\epsilon k^2}<\infty. \]

所以对每个固定 \(\epsilon>0\),该级数在任意 Sobolev 阶数收敛。取可数个整数 \(r\) 并用圆周上的 Sobolev 嵌入,可得 \(X_\epsilon\) 几乎必然是 \(C^\infty\) 函数。

现在固定任意 \(s<\alpha-1/2\),并让 \(X_\epsilon\) 与 \(X\) 使用同一组 \(\xi_k^c,\xi_k^s\)。则

\[ E\|X_\epsilon-X\|_{H^s}^2 =2\sum_{k\ge1}(1+k^2)^s k^{-2\alpha} \left(1-e^{-\epsilon k^2}\right)^2. \]

对每个固定 \(k\),括号随 \(\epsilon\downarrow0\) 趋于零;它又被 1 控制,而

\[ 2(1+k^2)^s k^{-2\alpha} \]

在所选 \(s\) 下可求和。离散支配收敛定理于是给出

\[ \boxed{ X_\epsilon\longrightarrow X \quad\text{于 }L^2(\Omega;H^s), \qquad s<\alpha-\frac12. } \]

这里同时写清了概率拓扑、空间拓扑和耦合方式。若每个 \(\epsilon\) 都独立重采样一个 \(X_\epsilon'\),则交叉项期望为零,反而有

\[ E\|X_\epsilon'-X\|_{H^s}^2 =E\|X_\epsilon'\|_{H^s}^2+E\|X\|_{H^s}^2 \longrightarrow2E\|X\|_{H^s}^2, \]

不会趋零。两边的分布律仍可能收敛,不等于这组独立样本在均方意义靠近。

原始随机Fourier场具有延伸到高频的谱,热核乘子e的负epsilon k平方压低高频并产生光滑场。下方Sobolev轴标出s小于alpha减二分之一时L2概率H s收敛,临界点倍增壳层趋于二倍log2;固定测试函数只观察频谱的一种投影。

热核让每个固定 \(\epsilon>0\) 的场变光滑;让 \(\epsilon\downarrow0\) 时,只有先选定 \(s<\alpha-1/2\) 并沿用同一噪声,才能得到这里证明的强 \(L^2(\Omega;H^s)\) 极限。

把收敛证明变成一条可以安排计算量的误差界

固定白噪声例 \(\alpha=0,s=-3/4\)。取整数 \(L\ge1\),令热时间 \(\epsilon=L^{-2}\),并在频率L处分开求和。低频利用 \(1-e^{-z}\le z\),以及 \((1+k^2)^{-3/4}\le k^{-3/2}\):

\[\begin{aligned} &2\sum_{k\le L}(1+k^2)^{-3/4} (1-e^{-k^2/L^2})^2\\ &\quad\le 2L^{-4}\sum_{k\le L}k^{5/2} \le2L^{-1/2}. \end{aligned}\]

高频直接用 \((1-e^{-z})^2\le1\) 和已经证明的尾界,得到至多 \(4L^{-1/2}\)。所以

\[\boxed{ E\|X_{L^{-2}}-X\|_{H^{-3/4}}^2 \le6L^{-1/2}.}\]

例如 \(L=100\)、\(\epsilon=10^{-4}\) 时,平方误差至多0.6;这组参数超出当前热时间滑块,适合作为离线精度安排的例子。它是保守上界,不能当作实际误差值。

若希望平方误差不超过 \(\delta>0\),这条充分条件是 \(L\ge(6/\delta)^2\)。若目标是 \(L^2(\Omega;H^{-3/4})\) 范数误差不超过 \(\eta\),则要先令 \(\delta=\eta^2\)。沿这组热时间,平方误差界按 \(\epsilon^{1/4}\) 衰减,范数界按 \(\epsilon^{1/8}\) 衰减。这样选参数,才同时控制了可见低频和全部未显示的高频。

7. 实验怎样读:频率滤波不是收敛证书

实验使用 \(\alpha\in[0,2]\)、\(s\in[-1.5,1]\)、截断 \(N\in\{8,16,32,64,128\}\) (截断滑块选择指数 \(n\),即 \(N=2^n\))与热时间 \(\epsilon\in[0.01,0.2]\);建议默认 \(\alpha=0\)、\(s=-3/4\)、\(N=16\)、\(\epsilon=0.05\)。它画出原始与热过滤后的模态贡献,并报告部分和、倍增壳层和理论尾界。

按下面顺序操作:

  1. 默认输入满足 \(s<\alpha-1/2\)。核对 \(E\|X_{16}\|_{H^{-3/4}}^2\approx3.267141087\),\(16<k\le32\) 的倍增壳层约为 \(0.282552038\)。
  2. 保持 \(\alpha=0,N=16\),把 \(s\) 改成临界值 \(-1/2\)。倍增壳层约为 \(1.354173150\);增大 \(N\) 时它趋近 \(2\log2\approx1.386294361\),不会趋零。
  3. 回到 \(s=-3/4\),减小 \(\epsilon\)。热滤波曲线逐步接近原谱,但越来越多高频重新出现;有限画布接近不等于已证明无穷级数收敛。
  4. 增大 \(\alpha\)。阈值 \(\alpha-1/2\) 向右移动,说明原场本身更平滑;这与固定 \(\epsilon>0\) 时热核带来的无限阶平滑是两件事。

默认收敛例有 \(p=3/2\)、\(C_s=1\),所以上一节尾界给

\[ E\|X-X_{16}\|_{H^{-3/4}}^2\le1. \]

这个界故意保守,但覆盖全部未显示频率。实验模态曲线和部分和只用于理解,阈值证明依赖完整级数与临界壳层。

无脚本后备:

参数 \(E\lVert X_N\rVert_{H^s}^2\) \(E\lVert X_{2N}-X_N\rVert_{H^s}^2\) 判定
\(\alpha=0,s=-3/4,N=16\) \(3.267141087\) \(0.282552038\) 位于收敛区,尾界为 1
\(\alpha=0,s=-1/2,N=16\) 不作为证书 \(1.354173150\) 临界;壳层趋于 \(2\log2\)
一般 \(s<\alpha-1/2\) 有有限极限 趋向 0 \(L^2(\Omega;H^s)\) Cauchy
一般 \(s=\alpha-1/2\) 发散 趋于 \(2\log2\) 非 Cauchy

8. 两道迁移题

题一。 取 \(\alpha=0\)。分别判断 \(X_N\) 在 \(H^{-3/4}\) 与 \(H^{-1/2}\) 中是否为 \(L^2(\Omega;H^s)\) Cauchy;对收敛情形给出平方尾误差界,对临界情形给出倍增反例。

查看答案:只差四分之一阶,结论完全不同

\(s=-3/4\) 时

\[ p=2\alpha-2s=\frac32>1, \]

所以序列收敛,且

\[ E\|X-X_N\|_{H^{-3/4}}^2 \le\frac{2}{1/2}N^{-1/2} =4N^{-1/2}. \]

\(s=-1/2\) 正好临界。此时

\[ E\|X_{2N}-X_N\|_{H^{-1/2}}^2 =2\sum_{k=N+1}^{2N}(1+k^2)^{-1/2} \longrightarrow2\log2. \]

因此它不是 Cauchy。白噪声属于每个 \(H^{-1/2-\delta}\)(以这里的 \(L^2(\Omega;H^s)\) 意义),却不属于临界 \(H^{-1/2}\)。

题二。 有人只用 \(\varphi(x)=\cos x\) 测试,发现 \(\langle X_N,\varphi\rangle\) 从 \(N=1\) 起就不变,于是宣称 \(X_N\) 在 \(H^{-1/2}\) 收敛。指出错误,并给出能直接反驳的量。

查看答案:一只低频探针看不见高频壳层

由正交性,

\[ \langle X_N,\cos x\rangle =\frac1{\sqrt2}\xi_1^c \]

对所有 \(N\ge1\) 都成立。这个测试只读取第一个余弦模,完全遗漏 \(k\ge2\)。

在 \(\alpha=0,s=-1/2\) 中,直接检查整个 Sobolev 壳层:

\[ E\|X_{2N}-X_N\|_{H^{-1/2}}^2\to2\log2\ne0. \]

因此 \(H^{-1/2}\) Cauchy 性失败。一个连续线性泛函上的收敛不能推出 Hilbert 空间范数收敛;需要控制全部方向的总能量。

9. 这一步离奇异 SPDE 还有多远?

本讲构造了一个线性 Gaussian 随机分布,证明 Fourier 截断和热核正则化在指定负 Sobolev 空间中收敛。它已经比单一测试函数配对更强,却还没有定义 \(X^2\)、Wick 幂或任何非线性方程。

进入奇异 SPDE时,困难会从“场在哪个空间”升级为“粗糙对象能否相乘、近似反项是否一致、解映射是否连续”。热核能制造光滑近似,但“每个 \(\epsilon\) 都可计算”本身不保证非线性结果随 \(\epsilon\downarrow0\) 收敛,更不保证不同正则化方案给出同一物理参数。

速查与资料

两份实 Fourier 模态 → \(H^s\) 二阶矩级数 → 阈值 \(s<\alpha-1/2\) → 临界壳层 \(2\log2\) → 热乘子 \(e^{-\epsilon k^2}\) → 同噪声下的支配收敛。回到Wick 平方时,应明确那里只证明了一个二次场的常数配对极限;本讲证明的是线性场自身在完整 Sobolev 范数中的极限。

Gaussian 随机变量、Hilbert 空间随机级数与正则化背景见 Hairer 的 Advanced Stochastic Analysis;热半群、空间白噪声与分布值 SPDE 的系统入口见作者讲义 An Introduction to Stochastic PDEs;圆周 Fourier 与 Sobolev 归一化可对照 Hairer 与 Kurniawan 的 Stochastic PDEs with multiscale structure。本讲阈值、临界壳层和尾界均由正文级数独立推出。资料核查:2026-09-08。

后续练习:分布乘积与重建从相位反例、负Sobolev单壳层极限走到准确的模型相容条件,区分重整化对象与重建定理的职责。