Skip to content

图像出现尖点以后,哪些结论还能交付 ​

本终点先交出两份经典解证书,再交出两份不可微解的选择与唯一性证书。拟线性特征线负责沿轨道求值,Hamilton–Jacobi构造还要核对梯度与作用量;黏性测试处理尖点,比较则负责唯一性和误差范围。

不要把“某条轨道可解”“函数在几乎处处满足方程”“同初值只有一个解”写成同一个验收项。它们需要不同的证据。

进入经典特征构造前,若还不熟悉行列式,先复习它对可逆性的判据;非特征条件用它检查初始面切方向与流向能否共同给出局部坐标。

任务一:变速输运的三个回代值 ​

求解

ut+xux+u=x,u(0,x)=x2.

沿标签 ξ 的曲线为

X(t,ξ)=etξ,U(t,ξ)=e−tξ2+ξ2(et−e−t).

位置导数 Xξ=et>0,所以每个有限时间都能反解 ξ=e−tx,得到

u(t,x)=e−3tx2+x2(1−e−2t).

交付初值和残差的恒等式。再在 t=log⁡2,x=2 处给出

u=54,ut=−1,ux=78,

并核验 −1+2(7/8)+5/4=2。公开程序把 E=e−t 作为正有理变量,用 ∂tEk=−kEk 检查整个多项式恒等式,数值点只是其中一个可读输出。

作为边界任务,对 ux=0 在 y=0 给数据:g(x)=x 为什么无解,g(x)=0 为什么允许任意 u=ψ(y)、ψ(0)=0?答案应指出初始曲线与流向相切,而不是将失败归咎于ODE求解器。

任务二:相空间曲线没有合并,位置投影却坍缩 ​

固定

H(x,p)=x2+p22,g(x)=0.

给出

X=ξcos⁡t,P=−ξsin⁡t,U=−ξ22sin⁡tcos⁡t.

首先核Hamilton方程与 U˙=PX˙−H。其次核标签恒等式

Uξ=−ξsin⁡tcos⁡t=PXξ.

最后在 |t|<π/2 反解标签,得到 u=−x2tan⁡t/2,直接检查 ux=P 和 ut+H(x,ux)=0。

取 cos⁡t=3/5,sin⁡t=4/5,ξ=5,得到

X=3,P=−4,U=−6,ut=−252,H=252.

这些数值全可用有理数复算。到 t=π/2,所有 X=0,但 P=−ξ 仍不同。因此失效的是标签到位置的可逆性;不能据此声称Hamilton ODE丧失唯一性,也不能继续把多值动量当作一个经典空间梯度。

任务三:同一个二次方程选择两种不同包络 ​

下面换成自主Hamilton量 H(p)=p2/2,不再含上一任务的 x2/2 势能。对初值 g+(x)=|x| 和 g−(x)=−|x|,分别计算Hopf–Lax最小值。

正初值给

U+(t,x)={x2/(2t),|x|≤t,|x|−t/2,|x|≥t;

其唯一最优标签为 y∗=sgn(x)max(|x|−t,0)。负初值给 U−(t,x)=−|x|−t/2;x>0 时标签 x+t,x<0 时标签 x−t,x=0 时恰有两个最优标签 ±t。

在 t=1/2、x=0,1/4,1,应得到

初值 三个函数值 对应最优标签
$ x $
$- x $

所有标签都满足 |x−y∗|≤t。对于正初值,接缝处两种表达式的函数值、时间导数和空间导数吻合;对于负初值,零点不可微,必须交付接触条件:上接触的 ϕt=−1/2、p∈[−1,1] 使残差非正,下接触不存在。

再拒绝一个貌似合理的答案 V=|x|−t/2。它在 x≠0 满足PDE,初值也正确,但 ϕ=−t/2 从下接触零点,残差 −1/2<0。因此它不是黏性超解。存在一份正确初值和几乎处处零残差,仍不足以通过验收。

用半群和扰动界做第二次核验 ​

对任意 s,t>0,两段代价的最优中间点是

z=ty+sxs+t,|z−y|22s+|x−z|22t=|x−y|22(s+t).

选总时间的最优标签 y,再检查 Ssg(z)+|x−z|2/(2t)=Ss+tg(x)。公开程序不仅检查正负绝对值,还用非凸分段仿射初值,逐段精确找出所有极小点再核该等式。

若给初值加入 esin⁡x,初始一致误差至多 |e|,所以两份精确包络在每个时间的差仍不超过 |e|。若一个同类近似函数的初始误差为 1/100、全时空黏性残差界为 1/200,则在 0≤t≤3 的误差至多

1100+3200=140.

这个数来自比较定理的条件化证书;不能只凭采样到的小残差宣称条件已满足。

任务四:认证最早到达时间,而不是只画一顶帐篷 ​

在 (−1,1) 上,单位速度的零边界到达时间是

d(x)=1−|x|.

应检查所有上接触斜率和零点不存在下接触的事实,并说明负帐篷如何被一条常数下测试拒绝。Eikonal距离定理的严格缩放比较随后给出整个连续闭包解类的唯一性;不是只在帐篷形候选中挑一个。

迁移到圆环 1<|x|<3,答案为 d(x)=min(|x|−1,3−|x|)。两侧边界均参与传播,|x|=2 时发生交汇。固定速度 c=2 时,到达时间变成 d/2;若速度依赖位置,除以终点速度不再是有效推导。

任务五:不等深双井,怎样选出竞争前沿 ​

取一个新的非凸初值

g(y)=min{|y+1|, |y−1|+1/2}.

左井在 y=−1 处深到零,右井在 y=1 处只深到 1/2。两份初始折线在 y=1/4 相交;完整分段为

g(y)={−y−1,y≤−1,y+1,−1≤y≤1/4,3/2−y,1/4≤y≤1,y−1/2,y≥1.

它仍是全局1-Lipschitz,但在两井之间有向下尖点。不能把两个井的最优标签平均成一个标签,也不能只选离观察点最近的井。

记单井包络为 Φt(a)=St|⋅|(a),其公式已在任务三求出。有限个分支与取下确界可以交换,故

(1)Stg(x)=infymin{|y+1|+(x−y)22t, |y−1|+12+(x−y)22t}=min{Φt(x+1), Φt(x−1)+1/2}.

每个分支内部有唯一极小标签

y−(t,x)=−1+sgn(x+1)max(|x+1|−t,0),y+(t,x)=1+sgn(x−1)max(|x−1|−t,0).

先比较式(1)的两个值,再保留胜出分支的标签;相等时取两份标签的并集。这确实列尽全局极小点:任何达到总最小值的原标签,至少在一个分支达到该分支的最小值,而每个分支又只有上述一个极小点。

固定半个时间单位,写出全部六段 ​

在 t=1/2,两包络在 x=1/4 等值。其差随 x 单调增加,等值点附近的表达式分别为 x+3/4 与 5/4−x,故这里只有一个分支切换位置。完整答案是

S1/2g(x)={−x−5/4,x≤−3/2,(x+1)2,−3/2≤x≤−1/2,x+3/4,−1/2≤x≤1/4,5/4−x,1/4≤x≤1/2,(x−1)2+1/2,1/2≤x≤3/2,x−3/4,x≥3/2.

前三区间的最优标签依次为 x+1/2,−1,x−1/2;后三个依次为 x+1/2,1,x−1/2。在抛物段与直线段的接点,相邻标签公式给出同一点;唯独分支切换点 x=1/4 有两个不同标签

argminy(g(y)+(1/4−y)21)={−1/4, 3/4}.

两标签的初值都是 3/4,移动代价都是 1/4,所以总值都是一。居中的 y=1/4 反而给 g(1/4)=5/4,不能取代这两个极小点。较浅右井的前沿被推向右侧,这就是井深参与选择的可算结果。

一条界面有两份半群证书 ​

把总时间 1/2 分成 s=t=1/4,仍取终点 x=1/4。总最优标签 y=−1/4 和 y=3/4 分别给中间点

z=y+x2=0或1/2.

由同一个双井公式计算

S1/4g(0)=S1/4g(1/2)=78,(x−z)22(1/4)=18.

两条路径都交付 7/8+1/8=1=S1/2g(1/4)。这不只是找到两条可行路径:任务三的两段配方给出任意中间点的总下界,因此两条路径都真正最优。对应的中间初始标签仍分别是 −1/4 与 3/4,分支没有在拼接时丢失。

公开程序的 piecewise 接口可直接输入断点 [-1,1/4,1]、四段斜率 [-1,1,-1,1],并指定首断点函数值为零。将所有分数写成 Fraction 后,hopf 返回精确最小值和全部并列标签;用它核对六段、切换点以及上面两份半群证书,不必将网格上的近似最低点冒充极小值。

下载与验收 ​

下载精确复算程序和结果记录。程序使用Python标准库的Fraction,不将浮点网格最小值作为精确极小值:分段仿射初值的每一段都通过投影二次极小点求解,并保留全部并列最优标签。

交付应包含经典输运恒等式、作用量/梯度恒等式、投影有效时间、两个单井包络及尖点拒绝、双井的全部竞争分支与并列标签、半群与扰动范围、零边界距离解与唯一性解类。程序的有限样例不能替代正式页的比较和存在性证明。

回到从特征曲线到唯一前沿路线。