Discussion :

[elglantosimpatico]

Bonjour,

Je voudrais intégrer la forme physique des polynômes d'Hermite dans mon code. Du coup, j'ai regardé la documentation de scipy à ce sujet (ici) mais pas moyen de trouver simplement le calcul d'un polynôme. Par exemple, "hermval(x, c[, tensor])" calcule la série d'Hermite avec un coéficient devant chaque polynôme mais moi, je voudrais simplement pouvoir appeler le polynôme d'ordre n et le calculer pour un vecteur de type sp.linspace. Du coup, suis-je bigleux ou bien n'y a-t-il pas le calcul tout simple inclus dans le package "numpy.polynomial.hermite" ?

[__dardanos__]

Salut,
Je ne comprends pas ce qui bloque. Il y a la méthode linspace :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
#!/usr/bin/python3
# -*- coding: utf-8 -*-
 
import numpy as np
import numpy.polynomial.hermite as ph
import matplotlib.pyplot as plt
 
if __name__ == '__main__':
    n = 6
    for i in range(1, n):
        c = np.zeros(n) ; c[i] = 1
        p = ph.Hermite(c)
        plt.plot(*p.linspace(), label=str(i))
    plt.legend()
    plt.show()

[elglantosimpatico]

Merci pour la réponse.

J'ai peut-être mal compris le fonctionnement de cette méthode mais n'y a-t-il pas un moyen plus simple ?
Plus précisément, j'ai une matrice x de nombre réels et je voudrais calculer les valeurs du polynôme pour toutes les valeurs de x avec un truc assez simple. Quelque chose comme cela :

Code :

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

1
2
3
4
axis = sp.linspace(-length/2,length/2,nb_points)
x, y = sp.meshgrid(axis,axis)
n = 3
result = hermite(n,x)

où n est l'ordre du polynôme et x contient les valeurs à partir desquelles calculer le polynôme.

Après, je suis peut-être encore un peu trop habitué à matlab (mais je me soigne) où cette fonction était directement intégrée comme cela.

[__dardanos__]

Ca n'existe pas tel quel, mais tu peux créer la fonction :

Code :

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

1
2
3
def hermite(n, x):
    c = np.zeros(n) ; c[-1] = 1
    return ph.hermval(x, c)

[elglantosimpatico]

Merci pour votre réponse.

En fait, j'ai finalement trouvé la fonction scipy.special.eval_hermite qui fait cela. Je ne cherchais donc pas au bon endroit.