MINISTERES DES ENSEIGNEMENTS TECHNIQUES ET
DE LA RECHERCHE SCIENTIFIQUE
DIRECTION GENERALE
DES ETUDES TECHNOLOGIQUES
INSTITUT SUPERIEUR DES ETUDES TECHNOLOGIQUES DE
SILIANA
DEPARTEMENT : GENIE ELECTRIQUE
Élaboré par : EMMANUEL AHUI
Professeur : Mme Zeinab
Classe : GE12
Année Universitaire : 2021-2022
ISET SILIANA
COMPE RENDU TP02 DE RESOLUTION NUMERIQUE :
I. Introduction
Dans ce TP, nous pouvons d’étudier les méthodes numériques pour le
calcul approché de :
Lorsqu’il s’agit d’une formule simple d’une fonction f(x), cette intégrale
peut se fait analytiquement et nous n’avons pas besoin d’utiliser des méthodes
numériques. Alors que dans le cas où la formule f(x) est compliquée ou lorsque
nous avons des mesures discrètes sans aucune formule mathématique qui relie
ces mesures, on fait recours aux méthodes numériques.
Fig 1. L’intégrale d’une fonction f
II. Méthode du point milieu
La formule classique du point milieu (ou du rectangle) est obtenue en
remplaçant f par sa valeur au milieu de l’intervalle [a, b].
Saisie sur Matlab :
Partie Editor:
Résolution numérique 2021-2022
ISET SILIANA
function [Ipmf]=Int_pt_milieu(f,a,b)
Ipmf=(b-a)*f((a+b)/2) ;
Partie Commande
>>Help inline;
>> f= inline (‘3*x^2+2*x’);
Entrer les variables a et b sur l’editor, puis saisie la fonction.
>> Ipmf=Int_pt_milieu(f,1,2) puis executer.
Ipmf=9,75 et ans=9,75
Fig 2. Formule du point milieu
III. Méthode du point milieu composite
La méthode du point milieu composite est obtenue en subdivisant
l’intervalle [a, b] en n sous-intervalles Ij=[x(j), x(j+1)], j=1,…,n.
Saisie sur Matlab :
Partie Editor :
function [Ipmcf]= pt_milieu_composite(f,a,b,n)
h=(b-a)/n;
Résolution numérique 2021-2022
ISET SILIANA
x=a:h:b;
Ipmcf=0;
for j=1:n
z(j)=(x(j+1)+x(j))/2;
Ipmcf= h*f(z(j))+Ipmcf;
end ;
Partie Command :
Entrer les variables a, b et n sur l’editor, puis saisie la fonction.
[Ipmcf]= pt_milieu_composite(f,a,b,n)
Exemple: a=1 ; b=2 et n=8 ;
[Ipmcf]= pt_milieu_composite(f,1,2,8) et appuyer sur la touche entrer.
Ans =*
9,98.
Pour n=6 ;
Ans =
9,99.
Pour n=10 ;
Ans =
9.98.
Pour n=6;
Ans =
9.99.
Résolution numérique 2021-2022
ISET SILIANA
Fig 3. Formule du point milieu composite
Remarque :
L’indice pm signifie “ point milieu “, et l’exposant c signifie “composite”.
IV. Méthode des trapèzes
Dans la methode du trapèze on joint f(x(j)) et f(x(j+1)) dans l’intervalle
[x1,x(n)]. Le calcul de l’intégrale dans ce cas revient au calcul de l’aire d’un
trapèze comme illustrer à la figure 4.
( petite base+ grande base )∗hauteur
S=
2
On donne la formule du trapèze ainsi :
n
h
I c ( f )= ∗∑ ( f ( x j ) + f ( x j+1 ) )
2 j=1
Saisie sur Matlab :
Partie Editor :
function [Itf]= Int_trapezes (f,a,b,n)
h=(b-a) /n ;
x=a:h:b;
Itf=0 ;
For j=1 : n
Résolution numérique 2021-2022
ISET SILIANA
Itf=(f(x(j)) +f(x(j+1))) *h/2+Itf ;
end ;
Partie Command :
>> help inline
>> f= inline (‘3*x^2+2*x’)
>> [Itf]= Int_trapezes (f,a,b,n)
Entrer les variables n, a et b, puis cliquez sue la touche entrer.
>> [Itf]= Int_trapezes (f,1,2,8)
Ans=
9.99
Fig 3. Méthode des trapèzes avec n=5
V. Conclusion
En somme, nous remarquons que chaque résultats obtenus avec les
différentes méthodes qui ont été effectués ont tous les même valeurs.
Résolution numérique 2021-2022
ISET SILIANA
Résolution numérique 2021-2022