MATLAB: How can i determine the analytical expression of the derived funtion “df(x)/dx”

matlab function

How can i determine the analytical expression of the derived funtion "df(x)/dx"?

x=linspace (0:40)

L=3;

sol(1)=2;

psi=1;

f (x) = @(x) cos ( (sol(1)/L) *x)- cosh( (sol(1)/L) *x)+ psi *( sin( (sol(1)/L) *x)- sinh( (sol(1))

Best Answer

Use the Symbolic Math Toolbox:

syms x 
L=sym(3);
sol(1)=sym(2);
psi=sym(1);
f (x) =  cos ( (sol(1)/L) *x)- cosh( (sol(1)/L) *x)+ psi *( sin( (sol(1)/L) *x)- sinh( (sol(1))));
df_dx = diff(f,x)
df_dx_fcn = matlabFunction(df_dx)
df_dx(x) =
(2*cos((2*x)/3))/3 - (2*sin((2*x)/3))/3 - (2*sinh((2*x)/3))/3
df_dx_fcn = @(x) cos(x.*(2.0./3.0)).*(2.0./3.0)-sin(x.*(2.0./3.0)).*(2.0./3.0)-sinh(x.*(2.0./3.0)).*(2.0./3.0)

Related Solutions

MATLAB: Problem with diff(f, diff())

Your first problem is that you did not declare a to be a symbolic function. diff() applied to a constant variable is going to yield 1, not a placeholder derivative.

syms a(x) adot
adot = diff(a);

but then you have

f(x, adot)=3*x+2*adot^2;

adot is now a function. It is not possible to define a function like f with a parameter that is a function. It is valid to define

f(x) = 3 * x + 2*adot(x)^2

Then you have

df=diff(f, adot)

which attempts take the derivative of f with respect to a function. diff() can only take derivatives with respect to variables. There is functionalDerivative() in newer versions of MATLAB, but the best it would be able to do would be to take the derivative with respect to the function a

>> functionalDerivative(f,a)
ans(x) =
-4*diff(a(x), x, x)

Perhaps you would prefer,

   syms AD
   f(x, AD) = subs(3 * x + 2*adot(x)^2, adot, AD)
   subs(diff(f, AD), AD, adot)

The problem with this is that it assumes that adot and x are independent of each other, which is not the case.

MATLAB: How can i display a funtion dependenet to x

This works:

L=3;
sol(1)=2;
psi=1;
f = @(x) cos ( (sol(1)/L) *x)- cosh( (sol(1)/L) *x)+ psi *( sin( (sol(1)/L) *x)- sinh( (sol(1)/L) *x));
x = linspace(-10, 10);
figure(1)
plot(x, f(x))
grid
xlabel('x')
ylabel('f(x)')

See the documentation for the linspace function to understand how to change the limits and number of elements in the ‘x’ vector.

Best Answer

Related Solutions

MATLAB: Problem with diff(f, diff())

MATLAB: How can i display a funtion dependenet to x

Related Question