Using symbolic matrix operations correctly?
4 vues (au cours des 30 derniers jours)
Afficher commentaires plus anciens
Eric Zhang
le 27 Juil 2016
Réponse apportée : Andrei Bobrov
le 27 Juil 2016
I am new to symbolic computations in MATLAB and am trying out matrix differentiations.
If we take derivative of trace(X*A) w.r.t. X, we should get A^T.
But MATLAB tells me
>> syms X A;
>> diff(trace(X*A), X)
ans =
A
Where am I wrong?
0 commentaires
Réponse acceptée
John D'Errico
le 27 Juil 2016
As you have defined them, X and A are SCALAR symbolic objects, not general matrices. So the trace operation is a no-op, essentially ignored as the trace of a scalar, and the differentiation does not see them as matrices, since A and X are indeed scalars.
However, if A and X are symbolic matrices, it appears that diff does not allow differentiation with respect to a symbolic matrix. So:
>> A = sym('A',[2,2])
A =
[ A1_1, A1_2]
[ A2_1, A2_2]
>> X = sym('X',[2,2])
X =
[ X1_1, X1_2]
[ X2_1, X2_2]
>> diff(trace(A*X),X)
Error using sym/diff (line 69)
The second argument must be a variable or a nonnegative integer specifying the number of differentiations.
So as you can see, this produces an error as you wish to do the operation. Yet, if you compute the derivatives wrt the scalar elements of X, we see the proper elements of the transpose of A are indeed produced.
>> diff(trace(A*X),X(1,2))
ans =
A2_1
>> diff(trace(A*X),X(2,1))
ans =
A1_2
0 commentaires
Plus de réponses (2)
Azzi Abdelmalek
le 27 Juil 2016
A is not defined as a symbolic matrix.
A=sym('a',5)
syms x
diff(trace(x*A),x)
0 commentaires
Voir également
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!