在Matlab上使用feval评估匿名函数时遇到问题

问题描述

我正在为bisection方法开发自己的Matlab代码，并定义了一个匿名函数，我试图为其寻找根。它会生成一个简单的图，我知道如果它能正确运行，它很容易找到根。

当我的代码开始以给定的x值评估函数时，它将返回：

Error using feval
Function to evaluate must be represented as a string scalar,character vector,or
function_handle object.

Error in bisection (line 3)
fa = feval(f,a);

我的完整代码是：

function m=bisection(f,a,b,imax,tol)

    fa = feval(f,a);
    fb = feval(f,b);
    i=0;

    if fa*fb>0
        disp('No root here: pick a new interval')
        return
    end
    while abs(b-a) >= tol
        i = i + 1;
        m=(a+b)/2;
        fm = feval(f,m);
        if fa*fb<0
            b = m;
        else
            a = m;
        end
        abs(fm);
    end
    % Show the last approximation considering the tolerance
    w = feval(f,m);
    fprintf('\n x = %f produces f(x) = %f \n %i iterations\n',m,fm,i-1);
    fprintf(' Approximation with tolerance = %f \n',tol);
end

我的功能是：

function [T] = freezing(x)

alpha = 0.138*10^-6;
Ti = 20;
Ts = -15;
t = 60*60*24*60;
tol = 10^-13;

%x = linspace(0,3);

T = @(x) (Ti-Ts)*erf(x./(2*sqrt(alpha*t)))+Ts;
dT = @(x) (Ti-Ts)*(1/sqrt(pi*alpha*t))*exp(-x.^2./(4*alpha*t));
T = T(x);
end

我真的不确定这是什么问题-在命令窗口中，我可以轻松输入x值并获取输出。任何建议都非常感谢！

解决方法

feval需要f才能成为function handle。当您调用bisection时，我怀疑您是将freezing的结果传递给bisection而不是函数句柄传递给freezing

您需要执行bisection(@freezing,...)而不是bisection(freezing(...),...)

@freezing创建一个传递给二等分的函数句柄。后者将冻结结果传递给bisection。