Discussion :

[space_int]

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 :

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

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

[parp1]

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.

[space_int]

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

[Guigui_]

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 :

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

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

et on a bien la bonne valeur.

[space_int]

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

[parp1]

Et en essayant :

Code :

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

special.jv(3, 4.0)

Est ce que ca peut fonctionner?

[Guigui_]

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

Code :

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

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

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

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.