How can I take the derivative of a symbolic matrix?
1 visualizzazione (ultimi 30 giorni)
Mostra commenti meno recenti
Benjamin Schuessler
il 7 Gen 2018
Modificato: Benjamin Schuessler
il 7 Gen 2018
Hello, I'm having a little trouble taking the derivative of a symbolic matrix.
Let's say I have a symbolic matrix 'f':
f =
[ 1, 1, 1, 1, 1]
[ cos(x)^2 - 1, cos(x)^2 - 1, cos(x)^2 - 1, 1, 1]
[ (cos(x)^2 - 1)^2, (cos(x)^2 - 1)^2, (cos(x)^2 - 1)^2, cos(x)^2 - 1, 1]
How would I go about taking the derivative of each element, where the order of the derivative is the value of the exponent in the function (i.e f(1,1) exponent = 0 (no derivative) f(1,2) exponent = 1 (first deriv.) and f(1,3) exponent = 2 (second deriv.)?
Thanks in advance for any advice!
EDIT* I started out with:
a =
0 0 0 1 2
1 1 1 2 3
2 2 2 3 4
then using:
syms x
f = ((cos(x)^2)-1).^a
f =
[ 1, 1, 1, cos(x)^2 - 1, (cos(x)^2 - 1)^2]
[ cos(x)^2 - 1, cos(x)^2 - 1, cos(x)^2 - 1, (cos(x)^2 - 1)^2, (cos(x)^2 - 1)^3]
[ (cos(x)^2 - 1)^2, (cos(x)^2 - 1)^2, (cos(x)^2 - 1)^2, (cos(x)^2 - 1)^3, (cos(x)^2 - 1)^4]
Trying this equation df = diff(f,(a(i,j))) results in:
The second argument must be a variable or a nonnegative integer specifying the number of differentiations.
0 Commenti
Risposta accettata
Walter Roberson
il 7 Gen 2018
Modificato: Walter Roberson
il 7 Gen 2018
temp = arrayfun(@(colidx) diff(f(:,colidx), colidx-1), 1:size(f,2), 'uniform', 0);
result = [temp{:}];
4 Commenti
Walter Roberson
il 7 Gen 2018
f(1,2) = 1, and f(1,3) = 1, but you defined the two as having different derivatives, so we have to use the column number minus one as the order of the derivative.
You have not consistently defined the derivative order to be used. The two expressions cos(x)^2 - 1, (cos(x)^2 - 1)^2 both have a term to power 2, so one might say the order should be 2; on the other hand, if you expand out (cos(x)^2 - 1)^2 then you will have a cos(x)^4 so perhaps it should be order 4.
Più risposte (0)
Vedere anche
Categorie
Scopri di più su Symbolic Math Toolbox in Help Center e File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!