Calcul d'intégrale Octave/MATLAB

Bonjour,

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

Code:

---

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

---

sous octave, j'obtiens :

Citation:

---

ans = 7.1853e-07

---

et sous matlab, j'obtiens :

Citation:

---

??? 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.

Citation:

---

Envoyé par help quad

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:

---

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...

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

Code:

---

hmin=eps(b-a)/1024

---

par

Code:

---

hmin=eps(b-a)/1e100

---

Ca permet d'obtenir un résultat avec

Code:

---

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é.

Citation:

---

Envoyé par Padri

sous octave, j'obtiens :

Citation:

---

ans = 7.1853e-07

---

---

Code:

---

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

---

Code:

---

1 2 3

ans =   7.1853e-007

---

Tiré de la documentation de QUAD :

Citation:

---

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.

---

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.

Effectivement, la fonction quadgk fonctionne bien.
Merci Dut.

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 ?

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

Citation:

---

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