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)

显然效果很好。

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

1 个答案: