- 积分
- 0
- 注册时间
- 2010-6-13
- 仿真币
-
- 最后登录
- 1970-1-1
|
本帖最后由 yutingting1986 于 2011-4-27 20:57 编辑
function F = myfun(x)
F =[x(1)^2/x(2)-43.5*T^1.5*(1-exp(-3340/T))*exp(-113165.6/T);
x(3)^2/x(4)-280.217*T^1.5*(1-exp(-2230/T))*exp(-59940.1/T);
x(5)^2/(x(4)*x(2))-5.4119*(1-exp(-2230/T))*(1-exp(-3340/T))*
(1-exp(-2690/T))^2*exp(-46788.55/T);
%空行中间是一个式子
x(6)*x(7)/x(5)-3.024e-7*T^2.5*exp(-107379.1/T);
2*x(2)+x(1)+x(5)+x(6)-7.054e-2;
2*x(4)+x(3)+x(5)+x(6)-1.875e-2;
x(6)-x(7)];
上面是函数m文件。
T=7000;
format long e;
x0 =[0.01;0.01;0.01;0.01;0.01;0.01;0.01];
options=optimset('Display','iter');
[x,fval] = fsolve(@myfun,x0,options);
这是我的求解程序,给定的X0值不同结果就不同,而且结果不对,0<x<0.04。不知道该怎么解决,急求求解方法!!!请各位高手帮忙! |
|