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

Approximation optimale en LpL^p d'ensembles convexes par des sous-ensembles convexes

Informations de base

  • ID de l'article : 2501.00928
  • Titre : Optimal LpL^p-approximation of convex sets by convex subsets
  • Auteurs : Zakaria Fattah, Ilias Ftouhi, Enrique Zuazua
  • Classification : math.OC (Optimisation et Contrôle)
  • Date de publication : 3 janvier 2025
  • Lien de l'article : https://arxiv.org/abs/2501.00928

Résumé

Cet article étudie un problème d'optimisation de forme consistant à trouver un sous-ensemble convexe ω⊂Ω\omega \subset \Omega de mesure donnée dans un ensemble convexe fixé Ω⊂Rn\Omega \subset \mathbb{R}^n, en minimisant la fonctionnelle de distance 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}} où 1≤p<∞1 \leq p < \infty, et hωh_\omega et hΩh_\Omega sont respectivement les fonctions d'appui de ω\omega et du conteneur fixe Ω\Omega. L'article établit l'existence de solutions et montre que lorsque p→+∞p \to +\infty, le problème de minimisation Γ\Gamma-converge vers le problème de recherche du sous-ensemble convexe minimisant la distance de Hausdorff à l'ensemble convexe Ω\Omega.

Contexte et motivation de la recherche

Contexte du problème

Cette recherche provient de problèmes de placement stratégique et de conception de capteurs et d'actionneurs, qui sont essentiels dans de nombreuses applications impliquant des modèles d'équations aux dérivées partielles ou des modèles purement géométriques. D'un point de vue mathématique, de nombreux problèmes intéressants peuvent être formulés dans le cadre de la conception optimale, visant à identifier les sous-domaines minimisant une certaine fonctionnelle d'énergie.

Motivation de la recherche

  1. Défis théoriques : La norme infinie de la distance de Hausdorff n'est pas différentiable, ce qui pose des défis pour l'analyse numérique et théorique
  2. Besoins pratiques : Nécessité de trouver des méthodes d'approximation à la fois rigoureuses théoriquement et numériquement viables
  3. Optimisation géométrique : Le problème d'approximation optimale d'ensembles convexes est fondamental en géométrie et en théorie de l'optimisation

Limitations des approches existantes

  • Le traitement direct de la norme infinie de la distance de Hausdorff est numériquement difficile
  • Les méthodes traditionnelles basées sur l'optimisation des coefficients de Fourier de la fonction d'appui ne peuvent pas capturer efficacement les segments de droite sur la frontière de la forme optimale
  • Manque de compréhension approfondie des caractéristiques structurelles de la forme optimale

Contributions principales

  1. Preuve d'existence : Établissement de l'existence de solutions au problème (Pp)(P_p)
  2. Théorie de Γ\Gamma-convergence : Établissement de la Γ\Gamma-convergence du problème d'approximation en LpL^p vers le problème de minimisation de la distance de Hausdorff lorsque p→+∞p \to +\infty
  3. Théorème de caractérisation structurelle : Preuve que dans le cas plan, la partie de frontière libre de la forme optimale est composée de lignes polygonales
  4. Cadre théorique général : Proposition du Théorème 3 fournissant des conditions suffisantes pour les problèmes connexes
  5. Innovation en méthodes numériques : Combinaison de l'analyse de Fourier et de nouveaux schémas numériques pour traiter efficacement les segments de frontière

Détails des méthodes

Définition du problème

Étant donné un ensemble convexe Ω⊂Rn\Omega \subset \mathbb{R}^n et une constante c∈[0,∣Ω∣]c \in [0, |\Omega|], résoudre : (Pp):σp:=inf⁡{Jp(ω)∣ω⊂Ω est convexe et ∣ω∣=c}(P_p): \quad \sigma_p := \inf\{\mathcal{J}_p(\omega) \mid \omega \subset \Omega \text{ est convexe et } |\omega| = c\}

Cadre théorique

Paramétrisation par fonction d'appui

Pour un ensemble convexe Ω⊂Rn\Omega \subset \mathbb{R}^n, sa fonction d'appui est définie par : hΩ:θ∈Sn−1→sup⁡{⟨θ,y⟩∣y∈Ω}h_\Omega: \theta \in \mathbb{S}^{n-1} \to \sup\{\langle\theta, y\rangle \mid y \in \Omega\}

Fonctionnelle de distance 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}}

Représentation de la distance de 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)|

Points techniques innovants

1. Analyse de Γ\Gamma-convergence

Théorème 1 : Preuve que la suite de fonctionnelles (Jp)(\mathcal{J}_p) Γ\Gamma-converge vers J∞\mathcal{J}_\infty lorsque p→+∞p \to +\infty, d'où :

  • lim⁡p→+∞σp=σ∞\lim_{p \to +\infty} \sigma_p = \sigma_\infty
  • Tout point d'accumulation des solutions du problème (Pp)(P_p) est une solution du problème (P∞)(P_\infty)

2. Caractérisation de la structure de frontière

Théorème 2 : Dans le cas plan, si ω∗\omega^* est une solution du problème (Pp)(P_p), alors sa partie de frontière libre ∂ω∗∖∂Ω\partial\omega^* \setminus \partial\Omega est une union de lignes polygonales.

3. Théorie d'équivalence générale

Théorème 3 : Fournit des conditions suffisantes pour l'équivalence entre fonctionnelles de forme, établissant des connexions entre différents problèmes d'optimisation.

Méthodes numériques

Méthode 1 : Optimisation des coefficients de Fourier

Représentation de la fonction d'appui comme série de 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))

Conditions de contrainte :

  • Contrainte d'inclusion : h≤hΩh \leq h_\Omega
  • Contrainte de convexité : h′′+h≥0h'' + h \geq 0
  • Contrainte d'aire : ∣ω∣=π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

Méthode 2 : Paramétrisation par convexité stricte

Utilisation de la méthode de Bogosel 4, adoptant des conditions de convexité discrètes : hj+1+hj−1−2hjcos⁡2πN≥0h_{j+1} + h_{j-1} - 2h_j \cos\frac{2\pi}{N} \geq 0

Approximation d'aire : ∣ω∣≈π/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})

Configuration expérimentale

Implémentation numérique

  • Comparaison de deux méthodes de paramétrisation différentes
  • Initialisations aléatoires multiples pour éviter les optima locaux
  • Sélection du résultat avec l'énergie minimale comme solution finale

Cas de test

  • Différentes formes de conteneur Ω\Omega
  • Différentes valeurs de pp : p∈{1,2,8}p \in \{1, 2, 8\}
  • Différents rapports d'aire : α∈{0.2,0.5,0.8}\alpha \in \{0.2, 0.5, 0.8\}

Critères d'évaluation

  • Comparaison des valeurs de la fonction objectif
  • Analyse de l'historique de convergence
  • Vérification des caractéristiques géométriques de la forme optimale

Résultats expérimentaux

Comparaison des méthodes

Les expériences montrent que la Méthode 2 (paramétrisation par convexité stricte) surpasse nettement la Méthode 1 (optimisation des coefficients de Fourier) pour traiter les formes optimales contenant des segments de frontière :

Cas de testValeur énergétique Méthode 1Valeur énergétique Méthode 2Amélioration
p=10,α=0.7p=10, \alpha=0.70.9420.9133.1%
p=4,α=0.4p=4, \alpha=0.41.1851.05311.1%

Analyse de convergence

  • La Méthode 2 montre une convergence plus stable dans les phases finales
  • Capacité supérieure à capturer les structures de segments de frontière
  • Vérification des caractéristiques de frontière polygonale prédites théoriquement

Vérification des caractéristiques de forme

Les résultats numériques confirment l'analyse théorique :

  • La frontière libre de la forme optimale présente effectivement des caractéristiques polygonales
  • Différentes valeurs de pp produisent différentes formes optimales
  • Le rapport d'aire α\alpha influence la complexité de la forme optimale

Travaux connexes

Domaine de l'optimisation de forme

  • Conception optimale d'actionneurs : 27, 28, 29 discutent des méthodes de placement d'actionneurs dans le cadre du contrôle optimal
  • Minimisation de distance moyenne : 7, 8, 22 étudient le problème classique de minimisation de distance moyenne dans les sous-ensembles
  • Optimisation de valeurs propres : 15 fournit un aperçu du placement de cavités pour optimiser les valeurs propres d'opérateurs différentiels

Théorie de la fonction d'appui

  • Métriques de distance pp : Vitale 34 et Florian 14 étudient les métriques classiques sur l'espace des corps convexes
  • Applications à l'optimisation de forme : Henrot et Harrel 18 étudient l'optimisation de forme impliquant des fonctionnelles de distance pp

Méthodes numériques

  • Méthodes de Fourier : Approche traditionnelle basée sur l'optimisation des coefficients de Fourier de la fonction d'appui
  • Convexité discrète : Conditions de convexité discrète stricte proposées par Bogosel 4

Conclusions et discussion

Conclusions principales

  1. Complétude théorique : Établissement d'un cadre théorique complet pour le problème d'approximation en LpL^p, incluant l'existence, la Γ\Gamma-convergence et les caractéristiques structurelles
  2. Intuitions géométriques : Révélation des caractéristiques polygonales de la frontière de la forme optimale, fournissant une preuve mathématique rigoureuse pour l'intuition géométrique
  3. Validité numérique : La méthode numérique hybride proposée peut traiter efficacement les problèmes de segments de frontière

Limitations

  1. Restriction dimensionnelle : Le théorème de caractérisation structurelle (Théorème 2) s'applique uniquement au cas plan
  2. Contrainte de convexité : L'analyse est limitée aux domaines convexes et sous-domaines convexes
  3. Complexité computationnelle : Le problème possède plusieurs optima locaux, nécessitant des initialisations aléatoires multiples

Directions futures

  1. Généralisation en dimension supérieure : Extension des résultats du cas plan aux dimensions supérieures
  2. Cas non-convexe : Étude du problème après suppression de la contrainte de convexité
  3. Autres contraintes : Considération d'autres contraintes géométriques telles que les contraintes de périmètre
  4. Approximation de Varadhan : Exploration des méthodes d'approximation par EDP pour les fonctions de distance

Évaluation approfondie

Points forts

  1. Profondeur théorique : Fourniture d'un cadre mathématique complet incluant l'existence, la convergence et les caractéristiques structurelles
  2. Innovation méthodologique : Combinaison ingénieuse de l'analyse fonctionnelle, de la géométrie et des méthodes numériques
  3. Valeur pratique : Résolution de problèmes réels en conception de capteurs et d'actionneurs
  4. Validation numérique : Les résultats théoriques sont suffisamment validés numériquement

Insuffisances

  1. Limitation dimensionnelle : Les résultats principaux sont limités au cas plan, la généralisation en dimension supérieure reste un problème ouvert
  2. Restriction de convexité : Les applications pratiques peuvent nécessiter de traiter des cas non-convexes
  3. Efficacité computationnelle : Pour les formes complexes, le coût computationnel des méthodes numériques peut être élevé

Impact

  1. Contribution théorique : Fourniture de nouveaux outils d'analyse et d'intuitions pour la théorie de l'optimisation de forme
  2. Perspectives d'application : Potentiel d'application large dans les réseaux de capteurs, la conception de matériaux, etc.
  3. Valeur méthodologique : Le cadre théorique proposé peut être généralisé à d'autres problèmes connexes

Scénarios d'application

  • Placement optimal de capteurs et d'actionneurs
  • Optimisation géométrique de conception de structures matérielles
  • Approximation de forme en traitement d'images
  • Problèmes de géométrie probabiliste et géométrie aléatoire

Références

L'article cite 35 références connexes, couvrant des travaux importants dans plusieurs domaines tels que l'optimisation de forme, la géométrie convexe et l'analyse numérique, fournissant une base théorique solide pour la recherche.

Évaluation globale : Cet article est un travail de recherche mathématique de haute qualité avec des contributions importantes tant en analyse théorique qu'en méthodes numériques. L'article résout un problème d'optimisation géométrique ayant une signification pratique, fournissant un cadre théorique complet et des méthodes numériques efficaces. Malgré certaines limitations telles que les restrictions dimensionnelles, il jette une base importante pour les recherches ultérieures dans les domaines connexes.