Anonymous Function Differentiation Problem?

127 visualizaciones (últimos 30 días)
Steve
Steve el 29 de Sept. de 2011
Respondida: Konstantin Ninidze el 24 de Abr. de 2021
Hello Experts,
I have the following function:
>>f = @(x) x^2-1.
I want to find a function that will do differentiation: g = f' = 2*x so that g will be anonymous too and it will be possible to evaluate for example: g(5) = 10.
Please help me,
Steve

Respuestas (5)

Richard Edmonds
Richard Edmonds el 5 de Nov. de 2015
Yes you can do it. Here is the code to acheive exactly what you want. (too late for you perhaps but this was the first link in my search so I came beack to drop in the best answer I found).
fdash=@(fun,x,epsilon)imag(fun(x(:)+1i*epsilon))/epsilon;
g=@(x)fdash(f,x,10^-20);

Matt Tearle
Matt Tearle el 29 de Sept. de 2011
If you know the independent variable is called 'x':
fprime = str2func(['@(x) ',char(diff(sym(regexprep(char(f),'^@\(x\)',''))))])
fprime(pi)
If not... well, you can do it, but it will be a bit messier.

Walter Roberson
Walter Roberson el 29 de Sept. de 2011
You cannot do that, not without converting the anonymous function in to a character string and doing symbolic differentiation on that, and constructing the appropriate anonymous function from the result.
If you attempt to do this numerically, then because you can only sample at a finite number of numeric points, there will be an infinite number of functions that match the numeric results, leaving you unable to select the correct function to differentiate. You will, for example, be unable to find removable discontinuities numerically (except by luck). And good luck doing a numeric reconstruction of tan(Pi*exp(x)) .

Doug Hull
Doug Hull el 29 de Sept. de 2011
If you know you will be dealing only with polynomials, can you make something like this work?

Konstantin Ninidze
Konstantin Ninidze el 24 de Abr. de 2021
you can try my function: handle_to_sym:
function Fs=handle_to_sym(F)
% input F function_handle
% output Fs symbolic expression
S=char(F);
S=S(5:end);
Fs = evalin(symengine,S);

Categorías

Más información sobre Loops and Conditional Statements en Help Center y File Exchange.

Productos

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by