查看: 935  |  回复: 6
【悬赏金币】回答本帖问题,作者自制奶酪将赠送您 100 个金币

自制奶酪

铜虫 (小有名气)

[求助] 用matlab求动力学参数已有3人参与

本人MATLAB小白,最近审稿人要求做个动力学模型,我不会用MATLAB拟合求速率常数,希望木虫的大神可以教我,谢谢。
(dC_A)/dt=-k1*CA-k7*CA
(dC_B)/dt=k1*CA-k2*CB- k5*CB
(dC_C)/dt=k2*CB-k3*CC- k6*CC
(dC_D)/dt=k3*CC-k4*CD
(dC_E)/dt= k4*CD
(dC_F)/dt=k5*CB+k6*CC+ k7*CA
此外还采用Nelder-Mead 方法通过最小化目标函数优化速率常数
F_obj=∑_m▒{∑_(n=1)^(n_p)▒[C_(m,n,calc)-C_(m,n,exp) ] ^2 }

C_(m,n,calc)是给定的时间n下组分的实验浓度;C_(m,n,exp)是给定的时间n下组分的计算浓度
第一组:
T= (0.5, 1, 1.5, 2, 3, 4, 5, 6)
A= (18.7, 6.5, 0.7, 2.1, 1.7, 1.5, 1.2, 0.2)
B= (48.9, 31.3, 15.3, 6.9, 2.2, 1.8, 1.7, 0.5)
C= (19.3, 34.7, 41.2, 40.9, 37.2, 31.7, 26.1, 19.4)
D= (9.7, 19.6, 28.7, 34.5, 38.4, 38.2, 36.3, 33.4)
E= (0.2, 2.6, 5.7, 9.3, 17.1, 24.8, 30.5, 36.7)

第二组:
T= (0.5, 1, 1.5, 2, 3, 4, 5, 6)
A= (23.7, 7.2, 1.7, 1.1, 0.7, 0.5, 0.2, 0.2)
B= (55.4, 44.3, 27.3, 11.9, 1.2, 1.1, 0.7, 0.5)
C= (14.3, 33.3, 43.2, 49.3, 51.5, 44.7,36.1, 30.4)
D= (3.7, 11.3, 19.6, 26.4, 33.4, 35.2, 34.1, 32.4)
E= (0.2, 1.2, 3.7, 7.5, 13.7, 18.3, 23.2, 27.3)

第三组:
T= (0.5, 1, 1.5, 2, 3, 4, 5, 6)
A= (33.7, 15.2, 3.7, 1.1, 0.7, 0.5, 0.2, 0.2)
B= (49.7, 53.4, 44.3, 28.3, 8.7, 1.1, 0.7, 0.5)
C= (9.3, 20.8, 32.2, 43.7, 53.5, 55.7, 51.4, 46.9)
D= (2.7, 7.9, 13.7, 18.4, 25.4, 28.2, 29.4, 30.4)
E= (0, 0.2, 2.1, 4.3, 8.7, 12.3, 15.7, 19.3)
回复此楼
简单快乐
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

独孤神宇

版主 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
数值计算
2楼2021-01-05 21:04:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
楼主是想3组数据求出一组共同的参数?还是每一组数据分别求出一组参数?自己用1stOpt试试。
3楼2021-01-05 21:21:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

自制奶酪

铜虫 (小有名气)

引用回帖:
3楼: Originally posted by dingd at 2021-01-05 21:21:18
楼主是想3组数据求出一组共同的参数?还是每一组数据分别求出一组参数?自己用1stOpt试试。

每一组分别求出一组参数k
简单快乐
4楼2021-01-05 21:41:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

自制奶酪

铜虫 (小有名气)

引用回帖:
2楼: Originally posted by 独孤神宇 at 2021-01-05 21:04:32
http://blog.sina.com.cn/s/blog_c0cb8ce60102ysqu.html

您好,链接打开加密了
简单快乐
5楼2021-01-05 21:42:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hzlhm

至尊木虫 (著名写手)

【答案】应助回帖

F_obj=∑_m▒{∑_(n=1)^(n_p)▒[C_(m,n,calc)-C_(m,n,exp) ] ^2 }这个关系式不明确,请重新书写。
QQ:2120156492
6楼2021-01-11 19:23:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

梅子炒米粉

新虫 (小有名气)

请问楼主问题解决了吗?方便交流交流吗?

发自小木虫Android客户端
7楼2021-04-24 21:52:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 自制奶酪 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 投票:  2022年8月10日,申请书的文件名来投票啦~~~~~~ +24 mmffdd26 2022-08-10 28/1400 2022-08-16 12:08 by derickofaze
[基金申请] 我也曾梦想仗剑走天涯,淡定,淡定! +3 了却无痕 2022-08-16 5/250 2022-08-16 12:04 by 了却无痕
[基金申请] i科研人 +9 Tmfds 2022-08-16 10/500 2022-08-16 11:41 by Tmfds
[基金申请] 【2022.8.16】刚打开JJW网站变成了400 Bad Request +21 kiyohawa 2022-08-16 30/1500 2022-08-16 10:53 by yaozizu
[找工作] 想咨询了解军队文职 +6 冀林昊 2022-08-12 13/650 2022-08-16 10:11 by 冀林昊
[基金申请] 国科办:不以承担科研项目和经费多少作为评价科研人员的标准 +42 X__X 2022-08-14 44/2200 2022-08-16 09:47 by daydreamren
[基金申请] 上会了,但是会评结束后大老板说给打听但是没消息,各位看看是怎么回事 24+4 好名字、 2022-08-14 22/1100 2022-08-16 09:11 by ljs_liu
[基金申请] 建议直接取消会评 +42 ffaaa 2022-08-11 48/2400 2022-08-16 08:53 by medcheng
[博后之家] 香港中文大学(深圳)2022年诚聘海内外人才 +3 ryutryi 2022-08-15 3/150 2022-08-15 23:39 by quwyrcb
[基金申请] 有收到联合基金答辩通知的吗? +4 weekend2016 2022-08-12 6/300 2022-08-15 14:21 by weekend2016
[基金申请] 写在国自然放榜前 +31 ZHANGxuefen 2022-08-13 43/2150 2022-08-15 11:44 by puluoman66
[有机交流] 脱甲氧基反应求助 +5 youwenhui 2022-08-13 6/300 2022-08-15 10:43 by 紫枫星域
[基金申请] 博后申请面上 +19 tswl 2022-08-11 41/2050 2022-08-15 10:42 by cjhaoren
[硕博家园] 邀人摇人啦!俺新建了『zotero互助学习群』 +5 Xiaokang2022 2022-08-13 6/300 2022-08-13 23:24 by Finnhu
[论文投稿] 论文版面费请教 +3 Fanninger 2022-08-13 5/250 2022-08-13 15:24 by 梧桐77
[论文投稿] 对自己科研能力不自信,结果出现奇迹 +8 Xzhang-23 2022-08-10 10/500 2022-08-12 20:04 by 选好路
[教师之家] 五部门联合发文要给青椒减负了 +11 wwwkkk83 2022-08-10 12/600 2022-08-12 20:04 by xiatc1981
[基金申请] 有时候是无趣的! +11 了却无痕 2022-08-12 17/850 2022-08-12 14:14 by 了却无痕
[基金申请] 请问公布中标课题的同时会不会同时公布基金号码 +10 爱生活爱拉登 2022-08-10 14/700 2022-08-11 10:40 by 爱生活爱拉登
[基金申请] +6 IEEEYLQ 2022-08-09 9/450 2022-08-10 11:30 by 了却无痕
信息提示
请填处理意见