Lemme de Bramble-Hilbert
En mathématiques , et en particulier en analyse numérique , le lemme de Bramble-Hilbert , qui porte les noms de James H. Bramble et Stephen Hilbert , donne une borne à l'erreur d'une approximation d'une fonction
u
{\displaystyle \textstyle u}
par un polynôme d'ordre au plus
m
−
1
{\displaystyle \textstyle m-1}
en fonction des dérivées de
u
{\displaystyle \textstyle u}
d'ordre
m
{\displaystyle \textstyle m}
. L'erreur de l'approximation et les dérivées de
u
{\displaystyle \textstyle u}
sont mesurées par des normes
L
p
{\displaystyle \textstyle L^{p}}
sur un domaine borné dans
R
n
{\displaystyle \textstyle \mathbb {R} ^{n}}
. Le lemme est proche d'un résultat classique en analyse numérique, qui indique, par exemple, que l'erreur d'une interpolation linéaire
u
{\displaystyle \textstyle u}
peut être bornée en utilisant la dérivée seconde de
u
{\displaystyle \textstyle u}
. Cependant, le lemme de Bramble-Hilbert s'applique pour un nombre quelconque de dimensions, et pas uniquement pour une dimension, et l'erreur d'approximation et les dérivées de
u
{\displaystyle \textstyle u}
sont mesurées par des normes plus générales utilisant des moyennes, et non juste la norme de la convergence uniforme .
Des hypothèse supplémentaires sur le domaine sont nécessaires pour la validité du lemme de Bramble-Hilbert. Principalement, la frontière du domaine doit être "raisonnable". Par exemple, les domaines qui ont une pointe ou une fente avec un angle nul sont exclus. Les domaines lipschitziens sont suffisamment raisonnables. Ils comprennent les domaines convexes et les domaines avec une frontière continûment différentiable .
Le lemme de Bramble-Hilbert est principalement utilisé pour trouver des bornes de l'erreur d'interpolation de fonction
u
{\displaystyle \textstyle u}
par un opérateur qui préserve les polynômes d'ordre au plus
m
−
1
{\displaystyle \textstyle m-1}
, en fonction des dérivées de
u
{\displaystyle \textstyle u}
d'ordre
m
{\displaystyle \textstyle m}
. C'est une étape essentielle dans l'estimation des erreurs de la méthode des éléments finis . Le lemme de Bramble-Hilbert est appliqué alors sur le domaine formé d'un seul élément.
Le cas unidimensionnel
À une dimension et pour une fonction
u
{\displaystyle \textstyle u}
qui a
m
{\displaystyle \textstyle m}
dérivées sur un intervalle
(
a
,
b
)
{\displaystyle \textstyle \left(a,b\right)}
, le lemme se réduit Ã
inf
v
∈
P
m
−
1
‖
u
(
k
)
−
v
(
k
)
‖
L
p
(
a
,
b
)
≤
C
(
m
)
(
b
−
a
)
m
−
k
‖
u
(
m
)
‖
L
p
(
a
,
b
)
,
{\displaystyle \inf _{v\in P_{m-1}}{\bigl \Vert }u^{\left(k\right)}-v^{\left(k\right)}{\bigr \Vert }_{L^{p}\left(a,b\right)}\leq C\left(m\right)\left(b-a\right)^{m-k}{\bigl \Vert }u^{\left(m\right)}{\bigr \Vert }_{L^{p}\left(a,b\right)},}
où
P
m
−
1
{\displaystyle \textstyle P_{m-1}}
est l'espace de tous les polynômes d'ordre au plus
m
−
1
{\displaystyle \textstyle m-1}
.
Dans le cas où
p
=
∞
{\displaystyle \textstyle p=\infty }
,
m
=
2
{\displaystyle \textstyle m=2}
,
k
=
0
{\displaystyle \textstyle k=0}
, et
u
{\displaystyle \textstyle u}
est deux fois différentiable, le lemme signifie qu'il existe un polynôme
v
{\displaystyle \textstyle v}
de degré un tel que pour tout
x
∈
(
a
,
b
)
{\displaystyle \textstyle x\in \left(a,b\right)}
,
|
u
(
x
)
−
v
(
x
)
|
≤
C
(
b
−
a
)
2
sup
(
a
,
b
)
|
u
′
′
|
.
{\displaystyle \left\vert u\left(x\right)-v\left(x\right)\right\vert \leq C\left(b-a\right)^{2}\sup _{\left(a,b\right)}\left\vert u^{\prime \prime }\right\vert .}
Cette inégalité alors résulte de manière bien connue de l'erreur estimée pour une interpolation linéaire en choisissant
v
{\displaystyle \textstyle v}
comme l'interpolant linéaire de
u
{\displaystyle \textstyle u}
.
Énoncé
Supposons que
Ω
{\displaystyle \textstyle \Omega }
soit un domaine borné dans
R
n
{\displaystyle \textstyle \mathbb {R} ^{n}}
,
n
≥
1
{\displaystyle \textstyle n\geq 1}
, avec une frontière
∂
Ω
{\displaystyle \textstyle \partial \Omega }
et un diamètre
d
{\displaystyle \textstyle d}
.
W
p
k
(
Ω
)
{\displaystyle \textstyle W_{p}^{k}(\Omega )}
est l'espace de Sobolev de toutes les fonctions
u
{\displaystyle \textstyle u}
sur
Ω
{\displaystyle \textstyle \Omega }
avec des dérivées partielles au sens faible
D
α
u
{\displaystyle \textstyle D^{\alpha }u}
d'ordre
|
α
|
{\displaystyle \textstyle \left\vert \alpha \right\vert }
jusqu'Ã
k
{\displaystyle \textstyle k}
dans
L
p
(
Ω
)
{\displaystyle \textstyle L^{p}(\Omega )}
. Ici,
α
=
(
α
1
,
α
2
,
…
,
α
n
)
{\displaystyle \textstyle \alpha =\left(\alpha _{1},\alpha _{2},\ldots ,\alpha _{n}\right)}
est un multi-indice ,
|
α
|
=
{\displaystyle \textstyle \left\vert \alpha \right\vert =}
α
1
+
α
2
+
⋯
+
α
n
{\displaystyle \textstyle \alpha _{1}+\alpha _{2}+\cdots +\alpha _{n}}
et
D
α
{\displaystyle \textstyle D^{\alpha }}
note les dérivées
α
1
{\displaystyle \textstyle \alpha _{1}}
fois par rapport Ã
x
1
{\displaystyle \textstyle x_{1}}
,
α
2
{\displaystyle \textstyle \alpha _{2}}
fois par rapport Ã
x
2
{\displaystyle \textstyle x_{2}}
, et ainsi de suite. La seminorme de Sobolev sur
W
p
m
(
Ω
)
{\displaystyle \textstyle W_{p}^{m}(\Omega )}
consiste en les
L
p
{\displaystyle \textstyle L^{p}}
normes des dérivées d'ordre le plus élevé,
|
u
|
W
p
m
(
Ω
)
=
(
∑
|
α
|
=
m
‖
D
α
u
‖
L
p
(
Ω
)
p
)
1
/
p
si
1
≤
p
<
∞
{\displaystyle \left\vert u\right\vert _{W_{p}^{m}(\Omega )}=\left(\sum _{\left\vert \alpha \right\vert =m}\left\Vert D^{\alpha }u\right\Vert _{L^{p}(\Omega )}^{p}\right)^{1/p}{\text{ si }}1\leq p<\infty }
et
|
u
|
W
∞
m
(
Ω
)
=
max
|
α
|
=
m
‖
D
α
u
‖
L
∞
(
Ω
)
{\displaystyle \left\vert u\right\vert _{W_{\infty }^{m}(\Omega )}=\max _{\left\vert \alpha \right\vert =m}\left\Vert D^{\alpha }u\right\Vert _{L^{\infty }(\Omega )}}
P
k
{\displaystyle \textstyle P_{k}}
est l'espace de tous les polynômes d'ordre au plus
k
{\displaystyle \textstyle k}
sur
R
n
{\displaystyle \textstyle \mathbb {R} ^{n}}
. Notons que
D
α
v
=
0
{\displaystyle \textstyle D^{\alpha }v=0}
pour tous
v
∈
P
m
−
1
{\displaystyle \textstyle v\in P_{m-1}}
et
|
α
|
=
m
{\displaystyle \textstyle \left\vert \alpha \right\vert =m}
. Ainsi
|
u
+
v
|
W
p
m
(
Ω
)
{\displaystyle \textstyle \left\vert u+v\right\vert _{W_{p}^{m}(\Omega )}}
a la même valeur pour tout
v
∈
P
k
−
1
{\displaystyle \textstyle v\in P_{k-1}}
.
Lemme (Bramble et Hilbert) Sous des hypothèse supplémentaires sur le domaine
Ω
{\displaystyle \textstyle \Omega }
, spécifiées plus bas, il existe une constante
C
=
C
(
m
,
Ω
)
{\displaystyle \textstyle C=C\left(m,\Omega \right)}
indépendante de
p
{\displaystyle \textstyle p}
et
u
{\displaystyle \textstyle u}
telle que pour tout
u
∈
W
p
k
(
Ω
)
{\displaystyle \textstyle u\in W_{p}^{k}(\Omega )}
, il existe un polynôme
v
∈
P
m
−
1
{\displaystyle \textstyle v\in P_{m-1}}
tel que pour tout
k
=
0
,
…
,
m
,
{\displaystyle \textstyle k=0,\ldots ,m,}
|
u
−
v
|
W
p
k
(
Ω
)
≤
C
d
m
−
k
|
u
|
W
p
m
(
Ω
)
.
{\displaystyle \left\vert u-v\right\vert _{W_{p}^{k}(\Omega )}\leq Cd^{m-k}\left\vert u\right\vert _{W_{p}^{m}(\Omega )}.}
Le résultat original
Le lemme a été prouvé par Bramble et Hilbert [1] sous l'hypothèse que
Ω
{\displaystyle \textstyle \Omega }
satisfasse la propriété forte du cône ; c'est-à -dire, qu'il existe un recouvrement ouvert fini
{
O
i
}
{\displaystyle \textstyle \left\{O_{i}\right\}}
de
∂
Ω
{\displaystyle \textstyle \partial \Omega }
et des cônes correspondants
{
C
i
}
{\displaystyle \textstyle \{C_{i}\}}
avec des sommets à l'origine tels que
x
+
C
i
{\displaystyle \textstyle x+C_{i}}
soit contenu dans
Ω
{\displaystyle \textstyle \Omega }
pour tout
x
{\displaystyle \textstyle x}
∈
Ω
∩
O
i
{\displaystyle \textstyle \in \Omega \cap O_{i}}
.
Notes et références
(en) J. H. Bramble et S. R. Hilbert , « Estimation of linear functionals on Sobolev spaces with application to Fourier transforms and spline interpolation » , SIAM J. Numer. Anal. , vol. 7,‎ 1970 , p. 112–124
Liens externes
Cet article est issu de
wikipedia . Text licence:
CC BY-SA 4.0 , Des conditions supplémentaires peuvent s’appliquer aux fichiers multimédias.