matlab求導

在matlab中求導要進行符號運算。

>>syms x;

>>y = x^cos(x);

>>ydot = diff(y, x, 1);%對x求一階導數

ydot = ? x^(cos(x) - 1)*cos(x) - x^cos(x)*log(x)*sin(x)??

>> y2dot = diff(y, x, 2)%求二階導數，求n階導數同理。??

y2dot = ? cos(x)*(x^(cos(x) - 2)*(cos(x) - 1) - x^(cos(x) - 1)*log(x)*sin(x)) - x^(cos(x) - 1)*sin(x) - log(x)*sin(x)*(x^(cos(x) - 1)*cos(x) - x^cos(x)*log(x)*sin(x)) - x^cos(x)*cos(x)*log(x) - (x^cos(x)*sin(x))/x ?

>> subs(ydot, x, pi)%將x = pi帶入ydot

ans =

? -0.101321183642338

?