Discussion :

[KR_Prog]

Bonjour,
j'aimerai saisir la formule* suivante avec SymPy:

*Oui ça se fait bien à la main, c'est plus un exemple "simple" avant de partir sur des équations plus complexes

mais les indices posent problème:

Code :

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

1
2
3
4
m, a, y, i = symbols('m a y i')
L = -(1/m) * Sum(Indexed('y',i)*log(Indexed('a',i)) + (1-Indexed('y',i))*log(1-Indexed('a',i)),(i,1,m))
Lp = diff(L, a)
print(Lp)

Résultat (faux):

Code :

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

-Sum(0, (i, 1, m))/m

Une idée?

[robinechuca]

Dans votre example, vous calculez un gradient, qui est un vecteur et non pas un scalaire tel que vous tentez de l'implementer.
Vous ecrivez

Code :

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

Lp = diff(L, a)

mais c'est un abus de langage, en toute rigueur vous voudriez ecrire

Code :

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

Lp = [diff(L, a0), diff(L, a1), ..., diff(L, an)]

Dans votre cas, les termes de la somme sont indépendant donc le plus simple c'est de ne pas ecrire la somme complète mais seulement un des membres.

Dans un cas plus general, vous devez a mon avis fixer i par une valeur numerique, et iterer a la main avec une boucle fort, et calculer independament la derivee pour chaque composante. Je ne suis pas sur qu'on puisse faire mieux.

PS: une astuce de syntaxe qui peux etre utile:

>>> sympy.symbols('i:5')
(i0, i1, i2, i3, i4)
>>>