| 查看: 1599 | 回复: 5 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
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 |
回复此楼» 猜你喜欢
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有47人回复
-
25届硕士申请，材料计算方向
已经有2人回复
-
物理学I论文润色/翻译怎么收费?
已经有139人回复
-
课题组招1-2名硕士
已经有2人回复
-
机时越来越贵之后，我把一部分 DFT 流程搬到了 AI 计算平台上
已经有0人回复
-
斯坦福大学与Elsevier世界前2%学者新排名查询
已经有1人回复
-
27申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有19人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
-
apl投稿状态
已经有2人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
-
日内瓦国际发明展申报完整流程（科研人员通用指南）
已经有10人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 1stOpt的函数作图不靠谱吗？
已经有4人回复
- 用origin或1stopt拟合一条隐函数方程曲线
已经有4人回复
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
4楼2012-10-11 15:06:32
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
2楼2012-10-11 08:46:51
3楼2012-10-11 11:35:09
5楼2012-10-12 20:50:33
