什么是适定问题?

计算科学 术语 调理
2021-12-12 21:46:10

我正在阅读Wikipedia 页面中的“适题问题”假设,如果一个问题是适定的,它必须满足以下条件:

  • 存在解决方案
  • 解决方案是独一无二的
  • 解的行为随着初始条件不断变化

我可以理解为什么前两个条件对于一个问题来说是可取的,但我无法理解最后一个条件。

例如,如果“解决方案的行为随着初始条件不断变化” ——这不是一件坏事吗?我们是否不希望“解决方案的行为随初始条件不断变化”?

如果解决方案的行为可以随着初始条件而改变,这是否会导致问题具有“混乱和不可预测”的能力,从而显示“不良行为”并因此 -不适定

3个回答

对初始条件的连续依赖意味着如果对于每个存在一个使得 我在这里使用了 给出的定义ϵ>0δ(ϵ)>0

x0x0~δ(ϵ)x(t,x0)x(t,x0~)ϵt0=0tT
t0=0

为什么这很重要的经典例子是洛伦兹系统/吸引子,它至少表现出对初始条件的敏感依赖性。该主题的一个很好的介绍是这篇博客文章,其中显示洛伦兹常微分方程对于初始条件中看似很小的偏差表现出非常不同的行为。

的变化,您也可能会出现混沌行为因此,您对参数系统适定性的定义可能会扩展到要求对参数的连续依赖。μf(x(t),t;μ)

为什么适定问题需要这些连续依赖标准?从计算科学的角度来看,您希望信任从模拟中获得的解决方案。问题是每台计算机只能以有限的精度表示初始条件 - 无法实现低于数据类型阈值的任何差异。但是,由于您也可以测量实际物理量,但只能达到一定的精度,因此实际上您永远不会准确地达到真正的物理起点。对于没有持续依赖初始数据的系统,您基本上可以不信任您的模拟!这就是为什么例如在天气模拟中,人们在不同的初始条件下运行多个模拟并获得平均值以预测可能发生的事情。

不适定问题的一些经典示例是从解的测量中推断 PDE 的系数。举一个具体的例子,考虑泊松问题

2u=f

在域上服从 Dirichlet 边界条件前向问题是在给定密度中找到这个 PDE的解 ,它也存在于这个空间中。前向问题是适定的,因为映射是一个连续线性算子,其中是对格林函数的积分。对于好的域,我们甚至可以计算的算子范数的显式界限。Ωu|Ω=0uL2(Ω)ffGfGG

现在假设我们有一组有限的线性泛函并且我们得到了一些测量值{μ1,,μN}

ξi=μi,u+ϵi

其中是不相关的正态随机变量,均值为 0,方差为请记住,我们不知道是什么,只知道这个向量我们的工作是估计的真实值。我们可以把它变成一个最小二乘问题:找到二次泛函ϵiσiuϵξff

J(f)=12ξMGfΣ12

其中到观测值的映射这个逆问题是不适定的,因为有许多可能的解决方案,并且对数据的估计值完全不同M:L2(Ω)RNuξξf

要理解为什么正向问题是适定的而逆向问题不是适定的,有助于理解拉普拉斯算子的谱特性。我之前说过求解泊松方程是一个连续的线性运算,但实际上我在那里掩盖了很多细节。解算子不仅是连续的——它是一个紧算子,输出比输入平滑得多。相比之下,您可以将拉普拉斯算子本身视为高通滤波器,因为它倾向于放大输入信号的高波数分量。解决上述的逆问题相当于将此高通滤波器应用于嘈杂的实验数据G,这只会进一步放大噪音。当试图在没有某种形式的正则化的情况下解决问题时,这种不适定性通常以极端振荡或嘈杂的估计具体来说维基百科文章中的定义,逆运算符作为从到自身的映射并不连续。到其对偶的映射是连续的,但这实际上对我们没有任何帮助。)fL2H1H1

摆脱这种困境的一种方法是使用贝叶斯方法。通常我们对并不完全无知,我们可能拥有的任何先验信息都可以用来规范问题。如果您有兴趣了解更多这方面的信息,我真的很喜欢Parameter Estimation and Inverse Problems一书——对 Dan Doe 的回答发表评论的 Brian Borchers 是第二作者。f

我怀疑您混淆了适定优化问题(无论初始猜测如何,您都希望有一个解决方案)和适定微分方程(您希望在初始条件下是连续的)