| 查看: 3276 | 回复: 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输出应为幂的形式 |
回复此楼» 猜你喜欢
科研孤儿太难了
已经有15人回复
-
让我中一个面上吧！
已经有10人回复
-
言之凿凿，明天官网会通知放榜？
已经有17人回复
-
时间戳又变了
已经有17人回复
-
今天系统多次维护，明天很可能放榜！
已经有11人回复
-
时间戳变了，能看出什么问题？
已经有23人回复
-
只有每年这种时候来逛逛小木虫
已经有22人回复
-
放榜前的不淡定
已经有12人回复
-
重要消息，中午系统在维护
已经有12人回复
-
今天基金会出结果吗？20260819
已经有15人回复
40