[SciPy] probleme scipy fonction bessel

bonjour,

j'ai un problême avec les fonctions de bessel, qui ne donnent pas le bon résultat.
mon code est le suivant:

Code:

---

1 2 3 4 5

from numpy import* from scipy import special   y=special.jv(3,4) print y

---

à l'exécution, le résultat est y=0.3971, alors que le bon résultat est y=0.4302

volia, je ne sais si j'ai fais une erreur dans mon programme ou s'agit il d'un bug

Merci

Tu peux me dire ce qu'est la fonction de bessel?

MErci, moi sur numpy j'avais aussi des petits soucis d'interprétations. Alors j'ai reprogrammé certain petits trucs.

Merci pour ta réponse,
les fonctions de Bessel sont très utilisées en modélisation électromagnétique, tu peux voir http://fr.wikipedia.org/wiki/Fonction_de_Bessel

et moi je les utilise dans la résolution numérique des équations de Maxwell .
et si ça donne de mauvais résultats coimme celui que j'ai cité, alors là c'est mal parti

Merci

Peut-être que c'est une valeur qu'il ne calcule pas et qu'il connaît par défaut (donc l'erreur serait là)
parce que j'ai testé:

Code:

---

1 2	>>> special.jv(3, 3.9999999) 0.430171469726

---

et on a bien la bonne valeur.

Merci Guigui,

peut-être? mais si on veut faire un algorithme qui fait des milliers de calculs, on ne peut pas à chaque faire un test pour savoir si c'est une valeur qu'il connait ou qu'il calcule. Je trouve qu'ici l'erreur est vraiment grande comme même.

Merci

Et en essayant :

Code:

---

special.jv(3, 4.0)

---

Est ce que ca peut fonctionner?

Ce serait pas plutôt la fonction jn que tu dois utiliser ?

Code:

---

1 2	>>> special.jn(3,4) 0.430171473876

---

En parcourant le fichier info.py du répertoire d'installation de scipy:

Citation:

---

Bessel Functions

jn -- Bessel function of integer order and real argument.
jv -- Bessel function of real-valued order and complex argument.
jve -- Exponentially scaled Bessel function.
yn -- Bessel function of second kind (integer order).
yv -- Bessel function of the second kind (real-valued order).
yve -- Exponentially scaled Bessel function of the second kind.
kn -- Modified Bessel function of the third kind (integer order).
kv -- Modified Bessel function of the third kind (real order).
kve -- Exponentially scaled modified Bessel function of the
third kind.
iv -- Modified Bessel function.
ive -- Exponentially scaled modified Bessel function.
hankel1 -- Hankel function of the first kind.
hankel1e -- Exponentially scaled Hankel function of the first kind.
hankel2 -- Hankel function of the second kind.
hankel2e -- Exponentially scaled Hankel function of the second kind.

---

Salut parp1,

Non, ça donne toujours le même résulat (ie 0.3671).

Merci Guigui,

Oui, je crois que la bonne fonction est special.jn, ça donne le bon résultat.
en faisant une petite comapraison entre special.jv et special.jn, je trouve qu'elles donnent les mêmes résultats sauf pour la valeur (3,4). alors je ne vois pas la différence entre jn et jv

Code:

---

1 2 3 4 5 6 7 8 9 10 11

from numpy import* from scipy import special from pylab import* x=arange(0,20,0.01) for k in range(1,10): y1=special.jn(k,x) y2=special.jv(k,x) y=y1-y2 plot(y)   show()

---

on voit un pic à x=4 et 0 ailleurs.

Quelle version de scipy ?

J'ai fait les tests avec numpy 1.0.3 et scipy 0.5.2

Bonjour à tous,

J'utilise python 2.5, numpy 1.0.1 et scipy 0.5.2,

en tous cas , special.jn fait bien ce que je veux

un grand Merci à tous:yaisse2: