0% ont trouvé ce document utile (0 vote)
113 vues4 pages

Calcul d'intégrales par méthodes numériques

Ce document décrit le principe du calcul d'intégrales définies sur un intervalle [a,b] en utilisant la méthode des points milieux. Il présente la formule du point milieu pour une fonction constante par morceaux, puis généralise la formule à n+1 points régulièrement répartis. La preuve de la formule est détaillée et un algorithme est proposé.

Transféré par

Pagnol Nouga
Copyright
© © All Rights Reserved
Nous prenons très au sérieux les droits relatifs au contenu. Si vous pensez qu’il s’agit de votre contenu, signalez une atteinte au droit d’auteur ici.
Formats disponibles
Téléchargez aux formats PDF, TXT ou lisez en ligne sur Scribd
0% ont trouvé ce document utile (0 vote)
113 vues4 pages

Calcul d'intégrales par méthodes numériques

Ce document décrit le principe du calcul d'intégrales définies sur un intervalle [a,b] en utilisant la méthode des points milieux. Il présente la formule du point milieu pour une fonction constante par morceaux, puis généralise la formule à n+1 points régulièrement répartis. La preuve de la formule est détaillée et un algorithme est proposé.

Transféré par

Pagnol Nouga
Copyright
© © All Rights Reserved
Nous prenons très au sérieux les droits relatifs au contenu. Si vous pensez qu’il s’agit de votre contenu, signalez une atteinte au droit d’auteur ici.
Formats disponibles
Téléchargez aux formats PDF, TXT ou lisez en ligne sur Scribd

1

GROUPE 1
NOMS DES MEMBRES :

NOUGA FAYO PAGNOL VANITE 20S49050


NOUGA FAYO PAGNOL VANITE 20S49050
NOUGA FAYO PAGNOL VANITE 20S49050
NOUGA FAYO PAGNOL VANITE 20S49050
NOUGA FAYO PAGNOL VANITE 20S49050
NOUGA FAYO PAGNOL VANITE 20S49050

EXAMINATRICE : DR TAPI MYRIAM

I-PRINCIPE
Il s’agit ici, en considerant une fonction définie sur [a,b] de
déterminer l’aire de la surface délimitée par les droites d’equation
x=0; x=b; y=0 et la courbe. Autrement dit, sous l’hypothèse
que f soit continue sur [a,b] donc RIEMANN-integrable sur cette
Rb
dernière calculer a f (x)dx . Il existe plusieurs formules qui perme-
ttent de calculer cette integrale notamment:
- Changement de variables
- Integration par parties,...

Nous pouvons noter des fonctions dont on sont calculées l’integrale


mais juste le fait éde chamger légèrement
R b l’expression de f suffirait
pour ne plus etre à meme de calculer a f (x)dx via les methodes
présentées. Cependant, meme si le calcul de cette integrale est
possible, il peut s’averer qu’il soit long. Il peut etre avantageux
de remplacer par un calcul approché à l’aide d’un certain nom-
bre de methodes notamment celle du point milieu utilisée lorsque
l’expression de la fonction à integrer n’est pas connue en certains
points de son domaine de définition.
2

0.1 Principe
0.1.1 En deux points
Soit f une fonction definie
 sur[a,b], supposons qu’ele soit constante alors on a :
Z b
a+b a+b
f (α)dx = (b−a)f où α = centre de l intervalle d’integration.
a 2 2
Z b  
a+b
On obtient alors f (x)dx = (b − a)f Posons h=b-a; on obtient
a 2
h
 
Z b 
2a + h
 a+
alors la formule du point milieu f (x)dx ≃ hf = hf 
 2=
a 2 2
h
 
b−
hf 
 2
2

Etant donnée que l’expression de cette surface(integrale) n’est qu’une ap-


Z b
proximation alors, ∆E l’erreur commise en calculant f (x)dx et supposons
a
′′
que f ∈ C 2 ([a, b]) et posons
 parla suite3 M2 = Sup|f (x)|x ∈ [a, b] on a donc
Z b
a+b h
∆E= | f (x)dx − hf | ≤ M2
a 2 24

0.1.2 En n+1 points régulièrement repartis


Supposons maintenant que notre fonction f soit définie sur [a,b] de grande am-
b−a
plitude. Soit h¿0 donné par h = n étant le nombre de sous-intervalles
n
de subdivision de meme amplitude h sous la forme [xn , xn+1 ]
Z b Z xn Z x1 Z x2 Z xn
f (x)dx = f (x)dx = f (x)dx + f (x)dx + ... + f (x)dx
a  x 0  x 0  x1 x n−1
 
x1 + x0 x1 + x2 xn−1 + xn
= (x1 −x0 )f + (x2 −x1 )f + ... + (xn −xn−1 )f
" n−1  2 # 2 2
  
X h h
=h f xn + +f a+
k=1
2 2
0.1. PRINCIPE 3

Pour estimer l’erreur de cette formule, supposons que f ∈ C 2 ([a, b]) et


posons par la suite M2 = Sup|f ′′ (x)|x ∈ [a, b] alors en appliquant la formule
de l’erreur pour n+1 points
" n−1donc n sous intervalles et en#sommant, on trouve
Z b    
X h h 1
: ∆E= | f (x)dx − h f xn + +f a+ | ≤ h3 M2
a k=1
2 2 24

0.1.3 Preuve :
En applicant la formule de Taylor-Lagrange à l’ordre 2 au voisinage de
xi+1 + xi
à un sous-intervalle de [a,b] étant [x+1 , xi ] on a :
2

     2
xi+1 + xi ′ xi+1 + xi xi+1 + xi 1 ′′ xi+1 + xi
f (x)dx = f +f x− + f (εi ) x −
2 2 2 2! 2
Z b n−1 xi+1 n−1 2
1 X xi+1 ′′
Z   Z  
X xi+1 + xi xi+1 + xi
=⇒ f (x)dx = f dx + f (εi ) x − dx
a i=0 x i
2 2! x i
2
Z xi+1   Z xi=0     xi+1
xi+1 + xi i+1
xi+1 + xi 1 xi+1 + xi
car x− dx=0 en effet, on a x− dx= x−
xi 2 xi 2 2 2 xi
0
Z b n−1   n−1 Z xi+1  2
X xi+1 + xi 1 X ′′ xi+1 + xi
=⇒ f (x)dx = h f dx + f (εi ) x− dx
a i=0
2 2! i=0 xi 2

en appliquant le theoreme de l’égalité de la moyenne, ∃αi ∈]a, b[ tel que


Z xi+1
 2 Z xi+1  2
xi+1 + xi ′′ xi+1 + xi
x− dx = f (αi ) x− dx
xi 2 xi 2

Z b n−1 n−1
"  3 #xi +1
X 1 X ′′ 1 xi+1 + xi
=⇒ f (x)dx = h f (xi )dx + f (αi ) x−
a i=0
2! i=0 3 2
xi
′′
Or αi in]a, b[⊆ [a, b] ⇒ αi ∈ [a,b] =⇒ f (αi ) est entre f”(a) et f”(b) et
n−1
X
′′ ′′
d’après le T.V.I, ∃α ∈ [a,b], f (αi ) = f (α) et on a f ′′ (αi ) = nf ′′ (α)
i=0
Z b n−1 3
X h
=⇒ f (x)dx = h f (xi )dx + nf ′′ (α)
a i=0
24
4

b n−1
h3 ′′ h3
Z X
=⇒ | f (x)dx - h f (xi )dx—= nf (α) ≤ nM2
a i=0
24 24

0.1.4 Exemple
exemple

0.2 ALGORITHME
1) Lire(a,b) les bornes de l’intervalle
2) Lire(n) le nombre de sous intervalles
b−a
3) h ←−
n
h
4)x ←− a +
2
h
5) I ←− f (a + )
2
6) Pour k allant de 1 à n-1 faire
∗x ←− x + h
∗I ←− I + f(x)
7)I ←− I ∗ h
8) Afficher I
9) fin pour

Vous aimerez peut-être aussi