SymPy et dérivée partielle avec index

Bonjour,
j'aimerai saisir la formule* suivante avec SymPy:
Pièce jointe 611144
*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:

---

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:

---

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

---

Une idée?

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:

---

Lp = diff(L, a)

---

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

Code:

---

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. :calim2:

PS: une astuce de syntaxe qui peux etre utile:

Citation:

---

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

---