| 查看: 1853 | 回复: 10 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
求助！这个程序错在哪里？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人回复
4楼2013-09-11 10:31:04
dingd
铁杆木虫 (职业作家)
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.5小时
- 虫号: 291104
- 注册: 2006-10-28
2楼2013-09-10 20:55:30
jjdg
版主 (知名作家)
- 程序强帖: 9
- 应助: 223 (大学生)
- 贵宾: 3.926
- 金币: 51003.4
- 红花: 238
- 帖子: 7227
- 在线: 474.2小时
- 虫号: 27000
- 注册: 2003-10-21
- 专业: 药物代谢与药物动力学
- 管辖: 程序语言
3楼2013-09-11 00:51:51
|
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; {实验数据，顺序与前面定义的输入数据+对比数据的顺序一一对应 例如： 输入数据分别为：温度、压力、入口摩尔流量、入口组成（H2,CO,CO2,N2）、出口N2组成 对比数据分别为：出口CO组成、出口CO2组成 以下数据分别为：温度、压力、入口摩尔流量、入口组成（H2,CO,CO2,N2）、出口N2组成，出口CO组成、出口CO2组成 } 201.3 8 0.68 0.7166 0.0998 0.1052 0.0784 0.083 0.0844 0.1034 215.6 8 0.602 0.7397 0.0862 0.0925 0.0816 0.0894 0.0592 0.0888 229.2 8.1 0.772 0.734 0.1028 0.0836 0.0796 0.0932 0.0581 0.0744 244.4 8.1 0.728 0.7229 0.1052 0.0866 0.0853 0.0992 0.064 0.0777 260.7 7.8 0.723 0.7831 0.0723 0.0622 0.0824 0.0899 0.0467 0.0547 199.3 7.1 0.695 0.7397 0.0862 0.0925 0.0816 0.0852 0.0743 0.0901 215 7.1 0.689 0.734 0.1028 0.0836 0.0796 0.0864 0.0772 0.0822 229.8 7.1 0.707 0.7229 0.1052 0.0866 0.0853 0.0986 0.062 0.0817 244.8 7.1 0.722 0.7831 0.0723 0.0622 0.0824 0.09 0.044 0.0564 259.7 7 0.695 0.7166 0.0998 0.1052 0.0784 0.0856 0.0735 0.1044 200.4 6.2 0.738 0.734 0.1028 0.0836 0.0796 0.0831 0.0907 0.0823 215.7 6 0.629 0.7229 0.1052 0.0866 0.0853 0.092 0.0801 0.087 230 6.1 0.697 0.7831 0.0723 0.0622 0.0824 0.09 0.0431 0.0576 245.2 6 0.981 0.7166 0.0998 0.1052 0.0784 0.0865 0.0717 0.1031 260.3 6.1 0.972 0.7397 0.0862 0.0925 0.0816 0.0862 0.0684 0.092 200.1 5 0.704 0.7229 0.1052 0.0866 0.0853 0.0879 0.0962 0.0862 215 5 0.655 0.7831 0.0723 0.0622 0.0824 0.086 0.0578 0.0611 229.5 5 0.765 0.7166 0.0998 0.1052 0.0784 0.0857 0.0711 0.1063 244.6 5 0.739 0.7397 0.0862 0.0925 0.0816 0.0879 0.0624 0.0914 259.8 5 0.648 0.734 0.1028 0.0836 0.0796 0.0841 0.0841 0.0841 |
5楼2013-09-11 10:31:39