24小时热门版块排行榜    

CyRhmU.jpeg
查看: 2229  |  回复: 14
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

thundernet8

木虫 (小有名气)

[求助] 动力学方程参数估计方法


如图片中所示,此类动力学方程如何进行准确的参数回归呢,其中一种方法:微分方法先拟合反应物浓度对时间曲线,再根据此拟合出反应速度对反应时间曲线,然后用最小二乘法求动力学参数,这个感觉不太准确,再者或者用一阶微分方程?有过动力学分析的友友请帮忙,先谢了
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

thundernet8

木虫 (小有名气)


fegg7502: 金币+1, 鼓励交流 2012-11-08 13:58:18
引用回帖:
6楼: Originally posted by dingd at 2012-11-07 08:40:00
把公式、数据等贴出来,看能否帮你跑下。

-r=-dC/dt=k0*e^(-Ea/RT)*(A^m)*(B^n)*(C^o)
单次实验T,A,B都是固定的,而C随时间变化并测量出来,所以T,A,B可以说是可变常数,即单个实验的数据就可以回归o值,但是这样弄一起回归好吗,所有实验数据如下
T(℃)       A       B(MPa)      t/min     C
260        0.05        4        0        3759.7
260        0.05        4        10        2772.45
260        0.05        4        20        2100
260        0.05        4        30        1679.03
260        0.05        4        40        1356.12
260        0.05        4        50        1122.72
260        0.05        4        60        818.6
260        0.05        4        70        643.28
260        0.05        4        80        485.26
260        0.05        4        90        361.16
260        0.05        4        100        290.17
260        0.1        4        0        3759.7
260        0.1        4        10        2450
260        0.1        4        20        1803.37
260        0.1        4        30        1246.87
260        0.1        4        40        895
260        0.1        4        50        496.14
260        0.1        4        60        410.13
260        0.1        4        70        179.72
260        0.1        4        80        142.38
260        0.1        4        90        94.73
260        0.1        4        100        93.51
260        0.15        4        0        3759.7
260        0.15        4        10        2300
260        0.15        4        20        1424.68
260        0.15        4        30        884.85
260        0.15        4        40        546.36
260        0.15        4        50        336.83
260        0.15        4        60        203.33
260        0.15        4        70        117.97
260        0.15        4        80        66.95
260        0.15        4        90        54.49
260        0.15        4        100        49.41
260        0.2        4        0        3759.7
260        0.2        4        10        2150
260        0.2        4        20        1200
260        0.2        4        30        650
260        0.2        4        40        380
260        0.2        4        50        230
260        0.2        4        60        130
260        0.2        4        70        80
260        0.2        4        80        65
260        0.2        4        90        50
260        0.2        4        100        45
200        0.1        4        0        3759.7
200        0.1        4        10        2893.56
200        0.1        4        20        2037.91
200        0.1        4        30        1319.34
200        0.1        4        40        1005.13
200        0.1        4        50        773.12
200        0.1        4        60        562.77
200        0.1        4        70        384.27
200        0.1        4        80        267.1
200        0.1        4        90        171.38
200        0.1        4        100        157.51
220        0.1        4        0        3759.7
220        0.1        4        10        2600
220        0.1        4        20        1691.77
220        0.1        4        30        1126.16
220        0.1        4        40        902.03
220        0.1        4        50        628.67
220        0.1        4        60        371.19
220        0.1        4        70        240.16
220        0.1        4        80        139.84
220        0.1        4        90        116.87
220        0.1        4        100        83.8
240        0.1        4        0        3759.7
240        0.1        4        10        2350.51
240        0.1        4        20        1450
240        0.1        4        30        1000
240        0.1        4        40        708.77
240        0.1        4        50        375.68
240        0.1        4        60        271.72
240        0.1        4        70        175.88
240        0.1        4        80        123.3
240        0.1        4        90        90
240        0.1        4        100        66.69
260        0.1        4        0        3759.7
260        0.1        4        10        2080
260        0.1        4        20        1200
260        0.1        4        30        884.85
260        0.1        4        40        546.36
260        0.1        4        50        336.83
260        0.1        4        60        203.33
260        0.1        4        70        117.97
260        0.1        4        80        66.95
260        0.1        4        90        54.49
260        0.1        4        100        49.41
280        0.1        4        0        3759.7
280        0.1        4        10        1935.79
280        0.1        4        20        1007.44
280        0.1        4        30        758.08
280        0.1        4        40        500
280        0.1        4        50        315.13
280        0.1        4        60        190.44
280        0.1        4        70        100
280        0.1        4        80        50
280        0.1        4        90        40
280        0.1        4        100        30
260        0.1        4        0        3759.7
260        0.1        4        10        2080
260        0.1        4        20        1200
260        0.1        4        30        749.04
260        0.1        4        40        500
260        0.1        4        50        280
260        0.1        4        60        203.33
260        0.1        4        70        117.97
260        0.1        4        80        66.95
260        0.1        4        90        54.49
260        0.1        4        100        49.41
260        0.1        3.5        0        3759.7
260        0.1        3.5        10        2244.58
260        0.1        3.5        20        1484.05
260        0.1        3.5        30        749.04
260        0.1        3.5        40        547.23
260        0.1        3.5        50        331.18
260        0.1        3.5        60        267.75
260        0.1        3.5        70        198.36
260        0.1        3.5        80        130.51
260        0.1        3.5        90        100.47
260        0.1        3.5        100        93.91
260        0.1        3        0        3759.7
260        0.1        3        10        2483.29
260        0.1        3        20        1669.27
260        0.1        3        30        1162.42
260        0.1        3        40        827.86
260        0.1        3        50        517.85
260        0.1        3        60        357.71
260        0.1        3        70        274.21
260        0.1        3        80        185.53
260        0.1        3        90        129.69
260        0.1        3        100        88.35
260        0.1        2.5        0        3759.7
260        0.1        2.5        10        2527.4
260        0.1        2.5        20        1857.92
260        0.1        2.5        30        1219.77
260        0.1        2.5        40        850.2
260        0.1        2.5        50        628.67
260        0.1        2.5        60        515.01
260        0.1        2.5        70        359.75
260        0.1        2.5        80        265.34
260        0.1        2.5        90        172.36
260        0.1        2.5        100        149.93
260        0.1        2        0        3759.7
260        0.1        2        10        2722.63
260        0.1        2        20        2089.13
260        0.1        2        30        1485.29
260        0.1        2        40        1062.05
260        0.1        2        50        831.48
260        0.1        2        60        766.41
260        0.1        2        70        547.1
260        0.1        2        80        461.41
260        0.1        2        90        347.28
260        0.1        2        100        244.06
7楼2012-11-07 15:15:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 15 个回答

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
csgt0: 金币+1, 多谢指导 2012-11-02 10:27:15
thundernet8: 金币+1, 有帮助, 感谢回复 2012-11-02 10:37:22
可归结为微分方程拟合问题,1stOpt可以很好地处理,推荐试试。
2楼2012-11-02 08:43:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

thundernet8

木虫 (小有名气)


fegg7502: 金币+1, 鼓励交流 2012-11-08 13:58:00
引用回帖:
2楼: Originally posted by dingd at 2012-11-02 08:43:11
可归结为微分方程拟合问题,1stOpt可以很好地处理,推荐试试。

1stopt用哪种优化算法呢,代码什么的能给个提示吗,谢谢
3楼2012-11-02 19:29:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zt_chem

木虫 (正式写手)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
thundernet8: 金币+1, 感谢参与 2012-11-05 11:25:46
csgt0: 金币+1, 欢迎交流 2012-11-08 10:13:07
接着楼上,现在1stopt1.0和1.5版本的不能拟合ode,能拟合ode的在网上找不到啊。我用matlab编程拟合,结果如楼主所说:误差很大!不知道该怎样解决,期待高手!
4楼2012-11-04 19:09:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见