查看: 933  |  回复: 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的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

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

独孤神宇

版主 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
数值计算
2楼2021-01-05 21:04:32
已阅   回复此楼   关注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的回帖
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 沈阳理工大学和德州学院?该选哪一个,纠结 +38 liu201404 2022-08-02 83/4150 2022-08-09 03:08 by liu201404
[教师之家] 大家把钱放在哪儿能稳定性理财? +24 zju2000 2022-08-02 28/1400 2022-08-09 01:00 by 若花而绽
[论文投稿] 论文小修返稿已经半个月了 +13 六月vvv 2022-08-04 25/1250 2022-08-09 00:19 by 六月vvv
[基金申请] 国家自然基金交流-F02分群 +8 lauren2018 2022-08-03 8/400 2022-08-08 20:04 by 流水不争
[高分子] 乙酸乙酯脱除 30+4 ljaxida 2022-08-04 6/300 2022-08-08 15:31 by 无视军团
[分析] ppm换算 5+3 spring69563 2022-08-07 7/350 2022-08-08 12:42 by 国平荣兴
[公派出国] 有准备明年申请的小伙伴吗? +3 nuogai1161 2022-08-03 4/200 2022-08-08 11:00 by shui9769
[基金申请] 基金过了 +62 zjjcn2001 2022-08-03 75/3750 2022-08-08 10:09 by zjjcn2001
[基金申请] 有做modelling的同仁么有? +7 tswl 2022-08-06 15/750 2022-08-08 00:22 by juta256
[基金申请] 科研影响力已有两次阅读 +7 jane_012 2022-08-07 8/400 2022-08-07 17:36 by lbjsd
[考博] 博导回复 +4 EmmaaaXie 2022-08-05 6/300 2022-08-07 15:08 by cwp0314
[论文投稿] JCR 2区 《Coatings》免费发表SCI论文 +8 dong_liu 2022-08-06 8/400 2022-08-07 10:00 by aasahr
[基金申请] 基金官宣前,科研处等小道消息准确性有多少? +15 风中劲草哦 2022-08-03 25/1250 2022-08-06 01:31 by itpfeng
[精细化工] 大佬们,请问有什么增稠剂能耐次氯酸钠呀? +4 domona 2022-08-02 5/250 2022-08-05 17:28 by wushuijin
[基金申请] 寒假写本子,累;暑假等结果,更累。 +11 wangjhoo 2022-08-03 14/700 2022-08-05 17:11 by deluck
[催化] 请问有人知道哪里能买到负载型的铂催化剂吗? +3 a379151515 2022-08-02 6/300 2022-08-05 15:17 by servicebio
[基金申请] 关于科研处的消息 +21 Vivien90 2022-08-03 23/1150 2022-08-05 13:48 by 视界有你
[硕博家园] 南通大学2022年招聘公告( 招聘快结束了,速报名!) +3 wusidun 2022-08-03 3/150 2022-08-04 09:20 by yzlxl88
[基金申请] 国家自然科学基金有没有第二申请代码中的 +3 zcfnm 2022-08-03 6/300 2022-08-03 21:54 by zcfnm
[基金申请] 第一年,没中 +12 噢噢噢00 2022-08-03 16/800 2022-08-03 15:59 by wsjing
信息提示
请填处理意见