Discussion :

[Padri]

Bonjour,

comment se fait-il que lorsque je tape la commande :

Code :

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

quad('exp(-x*x/2)',-inf,-5)

sous octave, j'obtiens :

ans = 7.1853e-07

et sous matlab, j'obtiens :

??? Error using ==> inlineeval at 15
Error in inline expression ==> exp(-x*x/2) ???
Error using ==> mtimes
Inner matrix dimensions must agree.
Error in ==> inline.subsref at 27
INLINE_OUT_ = inlineeval(INLINE_INPUTS_, INLINE_OBJ_.inputExpr, INLINE_OBJ_.expr);
Error in ==> quad at 77 y = f(x, varargin{:});

?

Je crois que ça ne vient pas de la fonction quad mais du calcul de exp(-x*x/2).

Merci d'avance.

[Caro-Line]

Use array operators .*, ./ and .^ in the definition of FUN
so that it can be evaluated with a vector argument.

Donc je dirais qu'il faut mettre :

Code :

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

quad('exp(-x.*x/2)',-inf,-5)

Ceci dit ça ne donne quand même pas le même résultat qu'avec Octave.
J'ai l'impression que MATLAB n'aime pas trop le -inf...

[Padri]

Bon, j'ai modifié le fichier quad.m. Dans la fonction quad, j'ai remplacé

Code :

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

hmin=eps(b-a)/1024

par

Code :

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

hmin=eps(b-a)/1e100

Ca permet d'obtenir un résultat avec

Code :

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

quad('exp(-x.*x/2)',-1e100,-5)

mais ça reste loin du résultat que j'obtiens avec octave.
Et donc le -inf ne marche pas sous matlab, ça va être très pratique!

En tout cas, merci de m'avoir aidé.

[Jerome Briot]

sous octave, j'obtiens :

ans = 7.1853e-07

Code :

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

quadgk(@(x)exp(-x.*x/2),-inf,-5)

Code :

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

1
2
3
ans =
 
  7.1853e-007

Tiré de la documentation de QUAD :

The quadgk function may be most efficient for high accuracies and oscillatory integrands. It supports infinite intervals and can handle moderate singularities at the endpoints. It also supports contour integration along piecewise linear paths.

[Padri]

J'ai résolu mon problème.

En fait, je voulais calculer la fonction de répartition de la loi normale centrée réduite N(0,1) et l'intégrale venait de la fonction erf (fonction d'erreur). Ce que je ne savais pas c'est que matlab et octave intègrent une fonction erf (qui a l'air d'être calculé analytiquement). Donc c'est bon, en plus les deux logiciels me donnent exactement le même résultat.

Merci.

[Padri]

Effectivement, la fonction quadgk fonctionne bien.
Merci Dut.

[Padri]

Bonjour,

juste une dernière question. La fonction quadgk n'existe pas dans les anciennes versions de matlab. Existe-t-il un équivalent pour les anciennes versions ?

[Padri]

C'est bon, j'ai trouvé. J'ai posté trop vite.

Matlab 7.5 introduced a vectorized quad ('quadgk.m') that is much much better. However, if you don't have version 7.5, you can download an alternative called 'quadva.m' written by Lawrence Shampine of SMU.

source : http://www.ldeo.columbia.edu/~spk/Co...nvolveWithGv.m