Skip to content

返回学习路线

局部平均与 Fourier 核:单元验收题及解答 ​

任务与阅读分支 ​

一台仪器记录一个周期方波,却只能保留有限频率,或读取平滑后的边界数据。你要选择能控制误差的重建方式,并说明“图看起来平滑”到底保证了什么。

主线依次阅读近似恒等核、Dirichlet 核、Fejér 求和与圆盘 Poisson 核,再由弱 Lp、最大不等式走到Lebesgue 微分。Riesz–Thorin 是算子范数的复插值分支;Marcinkiewicz 与 Hilbert 变换是奇异积分分支。前两项任务检验主线的重建能力,后三项检验几乎处处恢复与端点控制;全部完成后才算掌握整个单元。

设 f(x)=sign(sin⁡x),在跳点取零,并采用 dm=dx/(2π)。另在实线上设 g=1(0,1)。

  1. 写出 S3f、σ3f 与 Pr∗f,计算前两者在 x=π/6 的值,以及第三者在 r=3/5,x=π/6 的值。哪一种在这个点超出了输入范围?
  2. 分别说明这三种近似在连续点、跳点、L2 与一致范数中的行为。不得把方波当成连续函数应用 Fejér 的一致收敛结论。
  3. 给定 h∈L1(R) 与连续紧支撑近似 hm,证明缩小区间平均几乎处处恢复 h。明确写出坏点集合的测度界,不只说“由稠密性可得”。
  4. 对 g 算出 Hg 的完整超水平集长度,求其弱 L1 拟范数,并判断它属于哪些 Lp 空间。
  5. 由最大算子的弱 (1,1) 与 L∞ 界给出一个明确 L3 常数,再计算同一个 g 的真实范数比。解释上界为什么无需等于这个比值。

解答一:三个算子读的是同一份数据 ​

方波只有奇数正弦系数,故

S3f(x)=4π(sin⁡x+13sin⁡3x).

本单元的 Fejér 平均包含 S0,S1,S2,S3,所以

σ3f(x)=3πsin⁡x+13πsin⁡3x.

在 x=π/6,sin⁡x=1/2、sin⁡3x=1,得到

S3f(π/6)=103π≈1.061033,σ3f(π/6)=116π≈0.583568.

原函数在这个点等于一。普通部分和已超过一,而 Fejér 平均仍在 [−1,1] 中。原因是核的符号,不能归因于观测数据改变了。

Poisson 平均为

(Pr∗f)(x)=4π∑j=0∞r2j+1sin⁡((2j+1)x)2j+1.

对 |z|<1,由几何级数积分,

∑j=0∞z2j+12j+1=12log⁡1+z1−z.

令 z=reix 并取虚部。由于 r<1,商的实部为正,支路可连续选取,得到

(Pr∗f)(x)=2πarctan⁡2rsin⁡x1−r2.

在 r=3/5,x=π/6,分式为 15/16,所以结果为

2πarctan⁡1516≈0.479471.

这同样没有越出输入范围,因为 Poisson 核非负且质量为一。较低的值表示较强的平滑,不表示额外数据证明了平台温度更低。

解答二:把四种收敛声明分开 ​

方波分段光滑,普通 Fourier 部分和在连续点趋向 f(x),在跳点趋向左右极限的平均零;其 L2 误差也趋零,但跳点附近仍有 Gibbs 过冲。

Fejér 与 Poisson 的核都是偶的非负单位质量核。它们在连续点恢复函数,在有左右极限的跳点恢复两侧平均,也在 L2 中恢复方波。改变一个跳点的赋值不会改变任何系数或卷积。

三者都不能在整个圆周上一致趋向这条不连续方波。每一个近似都是连续函数,连续函数的一致极限必须连续。即使给跳点取平均值,左右邻域的跳跃仍存在。若把输入换成连续周期函数,Fejér 和 Poisson 则一致收敛;普通部分和仍需额外条件。

用核检查这项差别:DN 的有符号积分为一,但绝对质量按对数增长;FN 和 Pr 的绝对质量恰为一且向零角度集中。因此它们能把输入的一致连续性转换成统一的误差界。

解答三:真正控制所有小尺度的误差 ​

令

Dh(x)=lim supr↓012r∫x−rx+r|h(y)−h(x)|dy.

对任意连续紧支撑 hm,三角不等式和连续性给

Dh(x)≤M(h−hm)(x)+|h(x)−hm(x)|.

一维三倍覆盖给最大弱界常数 C1=3。因此对 η>0,

|{Dh>2η}|≤3+1η‖h−hm‖1.

令 m→∞,得到每个固定 η 的坏点集合测度为零;再对正有理数 η 取可数并,得到 Dh=0 几乎处处。绝对振荡趋零当然推出有符号平均恢复。

这里的上确界控制了所有半径,避免把逐个半径的近似误差误写成同时误差。如果仅有 L1 收敛而没有最大不等式,上面这个步骤就缺了依据。

解答四:一个可精算的弱端点 ​

已知

Hg(x)=1πlog⁡|xx−1|.

固定 λ>0,记 a=eπλ>1。正超水平集分成两段:

{Hg>λ}=(a1+a,1)∪(1,aa−1).

它们的长度为 1/(a+1) 与 1/(a−1)。变换满足 Hg(1−x)=−Hg(x),故负超水平集长度相同。于是

|{|Hg|>λ}|=2a+1+2a−1=2sinh⁡(πλ).

因为 sinh⁡u≥u,且 sinh⁡u/u→1,

‖Hg‖1,∞=supλ>02λsinh⁡(πλ)=2π.

所以这个输出的弱 L1 控制可以精确计算。它在零和一附近只有对数奇点,对每个有限 p 都局部 p 次可积;远处却像 1/(πx),因此全空间 Lp 积分恰在 p>1 时收敛。它不属于 L1,也因端点无界而不属于 L∞。由 L2 等距还可核验 ‖Hg‖22=‖g‖22=1。

解答五:从两个端点取得一个可用保证 ​

Marcinkiewicz 的无穷端点证明给

‖Mh‖pp≤p2pC1p−1‖h‖pp.

取 p=3,C1=3,得到

‖Mh‖3≤361/3‖h‖3.

对当前 g,区间内 Mg=1,两侧各有 1/(2x) 型尾,所以

‖Mg‖33=1+2∫1∞dx8x3=1+18=98.

‖g‖3=1,真实比值为 (9/8)1/3。这个单一输入没有使覆盖和分割步骤同时达到最坏情形,所以远小于统一上界是正常的。要声称最佳算子常数,需要对全部输入求上确界,而不是只测试一个区间。

验收标准 ​

  • 三个重建公式使用同一圆周归一化,Fejér 没有漏掉 S0 或错开一位
  • 能指出一次实际过冲,并由核的非负性证明另外两种输出保持取值范围
  • 把 L2、连续点、跳点和一致收敛分别判定;不以图像平滑代替证明
  • 写出坏点测度随 ‖h−hm‖1 趋零的界,说明最大算子在证明中的必要位置
  • 从两个实区间解出 Hilbert 超水平集,得到 2/sinh⁡(πλ),并解释弱与强端点的区别
  • 能从已知端点推出可复核的 L3 常数,同时不把一个样本函数的范数比误称为最佳常数