| 查看: 2369 | 回复: 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 ] |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有4人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有5人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有4人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助，用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.7小时
- 虫号: 291104
- 注册: 2006-10-28
2楼2013-09-10 20:55:30
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
