Question

我想了解以下代码中的错误？代码如下所示

function x0=secand2d(f,g,x,dx,tol)
% find solution of f(x,y)=0 and g(x,y)=0
for i=1:20
    x0=x;
    f0=[feval(f,x0),feval(g,x0)];
   x0=x+[dx(1),0];
   fx=[feval(f,x0),feval(g,x0)];
   fx=(fx-f0)/dx(1);
   x0=x+[0 dx(2)];
   fy=[feval(f,x0),feval(g,x0)];
   fy=(fy-f0)/dx(2);
   A=[fx;fy]+eps;
   x=x+dx;
   dx=(-A\f0'-dx')';
   if(norm(dx)<tol) ;return; end;
end;

 disp(norm(dx));
 pause;

end

表示二维割线方法，其中函数f和g由以下形式定义

f=@(x,y) (x-sin(x+y)); and   g=@(x,y) y-cos(x-y);

我还检查过，这些功能在它的参数上运行正常

f(3,3)

ans =

    3.2794

和

g(1,1)

ans =

     0

也

x=[0.9 0.9];
dx=[0.1 0.1];
tol=1e-5;

但是下面的代码产生了这种错误

secand2d(f,g,x,dx,tol)
Error using @(x,y)(x-sin(x+y))
Not enough input arguments.

Error in secand2d (line 5)
    f0=[feval(f,x0),feval(g,x0)];

请帮我澄清这个错误的原因

Answer 1

您收到该错误是因为此处仅将x0传递给f：

f0=[feval(f,x0),feval(g,x0)];

f和g期待两个输入。您可以像这样修复错误：

f0=[feval(f,x0,y0),feval(g,x0,y0)];

如果您使用feval没有特殊原因，可以这样做：

f0 = [f(x0, y0) g(x0, y0)];

二维secand方法中的错误

1 个答案: