查看: 934  |  回复: 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个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 公布时间,已排除8月15日 +15 风中劲草哦 2022-08-16 20/1000 2022-08-16 12:13 by Aa3493008
[基金申请] 投票:  2022年8月10日,申请书的文件名来投票啦~~~~~~ +24 mmffdd26 2022-08-10 28/1400 2022-08-16 12:08 by derickofaze
[基金申请] 基金结果得9月份才能公布了 +42 sunbiaoex 2022-08-15 57/2850 2022-08-16 11:17 by Doncic1989
[基金申请] 建议直接取消会评 +42 ffaaa 2022-08-11 48/2400 2022-08-16 08:53 by medcheng
[基金申请] 看到的,都很有可能中了! +12 紫色男 2022-08-15 13/650 2022-08-16 03:16 by joey052
[基金申请] 现在还是 没有消息就是最坏的消息的 观点吗? +4 wei_meng 2022-08-15 4/200 2022-08-15 23:36 by zp89
[基金申请] 科研之友心电图到底有用没 +21 creativekid 2022-08-14 29/1450 2022-08-15 15:05 by 秦苏辽
[基金申请] 验证了下 +15 410088 2022-08-14 21/1050 2022-08-15 14:49 by sjdbx
[基金申请] E04口青年基金函评有几个人? +8 foolishmani 2022-08-09 9/450 2022-08-15 13:15 by 黑白渡口
[论文投稿] 小修后需要多长时间才能有决定 +5 Eveleen 2022-08-11 10/500 2022-08-15 13:12 by 你真可爱呢
[找工作] 还有不坑的学校吗? +13 菜鸟夭夭 2022-08-12 13/650 2022-08-15 00:42 by shenqi2
[基金申请] 转发,一个第一次参加会评老师的抱怨 +44 babu2015 2022-08-10 48/2400 2022-08-14 14:47 by mahl05
[基金申请] 焦虑没日没夜 +14 Lydia8977 2022-08-13 14/700 2022-08-14 13:33 by zyh4014
[基金申请] 刚刚登录基金系统,啥变化都没有 +44 winnerche 2022-08-13 75/3750 2022-08-14 06:42 by 黑白渡口
[教师之家] offer求交流 +9 tomcat1234 2022-08-13 20/1000 2022-08-13 17:25 by tomcat1234
[论文投稿] EP编辑部是放暑假了吗 +10 Wichen 2022-08-11 12/600 2022-08-13 08:03 by Wichen
[基金申请] 这是要公布的前兆了吧 +14 day388383 2022-08-11 16/800 2022-08-12 14:27 by dafeiji2010
[有机交流] 如何去除溶液中剩余的溶剂 35+3 77777sss 2022-08-10 6/300 2022-08-11 21:04 by Nanamiwww
[基金申请] 那些花儿! +6 了却无痕 2022-08-10 16/800 2022-08-11 09:17 by 医学生的成长
[基金申请] 收到了调查问卷,收到了kyzy的成果认领,没有其它任何消息 +4 huozhisheng 2022-08-09 10/500 2022-08-09 21:42 by huozhisheng
信息提示
请填处理意见