从数字上看,弱收敛感觉如何?

计算科学 算法 收敛
2021-12-07 08:21:58

考虑一下,您在无限维 Hilbert 或 Banach 空间中遇到问题(想想 PDE 或此类空间中的优化问题),并且您有一个弱收敛到解决方案的算法。如果您将问题离散化并将相应的离散化算法应用于问题,那么弱收敛就是在每个坐标上的收敛,因此也是强的。我的问题是:

这种强收敛感觉或看起来与从原始无限算法的老式普通强收敛获得的收敛有什么不同?

或者,更具体的:

“离散化弱收敛方法”会发生什么样的不良行为?

当我只能证明弱收敛时,我自己通常不太高兴,但到目前为止,即使我将问题离散化到更高维度,我也无法观察到方法结果的一些问题。

请注意,我对“首先离散化而不是优化”与“首先优化而不是离散化”问题不感兴趣,并且我知道如果将算法应用于不与问题共享所有属性的离散化问题可能会出现问题该算法的设计目的。

更新:作为一个具体的例子,考虑一个变量的优化问题L2并使用(惯性)前向-后向分裂或其他仅收敛较弱的方法来解决它L2是已知的。对于离散化问题,您可以使用相同的方法,并且通过正确的离散化,如果您直接对算法进行离散化,您将获得相同的算法。提高离散化精度时会出现什么问题?

2个回答

确实,弱收敛在连续极限中是最关键的,因为h0(例如,无法观察到任何收敛速度)。至少在希尔伯特空间中,它也与极限的非唯一性密切相关,因此只有随后的收敛(例如,你可以在接近不同的极限点之间交替,再次破坏速率),并且很难区分两者就收敛。

专门针对弱收敛L2,你也有这样一个事实,即收敛不必是逐点的,你实际上可以在(足够精细的)离散化中观察到这一点。这是一系列最小化器的示例{uε}ε>0收敛为ε0

u(x)={1x<130x[13,23]1x>23
收敛较弱但不是逐点的[13,23](但几乎在其他任何地方都是逐点)。下图显示了序列中的三个代表性元素(对于ε已经很小了)。

弱收敛 1 弱收敛 2 弱收敛 3

这种现象在微分方程的 bang-bang 控制问题的近似中被称为“chittering”(即,解几乎处处都达到下限或上限的框约束问题)。

(这个具体的例子取自我们关于椭圆系统的 Multi-bang control 的论文,Ann. Henri Poincaré (C) 2014, 1109-1130, Remark 4.2。)

您提出的问题通常与实际问题无关,因为对于相同的解决方案序列,一个范数的弱收敛可能意味着另一个范数的强收敛。

举一个例子,假设我们用标准有限元在凸多边形域上用足够光滑的右手边求解拉普拉斯方程。然后解决方案uH2, 但当然是有限元解uh只在H1. 我们确实知道uhu强烈地在L2H1作为最大网格尺寸的规范h0因为我们有先验误差估计uuhL2Ch2uuhH1Ch.

但显然我们不能指望uhu强烈地H2因为uh只在H1. 但我们可能有uhu弱在H2(事实上​​,我认为这是成立的)。这可能意味着这样的陈述

2(uuh),2vo(1)vH2.

关键是弱收敛与强收敛的问题通常是您所查看的规范的问题,而不是您从方法中获得的解决方案序列的属性。