Question

请问这个matlab程序错在哪里？

function [x,y]=flow(f,a,b,e)
x1=a;
x2=b;
Q=fsxhesse(f,x1,x2);
x0=[x1 x2]';
fx1=diff(f,'x1');
fx2=diff(f,'x2');
g=[fx1 fx2]';
g1=subs(g);
d=-g1;
while (abs(norm(g1))>e)
t=(-d)'*d/((-d)'*Q*d);
x0=x0-t*g1;
v=x0;
a=[1 0]*x0;
b=[0 1]*x0;
x1=a;
x2=b;
g1=subs(g);
d=-g1;
end
x=v;
y=feval(f,v);

调用程序后显示
??? Error using ==> feval
Argument must contain a string or function_handle.

Error in ==> flow at 23
y=feval(f,v);
请问要怎么改？

Answer 1

怎么调用的？feval要求第一个参数是函数名称或者函数句柄，而不能用符号变量表达式。

追问

我想在最后输出的时候显示方程的最优解，就调用y=feval(f,v)，但matlab显示错误,请问要怎么改？

追答

feval要求第一个参数是函数名称或者函数句柄，而不能用符号变量表达式。

所以不能用feval，用subs就可以了。

追问

我把feval改成subs了,然后调用[x,y]=flow('x1^2+25*x2^2',2,2,10^(-6))，其中2,2是区间，10^(-6)是精度，就是用最速下降法求x1^2+25*x2^2，但是调用得出的结果是

怎么解决?

追答

很显然你还没给x2赋值呢，你只subs了x1吧？

追问

但是我程序里面不是已经有x2=b吗？虽然很不好意思,但是能麻烦你帮我完整得改一下吗?

追答

1. f是一个符号表达式，从程序中可见，f, x1, x2是三个独立的变量，无直接关联。

2. 给一个符号表达式中的符号赋值，应该用subs(f, x1, 1) % 相当于令x1=1，这里的x1不是flow
   函数里定义的变量x1，而是只符号表达式f中的x1。尽管看起来名字相同，这两个x1之间毫无关系。你可以把上面程序中的变量x1和x2改名字，应该完全不会有影响。