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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 动力学方程参数估计方法
北京石油化工学院2026年研究生招生接收调剂公告
151/2返回列表12下一页
查看: 2474  |  回复: 14
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

thundernet8

木虫 (小有名气)

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


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

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2012-11-01 22:24:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖置顶 ( 共有1个 )

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
csgt0: 金币+2, 谢谢 2012-11-08 10:13:40
thundernet8: 金币+8, ★★★★★最佳答案, 非常感谢^v^ 2012-11-08 11:05:48
fegg7502: 金币+1, 应助指数+1, 鼓励交流 2012-11-08 13:58:56
thundernet8: 回帖置顶 2012-11-10 15:13:58
1stOpt求解：
CODE:
Constant R=8.314;
VarConstant T1=[260,260,260,260,200,220,240,260,280,260,260,260,260,260],
            A=[0.05,0.1,0.15,0.2,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1],
            B=[4,4,4,4,4,4,4,4,4,4,3.5,3,2.5,2];
Variable t,C;
ODEFunction C'=-k0*exp(-Ea/(R*T1))*(A^m)*(B^n)*(C^o);
Data;
0        3759.7
10        2772.45
20        2100
30        1679.03
40        1356.12
50        1122.72
60        818.6
70        643.28
80        485.26
90        361.16
100        290.17
Data;
0        3759.7
10        2450
20        1803.37
30        1246.87
40        895
50        496.14
60        410.13
70        179.72
80        142.38
90        94.73
100        93.51
Data;
0        3759.7
10        2300
20        1424.68
30        884.85
40        546.36
50        336.83
60        203.33
70        117.97
80        66.95
90        54.49
100        49.41
Data;
0        3759.7
10        2150
20        1200
30        650
40        380
50        230
60        130
70        80
80        65
90        50
100        45
Data;
0        3759.7
10        2893.56
20        2037.91
30        1319.34
40        1005.13
50        773.12
60        562.77
70        384.27
80        267.1
90        171.38
100        157.51
Data;
0        3759.7
10        2600
20        1691.77
30        1126.16
40        902.03
50        628.67
60        371.19
70        240.16
80        139.84
90        116.87
100        83.8
Data;
0        3759.7
10        2350.51
20        1450
30        1000
40        708.77
50        375.68
60        271.72
70        175.88
80        123.3
90        90
100        66.69
Data;
0        3759.7
10        2080
20        1200
30        884.85
40        546.36
50        336.83
60        203.33
70        117.97
80        66.95
90        54.49
100        49.41
Data;
0        3759.7
10        1935.79
20        1007.44
30        758.08
40        500
50        315.13
60        190.44
70        100
80        50
90        40
100        30
Data;
0        3759.7
10        2080
20        1200
30        749.04
40        500
50        280
60        203.33
70        117.97
80        66.95
90        54.49
100        49.41
Data;
0        3759.7
10        2244.58
20        1484.05
30        749.04
40        547.23
50        331.18
60        267.75
70        198.36
80        130.51
90        100.47
100        93.91
Data;
0        3759.7
10        2483.29
20        1669.27
30        1162.42
40        827.86
50        517.85
60        357.71
70        274.21
80        185.53
90        129.69
100        88.35
Data;
0        3759.7
10        2527.4
20        1857.92
30        1219.77
40        850.2
50        628.67
60        515.01
70        359.75
80        265.34
90        172.36
100        149.93
Data;
0        3759.7
10        2722.63
20        2089.13
30        1485.29
40        1062.05
50        831.48
60        766.41
70        547.1
80        461.41
90        347.28
100        244.06

均方差(RMSE): 122.451713151013
残差平方和(SSE): 2099219.08750651
相关系数(R): 0.986538044771092
相关系数之平方(R^2): 0.973257313780769
决定系数(DC): 0.97324741144721
F统计(F-Statistic): 140.696085575023

参数                  最佳估算
--------------------        -------------
k0        0.126870111743644
ea        2526.23497262219
m        0.594482756783444
n        0.585283294110934
o        1.09436923590948

jv1.jpg
赞一下(1人) 回复此楼
10楼2012-11-08 09:25:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★
感谢参与，应助指数 +1
csgt0: 金币+1, 多谢指导 2012-11-02 10:27:15
thundernet8: 金币+1, ★有帮助, 感谢回复 2012-11-02 10:37:22
可归结为微分方程拟合问题，1stOpt可以很好地处理，推荐试试。
赞一下(1人) 回复此楼
2楼2012-11-02 08:43:11
已阅   回复此楼   关注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编程拟合，结果如楼主所说：误差很大！不知道该怎样解决，期待高手！
赞一下(1人) 回复此楼
4楼2012-11-04 19:09:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)
★
fegg7502: 金币+1, 鼓励交流 2012-11-08 13:57:53
均方差(RMSE): 122.115033528798
残差平方和(SSE): 2087691.39792352
相关系数(R): 0.986612411971849
相关系数之平方(R^2): 0.973404051456909
决定系数(DC): 0.973394321094808
F统计(F-Statistic): 140.194937516001

参数                  最佳估算
--------------------        -------------
k0        0.591058721457601
ea        11993.0119493646
m        0.594372666801216
n        0.584274660908677
o        1.09421452986405
赞一下(1人) 回复此楼
12楼2012-11-08 11:43:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

forcal

禁虫 (初入文坛)
本帖内容被屏蔽

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

thundernet8
14楼2012-11-10 09:32:07
已阅   回复此楼   关注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的回帖

thundernet8

木虫 (小有名气)
有高人指点吗
回复此楼
5楼2012-11-06 21:55:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

把公式、数据等贴出来，看能否帮你跑下。
赞一下 回复此楼
6楼2012-11-07 08:40:00
已阅   回复此楼   关注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的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★
fegg7502: 金币+1, 应助指数+1, 鼓励交流 2012-11-08 13:58:30
“-Ea/RT”等于“-Ea/R*T”还是“-Ea/(R*T)”？

Ea、R是待求参数还是常数？如果是常数的话是多少？

k0、m、n、o是待求参数吧？
赞一下 回复此楼
8楼2012-11-07 15:49:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

thundernet8

木虫 (小有名气)
引用回帖:
8楼: Originally posted by dingd at 2012-11-07 15:49:37
“-Ea/RT”等于“-Ea/R*T”还是“-Ea/(R*T)”？

Ea、R是待求参数还是常数？如果是常数的话是多少？

k0、m、n、o是待求参数吧？

R=8.3145 Ea,ko,m,n,o待求            是-Ea/(R*T)
赞一下 回复此楼
9楼2012-11-07 17:12:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 thundernet8 的主题更新
151/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 调剂 +5 不逢春 2026-04-05 6/300 2026-04-05 21:50 by 不逢春
[考研] 一志愿郑州大学085600求调剂 +10 吃的不少 2026-04-05 12/600 2026-04-05 21:38 by 醉翁wl
[考研] 材料0856 英一数二 323 求调剂 +14 袁sy 2026-04-01 14/700 2026-04-05 18:18 by cql1109
[考研] 一志愿上海海洋大学083200食品学硕，求调剂，接受其他专业083200 +4 what张 2026-04-04 5/250 2026-04-05 14:07 by chw1980_0
[考研] 生物学308分求调剂（一志愿华东师大） +8 相信必会光芒万� 2026-04-05 10/500 2026-04-05 12:19 by Hdyxbekcb
[考研] 环境科学与工程334分求调剂 +9 王一一依依 2026-03-30 12/600 2026-04-04 20:55 by dongzh2009
[考研] 359求调剂 +7 hhhhaaaa$ 2026-04-04 7/350 2026-04-04 18:49 by imissbao
[考研] 一志愿重庆大学085404，总分314分，求调剂 +4 zf83hn 2026-04-03 4/200 2026-04-03 21:25 by 啵啵啵0119
[考研] 301求调剂 +14 A_JiXing 2026-04-01 14/700 2026-04-03 18:31 by ls刘帅
[考研] 315分 085602 求调剂 +15 26考研上岸版26 2026-04-02 15/750 2026-04-03 12:45 by xingguangj
[考研] 08工科275分求调剂 +14 AaAa7420 2026-03-31 14/700 2026-04-03 11:13 by cocolv
[考研] 285求调剂 +7 AZMK 2026-04-02 9/450 2026-04-03 11:12 by wanwan00
[考研] 326求调剂 +10 崽崽仔 2026-04-02 10/500 2026-04-03 09:08 by 帕尔马拉特
[考研] 260求调剂 +3 朱芷琳 2026-04-02 3/150 2026-04-03 08:44 by yulian1987
[考研] 326求调剂 +3 9ahye 2026-04-02 4/200 2026-04-03 08:43 by Jaylen.
[考研] 食品学硕362求调剂 +3 xuanxianxian 2026-04-01 3/150 2026-04-01 21:05 by 啊李999
[考研] 307分求调剂 +14 (o~o) 2026-03-31 15/750 2026-04-01 20:43 by longlotian
[考研] 生物与医药考研调剂 +5 铁憨憨123425 2026-03-31 5/250 2026-04-01 18:01 by syh9288
[考研] 一志愿西交大080500材料学硕349 +6 jqx1258 2026-03-31 7/350 2026-03-31 21:08 by yuq
[考研] 求调剂 +8 11ggg 2026-03-30 8/400 2026-03-31 13:56 by nanaliuyun
信息提示
关闭
请填处理意见
关闭