| 查看: 2698 | 回复: 1 | |||
[求助]
fsolve计算结果问题
|
|
计算结果感觉误差好大，是由于初值取值造成的影响吗，还是要设置option来控制精度，第一次使用matlab一窍不通，请大家帮忙看下。 First-Order Norm of Iteration Func-count Residual optimality Lambda step 0 8 39.2026 1.13e+003 0.01 1 16 0.00103045 0.705 0.001 0.294456 2 24 3.55079e-005 0.0727 0.0001 0.237055 3 32 7.28141e-006 0.19 1e-005 0.144512 4 40 2.91808e-006 0.0123 1e-006 0.0248491 5 48 2.90763e-006 3.58e-007 1e-007 0.000338087 6 56 2.90763e-006 1.02e-009 1e-008 1.90914e-006 Equation solved, fsolve stalled. fsolve stopped because the relative size of the current step is less than the default value of the step size tolerance and the vector of function values is near zero as measured by the selected value of the function tolerance. y = -2.2413 2.2083 -0.1310 0.1318 1.0596 0.0211 -0.0272 fval = 0.0000 0.0008 -0.0008 0.0003 0.0000 0.0000 -0.0011 0.0006 下面是我要解的方程组 function F=myfun(x) b=1;a=1; F=[-0.0575*x(3)/a-x(7)*a*0.5-x(6); x(3)/a^2+x(4)+17.39*x(7)*log(a)+0.8695652174*b*log(a); -x(1)/(x(5)^2)-x(2)+x(3)/(x(5)^2)+x(4)-92.49880554*x(7)*log(x(5))+0.8695652174*b*log(x(5)); 0.0065*x(1)/x(5)-0.0026*x(2)*x(5)-0.0575*x(3)/x(5)+x(7)*x(5)*((0.7-0.4*log(x(5)))/1.4-0.9347826087); -2*x(1)/(x(5)^2)+109.8901099*x(7)-b; 2*x(1)/(1.5^2)-109.8901099*x(7)-b; 0.444444444444444*x(1)/a^2+x(2)+109.8901099*x(7)*log(1.5*a); -x(1)*log(1.5*a/x(5))+0.5*(1.5^2*a^2-x(5)^2)*x(2)-x(3)*log(x(5)/a)+0.5*(x(5)^2-a^2)*x(4)+17.39130435*x(7)*(0.5*x(5)^2*(log(x(5))+0.5)-0.5*a^2*(log(a)+0.5))+109.8901099*x(7)*(0.5*1.5*1.5*a^2*(0.5+log(1.5*a))-0.5*x(5)^2*(log(x(5))+0.5))+0.0209742*1.5^2*a^2*b+0.8695652174*b*(0.5*x(5)^2*(log(x(5))+0.5)-0.5*a^2*(log(a)+0.5))]; 选取的初值是x0=[-2,2,-0.5,0.5,1.05,0.01,-0.01] |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有8人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有8人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有11人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有11人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有8人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有10人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有12人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有8人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有10人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助 MATLAB-fsolve
已经有3人回复
- 使用fsolve求解非线性方程问题
已经有9人回复
- MATLAB求解代数方程组问题，请帮忙
已经有9人回复
- 各位朋友，如何用matlab编写这个方程，并求解？
已经有15人回复
- 求助 MATLAB解方程组-fslove
已经有7人回复
- vasp计算态密度结果出问题
已经有11人回复
- NBO计算结果问题
已经有8人回复
- 求助！matlab用fsolve函数求解非线性方程组的问题！
已经有19人回复
- Ansys载荷步计算结果的保存问题
已经有8人回复
- 各位大侠，请问下我在linux计算，结果老是出现下面的问题，请问怎么办呀。
已经有23人回复
- matlab 问题求助，是否计算结果错误，如何改正？
已经有3人回复
- 非线性方程组matlab拟合
已经有9人回复
- mathematica请教计算结果问题
已经有4人回复
- 为什么计算区域的大小对计算结果的精度影响这么大？颗粒扰流问题
已经有6人回复
- 气相色谱问题，出峰但是不计算结果
已经有5人回复
- Abinit 新旧版计算结果差别多大可以认为正常，没有安装错误之类的问题。
已经有4人回复
- 【求助】请高手帮我解答个比较疑惑的问题（UNIFAC估算）
已经有4人回复
- 【求助】vasp5 光学计算结果问题
已经有4人回复
- 【求助】请教一个关于解方程的问题
已经有8人回复
- 【求助】fluent的传热问题，计算结果有没有问题，请大家帮忙看看。【已解决】
已经有10人回复
nono2009
超级版主 (文学泰斗)
No gains, no pains.
-
专家经验: +21105 - 仿真EPI: 1
- 应助: 28684 (院士)
- 贵宾: 513.911
- 金币: 2555230
- 散金: 27828
- 红花: 2148
- 沙发: 66666
- 帖子: 1602255
- 在线: 65200.9小时
- 虫号: 827383
- 注册: 2009-08-13
- 性别: GG
- 专业: 工程热物理与能源利用
- 管辖: 科研家筹备委员会
2楼2013-01-07 14:23:41
