Discussion :

[garcia57]

Bonjour,

je suis étudiant en génie civil et je viens de débuter la programmation sous matlab.
Je cherche à écrire un programme me permettant de calculer une intégrale par la méthode du trapèze. Je cherche à résoudre uniquement des intégrales de fonction linéaire d'ordre 1 type ''x+a''. je dois écrire le programme grâce à une boucle for.

rappel : I=h/2*sum(f(xk)+f(xk+1)
avec :h=(b-a)/n et xk=a+k*h

Voici le programme que j'ai réalisé :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
function I = fonction_trapeze(a,b,n,f)
    h=(b-a)/n
    I=0
    for k=0:n-1
        x=a+(k*h)
        f1=f
        x=a+((k+1)*h)
        f2=f
        I=I+f1+f2
 
    end
    I=I*h/2
end

je trouve des résultats sous forme de matrice ligne lorsque j'utilise ce programme.
Que pourrai-je changer afin de pouvoir résoudre mes intégrales ?

merci d'avance.

[phryte]

Bonjour,

Quelques corrections :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
a=0;b=10;n=10;
f= @(x) x+1;
I=fonction_trapeze(a,b,n,f)
 
% Vérification
x=a:0.1:b;
I1=trapz(x,f(x))

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
function I = fonction_trapeze(a,b,n,f)
h=(b-a)/n;
I=0;
for k=0:n-1
    x=a+(k*h);
    f1=f(x);
    x=a+((k+1)*h);
    f2=f(x);
    I=I+f1+f2;
end
I=I*h/2;

[garcia57]

Merci beaucoup, effectivement avec ces corrections le programme fonctionne.