| 查看: 1601 | 回复: 4 | ||
[求助]
使用matlab 中的fsolve 老提示出错
|
|
我最近想用Matlab 中fsolve求解一组非线性方程。编写程序如下： function f= myfun(t) syms x a b t0 t4; a=5000;b=10000;t0=0;t4=1; t1=t(1);t2=t(2);t3=t(3); f(1)=-int(a+b*x,t1,t2)+(t1-t0)*(a+b*t1); f(2)=-int(a+b*x,t2,t3)+(t2-t1)*(a+b*t2); f(3)=-int(a+b*x,t3,t4)+(t3-t2)*(a+b*t3); end 但在调用fsolve(@myfun,[0.2 0.5 0.7]) 计算时，老报错，如下： Undefined function or method 'isfinite' for input arguments of type 'sym'. Error in ==> trustnleqn at 27 if any(~isfinite(Fvec)) Error in ==> fsolve at 366 [x,FVAL,JACOB,EXITFLAG,OUTPUT,msgData]=... 如果不用M函数文件，以上是有结果的，但不知编成M函数文件时始终提示不对。 望高手帮忙，万分感谢。 |
回复此楼» 本帖已获得的红花（最新10朵）
世纪乘风» 猜你喜欢
308求调剂
已经有4人回复
-
NSFC申报书里申请人简历中代表性论著还需要在申报书最后的附件里面再上传一遍吗
已经有14人回复
-
材料与化工一志愿南昌大学327求调剂推荐
已经有6人回复
-
化学调剂0703
已经有7人回复
-
327求调剂
已经有11人回复
-
调剂
已经有8人回复
-
梁成伟老师课题组欢迎你的加入
已经有7人回复
-
伙伴们，祝我生日快乐吧
已经有24人回复
-
中科院材料273求调剂
已经有3人回复
-
材料工程专硕274一志愿211求调剂
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助！matlab用fsolve函数求解非线性方程组的问题！
已经有19人回复
- 【求助】matlab中fsolve解的精度控制【已解决】
已经有4人回复
baobiao007
木虫 (职业作家)
中国特色
- 应助: 201 (大学生)
- 金币: 6482.7
- 散金: 557
- 红花: 40
- 帖子: 3050
- 在线: 1009.9小时
- 虫号: 505962
- 注册: 2008-02-18
- 专业: 应用地球物理学
【答案】应助回帖
★ ★
感谢参与，应助指数 +1
xiegangmai: 金币+2, 谢谢应助！ 2012-09-18 19:48:41
感谢参与，应助指数 +1
xiegangmai: 金币+2, 谢谢应助！ 2012-09-18 19:48:41
|
改成下面的形式: function f=myfun(t) syms x a b t0 t4 t1=t(1);t2=t(2);t3=t(3); f1=-int(a+b*x,t1,t2)+(t1-t0)*(a+b*t1); f2=-int(a+b*x,t2,t3)+(t2-t1)*(a+b*t2); f3=-int(a+b*x,t3,t4)+(t3-t2)*(a+b*t3); y1=subs(f1,{a,b,t0,t4},{5000 10000 0 1}); y2=subs(f2,{a,b,t0,t4},{5000 10000 0 1}); y3=subs(f3,{a,b,t0,t4},{5000 10000 0 1}); f=[y1,y2,y3]; end |
» 本帖已获得的红花（最新10朵）
2楼2012-09-18 13:53:47
3楼2012-09-18 19:06:18
4楼2012-09-18 19:08:14
5楼2015-12-10 12:01:37
20