24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

琚裕波

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 175.1
散金: 10
帖子: 51
在线: 20.7小时
虫号: 2094139
注册: 2012-10-29
性别: GG
专业: 能源化工

[求助] 求助！这个程序错在哪里？1stopt运行不了

Title "Kinetics";
Parameters k1_a,k1_b,k2_a,k2_b,KCO_a,KCO_b,KCO2_a,KCO2_b,KH2_a,KH2_b;  //需要拟合的参数
Variable TT,PP,Nin,yH2in,yCOin,yCO2in,yN2in,yN2out;                   //输入数据
Variable y1[Output],y2[Output];                                     //对比数据
StartProgram [VB];
Dim iter As Integer
dim TBA As Double
dim RCO As Double
dim RCO2 As Double
Dim KT1 As Double
dim KF1 As Double
dim KT2 As Double
dim KF2 As Double
dim KH2 As Double
dim KCO As Double
dim KCO2 As Double
Dim T As Double
dim P_Ca As Double
Dim fH2 As Double
dim fCO2 As Double
dim fCO As Double
dim fM As Double
dim fH2O As Double
Dim Beta1 As Double
dim Beta2 As Double
dim Nt as double
dim yCO as double
dim yCO2 as double
for iter=0 to Datalength-1
TBA=503.2             //实验平均温度
T=TT(iter)+273.15
P_Ca=PP(iter)/0.101325
//反应平衡常数计算
KF1=Exp(13.1652+9203.26/T-5.92839*Ln(T)-0.00352404*T+0.0000102264*T*T-0.00000000769446*T^3+2.38583E-12*T^4)*(0.101325^(-2))
KF2=Exp(1.6654+4553.34/T-2.72613*Ln(T)-0.01422914*T+0.000017206*T*T-0.00000001106294*T^3+3.19698E-12*T^4)*(0.101325^(-2))
//逸度系数计算
fH2=PP(iter)*yH2in(iter)*Exp((0.110785+35.3324/T-5005.47/(T*T*T)-19.6109*yH2in(iter)/T-20.9799*yH2in(iter)*yH2in(iter)/T)*P_Ca/T)
fCO2=PP(iter)*yCO2in(iter)*Exp((-0.343605+428.452/T-69217700/(T*T*T)-327.402*yCO2in(iter)/T-374.954*yCO2in(iter)*yCO2in(iter)/T)*P_Ca/T)
fCO=PP(iter)*yCOin(iter)*Exp((-0.093261+189.156/T-399940/(T*T*T)-181.527*yCOin(iter)/T+140.001*yCOin(iter)*yCOin(iter)/T)*P_Ca/T)
Beta1=fM/(KF1*fCO*fH2^2)
Beta2=fM*fH2O/(KF2*fCO2*fH2^3)
//吸附平衡常数计算
KCO=Exp(KCO_a+KCO_b*(1/T-1/TBA))
KCO2=Exp(KCO2_a+KCO2_b*(1/T-1/TBA))
KH2=Exp(KH2_a+KH2_b*(1/T-1/TBA))
//反应速率常数计算
KT1=k1_a*Exp(k1_b/8.314/T)
KT2=k2_a*Exp(k2_b/8.314/T)
//反应速率计算
RCO=KT1*fCO*(fH2^2)*(1-Beta1)/((1+KCO*fCO+KCO2*fCO2+KH2*fH2)^3)
RCO2=KT2*fCO2*(fH2^3)*(1-Beta2)/((1+KCO*fCO+KCO2*fCO2+KH2*fH2)^4)
//通过反应速率计算出口CO和CO2组成，其中1.9108为催化剂质量
Nt=Nin(iter)*yN2in(iter)/yN2out(iter)
yCO=(Nin(iter)*yCOin(iter)-RCO*1.9108)/Nt
yCO2=(Nin(iter)*yCO2in(iter)-RCO2*1.9108)/Nt
//计算值与输入值进行对比
y1(iter)=yCO
y2(iter)=yCO2
Next
EndProgram;
Data;

[ Last edited by nono2009 on 2013-9-10 at 20:21 ]

回复此楼

» 猜你喜欢

想换工作。大多数高校都是评职称时认可5年内在原单位取得的成果吗？已经有9人回复
求个博导看看已经有16人回复
上海工程技术大学张培磊教授团队招收博士生已经有4人回复
上海工程技术大学【激光智能制造】课题组招收硕士已经有5人回复
求助院士们，这个如何合成呀已经有4人回复
临港实验室与上科大联培博士招生1名已经有9人回复
需要合成515-64-0，50g，能接单的留言已经有4人回复
自荐读博已经有4人回复
写了一篇“相变储能技术在冷库中应用”的论文，论文内容以实验为主，投什么期刊合适？已经有6人回复
带资进组求博导收留已经有10人回复

» 本主题相关价值贴推荐，对您同样有帮助:

求助，用1stopt解方程组，两次得到的解不唯一？什么原因，大神帮忙已经有6人回复
sos!!1stopt解非线性方程组，一直出现常数定义错误的问题！求大神帮助！已经有12人回复
求大侠帮助跑一下1stopt 已经有12人回复
求助用1stOpt的sharedmodel 运行一组数据求参数，谢谢了啊，急需！已经有4人回复
哪位朋友有注册版1stOpt软件？急需已经有5人回复
1stOPT1.5版每次仿真结果不一样，求助！急！已经有3人回复
用1stopt拟合非线性方程结果与其他软件拟合结果差异大已经有5人回复
求高版本1stopt，解决一个问题，微分方程组参数拟合已经有3人回复
运算过程有复数出现，求高版本1stOpt 已经有6人回复
1stopt运行程序，总是出现字符函数错误的提示，请高手帮助已经有5人回复
大家在使用1stOpt进行数据拟合的时候有没有遇到过这种情况？已经有8人回复
1stopt4阶非线性常微分方程，帮小弟看一下，跪求。已经有7人回复
1stOpt高版本求助！求高版本的帮忙走一下谢了！已经有16人回复
1stopt或matlab如何求解以下的非线性方程并拟合出相应曲线？已经有13人回复
学1stopt要多久？已经有12人回复
哪位大哥能帮我用1stopt3.0以上版本跑一个拟合程序已经有11人回复
1stOpt的函数作图不靠谱吗？已经有4人回复
【求助】此程序没有错误，怎么运行不出来？已经有6人回复
用1stopt进行数据拟合时，无法运行sharedmodel ,求助高手已经有7人回复
【求助】安装ansys时，运行Keygen程序,但不会自动生成一个文件License.dat,为什么? 已经有12人回复

一切都会向着最好的方向发展！

1楼 2013-09-10 16:11:52

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

dingd

铁杆木虫 (职业作家)

应助: 1641 (讲师)
金币: 15037.3
散金: 101
红花: 234
帖子: 3410
在线: 1223.5小时
虫号: 291104
注册: 2006-10-28

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
jjdg: 金币+1, 辛苦了 2013-09-13 00:34:04
琚裕波: 回帖置顶 2013-09-13 14:28:13
琚裕波: 金币+20, ★★★很有帮助, 继续回答啊 2013-09-13 19:17:15

1.5破解版不行，5.0正版运行没问题，不用任何修改。

不知道你的参数有无范围限制，有无过拟合现象。

均方差(RMSE):0.00320856901844522
残差平方和(SSE):0.00041179660584506
相关系数(R): 0.999147398619087
相关系数之平方(R^2): 0.998295524167289
确定系数(DC): 0.998293391304997
F统计(F-Statistic): 16.4900938090552

参数最佳估算
-------------------- -------------
k1_a 16199.5274910722
k1_b -54640.7646512102
k2_a 194674.147413486
k2_b -70718.1386986521
kco_a -1.81845739744211
kco_b -23518.1117376316
kco2_a -1.95039686557898
kco2_b 15779.9670841516
kh2_a -0.851175550051731
kh2_b -780.551417096644