2025-11-15T13:40:12.030765

Optimal $L^p$-approximation of convex sets by convex subsets

Fattah, Ftouhi, Zuazua
Given a convex set $Ω$ of $\mathbb{R}^n$, we consider the shape optimization problem of finding a convex subset $ω\subset Ω$, of a given measure, minimizing the $p$-distance functional $$\mathcal{J}_p(ω) := \left(\int_{\mathbb{S}^{n-1}} |h_Ω-h_ω|^p d\mathcal{H}^{n-1}\right)^{\frac{1}{p}},$$ where $1 \le p <\infty$ and $h_ω$ and $h_Ω$ are the support functions of $ω$ and the fixed container $Ω$, respectively. We prove the existence of solutions and show that this minimization problem $Γ$-converges, when $p$ tends to $+\infty$, towards the problem of finding a convex subset $ω\subset Ω$, of a given measure, minimizing the Hausdorff distance to the convex $Ω$. In the planar case, we show that the free parts of the boundary of the optimal shapes, i.e., those that are in the interior of $Ω$, are given by polygonal lines. Still in the $2-d$ setting, from a computational perspective, the classical method based on optimizing Fourier coefficients of support functions is not efficient, as it is unable to efficiently capture the presence of segments on the boundary of optimal shapes. We subsequently propose a method combining Fourier analysis and a recent numerical scheme, allowing to obtain accurate results, as demonstrated through numerical experiments.
academic

Optimal LpL^p-approximation of convex sets by convex subsets

基本信息

  • 论文ID: 2501.00928
  • 标题: Optimal LpL^p-approximation of convex sets by convex subsets
  • 作者: Zakaria Fattah, Ilias Ftouhi, Enrique Zuazua
  • 分类: math.OC (Optimization and Control)
  • 发表时间: 2025年1月3日
  • 论文链接: https://arxiv.org/abs/2501.00928

摘要

本文研究在给定凸集 Ω⊂Rn\Omega \subset \mathbb{R}^n 中寻找具有给定测度的凸子集 ω⊂Ω\omega \subset \Omega 的形状优化问题，目标是最小化 pp-距离泛函： Jp(ω):=(∫Sn−1∣hΩ−hω∣pdHn−1)1p\mathcal{J}_p(\omega) := \left(\int_{\mathbb{S}^{n-1}} |h_\Omega-h_\omega|^p d\mathcal{H}^{n-1}\right)^{\frac{1}{p}} 其中 1≤p<∞1 \leq p < \infty，hωh_\omega 和 hΩh_\Omega 分别是 ω\omega 和固定容器 Ω\Omega 的支撑函数。论文证明了解的存在性，并显示当 p→+∞p \to +\infty 时，该最小化问题 Γ\Gamma-收敛到寻找最小化到凸集 Ω\Omega 的Hausdorff距离的凸子集问题。

研究背景与动机

问题背景

该研究源于传感器和执行器的战略放置与形状设计问题，这在涉及偏微分方程模型或纯几何模型的众多应用中至关重要。从数学角度看，可以在最优设计框架内制定大量有趣的问题，旨在识别最小化某种能量泛函的子域。

研究动机

  1. 理论挑战：Hausdorff距离的无穷范数是不可微的，这给数值和理论分析带来挑战
  2. 实际需求：需要找到既在理论上严谨又在计算上可行的近似方法
  3. 几何优化：凸集的最优近似问题在几何学和优化理论中具有基础性意义

现有方法局限性

  • 直接处理Hausdorff距离的无穷范数在数值上困难
  • 传统的基于支撑函数Fourier系数优化的方法无法有效捕捉最优形状边界上的线段
  • 缺乏对最优形状结构特性的深入理解

核心贡献

  1. 存在性证明：证明了问题 (Pp)(P_p) 解的存在性
  2. Γ\Gamma-收敛理论：建立了当 p→+∞p \to +\infty 时，LpL^p-近似问题向Hausdorff距离最小化问题的 Γ\Gamma-收敛
  3. 结构特征定理：在平面情况下，证明了最优形状的自由边界部分由多边形线组成
  4. 通用理论框架：提出了Theorem 3，为相关问题提供了充分条件
  5. 数值方法创新：结合Fourier分析和新的数值格式，有效处理边界线段问题

方法详解

任务定义

给定凸集 Ω⊂Rn\Omega \subset \mathbb{R}^n 和常数 c∈[0,∣Ω∣]c \in [0, |\Omega|]，求解： (Pp):σp:=inf⁡{Jp(ω)∣ω⊂Ω 是凸的且 ∣ω∣=c}(P_p): \quad \sigma_p := \inf\{\mathcal{J}_p(\omega) \mid \omega \subset \Omega \text{ 是凸的且 } |\omega| = c\}

理论框架

支撑函数参数化

对于凸集 Ω⊂Rn\Omega \subset \mathbb{R}^n，其支撑函数定义为： hΩ:θ∈Sn−1→sup⁡{⟨θ,y⟩∣y∈Ω}h_\Omega: \theta \in \mathbb{S}^{n-1} \to \sup\{\langle\theta, y\rangle \mid y \in \Omega\}

pp-距离泛函

Jp(ω):=∥hΩ−hω∥p=(∫Sn−1∣hΩ−hω∣pdHn−1)1p\mathcal{J}_p(\omega) := \|h_\Omega - h_\omega\|_p = \left(\int_{\mathbb{S}^{n-1}} |h_\Omega - h_\omega|^p d\mathcal{H}^{n-1}\right)^{\frac{1}{p}}

Hausdorff距离表示

J∞(ω):=dH(ω,Ω)=∥hΩ−hω∥∞=max⁡θ∈Sn−1∣hΩ(θ)−hω(θ)∣\mathcal{J}_\infty(\omega) := d_H(\omega, \Omega) = \|h_\Omega - h_\omega\|_\infty = \max_{\theta \in \mathbb{S}^{n-1}} |h_\Omega(\theta) - h_\omega(\theta)|

技术创新点

1. Γ\Gamma-收敛分析

Theorem 1：证明了泛函序列 (Jp)(\mathcal{J}_p) 当 p→+∞p \to +\infty 时 Γ\Gamma-收敛到 J∞\mathcal{J}_\infty，从而：

  • lim⁡p→+∞σp=σ∞\lim_{p \to +\infty} \sigma_p = \sigma_\infty
  • 问题 (Pp)(P_p) 解的任何聚集点都是问题 (P∞)(P_\infty) 的解

2. 边界结构刻画

Theorem 2：在平面情况下，如果 ω∗\omega^* 是问题 (Pp)(P_p) 的解，则其自由边界部分 ∂ω∗∖∂Ω\partial\omega^* \setminus \partial\Omega 是多边形线的并。

3. 通用等价性理论

Theorem 3：提供了形状泛函间等价性的充分条件，建立了不同优化问题间的联系。

数值方法

方法1：Fourier系数优化

将支撑函数表示为Fourier级数： h(θ)=a0+∑k=1N(akcos⁡(kθ)+bksin⁡(kθ))h(\theta) = a_0 + \sum_{k=1}^N (a_k \cos(k\theta) + b_k \sin(k\theta))

约束条件：

  • 包含约束：h≤hΩh \leq h_\Omega
  • 凸性约束：h′′+h≥0h'' + h \geq 0
  • 面积约束：∣ω∣=πa02+π2∑j=1N(1−j2)(aj2+bj2)=c|\omega| = \pi a_0^2 + \frac{\pi}{2}\sum_{j=1}^N (1-j^2)(a_j^2 + b_j^2) = c

方法2：严格凸性参数化

使用Bogosel 4的方法，采用离散凸性条件： hj+1+hj−1−2hjcos⁡2πN≥0h_{j+1} + h_{j-1} - 2h_j \cos\frac{2\pi}{N} \geq 0

面积近似： ∣ω∣≈π/N2−2cos⁡2πN∑j=1Nhj(hj+1+hj−1−2hjcos⁡2πN)|\omega| \approx \frac{\pi/N}{2-2\cos\frac{2\pi}{N}} \sum_{j=1}^N h_j(h_{j+1} + h_{j-1} - 2h_j\cos\frac{2\pi}{N})

实验设置

数值实现

  • 使用两种不同的参数化方法进行对比
  • 多次随机初始化以避免局部最优解
  • 选择能量最小的结果作为最终解

测试案例

  • 不同的容器形状 Ω\Omega
  • 不同的 pp 值：p∈{1,2,8}p \in \{1, 2, 8\}
  • 不同的面积比例：α∈{0.2,0.5,0.8}\alpha \in \{0.2, 0.5, 0.8\}

评价标准

  • 目标函数值的比较
  • 收敛历史分析
  • 最优形状的几何特征验证

实验结果

方法比较

实验显示方法2（严格凸性参数化）在处理包含边界线段的最优形状时明显优于方法1（Fourier系数优化）：

测试案例方法1能量值方法2能量值改进
p=10,α=0.7p=10, \alpha=0.70.9420.9133.1%
p=4,α=0.4p=4, \alpha=0.41.1851.05311.1%

收敛性分析

  • 方法2在后期收敛阶段表现更稳定
  • 能够更好地捕捉边界上的线段结构
  • 验证了理论预测的多边形边界特征

形状特征验证

数值结果证实了理论分析：

  • 最优形状的自由边界确实呈现多边形特征
  • 不同 pp 值产生不同的最优形状
  • 面积比例 α\alpha 影响最优形状的复杂度

相关工作

形状优化领域

  • 执行器最优设计：27, 28, 29 讨论了在最优控制框架下的执行器放置方法
  • 平均距离最小化：7, 8, 22 研究了子集内平均距离最小化的经典问题
  • 特征值优化：15 提供了关于腔体放置以优化微分算子特征值的综述

支撑函数理论

  • pp-距离度量：Vitale 34 和 Florian 14 研究了凸体空间上的经典度量
  • 形状优化应用：Henrot 和 Harrel 18 研究了涉及 pp-距离泛函的形状优化问题

数值方法

  • Fourier方法：传统的基于支撑函数Fourier系数优化的方法
  • 离散凸性：Bogosel 4 提出的严格离散凸性条件

结论与讨论

主要结论

  1. 理论完备性：建立了 LpL^p-近似问题的完整理论框架，包括存在性、Γ\Gamma-收敛性和结构特征
  2. 几何洞察：揭示了最优形状边界的多边形特征，为几何直觉提供了严格的数学证明
  3. 数值有效性：提出的混合数值方法能够有效处理边界线段问题

局限性

  1. 维度限制：结构特征定理（Theorem 2）仅适用于平面情况
  2. 凸性约束：分析局限于凸域和凸子域的情况
  3. 计算复杂性：问题存在多个局部最优解，需要多次随机初始化

未来方向

  1. 高维推广：将平面结果推广到高维情况
  2. 非凸情况：研究去除凸性约束后的问题
  3. 其他约束：考虑周长约束等其他几何约束
  4. Varadhan近似：探索距离函数的PDE近似方法

深度评价

优点

  1. 理论深度：提供了完整的数学理论框架，包括存在性、收敛性和结构特征
  2. 方法创新：巧妙结合了泛函分析、几何学和数值方法
  3. 实用价值：解决了传感器和执行器设计中的实际问题
  4. 数值验证：理论结果得到了充分的数值验证

不足

  1. 维度局限：主要结果局限于平面情况，高维推广仍是开放问题
  2. 凸性限制：实际应用中可能需要处理非凸情况
  3. 计算效率：对于复杂形状，数值方法的计算成本可能较高

影响力

  1. 理论贡献：为形状优化理论提供了新的分析工具和洞察
  2. 应用前景：在传感器网络、材料设计等领域有广泛应用潜力
  3. 方法论价值：提出的理论框架可推广到其他相关问题

适用场景

  • 传感器和执行器的最优布置
  • 材料结构的几何优化设计
  • 图像处理中的形状近似
  • 几何概率和随机几何问题

参考文献

论文引用了35篇相关文献，涵盖了形状优化、凸几何、数值分析等多个领域的重要工作，为研究提供了坚实的理论基础。

总体评价：这是一篇高质量的数学研究论文，在理论分析和数值方法方面都有重要贡献。论文解决了一个具有实际意义的几何优化问题，提供了完整的理论框架和有效的数值方法。尽管存在维度限制等局限性，但为相关领域的进一步研究奠定了重要基础。