| 查看: 1944 | 回复: 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朵）
世纪乘风» 猜你喜欢
让我中一个面上吧！
已经有10人回复
-
科研孤儿太难了
已经有14人回复
-
言之凿凿，明天官网会通知放榜？
已经有17人回复
-
时间戳又变了
已经有17人回复
-
今天系统多次维护，明天很可能放榜！
已经有11人回复
-
时间戳变了，能看出什么问题？
已经有23人回复
-
只有每年这种时候来逛逛小木虫
已经有22人回复
-
放榜前的不淡定
已经有12人回复
-
重要消息，中午系统在维护
已经有12人回复
-
今天基金会出结果吗？20260819
已经有15人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助！matlab用fsolve函数求解非线性方程组的问题！
已经有19人回复
- 【求助】matlab中fsolve解的精度控制【已解决】
已经有4人回复
baobiao007
木虫 (职业作家)
中国特色
- 应助: 201 (大学生)
- 金币: 6484.2
- 散金: 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
40