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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 动力学方程参数估计方法
151/2返回列表12下一页
查看: 2213  |  回复: 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下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭