24小时热门版块排行榜

返回列表

紫麒麟de

新虫 (小有名气)

应助: 1 (幼儿园)
金币: 2988.2
散金: 20
红花: 1
帖子: 130
在线: 128.4小时
虫号: 2222128
注册: 2013-01-04
性别: MM
专业: 金属材料的合金相、相变及

[求助] 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输出应为幂的形式

回复此楼

» 猜你喜欢

1楼 2014-06-06 10:48:48

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主紫麒麟de 的主题更新

返回列表