Iset Siliana 2222
Iset Siliana 2222
Iset Siliana 2222
Réaliser par :
Lakhdar mohamed moemen
Labidi aziz
1
Plan :
1) fonction matlab qui s’appelle fct calcules
des images de x note y en utilisant la
fonction :
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
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
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) 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 =
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