Iset Siliana 2222

Télécharger au format docx, pdf ou txt
Télécharger au format docx, pdf ou txt
Vous êtes sur la page 1sur 10

Institut supérieur des études technologiques de siliana

Tp2 : Les méthodes numériques


d’intégration(rectangle, trapèze)

Groupe génie électrique 13 grp 2

Réaliser par  :
Lakhdar mohamed moemen

Labidi aziz

Dirigé par : Mme Zaineb Rayouf

1
Plan :
1) fonction matlab qui s’appelle fct calcules
des images de x note y en utilisant la
fonction :

F(x)=3x^2+2X POUR x=0 :0.5 :5 avec graphe

2) méthode du point milieu

3) méthode du point milieu composite

4) méthode des trapèzes

5) application

2
Fonction matlab qui s’appelle fct calcules des
images :
function[y] = fct(x)
y= 3*x.^2+2*x
plot(x,y)

>> x=0:0.5:5

x=

Columns 1 through 9

0 0.5000 1.0000 1.5000 2.0000 2.5000 3.0000 3.5000 4.0000

Columns 10 through 11

4.5000 5.0000

>> y=fct(x)

y=

Columns 1 through 9

3
4
5
Méthode du point milieu
Function Ipm=point_milieu (a,b,f)
Ipm(f)= (b-a)f(a+b/2)
end

Méthode du point milieu composite


function[Ipmc]=point_milieu_composite(a,b,f,n)
h=(b-a)/n;
Ipmc=0;
x=a:h:b;
for j=1:n
xb=(x(j)+x(j+1))/2;
Ipmc=h*f(xb)+Ipmc;
end

Méthode des trapèzes


function[It]=trapeze(a,b,n,f)

h=(b-a)/n;

It=0;

x=a:h:b;

for j=1:n

It=h/2*(f((xj))+f((xj+1)))+It;

end

6
Application

1) Soit a intégrer la fonction f(x)= 3x^2+2x


dans l’intervalle [1,2].

Créer une fonction point milieu qui


calcule l’intégrale d’une fonction
functionIpm=point_milieu (a,b,f)
Ipm= (b-a)*f((a+b)/2)
End

1) a=1
2)
3) a =
4)
5) 1
6)
7) >> b=2
8)
9) b =
10)
11) 2
12) >> f= inline('3*x^2+2*x')
13)
14) f =
15)
16) Inlinefunction:
17) f(x) = 3*x^2+2*x
18)
19) >>Ipm=point_milieu (a,b,f)
20)
21) Ipm =
22)
23) 9.7500

7
2) Créer une fonction point_milieu
composite qui calcule l’integrale d’une
fonction . Pour N=4 , N+6 ,N=8 ,N=10

- Pour N=4
function[Ipmc]=point_milieu_composite(a,b,f,n)
h=(b-a)/n;
Ipmc=0;
x=a:h:b;
for j=1:n
xb=(x(j)+x(j+1))/2;
Ipmc=h*f(xb)+Ipmc;
end

a=1

a =

>> b=2

b =

>> f=inline('3*x^2+2*x')

f =

Inline function:
f(x) = 3*x^2+2*x

>> n=4

n =

>> [Ipmc]=point_milieu_composite(a,b,f,n)

Ipmc =

9.9844

8
3) Créer une fonction Trapèze qui calcule
l’intégrale d’une fonction .

function[It]=trapeze(a,b,n,f)

h=(b-a)/n;

It=0;

x=a:h:b;

for j=1:n

It=h/2*(f((xj))+f((xj+1)))+It;

end
>> f=inline('3*x^2+2*x')

f =

Inline function:
f(x) = 3*x^2+2*x

>> n=4

n =

>> [It]= trapeze (a,b,f,n)

9
4) comparer les resultats obtenus avec
differants methodes :

n 4 6 8 10
Methode 9,9844 9,9931 9,9961 9,9575
Point
milieu
composite
trapeze

10

Vous aimerez peut-être aussi