Yo !

Supposez X (nxn) symétrique et on definit la notation A : B = tr(A^TB) le produit de Frobenius de 2 matrices

Soit y un vecteur

Comment montrer que y^T X^{-1} y = yy^T : X ?

C'est ce qui est utilisé dans la réponse ici
https://math.stackexchange.com/questions/2992119/derivative-of-log-marginal-likelihood

Le 24 septembre 2022 à 20:08:01 :
Yo !

Supposez X (nxn) symétrique et on definit la notation A : B = tr(A^TB) le produit de Frobenius de 2 matrices

Soit y un vecteur

Comment montrer que y^T X^{-1} y = yy^T : X ?

C'est ce qui est utilisé dans la réponse ici
https://math.stackexchange.com/questions/2992119/derivative-of-log-marginal-likelihood

Euh t'es sûr ? si tu prends X = a*I_n alors y a un côté qui tend vers 0 et l'autre vers +inf quand a tend vers inf

Bon j'ai lu la page et je pense que t'as juste oublié (ou pas remarqué ) le -1.
Dans ce cas ça vient juste du fait que tr(AB) = tr(BA)

attends c'est pas tr(A^TB) = tr(B^TA) ? et tu peux permuter les vecteurs ? Je pensais que c'était que les matrices carrées

Du coup on aurait tr(y^TXy) = tr(yy^TX) = yy^T : X ? (Si j'omets le ^-1 par simplicité )

[06:01:53] <Sarazinade>
Du coup on aurait tr(y^TXy) = tr(yy^TX) = yy^T : X ? (Si j'omets le ^-1 par simplicité )

Oui

Le 25 septembre 2022 à 05:57:35 :
attends c'est pas tr(A^TB) = tr(B^TA) ? et tu peux permuter les vecteurs ? Je pensais que c'était que les matrices carrées

Non, enfin ça aussi c'est vrai, mais ça vient surtout du fait que Tr(M) = Tr(M^T), c'est une autre propriété de la trace, ce que je te dis n'a rien a voir

Mais du coup si on a z de dimension nxm , X de dimension mxp et y de dimension pxq alors on a clairement pas tr(zXy) = tr(yzX) suffit de regarder les dimensions de yzX. C'est quoi qui justifie que tr(y^TXy) = tr(yy^TX) ?

J'arrive pas à voir ce qui est A et B dans ton tr(AB) = tr(BA)

[00:30:47] <Sarazinade>
Mais du coup si on a z de dimension nxm , X de dimension mxp et y de dimension pxq alors on a clairement pas tr(zXy) = tr(yzX) suffit de regarder les dimensions de yzX. C'est quoi qui justifie que tr(y^TXy) = tr(yy^TX) ?

J'arrive pas à voir ce qui est A et B dans ton tr(AB) = tr(BA)

Dans ton exemple il faudrait que y soit de dim (p,n), car pour avoir tr(AB)=tr(BA) il faut que A soit de dim (n,p) et B soit de dim (p,n)

A=y^TX est de dim (1,n), B=y est de dim (n,1)

Aaaah oko kok ok ok o ko k ok ok ko ok ko ok ok ok ko ben oui la trace c'est sur une matrice carrée donc mon exemple marche pas, merci beaucoup l'ami je vois plus clair, j'avais fait un blocage