| 查看: 3163 | 回复: 0 | ||
紫麒麟de新虫 (小有名气)
|
[求助]
matlab运行后出现 function tolerance.<stopping criteria details>是什么原因
|
|
function F = myfun(x) global P; y=P; F = [x(8)^2 - y*2*(4.58794*10^(-19))*(0.1057-x(5))/(x(4)*x(2)/2-x(5))-0.1276*10^(-9); x(8)^2 - y*2*(1.49751*10^(-18))*(0.106-x(6))/(x(4)*(1-x(2))/2-x(6))-0.1276*10^(-9); x(8)^2 - y*2*(1.01476*10^(-11))*(0.0582-x(7))/(x(4)/2-x(7))-0.1276*10^(-9); x(5)*x(7)-x(2)*(2.71617*10^(-7))*exp(x(3)/x(8)); x(6)*x(7)-(1-x(2))*(5.71092*10^(-5))*exp(x(3)/x(8)); -x(1)+x(2)*1.24287*10^(-5)+(1-x(2))*1.13096*10^(-5); -x(3)+2*(x(2)*0.7586+(1-x(2))*0.5718)*x(1)/873.15; -x(4)+0.709165/(x(2)*1.24287+(1-x(2))*1.13096)]; clc clear x = [10^(-5) 0.5 10^(-8) 1 10^(-5) 10^(-5) 10^(-5) 10^(-9)] x0=x global P for i=1:1:10 P=10^(i-1) options = optimset('TolFun',1e-9); [x,fval] = fsolve(@myfun,x0,options) A(:,i)=x x0=x end 输出结果也不对，x5~x8输出应为幂的形式 |
回复此楼» 猜你喜欢
博士自荐
已经有6人回复
-
博士推荐
已经有4人回复
-
求环氧树脂研发1名
已经有10人回复
-
280求调剂
已经有5人回复
-
什么是人一生最重要的？
已经有10人回复
-
面上可以超过30页吧？
已经有13人回复
-
为什么中国大学工科教授们水了那么多所谓的顶会顶刊，但还是做不出宇树机器人？
已经有13人回复
-
版面费该交吗
已经有17人回复
-
【博士招生】太原理工大学2026化工博士
已经有8人回复