| 查看: 2659 | 回复: 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] |
回复此楼» 猜你喜欢
欢迎发来filecode的Mz6后的代码验证其规律
已经有74人回复
-
岩土工程学报什么时候才能终审完呐
已经有4人回复
-
时间戳又变了8-15
已经有25人回复
-
2027广东省杰青
已经有6人回复
-
快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。
已经有3人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有9人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有26人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有16人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有12人回复
-
招聘 金属平磨液，抛光液研发工程师
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助 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
25
