| 查看: 3143 | 回复: 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输出应为幂的形式 |
回复此楼» 猜你喜欢
溴的反应液脱色
已经有4人回复
-
参与限项
已经有5人回复
-
推荐一本书
已经有7人回复
-
有没有人能给点建议
已经有5人回复
-
假如你的研究生提出不合理要求
已经有12人回复
-
全日制（定向）博士
已经有5人回复
-
萌生出自己或许不适合搞科研的想法，现在跑or等等看？
已经有4人回复
-
Materials Today Chemistry审稿周期
已经有4人回复
-
对氯苯硼酸纯化
已经有3人回复
-
所感
已经有4人回复