C:/Documents and Settings/Denis/Mes documents/Tex-Txt/PLC1/Ecrit/Interpolation et quadrature/quadrature.dvi

Description

• leçon - matière potentielle : sur l' interpolation

Quadrature Denis Vekemans ∗ Les notations utilisées sont celles de la leçon sur l'interpolation. Le but est d'obtenir une valeur numérique approchée de l'intégrale définie I = ∫ b a f(x)ω(x)dx où a < b et ω est une fonction donnée qui ne s'annule pas sur ]a, b[ et qui est intégrable sur [a, b]. Les méthodes du type interpolation.

• polynôme vn satisfaisant aux conditions du théorème
• formule d'interpolation
• méthode de quadrature
• ζx ∈
• abscisses d'interpolation de façon optimale
• plc1 quadrature
• ︷︷ ︸
• vn
• théorèmes
• théorème
• conditions
• condition

Sujets

• Liouville
• La Grange
• Newton
• Condition
• Théorème
• Interpolation
• Quadrature
• Méthode

Informations

• Publié par: vosom
• Nombre de lectures: 26
• Langue: Français

Extrait

Quadrature
∗Denis Vekemans
Les notations utilisées sont celles de la leçon sur l’interpolation.
Le but est d’obtenir une valeur numérique approchée de l’intégrale définie
Z b
I = f(x)ω(x)dx
a
où a<b et ω est une fonction donnée qui ne s’annule pas sur ]a,b[ et qui est intégrable sur [a,b].
Les méthodes du type interpolation.
Dans ces méthodes, l’idée est d’approcher la fonctionf par son polynôme d’interpolationP en les pointsn
d’abscisses x , x , ..., x .0 1 n
On a f =P +E .n n
D’où Z Z Zb b b
f(x)ω(x)dx = P (x)ω(x)dx+ E (x)ω(x)dx.n n
a a a
En utilisant la formule d’interpolation de Lagrange, il vient
Z Zn b bX
(n)I = L (x)ω(x)dxf(x )+ E (x)ω(x)dx.i ni
a ai=0| {z }
(n)
=A
i
RP (n) bnOn prendra A f(x ) comme approximation de I et R (f) = E (x)ω(x)dx comme erreur dei n ni=0 i a
quadrature.
Etude de la stabilité et de la convergence de quadrature.
La méthode de quadrature est stable si
nX
(n)
∀(ε ,ε ,...,ε ),∃M indépendante de n, telle que | A ε |≤M max|ε |.0 1 n i ii
i
i=0
La méthode de quadrature est convergente si
lim R (f) = 0.n
n→∞
Théorème 1
1. S’il existe M indépendant de n tel que
nX
(n)
|A |≤M,i
i=0
alors la méthode de quadrature est stable.
∗Laboratoire de mathématiques pures et appliquées Joseph Liouville ; 50, rue Ferdinand Buisson BP 699 ; 62 228 Calais
cedex ; France
1PLC1 Quadrature 2008
2. Si la méthode de quadrature est convergente sur l’espace des polynômes et s’il existe M indépendant
de n tel que
nX
(n)
|A |≤M,i
i=0
alors la méthode de quadrature est convergente.
Démonstration
1. On a
n nX X
(n) (n)
| A ε |≤ |A |max|ε |, ∀n∈N,i ii i
i
i=0 i=0
| {z }
≤M
et par conséquent, la méthode est stable.
2. On a Z nb X
(n)
I = f(x)ω(x)dx = A f(x )+R (f), ∀n∈N,i ni
a i=0
et Z nb X
(n)
p (x)ω(x)dx = A p (x )+R (p ),n n i n ni | {z }a i=0
−→0
où p ∈P .n n
Par différence, on déduit
Z nb X
(n)(f(x)−p (x))ω(x)dx = A (f(x )−p (x ))+R (f)−R (p ).n i n i n n ni
a i=0
u
Le théorème de Weierstrass nous dit qu’il existe une suite de polynômes (p ) telle que p −→f.n n
Dès lors, on a
nX
(n)
A (f(x )−p (x ))−→ 0i n ii
i=0
car
nX
(n)
| A (f(x )−p (x ))|−→ 0i n ii
i=0| {z }
nX
(n)
|A |max |f(x )−p (x )|i i n ii
i=0
| {z }
≤M
et Z b
(f(x)−p (x))ω(x)dx−→ 0,n
a| {z }
R
b
≤||f−p || ω(x)dxn ∞ a
puis
R (f)−→ 0.n

–2/12– Mathématiques6
6
6
6
PLC1 Quadrature 2008
Formules de Newton-Côtes.
Il s’agit du cas où les abscisses d’interpolation sont équidistantes :
b−a
x =a+i , ∀i∈{0,1,...,n}.i
n|{z}
=h
De plus, on choisit ω(x) = 1, ∀x∈ [a,b].
On a vu que
nX f(x )i
P (x) =v (x) ,n n ′(x−x )v (x )i ini=0
avec
nY
v (x) = (x−x ).n j
j=0
Ainsi, Z b v (x)(n) n
A = dx.i ′(x−x )v (x )i ia n
x−aOn pose le changement de variable y = .
h
On calcule
n nY Y
′ n n n−iv (x ) = (x −x ) =h (i−j) =h (−1) i!(n−i)!i i jn
j=0 et j=i j=0 et j=i
et
nYv (x)n n=h (y−j).
x−xi
j=0 et j=i
Puis,
Z nn−i n Yb−a (−1)(n)A = (y−j)dy.
i n i!(n−i)! 0
j=0 et j=i
Remarques :
(1) (1) b−a b−a1. Pour n = 1, A =A = , et I = (f(a)+f(b))+R (f).10 1 2 2
(2) (1) (2)b−a b−a b−a a+b2. Pour n = 2, A =A = , A = 4 , et I = f(a)+4f( )+f(b) +R (f).20 2 16 6 6 2
Etude de l’erreur de quadrature de Newton-Côtes.
n+1On a vu que si f ∈C ([a,b]), alors
v (x)n (n+1)E (x) = (x−x )(x−x )...(x−x )⌈x,x ,x ,...,x ⌋ = f (ζ ),n 0 1 n 0 1 n xf| {z } (n+1)!
=v (x)n
avec ζ ∈ [a,b].x
On déduit l’expression de l’erreur de quadrature
Z Zb b v (x)n (n+1)R (f) = (x−x )(x−x )...(x−x )⌈x,x ,x ,...,x ⌋ dx = f (ζ )dx,n 0 1 n 0 1 n xf| {z } (n+1)!a a
=v (x)n
avec ζ ∈ [a,b].x
–3/12– MathématiquesPLC1 Quadrature 2008
Lemme 2
n+2Si f ∈C ([a,b]), alors
(n+2)d f (ζ )x⌈x,x ,x ,...,x ⌋ = ,0 1 n fdx (n+2)!
où ζ ∈ [a,b].x
Démonstration
n+1On a vu que si f ∈C ([a,b]), alors
(n+1)f (η )x
⌈x,x ,x ,...,x ⌋ = ,0 1 n f (n+1)!
avec η ∈ [a,b].x
D’où
(n+2)f (ζ )x
⌈x,x ,x ,...,x ,x+h⌋ = ,0 1 n f (n+2)!
avec ζ ∈ [a,b].x
Et, par définition des différences divisées,
=⌈x+h,x ,x ,...,x ⌋0 1 n f
z }| {
⌈x,x ,x ,...,x ⌋ −⌈x ,x ,...,x ,x+h⌋0 1 n 0 1 nf f
⌈x,x ,x ,...,x ,x+h⌋ = .0 1 n f (x+h)−x
| {z }
=h
Et, en faisant tendre h vers 0,
(n+2)d f (ζ )x
⌈x,x ,x ,...,x ⌋ = ,0 1 n fdx (n+2)!
où ζ ∈ [a,b].x

On donne la fonction V parn Z x
V (x) = v (t)dt.n n
a
Lemme 3
Pour n pair, on a V (a) =V (b) = 0 et V (x)> 0, ∀x∈ [a,b].n n n
Démonstration
a+b a+b a+bPuisque n est pair, v est antisymétrique par rapport à (i.e. v ( +x) =−v ( −x)).n n n2 2 2
Il s’ensuit V (a) =V (b) = 0.n n
a+bToujours suite à l’antisymétrie, il ne reste plus qu’à montrer que V (x)> 0, ∀x∈ [a, ].n 2
A l’aide du changement de variable t =a+uh,
Z x−a
h
n+2V (x) =h u(u−1)...(u−n)du.n
0
–4/12– MathématiquesPLC1 Quadrature 2008
– Pour u∈]0,1[, on a u(u−1)...(u−n)> 0 et donc V (x)> 0 pour x∈ [a,x ].n 1
– Pour u∈]1,2[, on a u(u−1)...(u−n)< 0 et donc V (x)>V (x ) pour x∈ [x ,x ].n n 2 1 2
Or
V (x )n 2
Z 2
n+2= h u(u−1)...(u−n)du
0 
Z Z 
1 2 
n+2 = h u(u−1)...(u−n)du+ v(v−1)...(v−n)dv 
0 1 | {z }
R1
= (u+1)u...(u+1−n)du avec u=v−1
0
Z 1
n+2= h (2u−n+1)u...(u+1−n)du .
0
Mais, pour u ∈]0,1[, on a (2u−n + 1)u(u− 1)...(u + 1−n) > 0 et donc V (x) > V (x ) > 0 pourn n 2
x∈ [x ,x ].1 2
– On suppose que V (x)> 0 pour x∈ [x ,x ] et que V (x)> 0 pour x∈ [x ,x ] et on veutn 2m 2m+1 n 2m+1 2m+2
montrer que V (x) > 0 pour x ∈ [x ,x ] et que V (x) > 0 pour x ∈ [x ,x ] lorsquen 2m+2 2m+3 n 2m+3 2m+4
n2m+4≤ .
2
– Pour u ∈]2m + 2,2m + 3[, on a u(u− 1)...(u− n) > 0 et donc V (x)− V (x ) > 0 pourn n 2m+2
x∈ [x ,x ], puis V (x)> 0 pour x∈ [x ,x ].2m+2 2m+3 n 2m+2 2m+3
– Pour u ∈]2m + 3,2m + 4[, on a u(u− 1)...(u− n) < 0 et donc V (x) > V (x ) pour x ∈n n 2m+4
[x ,x ].2m+3 2m+4
Or
V (x )n 2m+4
Z 2m+4
n+2= h u(u−1)...(u−n)du
0 Z 2m+2
n+2= h u(u−1)...(u−n)du
0 
 Z Z 2m+3 2m+4
 n+2+h  u(u−1)...(u−n)du+ v(v−1)...(v−n)dv 
 2m+2 2m+3 | {z }
R
2m+3
= (u+1)u...(u+1−n)du avec u=v−12m+2
Z 2m+3
n+2= V (x )+h (2u−n+1)u...(u+1−n)du .n 2m+2
2m+2
Mais, pour u ∈]2m + 2,2m + 3[, on a (2u−n + 1)u(u− 1)...(u + 1−n) > 0 et donc V (x) >n
V (x )> 0 pour x∈ [x ,x ].n 2m+4 2m+3 2m+4

–5/12– MathématiquesPLC1 Quadrature 2008
Théorème 4
n+2Si f ∈C ([a,b]), alors l’erreur de la formule de quadrature de Newton-Côtes pour n pair est
Zn+3 nh (n+2) 2R (f) = f (ζ) t (t−1)(t−2)...(t−n)dt,n
(n+2)! 0
où ζ ∈ [a,b].
Démonstration
On a Z b
R (f) = v (x)⌈x,x ,x ,...,x ⌋ dx.n n 0 1 n f
a
Une intégration par parties donne
Zh i bb d
R (f) = V (x)⌈x,x ,x ,...,x ⌋ − V (x) ⌈x,x ,x ,...,x ⌋ dx.n n 0 1 n n 0 1 nf f
a dxa| {z } | {z }
(n+2)f (ζ )=0, d’après le lemme 3. x= où ζ ∈[a,b], d’après le lemme 2.x(n+2)!
Puis, Z b(n+2)f (ζ)
R (f) = V (x)dx,n n(n+2)! a
où ζ ∈ [a,b], d’après le théorème de la moyenne.
Une intégration par parties donne
Z Zb b
bV (x)dx = [(x−a)V (x)] − (x−a)v (x)dx.n n na| {z }a a
=0, d’après le lemme 3.
En posant le changement de variable x =a+uh, on déduit
Z Z Zb n n
n+2 n+3 2xv (x)dx =h (uh)u(u−1)...(u−n)dx =h u (u−1)...(u−n)dx.n
a 0 0
Et, enfin, Zn+3 nh (n+2) 2R (f) = f (ζ) t (t−1)(t−2)...(t−n)dt,n
(n+2)! 0
où ζ ∈ [a,b].

Théorème 5
n+1Si f ∈C ([a,b]), alors l’erreur de la formule de quadrature de Newton-Côtes pour n impair est
Z nn+2h (n+1)R (f) = f (ζ) t(t−1)(t−2)...(t−n)dt,n
(n+1)! 0
où ζ ∈ [a,b].
–6/12– MathématiquesPLC1 Quadrature 2008
Démonstration
On a
Z Z Zb b−h b
R (f) = v (x)⌈x,x ,x ,...,x ⌋ dx = v (x)⌈x,x ,x ,...,x ⌋ dx+ v (x)⌈x,x ,x ,...,x ⌋ dx.n n 0 1 n n 0 1 n n 0 1 nf f f
a a b−h
v (x)< 0 sur [b−h,b] et donc ne change pas de signe sur [b−h,b]. Donc,n
Z Zb (n+1) bf (ζ )1
v (x)⌈x,x ,x ,...,x ⌋ dx = v (x)dx,n 0 1 n nf (n+1)!b−h b−h
avec ζ ∈ [b−h,b].1
D’autre part,
Z b−h
v (x) ⌈x,x ,x ,...,x ⌋ dxn 0 1 n f|{z} | {z }a
⌈x,x ,x ,...,x ⌋ −⌈x ,x ,...,x ⌋=(x−b) v (x) nn−1 0 1 n−1 0 1 ff
=| {z } x−b
d= V (x)n−1dx
Z Zb−h b−h dd
= V (x)⌈x,x ,x ,...,x ⌋ dx− V (x)⌈x ,x ,...,x ⌋ dxn−1 0 1 n−1 n−1 0 1 nf fdx dxa a| {z }
=0, d’après le lemme 3.
Et, une intégration par parties donne
Z b−h d
V (x)⌈x,x ,x ,...,x ⌋ dxn−1 0 1 n−1 fdxa
Zh i b−hb−h d
= V (x)⌈x,x ,x ,...,x ⌋ − V (x) ⌈x,x ,x ,...,x ⌋ dx.n−1 0 1 n−1 n−1 0 1 n−1f f
a dxa| {z }
=0, d’après le lemme 3.
Puis, d’après le théorème de la moyenne,
Z Zb−h b−h(n+1)d f (ζ )2
V (x) ⌈x,x ,x ,...,x