如何将1x1sym转换为函数/表达式？

Question

下面你可以找到实现的牛顿法。

function y = NewtonRoot(Fun, DFun, Xest,Err, imax)
%Fun - function
%DFun- derivative of F
%Xest - initial estimate of solution
%Err - maximum error
%y - solution

%EXAMPLE: NewtonRoot(@(x)x^2-4,@(x)2*x,1.3, 0.001, 100)

for i= 1: imax
    Xi = Xest - feval(Fun,Xest)/feval(DFun,Xest);
    if abs((Xi-Xest)/Xest) < Err
        y = Xi;
        break
    end
    Xest= Xi;
end

if i== imax
    fprint('Solution was not obtained in %i iterations.\n', imax)
    y=('No answer');
end

它正在工作：

NewtonRoot(@(x)x^2-4,@(x)2*x,1.3, 0.001, 100)

但实际上我想要计算一个更复杂的函数的导数。因此，我尝试使用diff函数，但它不起作用…你能帮帮我吗？

这是我的试探性的：

syms y(x) x
y=@(x)x^2-4
dy = diff(y,x)

NewtonRoot(y,@(x)diff(y,x),1.3, 0.001, 100)

Answer 1

您可以使用matlabFunction函数，该函数允许将符号表达式转换为函数句柄。因此，对于这个示例：

syms y(x) x
y=@(x)x^2-4;
dy = diff(y,x);

NewtonRoot(y, matlabFunction( diff(y,x)), 1.3, 0.001, 100)

这显然效果很好。