24小时热门版块排行榜    

查看: 2858  |  回复: 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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 时间戳变了,能看出什么问题? +14 基诺咪客 2026-08-17 18/900 2026-08-20 12:13 by saint2016
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +16 Ldrop2023 2026-08-13 20/1000 2026-08-20 11:07 by xskun
[基金申请] 言之凿凿,8月21号(明天)放榜 +12 医学老男孩 2026-08-20 14/700 2026-08-20 10:53 by ran2233
[基金申请] 只有每年这种时候来逛逛小木虫 +21 yaoyewhu2008 2026-08-20 21/1050 2026-08-20 10:37 by Vivilian
[基金申请] 今天基金会出结果吗?20260819 +14 kkkl_v 2026-08-19 15/750 2026-08-20 08:21 by sunzitan
[基金申请] filecode,4个jtjc了 +11 ziyangfang 2026-08-19 13/650 2026-08-19 21:28 by aasahr
[基金申请] 今天放榜没戏了吧 +4 yuleib84 2026-08-19 5/250 2026-08-19 19:52 by hhs666
[基金申请] filecode=后面第一个是大写字母 +10 wangze12014 2026-08-14 12/600 2026-08-19 16:56 by 苦难博士
[基金申请] 2027广东省杰青 +4 奶牛小黑 2026-08-15 10/500 2026-08-19 11:02 by wanfengnew
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +6 Tide man 2026-08-14 7/350 2026-08-19 09:56 by ZJTJZ
[基金申请] 什么时候开奖? +6 CrisMessi 2026-08-18 6/300 2026-08-19 08:06 by Equinoxhua
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 6/300 2026-08-18 15:36 by wwm09
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[基金申请] 今天系统多次维护,明天很可能放榜! +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
信息提示
请填处理意见