版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

返回列表

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

shuizl

木虫 (正式写手)

应助: 1 (幼儿园)
金币: 3840.6
散金: 499
红花: 5
帖子: 701
在线: 402.4小时
虫号: 792716
注册: 2009-06-12
性别: GG
专业: 化学反应工程

[求助] 1stopt 的运行结果之输出结果全部为0 求助

1stopt 的运行结果之输出结果全部为0

源代码：
Title "AC2O";
Parameters k_1[0,1e18],k_2[0,1e18],k_3[0,1e18],k_4[0,1e18],k_5[0,1e18],k_6[0,1e18],E1[2e4,2e5],E2[2e4,2e5],E3[2e4,2e5],E4[2e4,2e5],E5[2e4,2e5],E6[2e4,2e5];
Variable TT,PDMEin,CCH3Iin,CRhin,PCOin,CLiIin;
Variable C1out[Output],C2out[Output];
StartProgram [VB];
Dim iter As Integer
Dim RK As Integer
Dim j As Integer
Dim RK_i As Integer
Dim C_1 As Double
Dim C_2 As Double
Dim k_1 As Double
Dim k_2 As Double
Dim k_3 As Double
Dim k_4 As Double
Dim k_5 As Double
Dim k_6 As Double
Dim E1 As Double
Dim E2 As Double
Dim E3 As Double
Dim E4 As Double
Dim E5 As Double
Dim E6 As Double
Dim ti As Double
Dim T As Double
Dim TT As Double
Dim PDMEin As Double
Dim CCH3Iin As Double
Dim CRhin As Double
Dim PCOin As Double
Dim CLiIin As Double
Dim C1out As Double
Dim C2out As Double
Dim PDMEi As Double
Dim CCH3I As Double
Dim CRh As Double
Dim PCO As Double
Dim CLiI As Double
Dim ee1 As Double
Dim ee2 As Double
Dim ee3 As Double
Dim ee4 As Double
Dim ee5 As Double
Dim x_x1 As Double
Dim x_x2 As Double
Dim y_y1 As Double
Dim y_y2 As Double
Dim C1 As Double
Dim C2 As Double
Dim K1 As Double
Dim K2 As Double
Dim K3 As Double
Dim K4 As Double
Dim K5 As Double
Dim K6 As Double
Dim d_c1 As Double
Dim d_c2 As Double
Dim c As Double
For RK=0 to Datalength-1
  T=TT(iter)
  PDMEi=PDMEin(iter)
  CCH3I=CCH3Iin(iter)
  CRh=CRhin(iter)
  PCO=PCOin(iter)
  CLiI=CLiIin(iter)
  C_1=0;
  C_2=0;
  ti=6/20;
  For RK_i=1 to 20
      ee1=0.5
      ee2=0.5
      ee3=1
      ee4=1
      ee5=0.5
      x_x1=C_1
      x_x2=C_2
      y_y1=C_1
      y_y2=C_2
      For j=1 To 4
         C1=x_x1
         C2=x_x2
         K1=k_1*exp(-E1/T/8.314)
         K2=k_2*exp(-E2/T/8.314)
         K3=k_3*exp(-E3/T/8.314)
         K4=k_4*exp(-E4/T/8.314)
         K5=k_5*exp(-E5/T/8.314)
         K6=k_6*exp(-E6/T/8.314)
         d_c1=K1*(PDMEi-8.314*T*(C1+C2))^2*CCH3I^2*CRh*PCO/PDMEi^3-K2*(PDMEi-8.314*T*(C1+C2))*CCH3I*CRh*C1/PDMEi^2
         d_c2=K3*(PDMEi-8.314*T*(C1+C2))*CRh*PCO*C2/PDMEi/C1-K4*CRh*C2+K5*CCH3I*C1^2*PDMEi/(PDMEi-8.314*T*(C1+C2))-K6*CCH3I*C1*C2*CRh*PDMEi/CLiI/(PDMEi-8.314*T*(C1+C2))
         c=d_c1*ti
         if j=1 then
            x_x1=C_1+ee1*c
            y_y1=y_y1+ee2*c/3
         elseif j=2 then
            x_x1=C_1+ee2*c
            y_y1=y_y1+ee3*c/3
         elseif j=3 then
            x_x1=C_1+ee3*c
            y_y1=y_y1+ee4*c/3
         elseif j=4 then
            x_x1=C_1+ee4*c
            y_y1=y_y1+ee5*c/3
         end if
         c=d_c2*ti
         if j=1 then
            x_x2=C_2+ee1*c
            y_y2=y_y2+ee2*c/3
         elseif j=2 then
            x_x2=C_2+ee2*c
            y_y2=y_y2+ee3*c/3
         elseif j=3 then
            x_x2=C_2+ee3*c
            y_y2=y_y2+ee4*c/3
         elseif j=4 then
            x_x2=C_2+ee4*c
            y_y2=y_y2+ee5*c/3
         end if
      Next
      C_1=y_y1
      C_2=y_y2
  Next
  C1out(iter)=y_y1
  C2out(iter)=y_y2
Next
ENDPROGRAM;
Data;
453.00 150649.68 126.76 0.41 3049350.32 37.31 19.32 8.80
453.00 150649.68 126.76 0.83 3049350.32 37.31 13.00 19.28
453.00 150649.68 126.76 1.24 3049350.32 37.31 9.20 26.32
453.00 150649.68 126.76 1.66 3049350.32 37.31 8.72 29.52
453.00 150649.68 126.76 2.07 3049350.32 37.31 8.06 31.82
453.00 150649.68 126.76 2.48 3049350.32 37.31 6.95 30.41
453.00 150649.68 126.76 2.90 3049350.32 37.31 5.59 26.17
453.00 150649.68 42.25 2.07 3049350.32 37.31 6.85 13.43
453.00 150649.68 70.42 2.07 3049350.32 37.31 8.57 20.71
453.00 150649.68 84.51 2.07 3049350.32 37.31 10.04 26.91
453.00 150649.68 126.76 2.07 3049350.32 37.31 8.06 31.82
453.00 150649.68 169.01 2.07 3049350.32 37.31 6.47 29.89
453.00 150649.68 211.27 2.07 3049350.32 37.31 5.34 29.84
453.00 150649.68 253.52 2.07 3049350.32 37.31 4.40 29.96
453.00 150649.68 126.76 2.07 3049350.32 7.46 12.08 25.16
453.00 150649.68 126.76 2.07 3049350.32 14.93 10.42 27.91
453.00 150649.68 126.76 2.07 3049350.32 22.39 8.96 29.84
453.00 150649.68 126.76 2.07 3049350.32 29.85 8.45 30.95
453.00 150649.68 126.76 2.07 3049350.32 37.31 8.06 31.82
453.00 150649.68 126.76 2.07 3049350.32 44.78 7.75 32.12
453.00 150649.68 126.76 2.07 3049350.32 52.24 7.17 32.55
453.00 75324.84 126.76 2.07 3124675.16 37.31 16.12 3.88
453.00 112987.26 126.76 2.07 3087012.74 37.31 12.25 17.75
453.00 150649.68 126.76 2.07 3049350.32 37.31 8.06 31.82
453.00 188312.10 126.76 2.07 3011687.90 37.31 8.73 41.09
453.00 225974.52 126.76 2.07 2974025.48 37.31 9.48 50.40
453.00 263636.94 126.76 2.07 2936363.06 37.31 9.44 52.55
453.00 301299.36 126.76 2.07 2898700.64 37.31 7.60 43.20
453.00 150649.68 126.76 2.07 2349350.32 37.31 12.85 13.63
453.00 150649.68 126.76 2.07 2599350.32 37.31 12.80 19.71
453.00 150649.68 126.76 2.07 2849350.32 37.31 11.47 26.70
453.00 150649.68 126.76 2.07 3099350.32 37.31 7.93 31.95
453.00 150649.68 126.76 2.07 3349350.32 37.31 4.74 35.13
453.00 150649.68 126.76 2.07 3599350.32 37.31 4.02 35.82
453.00 150649.68 126.76 2.07 3849350.32 37.31 2.99 36.91
433.00 143998.48 126.76 2.07 3056001.52 37.31 13.52 18.40
438.00 145661.28 126.76 2.07 3104338.72 37.31 12.55 20.47
443.00 147324.08 126.76 2.07 3052675.92 37.31 11.55 22.17
448.00 148986.88 126.76 2.07 3101013.12 37.31 9.59 28.23
453.00 150649.68 126.76 2.07 3049350.32 37.31 8.06 31.82
458.00 152312.48 126.76 2.07 3097687.52 37.31 6.04 31.06
463.00 153975.28 126.76 2.07 3046024.72 37.31 4.40 27.24

运行结果：

"AC2O"

迭代数: 16
计算用时(时:分:秒:微秒): 00:00:00:322
优化算法: 准牛顿法(BFGS) + 通用全局优化法
计算结束原因: 达到收敛判定标准
均方差(RMSE): 22.1009039084701
残差平方和(SSE): 41029.7961
相关系数(R): 0
相关系数之平方(R^2): 0
决定系数(DC): 8.8666753259632E-16
F统计(F-Statistic): 0

参数       最佳估算
-------------------- -------------
k_1         1.5926178987138E17
k_2         5.34790054196492E17
k_3         5.79585216706619E17
k_4         5.18475014017895E17
k_5         5.73410185752437E17
k_6         1.99592327931896E17
E1          76320.0160069391
E2          34522.3072776571
E3          0.659441070165485
E4          1.32150601642206
E5          67320.5888597295
E6          20693.1685144082

====== 输出结果 =====

No 实测C1out 计算C1out 实测C2out 计算C2out
1 19.32 0 8.8 0
2 13 0 19.28 0
3 9.2 0 26.32 0
4 8.72 0 29.52 0
5 8.06 0 31.82 0
6 6.95 0 30.41 0
7 5.59 0 26.17 0
8 6.85 0 13.43 0
9 8.57 0 20.71 0
10 10.04 0 26.91 0
11 8.06 0 31.82 0
12 6.47 0 29.89 0
13 5.34 0 29.84 0
14 4.4 0 29.96 0
15 12.08 0 25.16 0
16 10.42 0 27.91 0
17 8.96 0 29.84 0
18 8.45 0 30.95 0
19 8.06 0 31.82 0
20 7.75 0 32.12 0
21 7.17 0 32.55 0
22 16.12 0 3.88 0
23 12.25 0 17.75 0
24 8.06 0 31.82 0
25 8.73 0 41.09 0
26 9.48 0 50.4 0
27 9.44 0 52.55 0
28 7.6 0 43.2 0
29 12.85 0 13.63 0
30 12.8 0 19.71 0
31 11.47 0 26.7 0
32 7.93 0 31.95 0
33 4.74 0 35.13 0
34 4.02 0 35.82 0
35 2.99 0 36.91 0
36 13.52 0 18.4 0
37 12.55 0 20.47 0
38 11.55 0 22.17 0
39 9.59 0 28.23 0
40 8.06 0 31.82 0
41 6.04 0 31.06 0
42 4.4 0 27.24 0

回复此楼

» 猜你喜欢

求标准粉末衍射卡号 ICDD 01-076-1802 已经有0人回复
新西兰Robinson研究所招收全奖PhD 已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有256人回复
石墨烯转移--二氧化硅衬底石墨烯已经有0人回复
笼目材料中量子自旋液体基态的证据已经有0人回复
数学教学论硕士可以读数学物理博士吗？已经有0人回复
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生已经有4人回复
澳门大学应用物理及材料工程研究院潘晖教授课题组诚招博士后已经有11人回复
求助NH4V4O10晶体的CIF文件已经有0人回复
英国全奖博士招聘-深度学习与量子物理已经有0人回复
间接带隙半导体有效质量求助已经有0人回复

somethingofeverythingandeverythingofsomething

1楼 2015-07-29 16:59:53

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

shuizl

木虫 (正式写手)

应助: 1 (幼儿园)
金币: 3840.6
散金: 499
红花: 5
帖子: 701
在线: 402.4小时
虫号: 792716
注册: 2009-06-12
性别: GG
专业: 化学反应工程

送红花一朵

引用回帖:

2楼: Originally posted by dingd at 2015-07-29 22:18:24
程序中： d_c2=K3*(PDMEi-8.314*T*(C1+C2))*CRh*PCO*C2/PDMEi/C1-.........

第一次循环时C1赋值是为0，上面C1出现在分母位置，被0除？

我看到了，谢谢你我已经修改了，但是还是运行结果为0.能否帮忙再看一下，谢谢！

赞一下

回复此楼

somethingofeverythingandeverythingofsomething

4楼2015-07-30 08:11:34

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

查看全部 5 个回答

dingd

铁杆木虫 (职业作家)

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

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
shuizl: 金币+30, ★★★★★最佳答案, 15 2015-07-30 09:23:46

程序中： d_c2=K3*(PDMEi-8.314*T*(C1+C2))*CRh*PCO*C2/PDMEi/C1-.........

第一次循环时C1赋值是为0，上面C1出现在分母位置，被0除？

赞一下(1人)

回复此楼

» 本帖已获得的红花（最新10朵）

shuizl

2楼2015-07-29 22:18:24

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

shuizl

木虫 (正式写手)

应助: 1 (幼儿园)
金币: 3840.6
散金: 499
红花: 5
帖子: 701
在线: 402.4小时
虫号: 792716
注册: 2009-06-12
性别: GG
专业: 化学反应工程

引用回帖:

Title "AC2O";
Parameters k_1[0,1e180],k_2[0,1e180],k_3[0,1e180],k_4[0,1e180],k_5[0,1e180],k_6[0,1e180],E1[2e4,2e5],E2[2e4,2e5],E3[2e4,2e5],E4[2e4,2e5],E5[2e4,2e5],E6[2e4,2e5];
Variable TT,PDMEin,CCH3Iin,CRhin,PCOin,CLiIin;
Variable C1out[Output],C2out[Output];
StartProgram [VB];
Dim iter As Integer
Dim RK As Integer
Dim j As Integer
Dim RK_i As Integer
Dim C_1 As Double
Dim C_2 As Double
Dim k_1 As Double
Dim k_2 As Double
Dim k_3 As Double
Dim k_4 As Double
Dim k_5 As Double
Dim k_6 As Double
Dim E1 As Double
Dim E2 As Double
Dim E3 As Double
Dim E4 As Double
Dim E5 As Double
Dim E6 As Double
Dim ti As Double
Dim T As Double
Dim TT As Double
Dim PDMEin As Double
Dim CCH3Iin As Double
Dim CRhin As Double
Dim PCOin As Double
Dim CLiIin As Double
Dim C1out As Double
Dim C2out As Double
Dim PDMEi As Double
Dim CCH3I As Double
Dim CRh As Double
Dim PCO As Double
Dim CLiI As Double
Dim ee1 As Double
Dim ee2 As Double
Dim ee3 As Double
Dim ee4 As Double
Dim ee5 As Double
Dim x_x1 As Double
Dim x_x2 As Double
Dim y_y1 As Double
Dim y_y2 As Double
Dim C1 As Double
Dim C2 As Double
Dim K1 As Double
Dim K2 As Double
Dim K3 As Double
Dim K4 As Double
Dim K5 As Double
Dim K6 As Double
Dim d_c1 As Double
Dim d_c2 As Double
Dim c As Double
For RK=0 to Datalength-1
  T=TT(iter)
  PDMEi=PDMEin(iter)
  CCH3I=CCH3Iin(iter)
  CRh=CRhin(iter)
  PCO=PCOin(iter)
  CLiI=CLiIin(iter)
  C_1=0;
  C_2=0;
  ti=6/20;
  For RK_i=1 to 20
      ee1=0.5
      ee2=0.5
      ee3=1
      ee4=1
      ee5=0.5
      x_x1=C_1
      x_x2=C_2
      y_y1=C_1
      y_y2=C_2
      For j=1 To 4
         C1=x_x1
         C2=x_x2
         K1=k_1*exp(-E1/T/8.314)
         K2=k_2*exp(-E2/T/8.314)
         K3=k_3*exp(-E3/T/8.314)
         K4=k_4*exp(-E4/T/8.314)
         K5=k_5*exp(-E5/T/8.314)
         K6=k_6*exp(-E6/T/8.314)
         d_c1=K1*(PDMEi-8.314*T*(C1+C2))^2*CCH3I^2*CRh*PCO/PDMEi^3-K2*(PDMEi-8.314*T*(C1+C2))*CCH3I*CRh*C1/PDMEi^2
         c=d_c1*ti
         if j=1 then
            x_x1=C_1+ee1*c
            y_y1=y_y1+ee2*c/3
         elseif j=2 then
            x_x1=C_1+ee2*c
            y_y1=y_y1+ee3*c/3
         elseif j=3 then
            x_x1=C_1+ee3*c
            y_y1=y_y1+ee4*c/3
         elseif j=4 then
            x_x1=C_1+ee4*c
            y_y1=y_y1+ee5*c/3
         end if
         C1=x_x1
         d_c2=K3*(PDMEi-8.314*T*(C1+C2))*CRh*PCO*C2/PDMEi/C1-K4*CRh*C2+K5*CCH3I*C1^2*PDMEi/(PDMEi-8.314*T*(C1+C2))-K6*CCH3I*C1*C2*CRh*PDMEi/CLiI/(PDMEi-8.314*T*(C1+C2))
         c=d_c2*ti
         if j=1 then
            x_x2=C_2+ee1*c
            y_y2=y_y2+ee2*c/3
         elseif j=2 then
            x_x2=C_2+ee2*c
            y_y2=y_y2+ee3*c/3
         elseif j=3 then
            x_x2=C_2+ee3*c
            y_y2=y_y2+ee4*c/3
         elseif j=4 then
            x_x2=C_2+ee4*c
            y_y2=y_y2+ee5*c/3
         end if
      Next
      C_1=y_y1
      C_2=y_y2
  Next
  C1out(iter)=y_y1
  C2out(iter)=y_y2
Next
ENDPROGRAM;
Data;
453.00 150649.68 126.76 0.41 3049350.32 37.31 19.32 8.80
453.00 150649.68 126.76 0.83 3049350.32 37.31 13.00 19.28
453.00 150649.68 126.76 1.24 3049350.32 37.31 9.20 26.32
453.00 150649.68 126.76 1.66 3049350.32 37.31 8.72 29.52
453.00 150649.68 126.76 2.07 3049350.32 37.31 8.06 31.82
453.00 150649.68 126.76 2.48 3049350.32 37.31 6.95 30.41
453.00 150649.68 126.76 2.90 3049350.32 37.31 5.59 26.17
453.00 150649.68 42.25 2.07 3049350.32 37.31 6.85 13.43
453.00 150649.68 70.42 2.07 3049350.32 37.31 8.57 20.71
453.00 150649.68 84.51 2.07 3049350.32 37.31 10.04 26.91
453.00 150649.68 126.76 2.07 3049350.32 37.31 8.06 31.82
453.00 150649.68 169.01 2.07 3049350.32 37.31 6.47 29.89
453.00 150649.68 211.27 2.07 3049350.32 37.31 5.34 29.84
453.00 150649.68 253.52 2.07 3049350.32 37.31 4.40 29.96
453.00 150649.68 126.76 2.07 3049350.32 7.46 12.08 25.16
453.00 150649.68 126.76 2.07 3049350.32 14.93 10.42 27.91
453.00 150649.68 126.76 2.07 3049350.32 22.39 8.96 29.84
453.00 150649.68 126.76 2.07 3049350.32 29.85 8.45 30.95
453.00 150649.68 126.76 2.07 3049350.32 37.31 8.06 31.82
453.00 150649.68 126.76 2.07 3049350.32 44.78 7.75 32.12
453.00 150649.68 126.76 2.07 3049350.32 52.24 7.17 32.55
453.00 75324.84 126.76 2.07 3124675.16 37.31 16.12 3.88
453.00 112987.26 126.76 2.07 3087012.74 37.31 12.25 17.75
453.00 150649.68 126.76 2.07 3049350.32 37.31 8.06 31.82
453.00 188312.10 126.76 2.07 3011687.90 37.31 8.73 41.09
453.00 225974.52 126.76 2.07 2974025.48 37.31 9.48 50.40
453.00 263636.94 126.76 2.07 2936363.06 37.31 9.44 52.55
453.00 301299.36 126.76 2.07 2898700.64 37.31 7.60 43.20
453.00 150649.68 126.76 2.07 2349350.32 37.31 12.85 13.63
453.00 150649.68 126.76 2.07 2599350.32 37.31 12.80 19.71
453.00 150649.68 126.76 2.07 2849350.32 37.31 11.47 26.70
453.00 150649.68 126.76 2.07 3099350.32 37.31 7.93 31.95
453.00 150649.68 126.76 2.07 3349350.32 37.31 4.74 35.13
453.00 150649.68 126.76 2.07 3599350.32 37.31 4.02 35.82
453.00 150649.68 126.76 2.07 3849350.32 37.31 2.99 36.91
433.00 143998.48 126.76 2.07 3056001.52 37.31 13.52 18.40
438.00 145661.28 126.76 2.07 3104338.72 37.31 12.55 20.47
443.00 147324.08 126.76 2.07 3052675.92 37.31 11.55 22.17
448.00 148986.88 126.76 2.07 3101013.12 37.31 9.59 28.23
453.00 150649.68 126.76 2.07 3049350.32 37.31 8.06 31.82
458.00 152312.48 126.76 2.07 3097687.52 37.31 6.04 31.06
463.00 153975.28 126.76 2.07 3046024.72 37.31 4.40 27.24

赞一下

回复此楼

somethingofeverythingandeverythingofsomething

3楼2015-07-30 08:09:40

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

dingd

铁杆木虫 (职业作家)

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

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
shuizl: 金币+20, ★★★★★最佳答案, 非常感谢，以后有问题还请多帮助。 2015-07-31 09:15:17

改动如下，可以正常运行，效果不好，不知模型逻辑构成方面是否正确，自己慢慢改动运算：

CODE:

Parameters k_1=[0,],k_2=[0,],k_3=[0,],k_4=[0,],k_5=[0,],k_6=[0,],E1=[2e4,2e5],E2=[2e4,2e5],E3=[2e4,2e5],E4=[2e4,2e5],E5=[2e4,2e5],E6=[2e4,2e5];

Variable TT,PDMEin,CCH3Iin,CRhin,PCOin,CLiIin;

Variable C1out[Output],C2out[Output];

StartProgram [Basic];

Dim iter As Integer

Dim RK As Integer

Dim j As Integer

Dim RK_i As Integer

Dim C_1 As Double

Dim C_2 As Double

Dim ti As Double

Dim T As Double

Dim PDMEi As Double

Dim CCH3I As Double

Dim CRh As Double

Dim PCO As Double

Dim CLiI As Double

Dim ee1 As Double

Dim ee2 As Double

Dim ee3 As Double

Dim ee4 As Double

Dim ee5 As Double

Dim x_x1 As Double

Dim x_x2 As Double

Dim y_y1 As Double

Dim y_y2 As Double

Dim C1 As Double

Dim C2 As Double

Dim K1 As Double

Dim K2 As Double

Dim K3 As Double

Dim K4 As Double

Dim K5 As Double

Dim K6 As Double

Dim d_c1 As Double

Dim d_c2 As Double

Dim c As Double

For RK=0 to Datalength-1

  iter=RK

  T=TT(iter)

  PDMEi=PDMEin(iter)

  CCH3I=CCH3Iin(iter)

  CRh=CRhin(iter)

  PCO=PCOin(iter)

  CLiI=CLiIin(iter)

  C_1=0

  C_2=0

  ti=6/20

  For RK_i=1 to 20

        ee1=0.5

        ee2=0.5

        ee3=1

        ee4=1

        ee5=0.5

        x_x1=C_1

        x_x2=C_2

        y_y1=C_1

        y_y2=C_2

        For j=1 To 4

            C1=x_x1

            C2=x_x2

            K1=k_1*exp(-E1/T/8.314)

            K2=k_2*exp(-E2/T/8.314)

            K3=k_3*exp(-E3/T/8.314)

            K4=k_4*exp(-E4/T/8.314)

            K5=k_5*exp(-E5/T/8.314)

            K6=k_6*exp(-E6/T/8.314)

            d_c1=K1*(PDMEi-8.314*T*(C1+C2))^2*CCH3I^2*CRh*PCO/PDMEi^3-K2*(PDMEi-8.314*T*(C1+C2))*CCH3I*CRh*C1/PDMEi^2

            c=d_c1*ti

            if j=1 then

                x_x1=C_1+ee1*c

                y_y1=y_y1+ee2*c/3

            elseif j=2 then

                x_x1=C_1+ee2*c

                y_y1=y_y1+ee3*c/3

            elseif j=3 then

                x_x1=C_1+ee3*c

                y_y1=y_y1+ee4*c/3

            elseif j=4 then

                x_x1=C_1+ee4*c

                y_y1=y_y1+ee5*c/3

            end if

            C1=x_x1

            d_c2=K3*(PDMEi-8.314*T*(C1+C2))*CRh*PCO*C2/PDMEi/C1-K4*CRh*C2+K5*CCH3I*C1^2*PDMEi/(PDMEi-8.314*T*(C1+C2))-K6*CCH3I*C1*C2*CRh*PDMEi/CLiI/(PDMEi-8.314*T*(C1+C2))

            c=d_c2*ti

            if j=1 then

                x_x2=C_2+ee1*c

                y_y2=y_y2+ee2*c/3

            elseif j=2 then

                x_x2=C_2+ee2*c

                y_y2=y_y2+ee3*c/3

            elseif j=3 then

                x_x2=C_2+ee3*c

                y_y2=y_y2+ee4*c/3

            elseif j=4 then

                x_x2=C_2+ee4*c

                y_y2=y_y2+ee5*c/3

            end if

        Next

        C_1=y_y1

        C_2=y_y2

  Next

  C1out(iter)=y_y1

  C2out(iter)=y_y2

Next

ENDPROGRAM;

Data;

453.00         150649.68         126.76         0.41         3049350.32         37.31         19.32         8.80

453.00         150649.68         126.76         0.83         3049350.32         37.31         13.00         19.28

453.00         150649.68         126.76         1.24         3049350.32         37.31         9.20         26.32

453.00         150649.68         126.76         1.66         3049350.32         37.31         8.72         29.52

453.00         150649.68         126.76         2.07         3049350.32         37.31         8.06         31.82

453.00         150649.68         126.76         2.48         3049350.32         37.31         6.95         30.41

453.00         150649.68         126.76         2.90         3049350.32         37.31         5.59         26.17

453.00         150649.68         42.25         2.07         3049350.32         37.31         6.85         13.43

453.00         150649.68         70.42         2.07         3049350.32         37.31         8.57         20.71

453.00         150649.68         84.51         2.07         3049350.32         37.31         10.04         26.91

453.00         150649.68         126.76         2.07         3049350.32         37.31         8.06         31.82

453.00         150649.68         169.01         2.07         3049350.32         37.31         6.47         29.89

453.00         150649.68         211.27         2.07         3049350.32         37.31         5.34         29.84

453.00         150649.68         253.52         2.07         3049350.32         37.31         4.40         29.96

453.00         150649.68         126.76         2.07         3049350.32         7.46         12.08         25.16

453.00         150649.68         126.76         2.07         3049350.32         14.93         10.42         27.91

453.00         150649.68         126.76         2.07         3049350.32         22.39         8.96         29.84

453.00         150649.68         126.76         2.07         3049350.32         29.85         8.45         30.95

453.00         150649.68         126.76         2.07         3049350.32         37.31         8.06         31.82

453.00         150649.68         126.76         2.07         3049350.32         44.78         7.75         32.12

453.00         150649.68         126.76         2.07         3049350.32         52.24         7.17         32.55

453.00         75324.84         126.76         2.07         3124675.16         37.31         16.12         3.88

453.00         112987.26         126.76         2.07         3087012.74         37.31         12.25         17.75

453.00         150649.68         126.76         2.07         3049350.32         37.31         8.06         31.82

453.00         188312.10         126.76         2.07         3011687.90         37.31         8.73         41.09

453.00         225974.52         126.76         2.07         2974025.48         37.31         9.48         50.40

453.00         263636.94         126.76         2.07         2936363.06         37.31         9.44         52.55

453.00         301299.36         126.76         2.07         2898700.64         37.31         7.60         43.20

453.00         150649.68         126.76         2.07         2349350.32         37.31         12.85         13.63

453.00         150649.68         126.76         2.07         2599350.32         37.31         12.80         19.71

453.00         150649.68         126.76         2.07         2849350.32         37.31         11.47         26.70

453.00         150649.68         126.76         2.07         3099350.32         37.31         7.93         31.95

453.00         150649.68         126.76         2.07         3349350.32         37.31         4.74         35.13

453.00         150649.68         126.76         2.07         3599350.32         37.31         4.02         35.82

453.00         150649.68         126.76         2.07         3849350.32         37.31         2.99         36.91

433.00         143998.48         126.76         2.07         3056001.52         37.31         13.52         18.40

438.00         145661.28         126.76         2.07         3104338.72         37.31         12.55         20.47

443.00         147324.08         126.76         2.07         3052675.92         37.31         11.55         22.17

448.00         148986.88         126.76         2.07         3101013.12         37.31         9.59         28.23

453.00         150649.68         126.76         2.07         3049350.32         37.31         8.06         31.82

458.00         152312.48         126.76         2.07         3097687.52         37.31         6.04         31.06

463.00         153975.28         126.76         2.07         3046024.72         37.31         4.40         27.24

Root of Mean Square Error (RMSE):11.5499416460334
Sum of Square Error:11205.6967702493
Correlation Coef. (R): -0.025193150477488
R-Square: 0.000634694830981356
Determination Coef. (DC): -3.6524379110701
F-Statistic: 0.379799714772977

Parameter Best Estimate
-------------------- -------------
k_1 0.0443887227013897
k_2 1.48988141904544
k_3 0.15806362988123
k_4 0.551727467459332
k_5 1.27971229562335
k_6 3.90080527521901
e1 34064.6206926927
e2 33558.1459707385
e3 46339.7668465127
e4 32936.9747130519
e5 33108.0993428242
e6 32934.5560026088

====== Output Results ======

File: Data File-1
No Target c1out Calculated c1out Target c2out Calculated c2out
1 19.32 3.72446132159712 8.8 1.45941820895285
2 13 6.53507135688917 19.28 5.43948599766212
3 9.2 8.3844585821214 26.32 10.9146624801039
4 8.72 9.57375357447746 29.52 17.4176359596032
5 8.06 10.2830886011808 31.82 24.8297659261248
6 6.95 10.8024412001823 30.41 27.2516499967275
7 5.59 11.3555309700269 26.17 39.278137463255
8 6.85 0.481409991498034 13.43 34.6582295922687
9 8.57 2.98065943110798 20.71 34.7680122723834
10 10.04 5.92518851451692 26.91 18.8384744557624
11 8.06 10.2830886011808 31.82 24.8297659261248
12 6.47 11.9449557638363 29.89 26.0025570402153
13 5.34 39.9686546702136 29.84 23.268772170038
14 4.4 25.6463940788978 29.96 36.9515358726961
15 12.08 11.2082725915566 25.16 12.5511124159842
16 10.42 10.6780692280548 27.91 17.9381468602351
17 8.96 10.4668996961607 29.84 20.9568431588697
18 8.45 10.3536541774117 30.95 23.0860941405943
19 8.06 10.2830886011808 31.82 24.8297659261248
20 7.75 10.2348920489271 32.12 26.4876454219494
21 7.17 10.2000642722205 32.55 28.5810076025036
22 16.12 10.2439751831954 3.88 6.01419821549014
23 12.25 9.95652354737421 17.75 15.1345857359885
24 8.06 10.2830886011808 31.82 24.8297659261248
25 8.73 10.1175326870607 41.09 20.3228972672854
26 9.48 9.44413185248538 50.4 19.1642620585202
27 9.44 8.58943190637761 52.55 19.8006045683615
28 7.6 7.69302617652497 43.2 22.6514804670537
29 12.85 9.45991413360167 13.63 15.8282959515173
30 12.8 9.81660871068178 19.71 18.6483271853985
31 11.47 10.0999291047028 26.7 21.769488046076
32 7.93 10.3237446946982 31.95 25.8185269212792
33 4.74 10.532154045967 35.13 35.2973893006523
34 4.02 10.7334779879877 35.82 44.776189116222
35 2.99 32.9322212693294 36.91 83.6274114464888
36 13.52 9.63412367023005 18.4 8.29143322070249
37 12.55 10.0548474744031 20.47 11.3109586906801
38 11.55 10.2064105002571 22.17 14.4464803808188
39 9.59 10.3619096378366 28.23 19.3200897771154
40 8.06 10.2830886011808 31.82 24.8297659261248
41 6.04 10.3746296215376 31.06 36.6337375772535
42 4.4 10.0909841087201 27.24 30.7266741146468

赞一下

回复此楼

5楼2015-07-30 14:54:35

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

查看全部 5 个回答