[论文解读] Convex integration and infinitely many weak solutions to the Perona-Malik equation in all dimensions

Seonghak Kim, Baisheng Yan|arXiv (Cornell University)|Mar 2, 2015

Numerical methods in inverse problems参考文献 15被引用 3

一句话总结

本文在所有维度的光滑有界凸域上，无论初始数据是否为亚临界、超临界或跨临界，均证明了Perona-Malik方程存在无穷多组Lipschitz弱解。作者采用一种新颖的凸积分方法，并结合Baire范畴法，证明解集在合适的函数空间中为残集，从而解决了图像处理PDE领域中长期存在的高维开放问题。

ABSTRACT

We prove that for all smooth nonconstant initial data the initial-Neumann boundary value problem for the Perona-Malik equation in image processing possesses infinitely many Lipschitz weak solutions on smooth bounded convex domains in all dimensions. Such existence results have not been known except for the one-dimensional problems. Our approach is motivated by reformulating the Perona-Malik equation as a nonhomogeneous partial differential inclusion with linear constraint and uncontrollable components of gradient. We establish a general existence result by a suitable Baire's category method under a pivotal density hypothesis. We finally fulfill this density hypothesis by convex integration based on certain approximations from an explicit formula of lamination convex hull of some matrix set involved.

研究动机与目标

• 解决高维（大于一维）Perona-Malik方程弱解存在性的长期开放问题。
• 将先前局限于一维或径向对称情形的结果推广至一般光滑有界凸域。
• 为所有光滑非恒定初始数据（包括跨临界情形）建立无穷多组Lipschitz弱解的存在性。
• 证明弱解的梯度最大值原理，表明梯度始终被初始梯度的可控倍数所控制。

提出的方法

• 将Perona-Malik方程重新表述为带有线性约束和不可控梯度分量的非齐次偏微分包含式。
• 在合适的函数空间中应用Baire范畴法，以建立解的残集集合的存在性。
• 利用凸积分技术，通过在小立方体和区间上拼接振荡校正项来构造近似解。
• 通过源自PDE包含式中涉及的矩阵集的显式公式推导出的层状凸包逼近，建立密度假设。
• 使用具有受控时间和空间导数的分片$ C^1 $函数构造扰动项，以满足弱形式。
• 验证所构造的解满足弱形式，并在$ L^1 $-范数下逼近Perona-Malik函数$ \sigma(p) = \frac{p}{1+|p|^2} $。

实验结果

研究问题

• RQ1 即使初始数据为跨临界情形，Perona-Malik方程在高维（大于一维）下是否仍可能存在无穷多组弱解？
• RQ2 能否在超越一维或径向对称设定的背景下，利用凸积分技术构造此类解？
• RQ3 尽管原始PDE存在不适定性，解集在弱拓扑下是否仍保持在合适函数空间中为稠密或残集？
• RQ4 对于此类弱解，是否仍可保持梯度最大值原理，即使经典解不存在？
• RQ5 对于解的存在性而言，区域的何种结构性质（如凸性）是本质性的？

主要发现

• 在任意$ \mathbb{R}^n $中的光滑有界凸域上，Perona-Malik方程的初值-Neumann边界问题对所有$ n \geq 1 $均存在无穷多组Lipschitz弱解。
• 该结果适用于所有光滑非恒定初始数据，包括具有亚临界与超临界区域的跨临界情形。
• 弱解满足几乎梯度最大值原理：对任意$ \lambda > 0 $，有$ \|Du\|_{L^\infty(\Omega_T)} \leq \|Du_0\|_{L^\infty(\Omega)} + \lambda $，从而保证梯度有界。
• 解集在合适的函数空间中为残集，意味着此类解在拓扑意义上具有普遍性。
• 该证明依赖于一种新型凸积分方案，通过在有限个小区集上拼接振荡校正项来构造解，从而克服了高维中梯度控制缺失的问题。
• 关键技术突破在于显式构造了涉及PDE包含式的矩阵集的层状凸包，该构造满足Baire范畴论证所需的密度假设。