24小时热门版块排行榜    

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

小小莓邮

铁虫 (小有名气)

[求助] 拜求!matlab拟合动力学参数 已有2人参与

实验数据如下
t(h)[RH]              [ROOH]
0        8.32016                0
1        7.37673                0.21115
2        6.74433                0.7379
3        5.92613                1.10562
4        5.68806                1.29586
5        5.01584                1.58327
6        6.09249                1.36404
8        5.34827                1.84065
10        4.9961                2.00986
动力学模型如图
其中[o2]浓度未测量但认为恒定,想把K4[o2]看成一个参数,并估计出其他的k1 k1‘ K4 K5 K8 f 值,matlab一窍不通,看了几天了还是不明所以,所以恳请各位大神给拟合一下,并给出代码,因为还要拟合其他的,所以想大神给出一个后面好跟着做。。。
拜求!matlab拟合动力学参数
捕获.JPG
回复此楼

» 猜你喜欢

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

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

dingd

铁杆木虫 (职业作家)

不知理解的对否,用1stOpt求解:
CODE:
Variable t, RH, ROOH;
ODEFunction RH'=-(k4/(k4+k5*RH)*k5*RH*sqrt(k1*ROOH/k8)+2*f*k1*ROOH);
Data;
//t(h)[RH]              [ROOH]
0        8.32016                0
1        7.37673                0.21115
2        6.74433                0.7379
3        5.92613                1.10562
4        5.68806                1.29586
5        5.01584                1.58327
6        6.09249                1.36404
8        5.34827                1.84065
10        4.9961                2.00986

目标函数值稳定但参数组值不唯一:
1:
均方差(RMSE):0.337408531472899
残差平方和(SSE):0.910756136885585
相关系数(R): 0.901867531322994
相关系数之平方(R^2): 0.813365044054631
确定系数(DC): 0.81334823805307
F统计(F-Statistic): 4.01817798116069

参数                  最佳估算
--------------------        -------------
k4        -0.00960366653113242
k5        0.00108370577772373
k1        -3.63159583413505
k8        -0.00354975409632072
f        0.0406346976895024

2:
均方差(RMSE):0.337408531472899
残差平方和(SSE):0.910756136885585
相关系数(R): 0.901867531322902
相关系数之平方(R^2): 0.813365044054465
确定系数(DC): 0.81334823805307
F统计(F-Statistic): 4.01817798116068

参数                  最佳估算
--------------------        -------------
k4        -0.000186819014387768
k5        2.10812031658941E-5
k1        -0.254411174415381
k8        -9.41033214347948E-8
f        0.580040555144359

有过拟合现象。其实微分方程可简化为:RH'=-(k8/(1+k5*RH)*RH*sqrt(ROOH)+f*ROOH);此时参数仅3个,结果稳定唯一:
均方差(RMSE):0.337408531472898
残差平方和(SSE):0.910756136885582
相关系数(R): 0.901867531324836
相关系数之平方(R^2): 0.813365044057954
确定系数(DC): 0.813348238053071
F统计(F-Statistic): 13.393926603869

参数                  最佳估算
--------------------        -------------
k8        0.0346625904471096
k5        -0.112842920106452
f        -0.295137602915853
拜求!matlab拟合动力学参数-1
f1.jpg

6楼2014-02-19 10:54:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 15 个回答

钟静

新虫 (初入文坛)

有几点不懂哈:
1.等式右边的是什么,-r_RH?怎么没给出来?直接把原模型放出来看看啦!在等式里面也没有体现时间,可见这个式子和数据,给出的数据和给出的模型等式不匹配啊兄弟!
2.[o2]浓度恒定,你居然不知道数值?测量困难?多一个参数对其他参数影响也很大呀!
3.最小二乘法啦,Matlab的lsqcurvefit函数可以满足你的计算要求!
2楼2014-02-18 16:54:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小小莓邮

铁虫 (小有名气)

引用回帖:
2楼: Originally posted by 钟静 at 2014-02-18 16:54:32
有几点不懂哈:
1.等式右边的是什么,-r_RH?怎么没给出来?直接把原模型放出来看看啦!在等式里面也没有体现时间,可见这个式子和数据,给出的数据和给出的模型等式不匹配啊兄弟!
2.浓度恒定,你居然不知道数值? ...

1、等式左边-rRH=-d[RH]/dt 速率表示
2、氧气浓度目前没法测,但认为其浓度恒定,主要是为了确定速率与几种产物的关系,并不是深究每个k的真实值,所以才认为k4[O2]可以看做一个常数
3、matlab真心不会,在坛子里看了好几天了,也看到了一下案例及代码,但还是整不明白,也写不出来,才求教各位大神,先看看能不能帮忙写出这个得代码,后面的自己再慢慢来!
小弟先在此谢过了!
3楼2014-02-18 19:00:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jerkwin

专家顾问 (正式写手)

你的模型没说全, 单凭你这几个数据不可能做出来
提问之前连最基本的背景都不说, 你以为别人都是神, 知道你要的是什么?
4楼2014-02-18 22:54:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 考研二轮调剂 +7 故人?? 2026-04-09 7/350 2026-04-10 07:49 by 314126402
[考研] 282,求调剂 +9 jggshjkkm 2026-04-09 10/500 2026-04-10 07:26 by swb0218
[考研] 086004 求调剂 309 +7 Yin DY 2026-04-08 7/350 2026-04-09 13:59 by Delta2012
[考研] 一志愿华南师范大学0702物理学305调剂 +4 念常安 2026-04-07 6/300 2026-04-08 22:53 by bljnqdcc
[考研] 求调剂 +13 柒luck 2026-04-07 13/650 2026-04-08 22:46 by 猪会飞
[考研] 材料与化工专硕306分找合适调剂 +27 沧海轻舟e 2026-04-06 28/1400 2026-04-08 22:06 by wdyheheeh
[考研] 调剂 +3 电气300求调剂不 2026-04-08 6/300 2026-04-08 09:39 by 电气300求调剂不
[考研] 生物医药调剂|SCI中科院三区一作+多项科研成果 +8 likangxing 2026-04-07 11/550 2026-04-08 00:02 by lys0704
[考研] 农学,求调剂,314分 +4 访客记录可爱 2026-04-04 4/200 2026-04-07 21:07 by 等岸
[考研] 一志愿西南090202求调剂 +4 在线求有学上 2026-04-07 4/200 2026-04-07 19:47 by biomichael
[考研] 一志愿太原理工大学计算机技术专硕348,求调剂指导 +3 nexious 2026-04-05 3/150 2026-04-07 08:19 by jp9609
[考研] 081200-11408-367学硕求调剂 +4 1_2_3111 2026-04-06 4/200 2026-04-07 08:13 by jp9609
[考研] 生物与医药求调剂 +7 heguanhua 2026-04-05 8/400 2026-04-06 18:41 by macy2011
[考研] 362求调剂一志愿中国石油大学 +4 我要考大 2026-04-06 6/300 2026-04-06 14:11 by 无际的草原
[考研] 一志愿武汉理工大学080200机械工程308分,求调剂 +4 终不似从前 2026-04-05 4/200 2026-04-06 11:46 by 考研学校招点人
[考研] 296求调剂 +3 汪!?! 2026-04-05 4/200 2026-04-05 20:13 by 啵啵啵0119
[考研] 313求调剂 +3 海日海日 2026-04-04 3/150 2026-04-05 07:48 by 544594351
[考研] 一志愿东北大学085901土木专硕345求调剂 +3 zxt11111 2026-04-04 3/150 2026-04-04 14:21 by 土木硕士招生
[考研] 考研调剂 +3 15615482637 2026-04-03 3/150 2026-04-03 22:50 by ms629
[考研] 一志愿重庆大学085404,总分314分,求调剂 +4 zf83hn 2026-04-03 4/200 2026-04-03 21:25 by 啵啵啵0119
信息提示
请填处理意见