0

我想了解以下代码中的错误是什么？代码如下

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)];

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

4

1 回答 1

2

您收到该错误是因为您只传递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)];
于 2013-05-09T13:18:59.293 回答