fsolve求方程组时的错误问题

zhangyong3725

function f=funsss(p)
t=sym('exp(0.1*(x-0.925656)/(0.99999-x))+exp(0.3*(y-0.925656)/(0.99999-y))-z*(x*y-0.95)');
dyf1=diff(t,'x'); dyf2=diff(t,'y'); dyf3=diff(t,'z');
x=p(1); y=p(2); z=p(3); f=[dyf1;dyf2;dyf3];
当输入
options.MaxIter=1e12; options.MaxFunEvals=1e8;
fsolve(@funsss,[0.986;0.963;0.5],options)
后，出错。
??? Undefined function or method 'norm' for input arguments of type 'sym'.
Error in ==> optim\private\trustnleqn at 138
normgradinf = norm(grad,inf);
Error in ==> fsolve at 295
    [x,FVAL,JACOB,EXITFLAG,OUTPUT,msg]=...
但是，当我把M文件中计算出来的偏导数手动带入到f中时，就不会报错。请帮忙指出来。谢谢。

[ 本帖最后由 ChaChing 于 2010-7-22 09:49 编辑 ]

rocwoods

楼主和我论坛生日一样哦，呵呵，你的函数用到符号运算了，f返回的是符号数值，需要转化下。关于MATLAB符号运算，最好用它来推导方程，得到符号函数表达式后用char转化成字符串然后用eval转化成匿名函数。如下：

1. 
   
2. function test
   
3. t=sym('exp(0.1*(x-0.925656)/(0.99999-x))+exp(0.3*(y-0.925656)/(0.99999-y))-z*(x*y-0.95)');
   
4. dyf1 = eval( ['@(x,y,z) ', char( diff(t,'x') ),';']);
   
5. dyf2 = eval( ['@(x,y,z) ', char( diff(t,'y') ),';']);
   
6. dyf3 = eval( ['@(x,y,z) ', char( diff(t,'z') ),';']);
   
7. function f = funsss(p)
   
8. x=p(1);
   
9. y=p(2);
   
10. z=p(3);
    
11. f=[dyf1(x,y,z);dyf2(x,y,z);dyf3(x,y,z)];
    
12. end
    
13. 
    
14. options.MaxIter=1e12;
    
15. options.MaxFunEvals=1e8;
    
16. fsolve(@funsss,[0.986;0.963;0.5],options)
    
17. 
    
18. end
    
19. 
    

复制代码

运行结果：

1. 
   
2. Optimization terminated: first-order optimality is less than options.TolFun.
   
3. ans =
   
4.     0.9818
   
5.     0.9676
   
6.    31.8249
   

复制代码

[ 本帖最后由 rocwoods 于 2010-7-16 16:10 编辑 ]

zhangyong3725

谢谢你了，兄弟。我的Matlab太烂了，如果可以的话，可以向兄弟你学习下。
这个问题困惑我一天了，太感谢你了。
如果不嫌我学识浅薄，可以交流下。

rocwoods

呵呵，客气了，在论坛里面有很多热心的人，欢迎你加入MATLAB学习者这个阵营，只要坚持每天看帖，时间长了自然水平就上去了！