24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 2464  |  回复: 14
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

thundernet8

木虫 (小有名气)

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


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

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

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

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

10楼2012-11-08 09:25:37
已阅   回复此楼   关注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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 考研调剂 +9 小蜡新笔 2026-03-29 10/500 2026-03-31 19:52 by Dyhoer
[考研] 已决定调剂院校 +8 JKSOIID 2026-03-26 8/400 2026-03-31 19:51 by mg1014
[硕博家园] 求调剂 有机化学考研356分 +11 Nadiums 2026-03-25 12/600 2026-03-31 18:24 by oooqiao
[考研] 本2一志愿C9-333分,材料科学与工程,求调剂 +9 升升不降 2026-03-31 9/450 2026-03-31 18:01 by 无际的草原
[考研] 070300化学354求调剂 +15 101次希望 2026-03-28 15/750 2026-03-31 17:58 by jp9609
[考研] 英一数一总分334求调剂 +4 陈阳坤 2026-03-31 4/200 2026-03-31 14:22 by 记事本2026
[考研] 调剂310 +13 温柔的晚安 2026-03-25 14/700 2026-03-31 13:03 by 记事本2026
[考研] 一志愿大连理工大学,机械工程学硕,341 +3 西瓜田的守望者 2026-03-30 3/150 2026-03-31 11:08 by asdfzly
[考研] 334分 一志愿武理 材料求调剂 +16 李李不服输 2026-03-26 16/800 2026-03-30 17:47 by wangjy2002
[考研] 0703本科郑州大学求调剂 +7 nhj_ 2026-03-25 7/350 2026-03-30 12:44 by fangnagu
[考研] 085602 化学工程专硕 340分求调剂 +4 qianbai11 2026-03-29 4/200 2026-03-30 11:34 by 唐沐儿
[考研] 318一志愿吉林大学生物与医药 求调剂 +5 笃行致远. 2026-03-28 5/250 2026-03-30 06:56 by ilovexiaobin
[考研] 总分293求调剂 +8 加一一九 2026-03-25 11/550 2026-03-29 19:53 by 无际的草原
[考研] 一志愿北京理工大学本科211材料工程294求调剂 +8 mikasa的围巾 2026-03-28 8/400 2026-03-29 12:48 by 无际的草原
[考研] 316求调剂 +7 江辞666 2026-03-26 7/350 2026-03-28 21:28 by sanrepian
[考研] 295求调剂 +5 1428151015 2026-03-27 6/300 2026-03-28 04:04 by fmesaito
[考研] 求调剂 +4 零八# 2026-03-27 4/200 2026-03-27 18:07 by yu221
[考研] 272求调剂 +7 脚滑的守法公民 2026-03-27 7/350 2026-03-27 17:23 by laoshidan
[考研] 305求调剂 +5 哇卢卡库 2026-03-26 5/250 2026-03-27 14:01 by laoshidan
[考研] 求调剂 +3 李李不服输 2026-03-25 3/150 2026-03-25 13:03 by cmz0325
信息提示
请填处理意见