Skip to content

第一类紧算子方程:何时有解,怎样抗噪恢复 ​

第二类方程 u−λTu=f 的高频分母趋于一;第一类方程 Tu=g 却要除以趋零的奇异值。这份练习沿用已经求好的 Green 谱,完成三个新的判断:哪些数据真的有原像,如何定义带噪恢复,以及哪条额外条件足以给出具体恢复速度。

最短入口与可选分支 ​

核心只需三站:紧算子的奇异系统 → Picard 判据与最小范数逆 → 确定性谱正则化。第一站给出完整左右坐标,第二站检查倒数级数是否合法,第三站控制同时变化的参数与噪声。

补课按缺口选择:不熟悉伴随和闭包,先看伴随的值域恒等式与投影定理;不熟悉紧性,先看紧算子与紧自伴谱定理。只想做 Green 计算,可直接复用已证明的正弦谱,不必重新推导 Green 核。

可选的结构迁移是下文 Volterra 题:它把同空间自伴特征展开改成非自伴的两端基。有限维求解保留 QR、数值秩和阈值算法的职责;核岭回归与早停则是统计风险分支。这里的噪声只有范数界,不需要给它指定概率分布。

复算一:Green 的每个平方可积数据都有原像吗 ​

设 H=L2(0,1),

Tf(x)=∫01{min(x,t)−xt}f(t)dt,en(x)=2sin⁡(nπx),Ten=σnen,σn=1π2n2.

给定 q>1/2 和 gq=∑n−qen。任务是判断精确解、写出它的系数,并单独处理 q=5/2。

答案从能量账本开始。‖gq‖2=∑n−2q<∞,所以数据合法。T=T∗ 且两核为零,余核条件自动满足。若 Tf=gq,每个系数必须是

fn=n−qσn=π2n2−q,‖f‖2=π4∑nn4−2q.

后一个和有限恰好要求 q>5/2。在这个区间,正交级数构成真解,T 的连续性允许将算子移入其范数极限,逐项乘回后得到 gq;单射又保证唯一。

q=5/2 给调和级数 π4∑1/n,无解;q=2 给常数输入系数 π2,也无解;q=3 则给

f3=π2∑n≥11nen∈L2.

对无解的 q,有限部分和 gq,N 仍在值域并收敛到 gq。因此最小残差的下确界是零,却没有取得它的输入。这回答的是值域与值域闭包的差别,不能用“有限截断都可解”略过。

复算二:同一高频噪声怎样经过三个逆 ​

现在固定 y=g3、x†=f3,观测为 yδ=y+δem。直接逆给

T†yδ=x†+π2m2δem.

令 m→∞、δm=1/(π2m2),观测误差趋零,恢复误差却始终为一。每一个观测都有原像,所以这个失败不能归因于“恰好遇到无解数据”;它直接暴露逆的不连续。

TSVD 取 τ>0,保留 σn≥τ,即 n≤Nτ=⌊1/(πτ)⌋。答案是

xτδ=∑n≤Nτ(π2n+π2n2δ1{n=m})en.

偏差只在删掉的尾部,噪声若被保留则只在头部,两者正交。因此

‖xτδ−x†‖2=π4∑n>Nτn−2+π4m4δ21{m≤Nτ}.

对任意位置的允许噪声,还可统一控制为 δ/τ 加上真解尾范数。取 τ=δ,噪声项趋零,Nτ→∞ 使平方可和尾部也趋零。这是证明,不依赖检查多少个数值频率。

Tikhonov 最小化 ‖Tx−yδ‖2+α‖x‖2。其答案为

⟨en,xαδ⟩=π2/n+π2n2δ1{n=m}1+απ4n4.

噪声放大至多 δ/(2α);偏差平方为

∑n≥1π4n2(απ4n41+απ4n4)2.

取 α=δ,前者趋零。对后者先用 ∑n>Nn−2≤1/N 控制尾部,再让有限头部的收缩残差趋零,便得到总误差趋零。两种方法的保证都是对固定 x†、整个噪声球统一成立;并未声称它们在每个有限噪声水平下都比直接逆误差小。

迁移:把 Green 换成非自伴 Volterra ​

将算子换为 Vf(x)=∫0xf(t)dt。它紧且没有非零普通特征值。任务是重新找到完整坐标,并判断常数数据一和线性数据 x 哪个有原像。

V∗g(t)=∫t1g(x)dx。对 h=V∗Vf,两次微分与端点给

h″=−f,h′(0)=0,h(1)=0.

若 h=λf、λ>0,积分正则性先使 f 成为经典解,再解 f″+λ−1f=0。混合端点迫使

ωn=(n−12)π,σn=1/ωn,vn(t)=2cos⁡(ωnt),νn(x)=2sin⁡(ωnx).

直接积分得到 Vvn=σnνn、V∗νn=σnvn,故左基为 un=νn。两个积分算子的核都为零,且该 ODE 已穷尽 V∗V 的正特征对;紧自伴谱定理于是保证两组基都覆盖 L2。不能仅从几对积分关系推断完备性。

对 y=1,左系数 yn=2/ωn,所以逆系数恒为 2,无解。对 y=x,左系数为 2(−1)n−1/ωn2,逆系数平方可和,且直接积分验证唯一解为 f=1。

结构迁移改变了左右基、奇异值衰减与边界条件,Picard 的两个条件和正则化证明却保留。与 Green 的二阶放大 π2n2 相比,Volterra 对相应左噪声方向的直接逆放大为 ωn;这项比较来自各自算出的奇异值,而不是普通特征值。

边界题:删掉哪条条件会失败 ​

一、只查加权级数,不查余核 ​

令 K(f,a)=(Tf,0),数据为 (g3,1)。加权级数有限,但最后一个输出坐标永远为零,所以精确无解。最小二乘解是 (f3,a),最小残差一,最小范数代表为 (f3,0)。因此精确逆的可解判据与最小二乘逆的定义域不能混写。

二、参数与噪声都趋零就足够吗 ​

不够。在任意无限秩奇异系统上取真解零、噪声 yδn=σnun、δn=σn→0。TSVD 若取 τn=σn,输出 vn;Tikhonov 若取 αn=σn2,输出 vn/2。两者均不收敛到零。错误在于 δn/τn、δn/αn 都恒为一。

三、有解能否自动换来 δ 的 2/3 次幂速率 ​

不能。要使用 32ρ1/3δ2/3,本题要求 x†=(K∗K)w、‖w‖≤ρ。仍取 Green,若 w=e1、ρ=1,则

x†=π−4e1,y=π−6e1.

此时偏差系数界给 ‖Qαy−x†‖≤α,噪声界给 δ/(2α)。选 α=δ2/3,相加得 32δ2/3。这是指定选择的有效上界,不是声称它是所有选择的最佳常数。

原来的 g3 则要求源向量系数 wn=π6n3,不平方可和,所以不能套这个速率。它仍有已证明的固定真解收敛,只是需要另一种尾部信息才能得到相应定量速度。

四、能否恢复核分量,或在所有单位真解上统一收敛 ​

两种恢复都没有核分量。如果真实输入为 x†+z、z∈ker⁡K,数据不变,极限仍为 x†,误差至少为 ‖z‖。此外,无限秩时可随阈值选择一个被删掉的单位 vn,所以 TSVD 的单位球偏差上确界始终为一;Tikhonov 的相应上确界也是 supnα/(σn2+α)=1。固定真解的尾部论证不能升级为整个单位球的算子范数收敛。

一个只做有限复算的程序 ​

下载 Python 标准库复算器。默认执行即可输出:指定 Green 高频噪声的直接逆误差、两种恢复的前若干系数、带解析尾界的误差范围,以及独立源条件例的实际误差和认证上界。可用 --delta、--mode、--terms 改变噪声、模态与有限求和长度;--self-test 运行包含非法输入的有限测试。

程序将未计算的 g3 尾部用 ∑n>Ln−2≤1/L 包住。TSVD 还使用下界 1/(L+1);Tikhonov 尾部只给零下界与同一上界。它不把有限求和长度当成模型维数,也不会输出“已数值证明无限维收敛”。阈值边界用浮点比较复核相邻整数,恰好相等的数学判定仍以正文为准;过小到不能分辨相邻整数的噪声参数会被显式拒绝。浮点显示有舍入误差,严格收敛、精确阈值以及界的来源都在前面的解析证明中。

参考与交卷检查 ​

主来源是 Christian Clason Regularization of Inverse Problems, arXiv:2001.00617v2,Theorems 3.9–3.10、Example 3.11、Theorem 5.6、Example 5.7 与 Theorem 6.5;三个概念页逐项标出准确印刷页码和本文补充的计算。

交卷时应能从头复算五件事:Green 的 q>5/2 与临界失败;高频噪声随 m 的放大;Volterra 的混合边界与完整左右基;参数同步失败的单位方向反例;以及源条件例中 2/3 指数与 3/2 常数。每个答案都要标明使用的空间、范数和量词,才构成可迁移的恢复保证。