[SciPy] probleme scipy fonction bessel

**space_int** · 13/06/2007, 14h54

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** · 13/06/2007, 15h15

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** · 13/06/2007, 15h40

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_** · 13/06/2007, 15h53

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** · 13/06/2007, 16h00

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** · 13/06/2007, 16h04

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_** · 13/06/2007, 16h05

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.

**space_int** · 13/06/2007, 16h07

Salut parp1,

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

**space_int** · 13/06/2007, 16h22

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 :

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

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.