Discussion :

[hugo_n6k]

Bonjour,

Alors mon pb est assez simple je n'arrive pas à calculer correctement des intégrales. J'ai défini une gaussienne avec un ecart type faible (donc proche d'un dirac) :

Code :

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

1
2
 
y = 1/(0.01*sqrt(2*pi))*exp( -((x-7).^2)/(2*0.01^2) ) ;

Ensuite qd je calcule l'integrale, j'ai :

Code :

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

1
2
3
4
 
>> quad(@distrib, 5, 9)
ans =
    1.0000

ce qui me va car sigma est tres faible.

Code :

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

1
2
3
4
 
>> quad(@distrib, 0, 10)
ans =
    9.2279e-175

et donc matlab zappe completement mon 'pic'.

comment puis calculer correctement cette integrale ? j'ai deja essaye quad, quadv, quadgk...

merci d'avance.

[Mat32]

Essaie la fonction trapz().h

Intégration selon la méthode des trapèzes avec un pas de h

[hugo_n6k]

ok effectivement ca marche mais il faut bien choisir le pas... c'est ca qui me pose pb mnt bon je vais voir ce que je peux faire. merci !