24小时热门版块排行榜

返回列表

ssmile37

铁虫 (小有名气)

应助: 0 (幼儿园)
金币: 14.4
帖子: 87
在线: 36.8小时
虫号: 1416175
注册: 2011-09-25
性别: MM
专业: 热力学

[求助] 1stopt运行程序，总是出现字符函数错误的提示，请高手帮助

不知是哪里出现了错误，点击运行就会出现错误提示，现将程序和数据发上来，请大家帮忙查找下错误
Title 2;
Parameters a,a1,b,b1,c,c1,d,d1,e,e1;
Variable x1,x2,x3,x4,n2,x5,T;
conststr tao145=a+a1*(T-298.15);
conststr tao514=b+b1*(T-298.15);
conststr tao135=c+c1*(T-298.15);
conststr tao513=d+d1*(T-298.15);
conststr tao2414=e+e1*(T-298.15);
conststr tao1424=-e-e1*(T-298.15);
conststr tao4232=-1.0268-9.008*10^-3*(T-298.15);
conststr tao3242=1.0268+9.008*10^-3*(T-298.15);
conststr tao235=-5.7508+5.37*10^-3*(T-298.15);
conststr tao523=11.0427-8.159*10^-3*(T-298.15);
conststr tao245=-3.9126;
conststr tao524=8.2457;
conststr g15=(x3*exp(-0.2*tao135)+2*x4*exp(-0.2*tao145))/(x3+2*x4);
conststr g25=(x3*exp(-0.2*tao235)+2*x4*exp(-0.2*tao245))/(x3+2*x4);
conststr g35=(x1*exp(-0.2*tao135)+x2*exp(-0.2*tao235))/(x1+x2);
conststr g45=(x1*exp(-0.2*tao145)+x2*exp(-0.2*tao245))/(x1+x2);
conststr tao15=-ln(g15)/0.2;
conststr tao25=-ln(g25)/0.2;
conststr tao35=-ln(g35)/0.2;
conststr tao45=-ln(g45)/0.2;
conststr tao5131=tao15-tao135+tao513;
conststr tao5141=tao15-tao145+tao514;
conststr tao5424=tao45-tao245+tao524;
conststr tao5313=tao35-tao135+tao513;
conststr tao5414=tao45-tao145+tao514;
conststr tao5232=tao25-tao235+tao523;
conststr tao5242=tao25-tao245+tao524;
conststr tao5414=tao45-tao145+tao514;
conststr tao5424=tao45-tao245+tao524;
conststr tao5323=tao35-tao235+tao523;
conststr lnr11=x3/(x3+2*x4)*x5*exp(-0.2*tao5131)*tao5131/(x3+x5*exp(-0.2*tao5131))+2*x4/(2*x4+x3)*x5*exp(-0.2*tao5141)/(2*x4+x5*exp(-0.2*tao5141))+x5*g15/(x1*g15+x2*g25+x3*g35+2*x4*g45+x5)*(tao15-(x1*g15*tao15+x2*g25*tao25+x3*g35*tao35+2*x4*g45*tao45)/(x1*g15+x2*g25+x3*g35+2*x4*g45+x5))-x1/(x1+x2)*x3/(x1+x5*exp(-0.2*tao5313))*x5*exp(-0.2*tao5313)*tao5313/(x1+x5*exp(-0.2*tao5313))-x1/(x1+x2)*2*x4/(x1+x2*exp(-0.2*tao2414)+x5*exp(-0.2*tao5414))*(x2*tao2414*exp(-0.2*tao2414)+x5*tao5414*exp(-0.2*tao5414))/(x1+x2*exp(-0.2*tao2414)+x5*exp(-0.2*tao5414))+x2/(x1+x2)*2*x4*exp(-0.2*tao1424)/(x1*exp(-0.2*tao1424)+x2+x5*exp(-0.2*tao5424))*(tao1424-(x1*tao1424*exp(-0.2*tao1424)+x5*tao5424*exp(-0.2*tao5424))/(x1*exp(-0.2*tao1424)+x2+x5*exp(-0.2*tao5424)));
conststr lnr21=x3/(x3+2*x4)*(2*x4*tao4232*exp(-0.2*tao4232)+x5*tao5232*exp(-0.2*tao5232))/(x3+2*x4*exp(-0.2*tao4232)+x5*exp(-0.2*tao5232))+2*x4/(x3+2*x4)*(x3*tao3242*exp(-0.2*tao3242)+x5*tao5242*exp(-0.2*tao5242))/(x3*exp(-0.2*tao3242)+2*x4+x5*exp(-0.2*tao5242))+x5*g25/(x1*g15+x2*g25+x3*g35+2*x4*g45+x5)*(tao25-(x1*tao15*g15+x2*g25*tao25+x3*g35*tao35+2*x4*g45*tao45)/(x1*g15+x2*g25+x3*g35+2*x4*g45+x5))-x2/(x1+x2)*x3/(x2+x5*exp(-0.2*tao5323))*x5*tao5323*exp(-0.2*tao5323)/(x2+x5*exp(-0.2*tao5323))+x1/(x1+x2)*2*x4*exp(-0.2*tao2414)/(x1+x2*exp(-0.2*tao2414)+x5*exp(-0.2*tao5414))*(tao2414-(x2*tao2414*exp(-0.2*tao2414)+x5*tao5414*exp(-0.2*tao5414))/(x1+x2*exp(-0.2*tao2414)+x5*exp(-0.2*tao5414)))-x2/(x1+x2)*2*x4/(x1*exp(-0.2*tao1424)+x2+x5*exp(-0.2*tao5424))*(x1*tao1424*exp(-0.2*tao1424)+x5*tao5424*exp(-0.2*tao5424))/(x1*exp(-0.2*tao1424)+x2+x5*exp(-0.2*tao5424));
conststr A=-61.44534*exp(T/273.15-1)+2.864468*(exp(T/273.15-1))^2+183.5379*ln(T/273.15)-0.6820223*(T-273.15)+0.0007875695*(T^2-(273.15)^2)+58.95788*(273.15/T);
conststr I=0.5*(x1+x2+x3+4*x4);
conststr M=(n2*98+x5*18)/(n2+x5);
conststr lnr12=-sqrt(1000/M)*A*(2/14.9*ln(1+14.9*sqrt(I))+(sqrt(I)-2*I^1.5)/(1+14.9*sqrt(I)));
conststr lnr22=-sqrt(1000/M)*A*(2/14.9*ln(1+14.9*sqrt(I))+(sqrt(I)-2*I^1.5)/(1+14.9*sqrt(I)));
conststr r1=exp(lnr11+lnr12);
conststr r2=exp(lnr21+lnr22);
Function 0.189=r1*x1/(r2*x2);
Data;
0.01073 0.003444 0.00477 0.01417 0.0009472 0.9798 289.10
0.01166 0.002988 0.004278 0.01465 0.0009463 0.9789 292.65
0.01249 0.002595 0.003825 0.01508 0.0009455 0.9781 296.55
0.01459 0.001654 0.002626 0.01624 0.0009435 0.9760 307.05
0.01621 0.0009806 0.001647 0.01719 0.0009419 0.9744 317.75
0.01712 0.0006286 0.001079 0.01774 0.0009411 0.9735 323.05
0.01804 0.0002739 0.00049 0.01831 0.0009402 0.9726 328.95
0.01924 5.184E-13 9.564E-13 0.01888 0.0009390 0.9714 336.65

回复此楼