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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 理论与计算方法 » 用matlab求动力学参数
71/1返回列表
查看: 1922  |  回复: 6
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
【悬赏金币】回答本帖问题，作者自制奶酪将赠送您 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)
回复此楼

» 猜你喜欢
简单快乐
1楼 2021-01-05 15:31:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

独孤神宇

版主 (知名作家)

【答案】应助回帖

感谢参与，应助指数 +1
https://blog.sina.com.cn/s/blog_c0cb8ce60102ysqu.html
赞一下 回复此楼
数值计算
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
赞一下(1人) 回复此楼
简单快乐
4楼2021-01-05 21:41:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

自制奶酪

铜虫 (小有名气)
引用回帖:
2楼: Originally posted by 独孤神宇 at 2021-01-05 21:04:32
https://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的回帖
相关版块跳转 我要订阅楼主 自制奶酪 的主题更新
71/1返回列表
不应助 确定回帖应助 (注意：应助才可能被奖励，但不允许灌水，必须填写15个字符以上)
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭